diff options
| author | Jaron Kent-Dobias <jaron@kent-dobias.com> | 2018-09-08 19:10:08 -0400 |
|---|---|---|
| committer | Jaron Kent-Dobias <jaron@kent-dobias.com> | 2018-09-08 19:10:08 -0400 |
| commit | 749d387b48456b78e448ab33516c01de20a2624b (patch) | |
| tree | fa75fa90578b69ec81cb4bff0a10c98602222dcd /lib | |
| parent | 42efc30ddcb2d0b8471120d56d2c0c64d6866aa9 (diff) | |
| download | c++-749d387b48456b78e448ab33516c01de20a2624b.tar.gz c++-749d387b48456b78e448ab33516c01de20a2624b.tar.bz2 c++-749d387b48456b78e448ab33516c01de20a2624b.zip | |
small perturb change
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 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); } |