33 void starMesh::createPolyCells()
39 cellPolys_.setSize(cellShapes_.size());
45 cell& curCell = cellPolys_[cellI];
47 curCell.
setSize(cellFaces_[cellI].size());
54 maxFaces += cellFaces_[cellI].
size();
57 Info <<
"Maximum possible number of faces in mesh: " << maxFaces <<
endl;
59 meshFaces_.setSize(maxFaces);
76 const faceList& curFaces = cellFaces_[cellI];
84 label nNeighbours = 0;
90 if (cellPolys_[cellI][faceI] >= 0)
continue;
94 const face& curFace = curFaces[faceI];
103 const labelList& curNeighbours = PointCells[curPoints[pointI]];
106 forAll(curNeighbours, neiI)
108 label curNei = curNeighbours[neiI];
115 const faceList& searchFaces = cellFaces_[curNei];
117 forAll(searchFaces, neiFaceI)
119 if (searchFaces[neiFaceI] == curFace)
125 neiCells[faceI] = curNei;
126 faceOfNeiCell[faceI] = neiFaceI;
141 for (label neiSearch = 0; neiSearch < nNeighbours; neiSearch++)
145 label minNei = cellPolys_.size();
149 if (neiCells[ncI] > -1 && neiCells[ncI] < minNei)
152 minNei = neiCells[ncI];
159 meshFaces_[nInternalFaces_] = curFaces[nextNei];
162 cellPolys_[cellI][nextNei] = nInternalFaces_;
165 cellPolys_[neiCells[nextNei]][faceOfNeiCell[nextNei]] =
169 neiCells[nextNei] = -1;
177 <<
"Error in internal face insertion"