51 Foam::rawTopoChangerFvMesh::rawTopoChangerFvMesh(
const IOobject& io)
68 Info<<
"rawTopoChangerFvMesh : Checking for topology changes..."
78 bool hasChanged = topoChangeMap.
valid();
82 Info<<
"rawTopoChangerFvMesh : Done topology changes..."
93 bitSet mappedFace(nFaces());
95 const label nOldInternal = topoChangeMap().oldPatchStarts()[0];
98 for (label facei = 0; facei < nInternalFaces(); facei++)
102 mappedFace.
set(facei);
105 for (label facei = nInternalFaces(); facei < nFaces(); facei++)
109 mappedFace.
set(facei);
113 const List<objectMap>& fromFaces = topoChangeMap().facesFromFacesMap();
117 mappedFace.
set(fromFaces[i].index());
120 const List<objectMap>& fromEdges = topoChangeMap().facesFromEdgesMap();
124 mappedFace.
set(fromEdges[i].index());
131 mappedFace.
set(fromPts[i].index());
135 Info<<
"rawTopoChangerFvMesh : zeroing unmapped boundary values."
137 zeroUnmappedValues<scalar, fvPatchField, volMesh>(mappedFace);
138 zeroUnmappedValues<vector, fvPatchField, volMesh>(mappedFace);
139 zeroUnmappedValues<sphericalTensor, fvPatchField, volMesh>(mappedFace);
140 zeroUnmappedValues<symmTensor, fvPatchField, volMesh>(mappedFace);
141 zeroUnmappedValues<tensor, fvPatchField, volMesh>(mappedFace);
144 Info<<
"rawTopoChangerFvMesh :"
145 <<
" recreating phi for unmapped boundary values." <<
endl;
158 if (topoChangeMap().hasMotionPoints())
160 pointField newPoints = topoChangeMap().preMotionPoints();
163 Info<<
"rawTopoChangerFvMesh :"
164 <<
" calling modifyMotionPoints." <<
endl;
165 topoChanger_.modifyMotionPoints(newPoints);
168 Info<<
"rawTopoChangerFvMesh :"
169 <<
" calling movePoints." <<
endl;
171 movePoints(newPoints);