From 10510b7e2a09c102c307165e0cd716e17b17ef84 Mon Sep 17 00:00:00 2001
From: Jaron Kent-Dobias <jaron@kent-dobias.com>
Date: Wed, 5 Sep 2018 22:45:39 -0400
Subject: fixed a bug in perturbative generation and set up epsilon to be
automatically set by T and H
---
lib/orthogonal.h | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
(limited to 'lib')
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
--
cgit v1.2.3-70-g09d2