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

1 files changed, 4 insertions, 0 deletions

diff --git a/lib/cluster_finite.h b/lib/cluster_finite.h
index b2d764e..117087e 100644
--- a/lib/cluster_finite.h
+++ b/lib/cluster_finite.h
@@ -45,6 +45,10 @@ typedef struct {
   q_t *R;
   v_t *B;
   v_t *M;
+  double *ReF;
+  double *ImF;
+  double *precomputed_cos;
+  double *precomputed_sin;
 } state_finite_t;
 
 v_t flip_cluster_finite(state_finite_t *s, v_t v0, R_t rot, gsl_rng *r);