From 1160baa61bad605cf8a1d583e8ae356a54a942df Mon Sep 17 00:00:00 2001
From: Jaron Kent-Dobias <jaron@kent-dobias.com>
Date: Thu, 26 Jul 2018 16:18:40 -0400
Subject: many changes, including new spin spaces and groups and cleaning up
 core library code

---
 lib/cluster.h | 16 +++++++++-------
 1 file changed, 9 insertions(+), 7 deletions(-)

(limited to 'lib/cluster.h')

diff --git a/lib/cluster.h b/lib/cluster.h
index c5f2be7..f948586 100644
--- a/lib/cluster.h
+++ b/lib/cluster.h
@@ -31,7 +31,9 @@ void flip_cluster(state_t <R_t, X_t> *state, v_t v0, R_t r, gsl_rng *rand) {
       R_old = state->R;
       marks[v] = true;
 
-      if (v == state->nv) {
+      bool v_is_ghost = (v == state->nv);
+
+      if (v_is_ghost) {
         R_new = r.act(R_old);
       } else {
         si_old = state->spins[v];
@@ -40,19 +42,18 @@ void flip_cluster(state_t <R_t, X_t> *state, v_t v0, R_t r, gsl_rng *rand) {
 
       for (const v_t &vn : state->g.v_adj[v]) {
         X_t sj;
+        bool vn_is_ghost = (vn == state->nv);
 
-        if (vn != state->nv) {
+        if (!vn_is_ghost) {
           sj = state->spins[vn];
         }
 
         double prob;
 
-        bool is_ext = (v == state->nv || vn == state->nv);
-
-        if (is_ext) {
+        if (v_is_ghost || vn_is_ghost) {
           X_t rs_old, rs_new;
           v_t non_ghost;
-          if (vn == state->nv) {
+          if (vn_is_ghost) {
             rs_old = R_old.act_inverse(si_old);
             rs_new = R_old.act_inverse(si_new);
             non_ghost = v;
@@ -61,6 +62,7 @@ void flip_cluster(state_t <R_t, X_t> *state, v_t v0, R_t r, gsl_rng *rand) {
             rs_new = R_new.act_inverse(sj);
             non_ghost = vn;
           }
+
           double dE = state->H(rs_old) - state->H(rs_new);
 #ifdef FINITE_STATES
           prob = H_probs[state_to_ind(rs_old)][state_to_ind(rs_new)];
@@ -97,7 +99,7 @@ void flip_cluster(state_t <R_t, X_t> *state, v_t v0, R_t r, gsl_rng *rand) {
         }
       }
 
-      if (v == state->nv) {
+      if (v_is_ghost) {
         state->R = R_new;
       } else {
         state->spins[v] = si_new;
-- 
cgit v1.2.3-70-g09d2