From 66aa9399a7d2f5ba32d78397f2345d0a044c48d6 Mon Sep 17 00:00:00 2001
From: Martin Diehl <mail@martin-diehl.net>
Date: Sun, 17 Mar 2019 18:06:41 +0100
Subject: [PATCH] don'r rely on halton sequence

not really trustworthy as we saw some pattern for higher dimensional
data
---
 src/math.f90 | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/src/math.f90 b/src/math.f90
index f77f81f8b..8c0020799 100644
--- a/src/math.f90
+++ b/src/math.f90
@@ -1895,7 +1895,7 @@ real(pReal) function math_sampleGaussVar(meanvalue, stddev, width)
    myWidth = merge(width,3.0_pReal,present(width))                                                  ! use +-3*sigma as default value for scatter if not given
   
    do
-     rnd = halton([6,2])
+     call random_number(rnd)
      scatter = myWidth * (2.0_pReal * rnd(1) - 1.0_pReal)
      if (rnd(2) <= exp(-0.5_pReal * scatter ** 2.0_pReal)) exit                                     ! test if scattered value is drawn
    enddo