115 cellCtrs[own[facei]],
145 cellCtrs[own[facei]],
192 const point& fc = fCtrs[facei];
193 const vector& fa = fAreas[facei];
195 scalar dOwn =
mag(fa & (fc-cellCtrs[own[facei]]));
196 scalar dNei =
mag(fa & (cellCtrs[nei[facei]]-fc));
198 weight[facei] =
min(dNei,dOwn)/(dNei+dOwn+vSmall);
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);
248 scalar volOwn = vol[own[facei]];
249 scalar volNei = vol[nei[facei]];
251 ratio[facei] =
min(volOwn,volNei)/(
max(volOwn, volNei)+vSmall);
270 scalar volOwn = vol[own[facei]];
271 scalar volNei = neiVol[bFacei];
273 ratio[facei] =
min(volOwn,volNei)/(
max(volOwn, volNei)+vSmall);
#define forAll(list, i)
Loop across all elements in list.
Mesh consisting of general polyhedral cells.
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.
virtual bool coupled() const
Return true if this patch is geometrically coupled (i.e. faces and.
label start() const
Return start label of this patch in the polyMesh face list.
label nInternalFaces() const
A class for managing temporary objects.
T & ref() const
Return non-const reference or generate a fatal error.
intWM_LABEL_SIZE_t label
A label is an int32_t or int64_t as specified by the pre-processor macro WM_LABEL_SIZE.
Field< scalar > scalarField
Specialisation of Field<T> for scalar.
layerAndWeight min(const layerAndWeight &a, const layerAndWeight &b)
dimensioned< scalar > mag(const dimensioned< Type > &)
layerAndWeight max(const layerAndWeight &a, const layerAndWeight &b)
dimensionedTensor skew(const dimensionedTensor &dt)