From 30d0fee3be1b899e93c5af7cf9de585071bacd44 Mon Sep 17 00:00:00 2001
From: Jaron Kent-Dobias <jaron@kent-dobias.com>
Date: Fri, 2 Jul 2021 17:10:44 +0200
Subject: Work on the Ciamarra model.

---
 distinguishable.cpp | 16 +++++++++++++---
 1 file changed, 13 insertions(+), 3 deletions(-)

(limited to 'distinguishable.cpp')

diff --git a/distinguishable.cpp b/distinguishable.cpp
index 22f228a..b739927 100644
--- a/distinguishable.cpp
+++ b/distinguishable.cpp
@@ -83,13 +83,23 @@ int main() {
     */
     /*
      */
-    s.setInitialPosition();
     std::cout << T << " ";
-    for (unsigned i = 0; i < 1e4; i++) {
+    s.setInitialPosition();
+    auto start = std::chrono::high_resolution_clock::now();
+    for (unsigned i = 0; i < 1e5; i++) {
+      for (unsigned j = 0; j < s.vertices.size(); j++) {
+        s.tryRandomSwap(T, r);
+      }
+      /*
       Matrix<D> m = r.pick(ms);
       Vertex<D, DistinguishableState>& v = r.pick(s.vertices);
       s.wolff(Transformation<D>(L, m, v.position - m * v.position), T, r);
-      std::cout << s.selfIntermediateScattering() << " ";
+      */
+      if (i % 10 == 0) {
+        auto stop = std::chrono::high_resolution_clock::now();
+        auto duration = duration_cast<std::chrono::microseconds>(stop - start);
+        std::cout << duration.count() << " " << s.selfIntermediateScattering(ms) << " ";
+      }
     }
     /*
     for (unsigned i = 0; i < 1e5; i++) {
-- 
cgit v1.2.3-54-g00ecf