1 files changed, 0 insertions, 120 deletions
diff --git a/lib/include/wolff/cluster.hpp b/lib/include/wolff/cluster.hpp
deleted file mode 100644
index 3912446..0000000
--- a/lib/include/wolff/cluster.hpp
+++ /dev/null
@@ -1,120 +0,0 @@
-
-#ifndef WOLFF_CLUSTER_H
-#define WOLFF_CLUSTER_H
-
-#include <random>
-#include <cmath>
-#include <vector>
-#include <queue>
-
-#include "system.hpp"
-
-namespace wolff {
-
-template <class R_t, class X_t>
-void system<R_t, X_t>::flip_cluster(v_t i0, const R_t& r,
-                  std::mt19937& rng, measurement<R_t, X_t>& A) {
-  std::uniform_real_distribution<double> dist(0.0, 1.0);
-
-  std::queue<v_t> queue;
-  queue.push(i0);
-
-  std::vector<bool> marks(G.nv, false);
-
-  while (!queue.empty()) {
-    v_t i = queue.front();
-    queue.pop();
-
-    if (!marks[i]) { // don't reprocess anyone we've already visited!
-      marks[i] = true;
-
-      X_t si_new;
-#ifndef WOLFF_NO_FIELD
-      R_t s0_new;
-
-      bool we_are_ghost = (i == nv);
-
-      if (we_are_ghost) {
-        s0_new = r.act(s0);
-      } else
-#endif
-      {
-        si_new = r.act(s[i]);
-      }
-
-      for (const v_t &j : G.adjacency[i]) {
-        double dE, p;
-
-#ifndef WOLFF_NO_FIELD
-        bool neighbor_is_ghost = (j == nv);
-
-        if (we_are_ghost || neighbor_is_ghost) {
-          X_t s0s_old, s0s_new;
-          v_t non_ghost;
-
-          if (neighbor_is_ghost) {
-            non_ghost = i;
-            s0s_old = s0.act_inverse(s[i]);
-            s0s_new = s0.act_inverse(si_new);
-          } else {
-            non_ghost = j;
-            s0s_old = s0.act_inverse(s[j]);
-            s0s_new = s0_new.act_inverse(s[j]);
-          }
-
-#ifdef WOLFF_SITE_DEPENDENCE
-          dE = B(non_ghost, s0s_old) - B(non_ghost, s0s_new);
-#else
-          dE = B(s0s_old) - B(s0s_new);
-#endif
-
-#ifdef WOLFF_USE_FINITE_STATES
-          p = Bp[s0s_old.enumerate()][s0s_new.enumerate()];
-#endif
-
-          // run measurement hooks for encountering a ghost bond
-          A.ghost_bond_visited(*this, non_ghost, s0s_old, s0s_new, dE);
-        } else // this is a perfectly normal bond!
-#endif
-        {
-#ifdef WOLFF_BOND_DEPENDENCE
-          dE = Z(i, s[i], j, s[j]) - Z(i, si_new, j, s[j]);
-#else
-          dE = Z(s[i], s[j]) - Z(si_new, s[j]);
-#endif
-
-#ifdef WOLFF_USE_FINITE_STATES
-          p = Zp[s[i].enumerate()][si_new.enumerate()][s[j].enumerate()];
-#endif
-
-          // run measurement hooks for encountering a plain bond
-          A.plain_bond_visited(*this, i, si_new, j, dE);
-        }
-
-#ifndef WOLFF_USE_FINITE_STATES
-        p = 1.0 - exp(-dE / T);
-#endif
-
-        if (dist(rng) < p) {
-          queue.push(j); // push the neighboring vertex to the queue 
-        }
-      }
-
-#ifndef WOLFF_NO_FIELD
-      if (we_are_ghost) {
-        A.ghost_site_transformed(*this, s0_new);
-        s0 = s0_new;
-      } else
-#endif
-      {
-        A.plain_site_transformed(*this, i, si_new);
-        s[i] = si_new;
-      }
-    }
-  }
-}
-
-}
-
-#endif
-