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

#2945 Flip normals if the average well path normals doesn't match the requested direction.

Browse Source
This commit is contained in:
Gaute Lindkvist 2018-05-24 13:57:36 +02:00
parent 88c963dd38
commit 70bcd25c86
1 changed files with 16 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

21
ApplicationCode/ReservoirDataModel/RigWellPathGeometryTools.cpp
View File

@ -38,14 +38,16 @@ std::vector<cvf::Vec3d> RigWellPathGeometryTools::calculateLineSegmentNormals(co
pointNormals.reserve(vertices.size());
cvf::Vec3d up(0, 0, 1);
const cvf::Vec3d up(0, 0, 1);
const cvf::Vec3d rotatedUp = up.getTransformedVector(cvf::Mat3d::fromRotation(cvf::Vec3d(0.0, 1.0, 0.0), planeAngle));
cvf::Vec3d dominantDirection = estimateDominantDirectionInXYPlane(vertices);
const cvf::Vec3d dominantDirection = estimateDominantDirectionInXYPlane(vertices);
const cvf::Vec3d projectionPlaneNormal = (up ^ dominantDirection).getNormalized();
CVF_ASSERT(projectionPlaneNormal * dominantDirection <= std::numeric_limits<double>::epsilon());
cvf::Vec3d lastNormal;
double sumDotWithRotatedUp = 0.0;
for (size_t i = 0; i < vertices.size() - 1; ++i)
{
cvf::Vec3d p1 = vertices[i];
@ -60,10 +62,19 @@ std::vector<cvf::Vec3d> RigWellPathGeometryTools::calculateLineSegmentNormals(co
normal = normal.getTransformedVector(cvf::Mat3d::fromRotation(tangent, planeAngle));
}
pointNormals.push_back(normal);
lastNormal = normal;
sumDotWithRotatedUp += normal * rotatedUp;
}
pointNormals.push_back(pointNormals.back());
if (sumDotWithRotatedUp < 0.0)
{
for (cvf::Vec3d& normal : pointNormals)
{
normal *= -1.0;
}
}
pointNormals.push_back(lastNormal);
return interpolateUndefinedNormals(up, pointNormals, vertices);
}