diff options
author | Jaron <jaron@kent-dobias.com> | 2016-11-08 10:18:15 -0500 |
---|---|---|
committer | Jaron <jaron@kent-dobias.com> | 2016-11-08 10:18:15 -0500 |
commit | 0cf3090b05aa4c35ba20f3773db4f95d4600b55e (patch) | |
tree | 620bbc5ca8b42715e76e16837b90894ac9909534 /src/net_fracture.c | |
parent | c5b44e1604456ea42899c7e6cb3bde14b5bb88da (diff) | |
download | fuse_networks-0cf3090b05aa4c35ba20f3773db4f95d4600b55e.tar.gz fuse_networks-0cf3090b05aa4c35ba20f3773db4f95d4600b55e.tar.bz2 fuse_networks-0cf3090b05aa4c35ba20f3773db4f95d4600b55e.zip |
restarted all data generation to correct for errors, made certain data (stress, energy, thresholds) longs
Diffstat (limited to 'src/net_fracture.c')
-rw-r--r-- | src/net_fracture.c | 4 |
1 files changed, 2 insertions, 2 deletions
diff --git a/src/net_fracture.c b/src/net_fracture.c index 3bfcd83..b5ab524 100644 --- a/src/net_fracture.c +++ b/src/net_fracture.c @@ -46,7 +46,7 @@ data_t *net_fracture(net_t *net, cholmod_common *c, double cutoff, uint_t refact uint_t last_broke = get_next_broken(net, currents, cutoff); - double sim_current; + long double sim_current; if (net->voltage_bound) { sim_current = conductivity; @@ -54,7 +54,7 @@ data_t *net_fracture(net_t *net, cholmod_common *c, double cutoff, uint_t refact sim_current = 1; } - data_update(data, last_broke, fabsl(sim_current * (net->thres)[last_broke] / currents[last_broke]), conductivity); + data_update(data, last_broke, fabsl(sim_current * (net->thres)[last_broke] / ((long double)currents[last_broke])), conductivity); free(voltages); free(currents); |