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/geometry.c | 55 +++++++++++++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 55 insertions(+)
 create mode 100644 lib/geometry.c

(limited to 'lib/geometry.c')

diff --git a/lib/geometry.c b/lib/geometry.c
new file mode 100644
index 0000000..ec788f1
--- /dev/null
+++ b/lib/geometry.c
@@ -0,0 +1,55 @@
+
+#include "fracture.h"
+
+int edge_to_verts(unsigned int width, bool periodic, unsigned int edge,
+									bool index) {
+	assert(edge < pow(width, 2));
+
+	int x = edge / width + 1;
+	int y = edge % width + 1;
+
+	if (periodic) {
+		return (((index ^ (x % 2)) + 2 * ((y + (index ^ (!(x % 2)))) / 2) - 1) %
+								width +
+						(x - index) * width) /
+					 2;
+	} else {
+		return ((index ^ (x % 2)) + 2 * ((y + (index ^ (!(x % 2)))) / 2) +
+						(x - index) * (width + 1) - 1) /
+					 2;
+	}
+}
+
+int dual_edge_to_verts(unsigned int width, bool periodic, unsigned int edge,
+											 bool index) {
+	assert(edge < pow(width, 2));
+
+	int x = edge / width + 1;
+	int y = edge % width + 1;
+
+	if (periodic) {
+		return (((index ^ (!(x % 2))) + 2 * ((y + (index ^ (x % 2))) / 2) - 1) %
+								width +
+						(x - index) * width) /
+					 2;
+	} else {
+		return ((index ^ (!(x % 2))) + 2 * ((y + (index ^ (x % 2))) / 2) +
+						(x - index) * (width + 1) - 1) /
+					 2;
+	}
+}
+
+double dual_vert_to_coord(unsigned int width, bool periodic, unsigned int vert,
+													bool index) {
+	if (periodic) {
+		if (index)
+			return (2 * vert) % width + (2 * vert / width) % 2;
+		else
+			return 2 * vert / width;
+	} else {
+		if (index)
+			return (2 * vert) % (width + 1);
+		else
+			return (2 * vert) / (width + 1);
+	}
+}
-- 
cgit v1.2.3-70-g09d2