From ef2adac0b99d8e87343a0e333d04d8afd51c4eaf Mon Sep 17 00:00:00 2001 From: laurence Date: Fri, 17 May 2013 11:44:37 +0100 Subject: [PATCH] BUG: Test for infinite vertex --- .../DelaunayMesh/DistributedDelaunayMesh.C | 22 +++++++++++++------ 1 file changed, 15 insertions(+), 7 deletions(-) diff --git a/applications/utilities/mesh/generation/foamyHexMesh/conformalVoronoiMesh/DelaunayMesh/DistributedDelaunayMesh.C b/applications/utilities/mesh/generation/foamyHexMesh/conformalVoronoiMesh/DelaunayMesh/DistributedDelaunayMesh.C index 8ba0f2d7c2..43e9737a6f 100644 --- a/applications/utilities/mesh/generation/foamyHexMesh/conformalVoronoiMesh/DelaunayMesh/DistributedDelaunayMesh.C +++ b/applications/utilities/mesh/generation/foamyHexMesh/conformalVoronoiMesh/DelaunayMesh/DistributedDelaunayMesh.C @@ -857,12 +857,12 @@ Foam::DistributedDelaunayMesh::rangeInsertReferredWithInfo << " nearest : " << nearV->info(); } } - else if (lt == Triangulation::OUTSIDE_CONVEX_HULL) - { - hint = this->insert(pointToInsert, hint); - - inserted = true; - } +// else if (lt == Triangulation::OUTSIDE_CONVEX_HULL) +// { +// hint = this->insert(pointToInsert, hint); +// +// inserted = true; +// } else if (lt == Triangulation::OUTSIDE_AFFINE_HULL) { WarningIn @@ -889,7 +889,15 @@ Foam::DistributedDelaunayMesh::rangeInsertReferredWithInfo for (size_t i = 0; i < V.size(); ++i) { - if (V[i]->real() || V[i]->hasFarPoint()) + if + ( + !Triangulation::is_infinite(V[i]) + && + ( + V[i]->real() + || V[i]->hasFarPoint() + ) + ) { hint = Triangulation::insert_in_hole (