diff --git a/applications/utilities/mesh/generation/cvMesh/conformalVoronoiMesh/conformationSurfaces/conformationSurfaces.H b/applications/utilities/mesh/generation/cvMesh/conformalVoronoiMesh/conformationSurfaces/conformationSurfaces.H index d8e5fd016616633358be21ac667efdac3eab22bf..ce56a977aa8e676c6f30f655a2730a90f2e4c3d8 100644 --- a/applications/utilities/mesh/generation/cvMesh/conformalVoronoiMesh/conformationSurfaces/conformationSurfaces.H +++ b/applications/utilities/mesh/generation/cvMesh/conformalVoronoiMesh/conformationSurfaces/conformationSurfaces.H @@ -285,7 +285,9 @@ public: ) const; //- Find the nearest points on each feature edge that is within - // a given distance from the sample point + // a given distance from the sample point. Will need to check for + // a hit or a miss because near edges may not have a nearest point + // on them which is perpendicular to the sample point. void findAllNearestEdges ( const point& sample, diff --git a/src/edgeMesh/extendedFeatureEdgeMesh/extendedFeatureEdgeMesh.C b/src/edgeMesh/extendedFeatureEdgeMesh/extendedFeatureEdgeMesh.C index 0ebb7b58ab6d0e54c7bc91092e02a51df9fadaa8..81f853aa417f1b3f2f0b264437bcc49ee47695d5 100644 --- a/src/edgeMesh/extendedFeatureEdgeMesh/extendedFeatureEdgeMesh.C +++ b/src/edgeMesh/extendedFeatureEdgeMesh/extendedFeatureEdgeMesh.C @@ -546,6 +546,49 @@ Foam::extendedFeatureEdgeMesh::extendedFeatureEdgeMesh } featurePointNormals_ = featurePointNormals; + + + // Create featurePointEdges_ + // For each feature point, stores a list of edges which are arranged in + // order according to their connectivity + +// featurePointEdges_.setSize(nonFeatureStart_); +// +// Info<< sFeatEds.size() << " " << surf.pointEdges().size() << " " +// << edges().size() << endl; +// +// forAll(sFeatEds, eI) +// { +// Info<< "Edge " << eI << " " << sFeatEds[eI] << " " +// << edges()[eI] << endl; +// } +// +// const edgeList& edges = eds; +// +// for (label i = 0; i < nonFeatureStart_; i++) +// { +// const labelList& ptEds = surf.pointEdges()[i]; +// +// DynamicList<label> tmpFtEdges; +// +// forAll(ptEds, eI) +// { +// const label edgeI = ptEds[eI]; +// const edge& e = sFeatEds[edgeI]; +// +// Info<< "Edges: " << edgeI << " " << e << endl; +// +// forAll(sFeatEds, fEdgeI) +// { +// if (edges[fEdgeI] == e) +// { +// tmpFtEdges.append(fEdgeI); +// } +// } +// } +// +// featurePointEdges_[i] = tmpFtEdges; +// } } @@ -821,7 +864,7 @@ void Foam::extendedFeatureEdgeMesh::allNearestFeatureEdges if (!hitPoint.hit()) { - nearHit = pointIndexHit(true, hitPoint.missPoint(), hitIndex); + nearHit = pointIndexHit(false, hitPoint.missPoint(), hitIndex); } else {