diff options
| -rw-r--r-- | lib/orthogonal.h | 2 | ||||
| -rw-r--r-- | 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"; |