src/sampling/include/buildPatch.H

   1
   2 IOdictionary planeToPatchDict
   3 (
   4     IOobject
   5     (
   6         "planeToPatchDict",
   7         runTime.system(),
   8         runTime,
   9         IOobject::MUST_READ,
  10         IOobject::NO_WRITE
  11     )
  12 );
  13
  14 Switch planeMapping
  15 (
  16     planeToPatchDict.lookup("mappingPlane")
  17 );
  18
  19 Info << "Constructing planePatch" << endl;
  20 cuttingPlane cut(plane(planeToPatchDict), mesh);
  21
  22 primitiveFacePatch planePatch
  23 (
  24     cut.faces(),
  25     cut.points()
  26 );
  27
  28 Info << "Finding target patch" << endl;
  29
  30 word toPatchName
  31 (
  32     planeToPatchDict.lookup("targetPatchName")
  33 );
  34
  35 const fvPatchList& patches = mesh.boundary();
  36 bool targetPatchNamePresent(false);
  37 label targetPatchNumber = -1;
  38
  39 forAll(patches, patchI)
  40 {
  41     if (patches[patchI].name() == toPatchName)
  42     {
  43         targetPatchNamePresent = true;
  44         targetPatchNumber = patchI;
  45     }
  46 }
  47
  48 const labelList& cutCells(cut.cells());
  49
  50 if (!targetPatchNamePresent)
  51 {
  52     targetPatchNumber = 1;
  53     Serr<< "Target patch not present. " << endl;
  54 }
  55 else
  56 {
  57     Serr << "Target patch name: " << patches[targetPatchNumber].name()
  58          << endl;
  59 }
  60
  61 patchToPatchInterpolation interPatch
  62 (
  63     planePatch,
  64     patches[targetPatchNumber].patch()
  65 );
  66
  67
  68 PrimitivePatchInterpolation
  69 <
  70     PrimitivePatch<face, List, const pointField&>
  71 > planePatchInter
  72 (
  73     planePatch
  74 );
  75
  76
  77 PrimitivePatchInterpolation
  78 <
  79     PrimitivePatch<face, List, const pointField&>
  80 > toPatchInter
  81 (
  82     patches[targetPatchNumber].patch()
  83 );