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

#1462 Add functionality to check if a given coordinate is within a cell

Browse Source
This commit is contained in:
Bjørnar Grip Fjær
2017-05-15 15:49:49 +02:00
parent 9bf361ed64
commit d0833adb4c
1 changed files with 23 additions and 0 deletions
Download Patch File Download Diff File Expand all files Collapse all files

23
ApplicationCode/ReservoirDataModel/RigWellLogExtractionTools.h
View File

@@ -21,6 +21,7 @@
#include "cvfBoundingBox.h" #include "cvfBoundingBox.h"
#include "cvfGeometryTools.h" #include "cvfGeometryTools.h"
#include "cvfStructGrid.h" #include "cvfStructGrid.h"
#include "cvfRay.h"
//================================================================================================== //==================================================================================================
/// Internal class for intersection point info /// Internal class for intersection point info
@@ -90,6 +91,28 @@ struct RigHexIntersector
return intersectionCount; return intersectionCount;
} }
static bool isPointInCell(const cvf::Vec3d point, const cvf::Vec3d hexCorners[8], const size_t hexIndex)
{
cvf::Ray ray;
ray.setOrigin(point);
for (int face = 0; face < 6; ++face)
{
cvf::ubyte faceVertexIndices[4];
cvf::StructGridInterface::cellFaceVertexIndices(static_cast<cvf::StructGridInterface::FaceType>(face), faceVertexIndices);
cvf::Vec3d faceCenter = cvf::GeometryTools::computeFaceCenter(hexCorners[faceVertexIndices[0]], hexCorners[faceVertexIndices[1]], hexCorners[faceVertexIndices[2]], hexCorners[faceVertexIndices[3]]);
for (int i = 0; i < 4; ++i)
{
int next = i < 3 ? i + 1 : 0;
if (ray.triangleIntersect(hexCorners[faceVertexIndices[i]], hexCorners[faceVertexIndices[next]], faceCenter))
{
return true;
}
}
}
return false;
}
static bool isEqualDepth(double d1, double d2) static bool isEqualDepth(double d1, double d2)
{ {
double depthDiff = d1 - d2; double depthDiff = d1 - d2;