From e3fbb92e68f0410f106285c9a49ecf8cd0a488a9 Mon Sep 17 00:00:00 2001
From: Jaron Kent-Dobias <jaron@kent-dobias.com>
Date: Sat, 21 Jul 2018 19:43:16 -0400
Subject: added visualization, and started potts

---
 lib/potts.h | 100 ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 100 insertions(+)
 create mode 100644 lib/potts.h

(limited to 'lib/potts.h')

diff --git a/lib/potts.h b/lib/potts.h
new file mode 100644
index 0000000..e7f0899
--- /dev/null
+++ b/lib/potts.h
@@ -0,0 +1,100 @@
+#pragma once
+
+#include <cmath>
+#include <stdio.h>
+
+#include "types.h"
+
+/* The following is the minimum definition of a spin class.
+ *
+ * The class must contain an M_t and an F_t for holding the sum of an
+ * integer number of spins and a double-weighted number of spins,
+ * respectively.
+ *
+ * void init(X_t *p);
+ * void free_spin(X_t p);
+ * void free_spin(M_t p);
+ * void free_spin(F_t p);
+ * X_t copy(X_t x);
+ * void add(M_t *x1, int factor, X_t x2);
+ * void add(F_t *x1, double factor, X_t x2);
+ * M_t scalar_multiple(int factor, X_t x);
+ * F_t scalar_multiple(double factor, X_t x);
+ * double norm_squared(F_t x);
+ * void write_magnetization(M_t M, FILE *outfile);
+ *
+ */
+
+template <q_t q>
+class potts_t {
+  public:
+    q_t x;
+
+    typedef int *M_t;
+    typedef double *F_t;
+};
+
+template <q_t q>
+void init(potts_t <q> *p) {
+  p->x = 0;
+}
+
+template <q_t q>
+void free_spin(potts_t <q> s) {
+  // do nothing!
+}
+
+template <q_t q>
+void free_spin(typename potts_t<q>::M_t s) {
+  free(s);
+}
+
+template <q_t q>
+void free_spin(typename potts_t<q>::F_t s) {
+  free(s);
+}
+
+template <q_t q>
+potts_t <q> copy(potts_t <q> s) {
+  return s;
+}
+
+template <q_t q>
+void add(typename potts_t<q>::M_t s1, int a, potts_t <q> s2) {
+  s1[s2.x] += a;
+}
+
+template <q_t q>
+void add(typename potts_t<q>::F_t s1, double a, potts_t <q> s2) {
+  s1[s2.x] += a;
+}
+
+template <q_t q>
+typename potts_t<q>::M_t scalar_multiple(int factor, potts_t <q> s) {
+  int *M = (int *)calloc(q, sizeof(int));
+  M[s.x] += factor;
+  return M;
+}
+
+template <q_t q>
+typename potts_t<q>::F_t scalar_multiple(double factor, potts_t <q> s) {
+  int *F = (double *)calloc(q, sizeof(double));
+  M[s.x] += factor;
+  return M;
+}
+
+template <q_t q>
+double norm_squared(typename potts<q>::F_t s) {
+  double total = 0;
+  for (q_t i = 0; i < q; i++) {
+    total += pow(s[i], 2);
+  }
+
+  return total * (double)q / ((double)q - 1.0);
+}
+
+template <q_t q>
+void write_magnetization(typename potts_t<q>::M_t M, FILE *outfile) {
+  fwrite(&M, sizeof(int), q, outfile);
+}
+
-- 
cgit v1.2.3-70-g09d2