68 if (warn && nWarnings_ < maxNWarnings)
72 <<
"No base point for face " <<
face() <<
", " <<
f
73 <<
", produces a valid tet decomposition." <<
endl;
74 if (nWarnings_ == maxNWarnings)
77 <<
"Suppressing any further warnings." <<
endl;
82 label facePtI = (tetPti_ + faceBasePtI) %
f.
size();
83 label faceOtherPtI =
f.
fcIndex(facePtI);
87 std::swap(facePtI, faceOtherPtI);
90 return triFace(
f[faceBasePtI],
f[facePtI],
f[faceOtherPtI]);
108 if (warn && nWarnings_ < maxNWarnings)
112 <<
"No base point for face " <<
face() <<
", " <<
f
113 <<
", produces a valid tet decomposition." <<
endl;
114 if (nWarnings_ == maxNWarnings)
117 <<
"Suppressing any further warnings." <<
endl;
122 label facePtI = (tetPti_ + faceBasePtI) %
f.
size();
123 label faceOtherPtI =
f.
fcIndex(facePtI);
127 std::swap(facePtI, faceOtherPtI);
130 return triFace(faceBasePtI, facePtI, faceOtherPtI);
173 meshOldPoints[tri[0]],
174 meshOldPoints[tri[1]],
175 meshOldPoints[tri[2]]
189 ((
diff = (a.cell() -
b.cell())) != 0) ?
diff
190 : ((
diff = (a.face() -
b.face())) != 0) ?
diff
191 : (a.tetPt() -
b.tetPt())
204 && a.face() ==
b.face()
205 && a.tetPt() ==
b.tetPt()
void size(const label n)
Older name for setAddressableSize.
label fcIndex(const label i) const noexcept
A cell is defined as a list of faces with extra functionality.
A face is a list of labels corresponding to mesh vertices.
Mesh consisting of general polyhedral cells.
virtual const faceList & faces() const
Return raw faces.
virtual const labelList & faceOwner() const
Return face owner.
const labelIOList & tetBasePtIs() const
Return the tetBasePtIs.
virtual const pointField & oldPoints() const
Return old points (mesh motion)
virtual const pointField & points() const
Return raw points.
const vectorField & cellCentres() const
Storage and named access for the indices of a tet which is part of the decomposition of a cell.
triFace triIs(const polyMesh &mesh, const bool warn=true) const
static label compare(const tetIndices &a, const tetIndices &b) noexcept
triFace faceTriIs(const polyMesh &mesh, const bool warn=true) const
Return the indices corresponding to the tri on the face for.
tetPointRef tet(const polyMesh &mesh) const
Return the geometry corresponding to this tet.
triPointRef faceTri(const polyMesh &mesh) const
constexpr tetIndices() noexcept
Default construct, with invalid labels (-1)
triPointRef oldFaceTri(const polyMesh &mesh) const
A triangular face using a FixedList of labels corresponding to mesh vertices.
A triangle primitive used to calculate face normals and swept volumes.
#define WarningInFunction
Report a warning using Foam::Warning.
bool operator!=(const eddy &a, const eddy &b)
tetrahedron< point, const point & > tetPointRef
A tetrahedron using referred points.
Ostream & endl(Ostream &os)
Add newline and flush stream.
tmp< faMatrix< Type > > operator==(const faMatrix< Type > &, const faMatrix< Type > &)
triangle< point, const point & > triPointRef
A triangle using referred points.
scalar diff(const triad &A, const triad &B)
Return a quantity of the difference between two triads.