diff options
| author | Jaron Kent-Dobias <jaron@kent-dobias.com> | 2018-09-05 22:45:39 -0400 |
|---|---|---|
| committer | Jaron Kent-Dobias <jaron@kent-dobias.com> | 2018-09-05 22:45:39 -0400 |
| commit | 10510b7e2a09c102c307165e0cd716e17b17ef84 (patch) | |
| tree | de7838a18d015245e70ac2ce4c666ced2edd2ed5 /lib | |
| parent | 8a4d4521758f23c960986a0f59d5925eb1a4292b (diff) | |
| download | c++-10510b7e2a09c102c307165e0cd716e17b17ef84.tar.gz c++-10510b7e2a09c102c307165e0cd716e17b17ef84.tar.bz2 c++-10510b7e2a09c102c307165e0cd716e17b17ef84.zip | |
fixed a bug in perturbative generation and set up epsilon to be automatically set by T and H
Diffstat (limited to 'lib')
| -rw-r--r-- | lib/orthogonal.h | 2 |
1 files changed, 1 insertions, 1 deletions
diff --git a/lib/orthogonal.h b/lib/orthogonal.h index 829df10..4ae3ca6 100644 --- a/lib/orthogonal.h +++ b/lib/orthogonal.h @@ -182,7 +182,7 @@ orthogonal_t <q, double> generate_rotation_perturbation (gsl_rng *r, const vecto v2 = 0; - double factor = epsilon * gsl_ran_ugaussian(r); + 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 |