Transmissibility equations: guard against NaN in very small intersect…

…ions.
OPM · Mar 15, 2024 · 1c56b9e · 1c56b9e
1 parent 5446a9a
commit 1c56b9e
Showing 1 changed file with 3 additions and 0 deletions.
diff --git a/...ode/ReservoirDataModel/Completions/RigEclipseToStimPlanCellTransmissibilityCalculator.cpp b/...ode/ReservoirDataModel/Completions/RigEclipseToStimPlanCellTransmissibilityCalculator.cpp
@@ -263,6 +263,9 @@ void RigEclipseToStimPlanCellTransmissibilityCalculator::calculateStimPlanCellsM
         for ( double lengtXarea : lengthXareaOfFractureParts )
             totalAreaXLength += lengtXarea;
 
+        // Guard against numerical issues for very small intersections
+        if ( std::isnan( totalAreaXLength ) || std::isnan( fractureArea ) ) continue;
+
         double fractureAreaWeightedlength = totalAreaXLength / fractureArea;
 
         // Transmissibility for inactive cells is set to zero