From c48fd16fe1554c88c79a1f0d50e81c803da8f61f Mon Sep 17 00:00:00 2001
From: Jaron Kent-Dobias <jaron@kent-dobias.com>
Date: Tue, 24 Jul 2018 13:12:35 -0400
Subject: implemented updating the first fourier moment in wolff_finite, but
 also make wolff_finite obselete by adding a hacky preprocessor method for
 making wolff.h as efficient

---
 lib/wolff.h | 4 ++++
 1 file changed, 4 insertions(+)

(limited to 'lib/wolff.h')

diff --git a/lib/wolff.h b/lib/wolff.h
index 8286024..da4b7b6 100644
--- a/lib/wolff.h
+++ b/lib/wolff.h
@@ -5,6 +5,10 @@
 template <class R_t, class X_t>
 void wolff(count_t N, state_t <R_t, X_t> *s, std::function <R_t(gsl_rng *, X_t s0)> gen_R, std::function <void(const state_t <R_t, X_t> *)> measurements, gsl_rng *r, bool silent) {
 
+#ifdef FINITE_STATES
+  initialize_probs(s->J, s->H, s->T);
+#endif
+
   if (!silent) printf("\n");
   for (count_t steps = 0; steps < N; steps++) {
     if (!silent) printf("\033[F\033[JWOLFF: step %" PRIu64 " / %" PRIu64 ": E = %.2f, S = %" PRIv "\n", steps, N, s->E, s->last_cluster_size);
-- 
cgit v1.2.3-70-g09d2