summaryrefslogtreecommitdiff
diff options
context:
space:
mode:
-rw-r--r--src/wolff_On.cpp9
1 files changed, 3 insertions, 6 deletions
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";