1 files changed, 0 insertions, 68 deletions
diff --git a/lib/net_fracture.c b/lib/net_fracture.c
deleted file mode 100644
index 65ede9b..0000000
--- a/lib/net_fracture.c
+++ /dev/null
@@ -1,68 +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;
-}