From cd85d3696be4a7ee65b0cf6669fc62dc4841aef3 Mon Sep 17 00:00:00 2001
From: Jaron Kent-Dobias <jaron@kent-dobias.com>
Date: Wed, 18 Jul 2018 15:37:27 -0400
Subject:  - added support for computing spatial fourier transforms  -
 measurements now are functions passed to wolff in array

---
 lib/cluster.h | 8 ++++----
 1 file changed, 4 insertions(+), 4 deletions(-)

(limited to 'lib/cluster.h')

diff --git a/lib/cluster.h b/lib/cluster.h
index 5f770ad..8061d34 100644
--- a/lib/cluster.h
+++ b/lib/cluster.h
@@ -52,7 +52,7 @@ template <class T>
 T subtract(T, T);
 
 template <class R_t, class X_t>
-v_t flip_cluster(state_t <R_t, X_t> *state, v_t v0, R_t r, gsl_rng *rand) {
+void flip_cluster(state_t <R_t, X_t> *state, v_t v0, R_t r, gsl_rng *rand) {
   v_t nv = 0;
 
   ll_t *stack = NULL;     // create a new stack
@@ -104,8 +104,8 @@ v_t flip_cluster(state_t <R_t, X_t> *state, v_t v0, R_t r, gsl_rng *rand) {
           double dE = state->H(rs_old) - state->H(rs_new);
           prob = 1.0 - exp(-dE / state->T);
 
-          subtract (state->M, rs_old);
-          add (state->M, rs_new);
+          subtract (&(state->M), rs_old);
+          add (&(state->M), rs_new);
           state->E += dE;
 
           free_spin (rs_old);
@@ -137,6 +137,6 @@ v_t flip_cluster(state_t <R_t, X_t> *state, v_t v0, R_t r, gsl_rng *rand) {
 
   free(marks);
 
-  return nv;
+  state->last_cluster_size = nv;
 }
 
-- 
cgit v1.2.3-70-g09d2