78 int main(
int argc,
char *argv[])
88 runTime.functionObjects().off();
103 <<
"Cannot find patch " << patchName <<
endl
116 forAll(meshPoints, pointI)
118 label meshPointI = meshPoints[pointI];
124 cutCells.insert(pCells[pCellI]);
128 Info<<
"Selected " << cutCells.size()
129 <<
" cells connected to patch " << pp.
name() <<
endl <<
endl;
141 Info<<
"Subsetting cells to cut based on cellSet" << setName << endl
157 cutCells.erase(iter.key());
159 Info<<
"Removed from cells to cut all the ones not in set " << setName
167 forAll(meshPoints, pointI)
169 label meshPointI = meshPoints[pointI];
171 vertOnPatch[meshPointI] =
true;
181 forAll(meshPoints, pointI)
183 label meshPointI = meshPoints[pointI];
189 label edgeI = pEdges[pEdgeI];
195 if (!vertOnPatch[otherPointI])
197 allCutEdges.append(edgeI);
199 if (e.
start() == meshPointI)
201 allCutEdgeWeights.append(weight);
205 allCutEdgeWeights.append(1 - weight);
211 allCutEdges.shrink();
212 allCutEdgeWeights.shrink();
214 Info<<
"Cutting:" << endl
215 <<
" cells:" << cutCells.size() << endl
216 <<
" edges:" << allCutEdges.size() << endl
221 cutEdgeWeights.
transfer(allCutEdgeWeights);
222 allCutEdgeWeights.clear();
241 cutter.setRefinement(cuts, meshMod);
253 if (morphMap().hasMotionPoints())
259 cutter.updateMesh(morphMap());
267 Info <<
"Writing refined morphMesh to time " << runTime.timeName() <<
endl;