From 1e1fdfc2e3892667bccaf317a01defd8832041c7 Mon Sep 17 00:00:00 2001 From: Jaron Kent-Dobias Date: Mon, 16 Jan 2017 01:31:10 -0500 Subject: fixed voltage and torus conditions, current and free boundaries and broken right now --- src/fracture.c | 16 ---------------- 1 file changed, 16 deletions(-) (limited to 'src/fracture.c') diff --git a/src/fracture.c b/src/fracture.c index 1ed594f..6ad2f26 100644 --- a/src/fracture.c +++ b/src/fracture.c @@ -263,13 +263,6 @@ int main(int argc, char *argv[]) { graph_t *g = graph_create(lattice, boundary, L, use_dual, &c); net_t *net = net_create(g, inf, beta, crack_len, use_voltage_boundaries, &c); - if (net->marks[g->nv] != net->marks[g->nv+1]) { - if (save_crit_stress) crit_stress[i] = 0; - if (save_conductivity) conductivity[i] = 0; - if (save_damage) damage[net->num_broken]++; - if (save_energy) energy[i] = 0; - if (save_threshold) thresholds[i] = 0; - } else { net_t *tmp_net = net_copy(net, &c); data_t *data = net_fracture(tmp_net, &c, cutoff); net_free(tmp_net, &c); @@ -325,14 +318,6 @@ int main(int argc, char *argv[]) { } } - FILE *testout = fopen("test.txt", "w"); - double *tmp_voltage = net_voltages(net, &c); - for (uint_t j = 0; j < g->nv; j++) { - fprintf(testout, "%g ", tmp_voltage[j]); - } - fclose(testout); - free(tmp_voltage); - if (save_damage) { uint_t would_break = 0; double *tmp_voltage = net_voltages(net, &c); @@ -382,7 +367,6 @@ int main(int argc, char *argv[]) { } data_free(data); - } if (save_network) { FILE *net_out = fopen("network.txt", "w"); for (uint_t j = 0; j < g->nv; j++) { -- cgit v1.2.3-70-g09d2