- added support for computing spatial fourier transforms

 - measurements now are functions passed to wolff in array

1 files changed, 2 insertions, 1 deletions

diff --git a/lib/orthogonal.h b/lib/orthogonal.h
index d8ad33d..340ee2c 100644
--- a/lib/orthogonal.h
+++ b/lib/orthogonal.h
@@ -184,9 +184,10 @@ orthogonal_t <q, double> generate_rotation_perturbation (gsl_rng *r, const state
   }
 
   double v2 = 0;
+  double factor = gsl_ran_ugaussian(r);
 
   for (q_t i = 0; i < q; i++) {
-    tmp_v.x[i] = (tmp_v.x[i] - tmpM * s->M.x[i] / M2) + epsilon * gsl_ran_ugaussian(r);
+    tmp_v.x[i] = (tmp_v.x[i] - tmpM * s->M.x[i] / M2) + epsilon * factor * s->M.x[i] / sqrt(M2);
     v2 += tmp_v.x[i] * tmp_v.x[i];
   }