38 bool Foam::layerAdditionRemoval::setLayerPairing()
const
58 mesh.faceZones()[faceZoneID_.
index()].masterCells();
63 mesh.faceZones()[faceZoneID_.
index()].flipMap();
65 const faceList& faces = mesh.faces();
70 mesh.faceZones()[faceZoneID_.
index()]().localFaces();
73 mesh.faceZones()[faceZoneID_.
index()]().meshPoints();
77 if (pointsPairingPtr_ || facesPairingPtr_)
81 "void Foam::layerAdditionRemoval::setLayerPairing() const"
82 ) <<
"Problem with layer pairing data"
86 pointsPairingPtr_ =
new labelList(meshPoints.size(), -1);
89 facesPairingPtr_ =
new labelList(mf.size(), -1);
95 label nPointErrors = 0;
96 label nFaceErrors = 0;
101 face curLocalFace = mlf[faceI];
106 curLocalFace = curLocalFace.reverseFace();
110 oppositeFace lidFace =
111 cells[mc[faceI]].opposingFace(mf[faceI], faces);
113 if (!lidFace.found())
128 ftc[faceI] = lidFace.oppositeIndex();
131 forAll (curLocalFace, pointI)
133 const label clp = curLocalFace[pointI];
138 ptc[clp] = lidFace[pointI];
143 if (ptc[clp] != lidFace[pointI])
159 reduce(nPointErrors, sumOp<label>());
160 reduce(nFaceErrors, sumOp<label>());
162 if (nPointErrors > 0 || nFaceErrors > 0)
176 const Foam::labelList& Foam::layerAdditionRemoval::pointsPairing()
const
178 if (!pointsPairingPtr_)
182 "const labelList& layerAdditionRemoval::pointsPairing() const"
183 ) <<
"Problem with layer pairing data for object " <<
name()
187 return *pointsPairingPtr_;
190 const Foam::labelList& Foam::layerAdditionRemoval::facesPairing()
const
192 if (!facesPairingPtr_)
196 "const labelList& layerAdditionRemoval::facesPairing() const"
197 ) <<
"Problem with layer pairing data for object " <<
name()
201 return *facesPairingPtr_;
214 Pout<<
"void layerAdditionRemoval::modifyMotionPoints("
215 <<
"pointField& motionPoints) const for object "
221 Pout <<
"No motion point adjustment" <<
endl;