applications/solvers/surfaceTracking/freeSurface/boundaryProcessorFaPatchPoints.H

   1 // Boundary processor patch points
   2 {
   3     const labelList& curPointEdges = pointEdges[curPoint];
   4
   5     label patchID = -1;
   6     label edgeID = -1;
   7
   8     forAll (curPointEdges, edgeI)
   9     {
  10         label curEdge = curPointEdges[edgeI];
  11
  12         if (edgeFaces[curEdge].size() == 1)
  13         {
  14             forAll(aMesh().boundary(), patchI)
  15             {
  16                 const labelList& curEdges = aMesh().boundary()[patchI];
  17
  18                 label index = findIndex(curEdges, curEdge);
  19
  20                 if (index != -1)
  21                 {
  22                     if
  23                     (
  24                         aMesh().boundary()[patchI].type()
  25                      != processorFaPatch::typeName
  26                     )
  27                     {
  28                         patchID = patchI;
  29                         edgeID = index;
  30                         break;
  31                     }
  32                 }
  33             }
  34         }
  35     }
  36
  37     if (patchID != -1)
  38     {
  39         vector mirrorCtrlPoint =
  40             patchMirrorPoints[patchID][edgeID];
  41
  42         lsPoints[curPatchPoint].setSize(lsPoints[curPatchPoint].size() + 1);
  43         lsPoints[curPatchPoint][lsPoints[curPatchPoint].size() - 1] =
  44             mirrorCtrlPoint;
  45     }
  46 }