96 label edgeI = pEdges[pEdgeI];
108 "findEdge(const primitiveMesh&, const label, const label)"
109 ) <<
"Cannot find edge between mesh points " << v0 <<
" and " << v1
123 FatalErrorIn(
"checkPatch(const polyBoundaryMesh&, const word&)")
124 <<
"Cannot find patch " << name <<
endl
125 <<
"It should be present but of zero size" <<
endl
126 <<
"Valid patches are " << bMesh.
names()
130 if (bMesh[patchI].size())
132 FatalErrorIn(
"checkPatch(const polyBoundaryMesh&, const word&)")
133 <<
"Patch " << name <<
" is present but non-zero size"
141 int main(
int argc,
char *argv[])
152 runTime.functionObjects().off();
154 const word oldInstance = mesh.pointsInstance();
162 faceSet facesSet(mesh, setName);
164 Info<<
"Read " << facesSet.size() <<
" faces to split" <<
endl <<
endl;
176 <<
"Face " << faces[i] <<
" in faceSet " << setName
177 <<
" is not an internal face."
184 checkPatch(mesh.boundaryMesh(), masterPatch);
185 checkPatch(mesh.boundaryMesh(), slavePatch);
204 const labelList& meshPoints = fPatch.meshPoints();
212 forAll(allEdgeFaces, patchEdgeI)
214 if (allEdgeFaces[patchEdgeI].size() == 1)
216 const edge& e = fPatch.edges()[patchEdgeI];
222 meshPoints[e.
start()],
226 fenceEdges.insert(edgeI);
231 label startFaceI = faces[0];
247 zoneFlip[i] = !regionInfo.sideOwner().found(faces[i]);
266 Info <<
"Adding point and face zones" <<
endl;
267 mesh.addZones(pz, fz, cz);
289 Info<<
nl <<
"Constructed topologyModifier:" <<
endl;
290 splitter[0].writeDict(
Info);
301 mesh.setInstance(oldInstance);
304 Info<<
"Writing mesh to " << runTime.timeName() <<
endl;
308 <<
"Failed writing polyMesh."