Commit 6d224df7 authored by laurence's avatar laurence
Browse files

ENH: foamyHexMesh: change inside/outside test to check for nearest intersection

parent 5e3a0579
......@@ -683,14 +683,18 @@ Foam::Field<bool> Foam::conformationSurfaces::wellInOutSide
vector hitDir = info[0].rawPoint() - samplePts[i];
hitDir /= mag(hitDir) + SMALL;
if
pointIndexHit surfHit;
label hitSurface;
findSurfaceNearestIntersection
(
findSurfaceAnyIntersection
(
samplePts[i],
info[0].rawPoint() - 1e-3*mag(hitDir)*(hitDir)
)
)
samplePts[i],
info[0].rawPoint(),
surfHit,
hitSurface
);
if (surfHit.hit() && hitSurface != surfaces_[s])
{
continue;
}
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment