diff options
author | Jaron Kent-Dobias <jaron@kent-dobias.com> | 2018-07-18 15:37:27 -0400 |
---|---|---|
committer | Jaron Kent-Dobias <jaron@kent-dobias.com> | 2018-07-18 15:37:27 -0400 |
commit | cd85d3696be4a7ee65b0cf6669fc62dc4841aef3 (patch) | |
tree | 6429963091115b7b1ff0a0684ff8f01274a32042 /lib/cluster.h | |
parent | 722bc71ed8d4e1ae5616c5c8284fbffe21c4ffa4 (diff) | |
download | c++-cd85d3696be4a7ee65b0cf6669fc62dc4841aef3.tar.gz c++-cd85d3696be4a7ee65b0cf6669fc62dc4841aef3.tar.bz2 c++-cd85d3696be4a7ee65b0cf6669fc62dc4841aef3.zip |
- added support for computing spatial fourier transforms
- measurements now are functions passed to wolff in array
Diffstat (limited to 'lib/cluster.h')
-rw-r--r-- | lib/cluster.h | 8 |
1 files changed, 4 insertions, 4 deletions
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; } |