68 cellCentresPtr_(NULL),
69 faceCentresPtr_(NULL),
70 cellVolumesPtr_(NULL),
80 const label nInternalFaces,
88 nInternalFaces_(nInternalFaces),
110 cellCentresPtr_(NULL),
111 faceCentresPtr_(NULL),
112 cellVolumesPtr_(NULL),
129 label& nInternalPoints,
132 const label nInternalFaces,
146 label nBoundaryPoints = 0;
147 for (label faceI = nInternalFaces; faceI < faces.
size(); faceI++)
149 const face&
f = faces[faceI];
153 label pointI = f[fp];
155 if (oldToNew[pointI] == -1)
157 oldToNew[pointI] = nBoundaryPoints++;
164 nInternalPoints = nPoints - nBoundaryPoints;
169 if (oldToNew[pointI] != -1)
171 oldToNew[pointI] += nInternalPoints;
179 label internalPointI = 0;
183 for (label faceI = 0; faceI < nInternalFaces; faceI++)
185 const face&
f = faces[faceI];
189 label pointI = f[fp];
191 if (oldToNew[pointI] == -1)
193 if (pointI >= nInternalPoints)
197 oldToNew[pointI] = internalPointI++;
209 const label nInternalFaces,
218 nInternal0Edges_ = -1;
219 nInternal1Edges_ = -1;
220 nInternalEdges_ = -1;
222 nInternalFaces_ = nInternalFaces;
227 label nInternalPoints;
230 bool isOrdered = calcPointOrder
241 nInternalPoints_ = nInternalPoints;
245 nInternalPoints_ = -1;
250 Pout<<
"primitiveMesh::reset : mesh reset to"
251 <<
" nInternalPoints:" << nInternalPoints_
252 <<
" nPoints:" << nPoints_
253 <<
" nEdges:" << nEdges_
254 <<
" nInternalFaces:" << nInternalFaces_
255 <<
" nFaces:" << nFaces_
256 <<
" nCells:" << nCells_
265 const label nInternalFaces,
286 const label nInternalFaces,
289 const Xfer<cellList>& clst
314 "primitiveMesh::movePoints(const pointField& newPoints, "
315 "const pointField& oldPoints)"
316 ) <<
"Cannot move points: size of given point list smaller "
317 <<
"than the number of active points"
329 sweptVols[faceI] = f[faceI].sweptVol(oldPoints, newPoints);
346 return *cellShapesPtr_;