From 0cf3090b05aa4c35ba20f3773db4f95d4600b55e Mon Sep 17 00:00:00 2001 From: Jaron Date: Tue, 8 Nov 2016 10:18:15 -0500 Subject: restarted all data generation to correct for errors, made certain data (stress, energy, thresholds) longs --- src/net_fracture.c | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) (limited to 'src/net_fracture.c') 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); -- cgit v1.2.3-70-g09d2