From c584052a004f6c7a2932b5290f57ea62bb722a26 Mon Sep 17 00:00:00 2001 From: Jaron Kent-Dobias Date: Thu, 19 Oct 2017 12:38:50 -0400 Subject: Revert "fixed a mistake in the way wolff cluster tries were counted" This reverts commit 020266dc3745afc4f9d21a15b40e3a83adbd92ff --- lib/wolff_tools.c | 5 ++--- 1 file changed, 2 insertions(+), 3 deletions(-) (limited to 'lib/wolff_tools.c') diff --git a/lib/wolff_tools.c b/lib/wolff_tools.c index 629da4e..bab0908 100644 --- a/lib/wolff_tools.c +++ b/lib/wolff_tools.c @@ -113,7 +113,6 @@ cluster_t *flip_cluster(const graph_t *g, const double *ps, bool *x, bool stop_o if (x[vn] == x0) { // if the neighboring site matches the flipping cluster... (*bond_counter)++; - c->nt++; if (gsl_rng_uniform(r) < prob) { // and with probability ps[e]... if (is_ext && stop_on_ghost) { @@ -193,7 +192,7 @@ uint32_t wolff_step(double T, double H, ising_state_t *s, sim_t sim, gsl_rng *r, cluster_t *c = flip_cluster(s->g, ps, s->spins, false, r); s->M += - sign(H) * 2 * c->dHb; s->H += 2 * (c->dJb + sign (H) * H * c->dHb); - n_flips = c->nt; + n_flips = c->nv; free(c); } @@ -214,7 +213,7 @@ uint32_t wolff_step(double T, double H, ising_state_t *s, sim_t sim, gsl_rng *r, s->H += 2 * (c->dJb + sign (H) * H * c->dHb); } - n_flips = c->nt; + n_flips = c->nv; free(c); } -- cgit v1.2.3-70-g09d2