34 const Foam::label Foam::tetWedgeMatcher::vertPerCell = 5;
35 const Foam::label Foam::tetWedgeMatcher::facePerCell = 4;
36 const Foam::label Foam::tetWedgeMatcher::maxVertPerFace = 4;
71 if (!faceSizeMatch(faces, myFaces))
83 label numVert = calcLocalFaces(faces, myFaces);
85 if (numVert != vertPerCell)
91 calcEdgeAddressing(numVert);
97 vertLabels_.setSize(vertPerCell);
98 faceLabels_.setSize(facePerCell);
108 if (faceSize_[faceI] == 3)
115 const face& face0 = localFaces_[face0I];
118 for(label face0vert0 = 0; face0vert0 < faceSize_[face0I]; face0vert0++)
125 vertLabels_[0] = pointMap_[face0[face0vert0]];
126 faceLabels_[0] = faceMap_[face0I];
134 !(owner[faceMap_[face0I]] == cellI)
136 vertLabels_[1] = pointMap_[face0[face0vert1]];
148 if (faceSize_[face1I] != 3)
152 faceLabels_[1] = faceMap_[face1I];
163 !(owner[faceMap_[face0I]] == cellI)
165 vertLabels_[2] = pointMap_[face0[face0vert2]];
176 faceLabels_[3] = faceMap_[face3I];
187 faceLabels_[2] = faceMap_[face2I];
190 label face3vert2 = pointFaceIndex_[face0[face0vert2]][face3I];
198 (owner[faceMap_[face3I]] == cellI)
201 const face& face3 = localFaces_[face3I];
203 vertLabels_[4] = pointMap_[face3[face3vert4]];
211 (owner[faceMap_[face3I]] == cellI)
213 vertLabels_[3] = pointMap_[face3[face3vert3]];
235 if (myFaces.
size() != 4)
245 label size = faces[myFaces[myFaceI]].
size();
260 if ((nTris == 2) && (nQuads == 2))
293 makeIdentity(faces.
size())
317 shape =
cellShape(model(), vertLabels());