From eb65ff8ec89ab68a6a1bd062476244ab29bf39c3 Mon Sep 17 00:00:00 2001
From: h_vogt <h_vogt>
Date: Thu, 4 Aug 2016 19:30:47 +0200
Subject: [PATCH] inp.c, xpressn.c, cmath2.c, apply gauss1() instead of
 gauss0() -> reproducible start value if 'set rndseed=val' is given.

---
 src/frontend/inp.c              | 4 ++--
 src/frontend/numparam/xpressn.c | 4 ++--
 src/maths/cmaths/cmath2.c       | 2 +-
 3 files changed, 5 insertions(+), 5 deletions(-)

diff --git a/src/frontend/inp.c b/src/frontend/inp.c
index 8676effc3..923dda44c 100644
--- a/src/frontend/inp.c
+++ b/src/frontend/inp.c
@@ -1897,7 +1897,7 @@ agauss(double nominal_val, double abs_variation, double sigma)
 {
     double stdvar;
     stdvar = abs_variation / sigma;
-    return (nominal_val + stdvar * gauss0());
+    return (nominal_val + stdvar * gauss1());
 }
 
 
@@ -1906,7 +1906,7 @@ gauss(double nominal_val, double rel_variation, double sigma)
 {
     double stdvar;
     stdvar = nominal_val * rel_variation / sigma;
-    return (nominal_val + stdvar * gauss0());
+    return (nominal_val + stdvar * gauss1());
 }
 
 
diff --git a/src/frontend/numparam/xpressn.c b/src/frontend/numparam/xpressn.c
index a27bed2ec..37bdf5aa8 100644
--- a/src/frontend/numparam/xpressn.c
+++ b/src/frontend/numparam/xpressn.c
@@ -47,7 +47,7 @@ agauss(double nominal_val, double abs_variation, double sigma)
 {
     double stdvar;
     stdvar = abs_variation / sigma;
-    return (nominal_val + stdvar * gauss0());
+    return (nominal_val + stdvar * gauss1());
 }
 
 
@@ -56,7 +56,7 @@ gauss(double nominal_val, double rel_variation, double sigma)
 {
     double stdvar;
     stdvar = nominal_val * rel_variation / sigma;
-    return (nominal_val + stdvar * gauss0());
+    return (nominal_val + stdvar * gauss1());
 }
 
 
diff --git a/src/maths/cmaths/cmath2.c b/src/maths/cmaths/cmath2.c
index 85a1ea8f0..4b35335df 100644
--- a/src/maths/cmaths/cmath2.c
+++ b/src/maths/cmaths/cmath2.c
@@ -305,7 +305,7 @@ cx_sgauss(void *data, short int type, int length, int *newlength, short int *new
         d = alloc_d(length);
         *newtype = VF_REAL;
         for (i = 0; i < length; i++) {
-            d[i] = gauss0();
+            d[i] = gauss1();
         }
         return ((void *) d);
     }