diff options
| author | Jaron Kent-Dobias <jaron@kent-dobias.com> | 2018-09-06 23:41:22 -0400 |
|---|---|---|
| committer | Jaron Kent-Dobias <jaron@kent-dobias.com> | 2018-09-06 23:41:22 -0400 |
| commit | 42efc30ddcb2d0b8471120d56d2c0c64d6866aa9 (patch) | |
| tree | 33008ea0aebf84b29f0c56eb0d755f144aedbe4c | |
| parent | baca97666aebc1926861338af312c0560bf923f7 (diff) | |
| download | c++-42efc30ddcb2d0b8471120d56d2c0c64d6866aa9.tar.gz c++-42efc30ddcb2d0b8471120d56d2c0c64d6866aa9.tar.bz2 c++-42efc30ddcb2d0b8471120d56d2c0c64d6866aa9.zip | |
small changed to perturb again
| -rw-r--r-- | src/wolff_On.cpp | 4 |
1 files changed, 2 insertions, 2 deletions
diff --git a/src/wolff_On.cpp b/src/wolff_On.cpp index fada9bf..0de9d42 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 * T / (D + Hish / 2)); + epsilon = sqrt((N_COMP - 1) * T / (D + Hish / 2)); gen_R = std::bind(generate_rotation_perturbation <N_COMP>, std::placeholders::_1, std::placeholders::_2, epsilon, order); - pert_type = "PERTURB3"; + pert_type = "PERTURB4"; } else { gen_R = generate_rotation_uniform <N_COMP>; pert_type = "UNIFORM"; |