From 749d387b48456b78e448ab33516c01de20a2624b Mon Sep 17 00:00:00 2001
From: Jaron Kent-Dobias <jaron@kent-dobias.com>
Date: Sat, 8 Sep 2018 19:10:08 -0400
Subject: small perturb change

---
 lib/orthogonal.h | 2 +-
 src/wolff_On.cpp | 4 ++--
 2 files changed, 3 insertions(+), 3 deletions(-)

diff --git a/lib/orthogonal.h b/lib/orthogonal.h
index 4ae3ca6..e1bf33a 100644
--- a/lib/orthogonal.h
+++ b/lib/orthogonal.h
@@ -185,7 +185,7 @@ orthogonal_t <q, double> generate_rotation_perturbation (gsl_rng *r, const vecto
   double factor = gsl_ran_gaussian(r, epsilon);
 
   for (q_t i = 0; i < q; i++) {
-    m[0][i] = m[0][i] + factor * v[i]; // perturb orthogonal vector in original direction
+    m[0][i] += factor * v[i]; // perturb orthogonal vector in original direction
     v2 += pow(m[0][i], 2);
   }
 
diff --git a/src/wolff_On.cpp b/src/wolff_On.cpp
index 0de9d42..2356875 100644
--- a/src/wolff_On.cpp
+++ b/src/wolff_On.cpp
@@ -153,10 +153,10 @@ int main(int argc, char *argv[]) {
       Hish = sqrt(H2);
     }
 
-    epsilon = sqrt((N_COMP - 1) * T / (D + Hish / 2));
+    epsilon = sqrt((N_COMP - 1) * T / (D + Hish / 2)) / 2;
 
     gen_R = std::bind(generate_rotation_perturbation <N_COMP>, std::placeholders::_1, std::placeholders::_2, epsilon, order);
-    pert_type = "PERTURB4";
+    pert_type = "PERTURB5";
   } else {
     gen_R = generate_rotation_uniform <N_COMP>;
     pert_type = "UNIFORM";
-- 
cgit v1.2.3-70-g09d2