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);