diff options
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); |