void Foam::pointPatchDist::correct()
{
const pointBoundaryMesh& pbm = mesh().boundary();
label nPoints = 0;
forAllConstIter(labelHashSet, patchIDs_, iter)
{
label patchI = iter.key();
nPoints += pbm[patchI].meshPoints().size();
}
pointField points(points_.size());
points_.copyInto(points.begin());
externalPointEdgePoint::trackingData td(points);
// Set initial changed points to all the patch points(if patch present)
List<externalPointEdgePoint> wallInfo(nPoints);
labelList wallPoints(nPoints);
nPoints = 0;
forAllConstIter(labelHashSet, patchIDs_, iter)
{
label patchI = iter.key();
// Retrieve the patch now we have its index in patches.
const labelList& mp = pbm[patchI].meshPoints();
forAll(mp, ppI)
