diff options
author | Jaron <jaron@kent-dobias.com> | 2016-10-17 08:52:20 -0400 |
---|---|---|
committer | Jaron <jaron@kent-dobias.com> | 2016-10-17 08:52:20 -0400 |
commit | 381d78796e69624635ecb35ecca833d5c6ece151 (patch) | |
tree | 7a61a9c5cc67fe363502ace2e05f4e5cc09d81c0 /src/fracture.c | |
parent | 52e79c6f3ba3360501a7849df3bef3c6292dca08 (diff) | |
download | fuse_networks-381d78796e69624635ecb35ecca833d5c6ece151.tar.gz fuse_networks-381d78796e69624635ecb35ecca833d5c6ece151.tar.bz2 fuse_networks-381d78796e69624635ecb35ecca833d5c6ece151.zip |
added processed data to repo, changed damage to include dangling bonds
Diffstat (limited to 'src/fracture.c')
-rw-r--r-- | src/fracture.c | 14 |
1 files changed, 13 insertions, 1 deletions
diff --git a/src/fracture.c b/src/fracture.c index 79558c2..74ed676 100644 --- a/src/fracture.c +++ b/src/fracture.c @@ -361,7 +361,19 @@ int main(int argc, char *argv[]) { if (save_conductivity) conductivity[i] = data->conductivity[max_pos]; - if (save_damage) damage[max_pos]++; + if (save_damage) { + uint_t num_dead = 0; + double *tmp_voltages = net_voltages(net, &c); + double *tmp_currents = net_currents(net, tmp_voltages, &c); + for (uint_t i = 0; i < g->ne; i++) { + if (!net->fuses[i] && fabs(tmp_currents[i]) < cutoff) { + num_dead++; + } + } + damage[max_pos + num_dead]++; + free(tmp_voltages); + free(tmp_currents); + } if (save_stress_field || save_voltage_field) { double *tmp_voltages = net_voltages(net, &c); |