diff options
| -rw-r--r-- | src/wolff_On.cpp | 9 |
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"; |