From fa18734e2b1eb6f0de21f31240684519bc0c621a Mon Sep 17 00:00:00 2001 From: mattijs Date: Fri, 5 Sep 2008 11:07:12 +0100 Subject: [PATCH] removal of hanging points --- .../polyTopoChange/polyTopoChange/removeFaces.C | 20 ++++++++------------ 1 file changed, 8 insertions(+), 12 deletions(-) diff --git a/src/dynamicMesh/polyTopoChange/polyTopoChange/removeFaces.C b/src/dynamicMesh/polyTopoChange/polyTopoChange/removeFaces.C index 28fe2a2..fab2b74 100644 --- a/src/dynamicMesh/polyTopoChange/polyTopoChange/removeFaces.C +++ b/src/dynamicMesh/polyTopoChange/polyTopoChange/removeFaces.C @@ -807,7 +807,7 @@ void Foam::removeFaces::setRefinement // Edges to remove labelHashSet edgesToRemove(faceLabels.size()); - // Per face the region it is. -1 for removed faces, -2 for regions + // Per face the region it is in. -1 for removed faces, -2 for regions // consisting of single face only. labelList faceRegion(mesh_.nFaces(), -1); @@ -1258,10 +1258,15 @@ void Foam::removeFaces::setRefinement // are only used by 2 unremoved edges. { // Usage of points by non-removed edges. - labelList nEdgesPerPoint(mesh_.nPoints(), labelMax); + labelList nEdgesPerPoint(mesh_.nPoints()); const labelListList& pointEdges = mesh_.pointEdges(); + forAll(pointEdges, pointI) + { + nEdgesPerPoint[pointI] = pointEdges[pointI].size(); + } + forAllConstIter(labelHashSet, edgesToRemove, iter) { // Edge will get removed. @@ -1269,16 +1274,7 @@ void Foam::removeFaces::setRefinement forAll(e, i) { - label pointI = e[i]; - - if (nEdgesPerPoint[pointI] == labelMax) - { - nEdgesPerPoint[pointI] = pointEdges[pointI].size()-1; - } - else - { - nEdgesPerPoint[pointI]--; - } + nEdgesPerPoint[e[i]]--; } } -- 2.11.4.GIT