fix: Set probability density to zero for non-positive arguments

galacticusorg · Dec 11, 2024 · 76f43b6 · 76f43b6
1 parent c41896d
commit 76f43b6
Showing 1 changed file with 6 additions and 2 deletions.
diff --git a/source/statistics.distributions.lognormal.F90 b/source/statistics.distributions.lognormal.F90
@@ -212,8 +212,12 @@ double precision function logNormalDensity(self,x)
     class           (distributionFunction1DLogNormal), intent(inout) :: self
     double precision                                 , intent(in   ) :: x
 
-    logNormalDensity=+self%distributionFunction1DNormal%density(log(x)) &
-         &           /x
+    if (x > 0.0d0) then
+       logNormalDensity=+self%distributionFunction1DNormal%density(log(x)) &
+            &           /x
+    else
+       logNormalDensity=+0.0d0
+    end if
     return
   end function logNormalDensity