From c3d77467e1d6e7664764a9cb6263ff31ca0a4688 Mon Sep 17 00:00:00 2001
From: laurence <laurence>
Date: Thu, 9 Feb 2012 12:26:21 +0000
Subject: [PATCH] ENH: cvMesh: Change wellInside call to findSurfaceNearest

Do not need need to know whether point is inside or outside.
---
 .../conformalVoronoiMeshConformToSurface.C      | 17 +++++++++++++----
 1 file changed, 13 insertions(+), 4 deletions(-)

diff --git a/applications/utilities/mesh/generation/cvMesh/conformalVoronoiMesh/conformalVoronoiMesh/conformalVoronoiMeshConformToSurface.C b/applications/utilities/mesh/generation/cvMesh/conformalVoronoiMesh/conformalVoronoiMesh/conformalVoronoiMeshConformToSurface.C
index 79ca95a03ac..349a7299daf 100644
--- a/applications/utilities/mesh/generation/cvMesh/conformalVoronoiMesh/conformalVoronoiMesh/conformalVoronoiMeshConformToSurface.C
+++ b/applications/utilities/mesh/generation/cvMesh/conformalVoronoiMesh/conformalVoronoiMesh/conformalVoronoiMeshConformToSurface.C
@@ -184,14 +184,23 @@ void Foam::conformalVoronoiMesh::buildSurfaceConformation
         vit++
     )
     {
-        if (vit->internalPoint())
+        if (vit->internalPoint() && !vit->nearBoundary())
         {
             const Foam::point& pt = topoint(vit->point());
             const scalar range = sqr(2.0*targetCellSize(pt));
 
-            bool closeToBoundary = geometryToConformTo_.wellInside(pt, range);
+            pointIndexHit pHit;
+            label hitSurface;
 
-            if (!closeToBoundary)
+            geometryToConformTo_.findSurfaceNearest
+            (
+                pt,
+                range,
+                pHit,
+                hitSurface
+            );
+
+            if (pHit.hit())
             {
                 vit->setNearBoundary();
                 countNearBoundaryVertices++;
@@ -2179,7 +2188,7 @@ void Foam::conformalVoronoiMesh::buildSizeAndAlignmentTree() const
             treeDataPoint(sizeAndAlignmentLocations_),
             overallBb,  // overall search domain
             10,         // max levels
-            10.0,       // maximum ratio of cubes v.s. cells
+            20.0,       // maximum ratio of cubes v.s. cells
             100.0       // max. duplicity; n/a since no bounding boxes.
         )
     );
-- 
GitLab