applications/solvers/solidMechanics/elasticPlasticNonLinULSolidFoam/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);