Go to the documentation of this file.
44 for (
const label facei : cFaces)
46 nVerts += meshFaces[facei].
size();
52 const labelList& firstFace = meshFaces[cFaces[0]];
54 std::copy(firstFace.cbegin(), firstFace.cend(),
pointLabels.begin());
57 nVerts = firstFace.size();
63 for (label facei = 1; facei < cFaces.size(); ++facei)
65 for (
const label curPoint : meshFaces[cFaces[facei]])
67 bool pointFound =
false;
69 for (label checki = 0; checki < nVerts; ++checki)
116 for (
const label facei : cFaces)
118 nEdges += meshFaces[facei].nEdges();
127 for (
const edge& curEdge : meshFaces[cFaces[facei]].edges())
129 bool edgeFound =
false;
131 for (label checki = 0; checki < nEdges; ++checki)
133 if (curEdge == allEdges[checki])
142 allEdges[nEdges] = curEdge;
182 for (
const label facei : cFaces)
184 const scalar magArea = meshFaces[facei].mag(meshPoints);
185 ctr += meshFaces[facei].centre(meshPoints)*magArea;
189 ctr /= sumArea + VSMALL;
197 for (
const label facei : cFaces)
199 const face&
f = meshFaces[facei];
214 return sumVc/(sumV + VSMALL);
237 for (
const label facei : cFaces)
239 ctr += meshFaces[facei].centre(meshPoints);
241 ctr /= cFaces.
size();
246 for (
const label facei : cFaces)
248 const face&
f = meshFaces[facei];
262 if (a.size() !=
b.size())
269 for (
const label curLabel :
b)
275 if (a[ai] == curLabel)
void resize(const label len)
Adjust allocated size of list.
static constexpr const zero Zero
Global zero (0)
pointField points(const faceUList &meshFaces, const UList< point > &meshPoints) const
Return the cell vertices given the list of faces and mesh points.
static const char *const typeName
An edge is a list of two point labels. The functionality it provides supports the discretisation on a...
#define forAll(list, i)
Loop across all elements in list.
scalar mag(const UList< point > &meshPoints, const faceUList &meshFaces) const
Returns cell volume.
labelList labels(const faceUList &meshFaces) const
Return unordered list of cell vertices given the list of faces.
const dimensionedScalar b
Wien displacement law constant: default SI units: [m.K].
tmp< faMatrix< Type > > operator==(const faMatrix< Type > &, const faMatrix< Type > &)
edgeList edges(const faceUList &meshFaces) const
Return cell edges.
point centre(const UList< point > &meshPoints, const faceUList &meshFaces) const
Returns cell centre.
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
A face is a list of labels corresponding to mesh vertices.
void size(const label n)
Older name for setAddressableSize.
static constexpr direction size() noexcept
The number of elements in the VectorSpace = Ncmpts.
pyramid< point, const point &, const face & > pyramidPointFaceRef
A pyramid using referred points and faces.
A cell is defined as a list of faces with extra functionality.
labelList pointLabels(nPoints, -1)