OilfieldToolsNet/ResInsight
4
0
Fork 0
mirror of https://github.com/OPM/ResInsight.git synced 2025-01-10 08:03:05 -06:00
Code Issues Packages Projects Releases Wiki Activity

#2980 Include interpolated previous normal for determining direction of current normal

Browse Source
This commit is contained in:
Gaute Lindkvist 2018-05-25 12:38:25 +02:00
parent 70bcd25c86
commit cf64f5f70d
1 changed files with 3 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
ApplicationCode/ReservoirDataModel/RigWellPathGeometryTools.cpp
View File

@ -85,6 +85,7 @@ std::vector<cvf::Vec3d> RigWellPathGeometryTools::interpolateUndefinedNormals(co
{
std::vector<cvf::Vec3d> interpolated(normals);
cvf::Vec3d lastNormal(0, 0, 0);
cvf::Vec3d lastNormalIncludingInterpolated(0, 0, 0);
double distanceFromLast = 0.0;
for (size_t i = 0; i < normals.size(); ++i)
@ -123,7 +124,7 @@ std::vector<cvf::Vec3d> RigWellPathGeometryTools::interpolateUndefinedNormals(co
currentNormal = nextNormal;
}
}
if (i > 0 && currentNormal * lastNormal < -std::numeric_limits<double>::epsilon())
if (i > 0 && currentNormal * lastNormalIncludingInterpolated < -std::numeric_limits<double>::epsilon())
{
currentNormal *= -1.0;
}
@ -132,6 +133,7 @@ std::vector<cvf::Vec3d> RigWellPathGeometryTools::interpolateUndefinedNormals(co
lastNormal = currentNormal;
distanceFromLast = 0.0; // Reset distance
}
lastNormalIncludingInterpolated = currentNormal;
interpolated[i] = currentNormal;
}
return interpolated;