64 cEst_d += fCtrs_d[facei];
67 cEst /= fAreas.
size();
68 cEst_d /= fAreas.
size();
74 mag(fAreas[facei] & (fCtrs[facei] - cEst));
77 (fAreas[facei] & (fCtrs[facei] - cEst))
79 ((fCtrs[facei] - cEst) & fAreas_d[facei])
81 + (fAreas[facei] & (fCtrs_d[facei] - cEst_d))
85 vector pc = (3.0/4.0)*fCtrs[facei] + (1.0/4.0)*cEst;
86 pT pc_d = (3.0/4.0)*fCtrs_d[facei] + (1.0/4.0)*cEst_d;
89 cellCtrs += pyr3Vol*pc;
92 cellCtrs_d += (pc*pyr3Vol_d + pyr3Vol*pc_d);
96 cellVols_d += pyr3Vol_d;
102 cellVols_d *= (1.0/3.0);
void size(const label n)
Older name for setAddressableSize.
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
A traits class, which is primarily used for primitives and vector-space.
#define forAll(list, i)
Loop across all elements in list.
A Vector of values with scalar precision, where scalar is float/double depending on the compilation f...
const scalarField & cellVols
pT makeCellCentres_d(const vectorField &fAreas, const vectorField &fCtrs, const Field< pT > &fAreas_d, const Field< pT > &fCtrs_d)
Compute cell center variation wrt given face movement or derivative.
typeOfRank< typename pTraits< arg1 >::cmptType, direction(pTraits< arg1 >::rank)+direction(pTraits< arg2 >::rank) - 2 >::type type
static constexpr const zero Zero
Global zero (0)