Commit 4d428e8f authored by mattijs's avatar mattijs
Browse files

BUG: meshSearch : extraneous loop

parent 5004aa10
...@@ -617,30 +617,25 @@ bool Foam::meshSearch::pointInCell(const point& p, label cellI) const ...@@ -617,30 +617,25 @@ bool Foam::meshSearch::pointInCell(const point& p, label cellI) const
{ {
label faceI = cFaces[i]; label faceI = cFaces[i];
const face& f = mesh_.faces()[faceI]; pointHit inter = mesh_.faces()[faceI].ray
(
ctr,
dir,
mesh_.points(),
intersection::HALF_RAY,
intersection::VECTOR
);
forAll(f, fp) if (inter.hit())
{ {
pointHit inter = f.ray scalar dist = inter.distance();
(
ctr,
dir,
mesh_.points(),
intersection::HALF_RAY,
intersection::VECTOR
);
if (inter.hit()) if (dist < magDir)
{ {
scalar dist = inter.distance(); // Valid hit. Hit face so point is not in cell.
intersection::setPlanarTol(oldTol);
if (dist < magDir) return false;
{
// Valid hit. Hit face so point is not in cell.
intersection::setPlanarTol(oldTol);
return false;
}
} }
} }
} }
......
Markdown is supported
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