changed code to rely on jst

1 files changed, 0 insertions, 67 deletions

diff --git a/src/net_fracture.c b/src/net_fracture.c
deleted file mode 100644
index e7f18fc..0000000
--- a/src/net_fracture.c
+++ /dev/null
@@ -1,67 +0,0 @@
-
-#include "fracture.h"
-
-uint_t get_next_broken(net_t *net, double *currents, double cutoff) {
-	uint_t max_pos = UINT_MAX;
-	long double max_val = 0;
-
-	for (uint_t i = 0; i < net->graph->ne; i++) {
-		long double current = fabs(currents[i]);
-		bool broken = net->fuses[i];
-
-		if (!broken && current > cutoff) {
-			long double val = current / net->thres[i];
-
-			if (val > max_val) {
-				max_val = val;
-				max_pos = i;
-			}
-		}
-	}
-
-	if (max_pos == UINT_MAX) {
-		printf("GET_NEXT_BROKEN: currents is zero or NaN, no max_val found\n");
-		exit(EXIT_FAILURE);
-	}
-
-	return max_pos;
-}
-
-data_t *net_fracture(net_t *net, cholmod_common *c, double cutoff) {
-	data_t *data = data_create(net->graph->ne);
-
-	uint_t n = 0;
-	while (true) {
-		n++;
-		double *voltages = net_voltages(net, c);
-		double *currents = net_currents(net, voltages, c);
-
-		double conductivity = net_conductivity(net, voltages);
-
-		if (conductivity < cutoff) {
-			free(voltages);
-			free(currents);
-			break;
-		}
-
-		uint_t last_broke = get_next_broken(net, currents, cutoff);
-
-		long double sim_current;
-
-		if (net->voltage_bound) {
-			sim_current = conductivity;
-		} else {
-			sim_current = 1;
-		}
-
-		data_update(data, last_broke, fabsl(sim_current * (net->thres)[last_broke] / ((long double)currents[last_broke])), conductivity);
-
-		free(voltages);
-		free(currents);
-
-		break_edge(net, last_broke, c);
-	}
-
-	return data;
-}
-