From 901b9f16494f37890be17ef4bb66e6efc6873340 Mon Sep 17 00:00:00 2001
From: Jaron Kent-Dobias <jaron@kent-dobias.com>
Date: Fri, 10 Feb 2017 12:18:11 -0500
Subject: changed code to rely on jst

---
 lib/correlations.c | 47 +++++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 47 insertions(+)
 create mode 100644 lib/correlations.c

(limited to 'lib/correlations.c')

diff --git a/lib/correlations.c b/lib/correlations.c
new file mode 100644
index 0000000..63532ec
--- /dev/null
+++ b/lib/correlations.c
@@ -0,0 +1,47 @@
+
+#include "fracture.h"
+
+double *get_corr(net_t *instance, uint_t **dists, cholmod_common *c) {
+	uint_t nv = instance->graph->dnv;
+	uint_t ne = instance->graph->ne;
+	uint_t *ev = instance->graph->dev;
+	bool nulldists = false;
+	if (dists == NULL) {
+		dists = graph_dists(instance->graph, false);
+		nulldists = true;
+	}
+	double *corr = calloc(nv, sizeof(double));
+	uint_t *numat = calloc(nv, sizeof(uint_t));
+
+	for (uint_t i = 0; i < ne; i++) {
+		uint_t v1 = ev[2 * i];
+		uint_t v2 = ev[2 * i + 1];
+		for (uint_t j = 0; j < ne; j++) {
+			uint_t v3 = ev[2 * j];
+			uint_t v4 = ev[2 * j + 1];
+			uint_t dist1 = dists[v1][v3];
+			uint_t dist2 = dists[v1][v4];
+			uint_t dist3 = dists[v2][v3];
+			uint_t dist4 = dists[v2][v4];
+			uint_t dist = (dist1 + dist2 + dist3 + dist4) / 4;
+			corr[dist] += instance->fuses[i] && instance->fuses[j];
+			numat[dist]++;
+		}
+	}
+
+	for (uint_t i = 0; i < nv; i++) {
+		if (numat[i] > 0) {
+			corr[i] /= numat[i];
+		}
+	}
+
+	if (nulldists) {
+		for (int i = 0; i < nv; i++) {
+			free(dists[i]);
+		}
+		free(dists);
+	}
+
+	return corr;
+}
+
-- 
cgit v1.2.3-70-g09d2