From f72890584fbf7a8957b5aced448cb008ca408913 Mon Sep 17 00:00:00 2001
From: Kristian Bendiksen <kristian.bendiksen@gmail.com>
Date: Wed, 19 May 2021 15:02:47 +0200
Subject: [PATCH] StimPlan Model Calculator: guard against inf/nan in stress
 calculation.

---
 .../StimPlanModel/RimStimPlanModelCalculator.cpp             | 5 +++++
 1 file changed, 5 insertions(+)

diff --git a/ApplicationLibCode/ProjectDataModel/StimPlanModel/RimStimPlanModelCalculator.cpp b/ApplicationLibCode/ProjectDataModel/StimPlanModel/RimStimPlanModelCalculator.cpp
index 89e24b4193..7b8ef268af 100644
--- a/ApplicationLibCode/ProjectDataModel/StimPlanModel/RimStimPlanModelCalculator.cpp
+++ b/ApplicationLibCode/ProjectDataModel/StimPlanModel/RimStimPlanModelCalculator.cpp
@@ -436,6 +436,11 @@ bool RimStimPlanModelCalculator::calculateStressWithGradients( std::vector<doubl
         double Sv_diff               = 0.0;
         double deltaHorizontalStress = poissonsRatio / ( 1.0 - poissonsRatio ) * ( Sv_diff - biot * pressureDiff ) +
                                        ( biot * pressureDiff );
+        if ( std::isnan( deltaHorizontalStress ) || std::isinf( deltaHorizontalStress ) )
+        {
+            RiaLogging::warning( "Invalid horizontal stress delta calculated. Setting to zero." );
+            deltaHorizontalStress = 0.0;
+        }
 
         double depletionStress = Sh_init + deltaHorizontalStress;
         stress.push_back( RiaEclipseUnitTools::barToPsi( depletionStress ) );