38 vector cEst = vector::zero;
43 cEst += points[pointLabels[i]];
47 cEst /= scalar(pointLabels.
size());
53 vector sumVc = vector::zero;
55 const faceList cellFaces = faces(pointLabels);
59 const face& curFace = cellFaces[i];
70 WarningIn(
"cellModel::centre(const labelList&, const pointField&)")
71 <<
"zero or negative pyramid volume: " << -pyrVol
84 return sumVc/(sumV + VSMALL);
95 vector cEst = vector::zero;
100 cEst += points[pointLabels[i]];
104 cEst /= scalar(pointLabels.
size());
113 const faceList cellFaces = faces(pointLabels);
117 const face& curFace =cellFaces[i];
128 WarningIn(
"cellModel::mag(const labelList&, const pointField&)")
129 <<
"zero or negative pyramid volume: " << -pyrVol