fixed memory leak due to not emptying stack in the case where the cluster flip is accepted in the hybrid method

1 files changed, 3 insertions, 0 deletions

diff --git a/lib/wolff_tools.c b/lib/wolff_tools.c
index 6b1e006..bab0908 100644
--- a/lib/wolff_tools.c
+++ b/lib/wolff_tools.c
@@ -206,6 +206,9 @@ uint32_t wolff_step(double T, double H, ising_state_t *s, sim_t sim, gsl_rng *r,
           s->spins[v] = !s->spins[v];
         }
       } else {
+        while (c->spins != NULL) {
+          stack_pop(&(c->spins));
+        }
         s->M += - sign(H) * 2 * c->dHb;
         s->H += 2 * (c->dJb + sign (H) * H * c->dHb);
       }