35 template<
class>
class FaceList,
56 const labelList& faceNbs = edgeFaces[edgeI];
58 if (faceNbs.
size() > 2)
62 const edge&
e = edges[edgeI];
67 e2 /=
mag(e2) + VSMALL;
70 const Face&
f = localFaces[faceNbs[0]];
73 label fp1 = f.fcIndex(fp0);
74 label vertI = (f[fp1] != e[1] ? f[fp1] : f.
fcIndex(fp1));
78 vector e0 = e2 ^ (localPoints[vertI] - edgePt);
79 e0 /=
mag(e0) + VSMALL;
89 for (label nbI = 1; nbI < faceNbs.size(); nbI++)
92 const Face& f = localFaces[faceNbs[nbI]];
94 label fp1 = f.fcIndex(fp0);
95 label vertI = (f[fp1] != e[1] ? f[fp1] : f.
fcIndex(fp1));
97 vector vec = e2 ^ (localPoints[vertI] - edgePt);
98 vec /=
mag(vec) + VSMALL;
119 sortedEdgeFaces[edgeI] = faceNbs;
123 return sortedEdgeFaces;