From 90ae915b5a7961a36e6a33509b16229244c6615a Mon Sep 17 00:00:00 2001
From: Jaron Kent-Dobias <jaron@kent-dobias.com>
Date: Sun, 1 Jul 2018 00:52:31 -0400
Subject: fixed both the system for determining bond energy and how global
 state is tracked

---
 lib/symmetric.c | 43 +++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 43 insertions(+)

(limited to 'lib/symmetric.c')

diff --git a/lib/symmetric.c b/lib/symmetric.c
index 729b38c..4487538 100644
--- a/lib/symmetric.c
+++ b/lib/symmetric.c
@@ -25,3 +25,46 @@ q_t *symmetric_invert(q_t q, const q_t *g) {
   return g_inv;
 }
 
+void swap(q_t *q1, q_t *q2) {
+  q_t temp = *q1;
+  *q1 = *q2;
+  *q2 = temp;
+}
+
+R_t factorial(q_t q) {
+  if (q == 0) {
+    return 1;
+  } else {
+    return q * factorial(q - 1);
+  }
+}
+
+void permute(q_t *a, q_t l, q_t r, R_t pos, q_t *transformations) {
+  if (l == r - 1) {
+    for (q_t i = 0; i < r; i++) {
+      transformations[r * pos + i] = a[i];
+    }
+  } else {
+    for (q_t i = l; i < r; i++) {
+      swap((a+l), (a+i));
+      permute(a, l+1, r, pos + (i - l) * factorial(r - l - 1), transformations);
+      swap((a+l), (a+i));
+    }
+  }
+}
+
+q_t *symmetric_gen_transformations(q_t q) {
+  q_t *transformations = (q_t *)malloc(q * factorial(q) * sizeof(q_t));
+  q_t *tmp = (q_t *)malloc(q * sizeof(q_t));
+
+  for (q_t i = 0; i < q; i++) {
+    tmp[i] = i;
+  }
+
+  permute(tmp, 0, q, 0, transformations);
+
+  free(tmp);
+
+  return transformations;
+}
+
-- 
cgit v1.2.3-70-g09d2