Fixed Ciamarra model.

1 files changed, 5 insertions, 2 deletions

diff --git a/ciamarra.cpp b/ciamarra.cpp
index 8449cbb..6b932bf 100644
--- a/ciamarra.cpp
+++ b/ciamarra.cpp
@@ -58,6 +58,9 @@ int main() {
   } 
 
   std::cerr << "Found state with appropriate density." << std::endl;
+    if (!s.compatible()) {
+      std::cerr <<"whoops!" << std::endl;
+    }
 
   CiamarraSystem<D> s0 = s;
 
@@ -69,9 +72,9 @@ int main() {
       n++;
       std::cout << i << " " << s.overlap(s0) << std::endl;
     }
-    s.sweepLocal(r);
+//    s.sweepLocal(r);
 //    s.sweepSwap(r);
-//    s.swendsenWang(Transformation<D>(L, ms, r), r);
+    s.swendsenWang(Transformation<D>(L, ms, r), r);
   }
 
   return 0;