From baca97666aebc1926861338af312c0560bf923f7 Mon Sep 17 00:00:00 2001
From: Jaron Kent-Dobias <jaron@kent-dobias.com>
Date: Thu, 6 Sep 2018 23:36:19 -0400
Subject: tweaked the perturbation yet again

---
 src/wolff_On.cpp | 9 +++------
 1 file changed, 3 insertions(+), 6 deletions(-)

(limited to 'src')

diff --git a/src/wolff_On.cpp b/src/wolff_On.cpp
index 8a7b1a6..fada9bf 100644
--- a/src/wolff_On.cpp
+++ b/src/wolff_On.cpp
@@ -153,13 +153,10 @@ int main(int argc, char *argv[]) {
       Hish = sqrt(H2);
     }
 
-    if (Hish > 1.0) {
-      epsilon = sqrt(T / Hish);
-    } else {
-      epsilon = sqrt(T);
-    }
+    epsilon = sqrt(N_COMP * T / (D + Hish / 2));
+
     gen_R = std::bind(generate_rotation_perturbation <N_COMP>, std::placeholders::_1, std::placeholders::_2, epsilon, order);
-    pert_type = "PERTURB2";
+    pert_type = "PERTURB3";
   } else {
     gen_R = generate_rotation_uniform <N_COMP>;
     pert_type = "UNIFORM";
-- 
cgit v1.2.3-70-g09d2