applications/solvers/solidMechanics/elasticNonLinULSolidFoam/findBoundaryPoints.H

   1 //void volPointInterpolation::makeBoundaryAddressing() const
   2 //{
   3
   4 // Go through all patches and mark up the external edge points
   5 labelHashSet pointsCorrectionMap(max(mesh.nPoints()/10, 100));
   6
   7 const fvBoundaryMesh& bm = mesh.boundary();
   8
   9 forAll (bm, patchI)
  10 {
  11   // If the patch is empty, skip it
  12   // If the patch is coupled, and there are no cyclic parallels, skip it
  13   if
  14     (
  15      !isA<emptyFvPatch>(bm[patchI])
  16      && !(
  17           bm[patchI].coupled()
  18           //&& Pstream::parRun()
  19           //&& !mesh.parallelData().cyclicParallel()
  20           )
  21      )
  22     {
  23       const labelList& bp = bm[patchI].patch().boundaryPoints();
  24
  25       const labelList& meshPoints = bm[patchI].patch().meshPoints();
  26
  27       forAll (bp, pointI)
  28         {
  29           pointsCorrectionMap.insert(meshPoints[bp[pointI]]);
  30         }
  31     }
  32 }
  33
  34 // Grab the points to correct
  35 //boundaryPointsPtr_ = new labelList(pointsCorrectionMap.toc());
  36 //labelList& ptc = *boundaryPointsPtr_;
  37 labelList ptc(pointsCorrectionMap.toc());
  38
  39 sort(ptc);