35 bool Foam::directionInfo::equal(
const edge&
e,
const label v0,
const label
v1)
38 (e.start() == v0 && e.end() ==
v1)
39 || (e.start() == v1 && e.end() == v0);
45 const primitiveMesh&
mesh,
49 const edge& e = mesh.edges()[edgeI];
53 * (mesh.points()[e.start()] + mesh.points()[e.end()]);
60 const primitiveMesh& mesh,
66 forAll(edgeLabels, edgeLabelI)
68 label edgeI = edgeLabels[edgeLabelI];
70 const edge& e = mesh.edges()[edgeI];
79 <<
"Cannot find an edge among " << edgeLabels <<
endl
80 <<
"that uses vertices " << v0
88 Foam::label Foam::directionInfo::lowest
96 label a1 = (a + 1) % size;
104 label b1 = (b + 1) % size;
109 <<
"Problem : a:" << a <<
" b:" << b <<
" size:" << size
128 if ((edgeI < 0) || (edgeI >= mesh.
nEdges()))
131 <<
"Illegal edge label:" << edgeI
132 <<
" when projecting cut edge from cell " << cellI
133 <<
" to face " << faceI
153 return lowest(f.
size(), fpA, fpB);
192 if ((fpA != -1) && (fpB != -1))
194 return lowest(f.
size(), fpA, fpB);
210 if ((fpA != -1) && (fpB != -1))
212 return lowest(f.
size(), fpA, fpB);
216 <<
"Found connected faces " << mesh.
faces()[f0I] <<
" and "
217 << mesh.
faces()[f1I] <<
" sharing edge " << edgeI <<
endl
218 <<
"But none seems to be connected to face " << faceI
234 const label thisCellI,
235 const label neighbourFaceI,
250 if (neighbourInfo.
index() == -2)
255 else if (neighbourInfo.
index() == -1)
261 label edgeI = mesh.
faceEdges()[neighbourFaceI][0];
290 label v0 = f[neighbourInfo.
index()];
291 label v1 = f[(neighbourInfo.
index() + 1) % f.
size()];
303 n_ = neighbourInfo.
n();
313 const label thisFaceI,
314 const label neighbourCellI,
330 if (neighbourInfo.
index() >= 0)
340 neighbourInfo.
index()
350 n_ = neighbourInfo.
n();
372 index_ = neighbourInfo.
index();
374 n_ = neighbourInfo.
n();
389 return os << wDist.index_ << wDist.n_;
395 return is >> wDist.index_ >> wDist.n_;