From a43ff1f98e9b9814f858bccb11c174b418458491 Mon Sep 17 00:00:00 2001
From: Jaron Kent-Dobias <jaron@kent-dobias.com>
Date: Wed, 10 Oct 2018 21:45:32 -0400
Subject: big rearrangement of files to make libraries and example (research)
 files clearer, and changed to c++ std lib random numbers

---
 examples/src/models/potts/symmetric.hpp | 51 +++++++++++++++++++++++++++++++++
 1 file changed, 51 insertions(+)
 create mode 100644 examples/src/models/potts/symmetric.hpp

(limited to 'examples/src/models/potts/symmetric.hpp')

diff --git a/examples/src/models/potts/symmetric.hpp b/examples/src/models/potts/symmetric.hpp
new file mode 100644
index 0000000..8636f15
--- /dev/null
+++ b/examples/src/models/potts/symmetric.hpp
@@ -0,0 +1,51 @@
+
+#pragma once
+
+#include <stdlib.h>
+#include <array>
+#include <wolff/types.h>
+#include "potts.hpp"
+
+template <q_t q>
+class symmetric_t : public std::array<q_t, q> {
+  public:
+
+    symmetric_t() {
+      for (q_t i = 0; i < q; i++) {
+        (*this)[i] = i;
+      }
+    }
+
+    potts_t<q> act(const potts_t<q> &s) const {
+      return potts_t<q>((*this)[s.x]);
+    }
+
+    symmetric_t<q> act(const symmetric_t<q>& r) const {
+      symmetric_t<q> r_rot;
+      for (q_t i = 0; i < q; i++) {
+        r_rot[i] = (*this)[r[i]];
+      }
+
+      return r_rot;
+    }
+
+    potts_t<q> act_inverse(const potts_t<q>& s) const {
+      for (q_t i = 0; i < q; i++) {
+        if ((*this)[i] == s.x) {
+          return potts_t<q>(i);
+        }
+      }
+
+      exit(EXIT_FAILURE);
+    }
+
+    symmetric_t<q> act_inverse(const symmetric_t<q>& r) const {
+      symmetric_t<q> r_rot;
+      for (q_t i = 0; i < q; i++) {
+        r_rot[(*this)[i]] = r[i];
+      }
+
+      return r_rot;
+    }
+};
+
-- 
cgit v1.2.3-70-g09d2