OilfieldToolsNet/ResInsight
4
0
Fork 0
mirror of https://github.com/OPM/ResInsight.git synced 2025-02-14 17:44:20 -06:00
Code Issues Packages Projects Releases Wiki Activity

#5899 Fracture: Use robust way to find diff between two angles

Browse Source
This commit is contained in:
Magne Sjaastad 2020-05-25 20:43:59 +02:00
parent f11d7d3832
commit 6bef85ba4c
1 changed files with 21 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

23
ApplicationCode/ProjectDataModel/Completions/RimFracture.cpp
View File

@ -381,13 +381,29 @@ double RimFracture::endMD() const
}
}
//--------------------------------------------------------------------------------------------------
/// https://stackoverflow.com/a/52432897
//--------------------------------------------------------------------------------------------------
double getAbsoluteDiff2Angles( const double x, const double y, const double c )
{
// c can be PI (for radians) or 180.0 (for degrees);
return c - fabs( fmod( fabs( x - y ), 2 * c ) - c );
}
//--------------------------------------------------------------------------------------------------
///
//--------------------------------------------------------------------------------------------------
double RimFracture::wellFractureAzimuthDiff() const
{
double wellDifference = fabs( wellAzimuthAtFracturePosition() - m_azimuth );
return wellDifference;
// Compute the relative difference between two lines
// See https://github.com/OPM/ResInsight/issues/5899
double angle1 = wellAzimuthAtFracturePosition();
double angle2 = m_azimuth;
double smallestDiffDegrees = getAbsoluteDiff2Angles( angle1, angle2, 180.0 );
return smallestDiffDegrees;
}
//--------------------------------------------------------------------------------------------------
@ -399,6 +415,9 @@ QString RimFracture::wellFractureAzimuthDiffText() const
return QString::number( wellDifference, 'f', 2 );
}
//--------------------------------------------------------------------------------------------------
///
//--------------------------------------------------------------------------------------------------
QString RimFracture::wellAzimuthAtFracturePositionText() const
{
double wellAzimuth = wellAzimuthAtFracturePosition();