more sensible cutoff for scatter

2018-02-15 11:26:19 +01:00 · 2018-02-15 11:26:19 +01:00 · 9b1e72e7bc
parent 230ffe8f35
commit 9b1e72e7bc
1 changed files with 1 additions and 1 deletions
--- a/src/math.f90
+++ b/src/math.f90
@ -1783,7 +1783,7 @@ function math_sampleGaussOri(center,noise)
 real(pReal), dimension(3), parameter :: ORIGIN = 0.0_pReal
 real(pReal), dimension(5) :: rnd

- noScatter: if (abs(noise) < tol_math_check) then
+ noScatter: if (noise < 0.5_pReal*INRAD) then
   math_sampleGaussOri = center
 else noScatter
  ! Helming uses different distribution with Bessel functions