OilfieldToolsNet/ResInsight
4
0
Fork 0
mirror of https://github.com/OPM/ResInsight.git synced 2025-02-25 18:55:39 -06:00
Code Issues Packages Projects Releases Wiki Activity

#6544 Improve divide-by-zero check.

Browse Source
This commit is contained in:
Kristian Bendiksen 2020-09-30 10:33:22 +02:00
parent 50403e559b
commit 8086a1eb15
1 changed files with 4 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

7
ApplicationCode/GeoMech/GeoMechDataModel/RigFemPartResultCalculatorStressAnisotropy.cpp
View File

@ -218,6 +218,7 @@ RigFemScalarResultFrames*
s13.resize( valCount, 0.0 );
s23.resize( valCount, 0.0 );
double epsilon = 0.0000001;
#pragma omp parallel for schedule( dynamic )
for ( long vIdx = 0; vIdx < static_cast<long>( valCount ); ++vIdx )
{
@ -226,17 +227,17 @@ RigFemScalarResultFrames*
double diffS1 = s1t[vIdx] - s1b[vIdx];
double diffS2 = s2t[vIdx] - s2b[vIdx];
double diffS3 = s3t[vIdx] - s3b[vIdx];
if ( diffS1 + diffS2 != 0.0 )
if ( std::abs( diffS1 + diffS2 ) > epsilon )
s12[vIdx] = 2.0 * ( diffS1 - diffS2 ) / ( diffS1 + diffS2 );
else
s12[vIdx] = inf;
if ( diffS1 + diffS3 != 0.0 )
if ( std::abs( diffS1 + diffS3 ) > epsilon )
s13[vIdx] = 2.0 * ( diffS1 - diffS3 ) / ( diffS1 + diffS3 );
else
s13[vIdx] = inf;
if ( diffS2 + diffS3 != 0.0 )
if ( std::abs( diffS2 + diffS3 ) > epsilon )
s23[vIdx] = 2.0 * ( diffS2 - diffS3 ) / ( diffS2 + diffS3 );
else
s23[vIdx] = inf;