summaryrefslogtreecommitdiff
path: root/src/fracture.c
diff options
context:
space:
mode:
authorJaron Kent-Dobias <jaron@kent-dobias.com>2017-01-16 01:31:10 -0500
committerJaron Kent-Dobias <jaron@kent-dobias.com>2017-01-16 01:31:10 -0500
commit1e1fdfc2e3892667bccaf317a01defd8832041c7 (patch)
treecc5ef9adbfe4a8f11744f4b7afd23a37cfdd74d4 /src/fracture.c
parent57857b9ebfb2c0a78c2eb1128d3fb4ed8d597ec4 (diff)
downloadfuse_networks-1e1fdfc2e3892667bccaf317a01defd8832041c7.tar.gz
fuse_networks-1e1fdfc2e3892667bccaf317a01defd8832041c7.tar.bz2
fuse_networks-1e1fdfc2e3892667bccaf317a01defd8832041c7.zip
fixed voltage and torus conditions, current and free boundaries and broken right now
Diffstat (limited to 'src/fracture.c')
-rw-r--r--src/fracture.c16
1 files changed, 0 insertions, 16 deletions
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++) {