diff --git a/src/mesh/conformalVoronoiMesh/conformalVoronoiMesh/conformalVoronoiMesh.C b/src/mesh/conformalVoronoiMesh/conformalVoronoiMesh/conformalVoronoiMesh.C index 2526d419dd..1a687c266b 100644 --- a/src/mesh/conformalVoronoiMesh/conformalVoronoiMesh/conformalVoronoiMesh.C +++ b/src/mesh/conformalVoronoiMesh/conformalVoronoiMesh/conformalVoronoiMesh.C @@ -775,6 +775,8 @@ void Foam::conformalVoronoiMesh::insertFeaturePoints() Info<< " Inserted " << nFeatureVertices << " vertices" << endl; } + Info<< "SORT OUT FEATURE POINT DISTRIBUTION AND STORAGE" << endl; + featureVertices_.setSize(number_of_vertices()); label featPtI = 0; diff --git a/src/mesh/conformalVoronoiMesh/conformalVoronoiMesh/conformalVoronoiMeshConformToSurface.C b/src/mesh/conformalVoronoiMesh/conformalVoronoiMesh/conformalVoronoiMeshConformToSurface.C index 5e75c90b04..af791a589b 100644 --- a/src/mesh/conformalVoronoiMesh/conformalVoronoiMesh/conformalVoronoiMeshConformToSurface.C +++ b/src/mesh/conformalVoronoiMesh/conformalVoronoiMesh/conformalVoronoiMeshConformToSurface.C @@ -1632,6 +1632,8 @@ void Foam::conformalVoronoiMesh::addSurfaceAndEdgeHits { keepSurfacePoint = false; + // If the triggering Vertex is part of a feature point, allow it to + // conform to the surface if (vit->index() < startOfInternalPoints_) { surfaceHits.append(surfHit);