diff options
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); } |