From d19b858ec2376cefb38cce4849f8b3616b766c00 Mon Sep 17 00:00:00 2001
From: Will Bainbridge <http://cfd.direct>
Date: Thu, 24 Oct 2024 13:00:01 +0100
Subject: [PATCH] distributions::standardNormal: Prevent call of log with zero
 value

Resolves bug report https://bugs.openfoam.org/view.php?id=4167
---
 src/OpenFOAM/distributions/standardNormal/standardNormal.C | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/src/OpenFOAM/distributions/standardNormal/standardNormal.C b/src/OpenFOAM/distributions/standardNormal/standardNormal.C
index d2f2e611ea..fa5b0893a1 100644
--- a/src/OpenFOAM/distributions/standardNormal/standardNormal.C
+++ b/src/OpenFOAM/distributions/standardNormal/standardNormal.C
@@ -57,7 +57,7 @@ Foam::tmp<Foam::scalarField> Foam::distributions::standardNormal::approxErf
 
 Foam::scalar Foam::distributions::standardNormal::approxErfInv(const scalar y)
 {
-    const scalar l = log(1 - y*y), b = 2/(pi*a_) + l/2;
+    const scalar l = log(Foam::max(1 - y*y, small/2)), b = 2/(pi*a_) + l/2;
     return sign(y)*sqrt(-b + sqrt(b*b - l/a_));
 }