math_sampleGaussOri did not treat the special case of zero noise (scatter)

2007-04-26 11:46:58 +00:00 · 2007-04-26 11:46:58 +00:00 · 11cc6d1ee1
parent 6fd473eaa7
commit 11cc6d1ee1
1 changed files with 5 additions and 0 deletions
--- a/trunk/math.f90
+++ b/trunk/math.f90
@ -725,6 +725,11 @@
 real(pReal), dimension(5) :: rnd
 real(pReal) noise,scatter,cosScatter
 integer(pInt) i
+
+if (noise<=0.01) then
+    math_sampleGaussOri = center
+    return
+endif

 ! Helming uses different distribution with Bessel functions
 ! therefore the gauss scatter width has to be scaled differently