48 auto& ortho = tortho.ref();
73 label facei = pp.start() + i;
105 auto&
skew = tskew.ref();
117 cellCtrs[own[facei]],
135 label facei = pp.start() + i;
146 cellCtrs[own[facei]],
155 label facei = pp.start() + i;
188 auto& weight = tweight.ref();
193 const point& fc = fCtrs[facei];
194 const vector& fa = fAreas[facei];
196 scalar dOwn =
mag(fa & (fc-cellCtrs[own[facei]]));
197 scalar dNei =
mag(fa & (cellCtrs[nei[facei]]-fc));
199 weight[facei] =
min(dNei,dOwn)/(dNei+dOwn+VSMALL);
214 label facei = pp.start() + i;
217 const point& fc = fCtrs[facei];
218 const vector& fa = fAreas[facei];
220 scalar dOwn =
mag(fa & (fc-cellCtrs[own[facei]]));
221 scalar dNei =
mag(fa & (neiCc[bFacei]-fc));
223 weight[facei] =
min(dNei,dOwn)/(dNei+dOwn+VSMALL);
243 auto& ratio = tratio.ref();
248 scalar volOwn = vol[own[facei]];
249 scalar volNei = vol[nei[facei]];
251 ratio[facei] =
min(volOwn,volNei)/(
max(volOwn, volNei)+VSMALL);
266 label facei = pp.start() + i;
269 scalar volOwn = vol[own[facei]];
270 scalar volNei = neiVol[bFacei];
272 ratio[facei] =
min(volOwn,volNei)/(
max(volOwn, volNei)+VSMALL);
static autoPtr< Time > New()
Construct (dummy) Time - no functionObjects or libraries.
tmp< scalarField > faceSkewness() const
Return face skewness.
A polyBoundaryMesh is a polyPatch list with additional search methods and registered IO.
Mesh consisting of general polyhedral cells.
virtual const faceList & faces() const
Return raw faces.
readUpdateState
Enumeration defining the state of the mesh after a read update.
virtual const labelList & faceOwner() const
Return face owner.
const polyBoundaryMesh & boundaryMesh() const
Return boundary mesh.
virtual const labelList & faceNeighbour() const
Return face neighbour.
A patch is a list of labels that address the faces in the global face list.
label nInternalFaces() const noexcept
Number of internal faces.
label nFaces() const noexcept
Number of mesh faces.
A class for managing temporary objects.
label max(const labelHashSet &set, label maxValue=labelMin)
Find the max value in labelHashSet, optionally limited by second argument.
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
label min(const labelHashSet &set, label minValue=labelMax)
Find the min value in labelHashSet, optionally limited by second argument.
dimensionedTensor skew(const dimensionedTensor &dt)
#define forAll(list, i)
Loop across all elements in list.