69 const cell&
c = mesh_.cells()[cellI_];
71 bool fullyBelow =
true;
72 bool fullyAbove =
true;
74 label nFaceBelowInterface = 0;
77 for (
const label facei :
c)
79 const label faceStatus = cutFace_.calcSubFace(facei, cutValue_);
83 cutFaceCentres_.append(cutFace_.subFaceCentre());
84 cutFaceAreas_.append(cutFace_.subFaceArea());
85 isoFaceEdges_.append(cutFace_.surfacePoints());
89 else if (faceStatus == -1)
91 cutFaceCentres_.append(cutFace_.subFaceCentre());
92 cutFaceAreas_.append(cutFace_.subFaceArea());
94 nFaceBelowInterface++;
102 if (!fullyBelow && !fullyAbove)
118 if (
mag(faceArea_) < ROOTVSMALL)
120 if (nFaceBelowInterface == 0)
124 subCellCentre_ =
Zero;
133 subCellCentre_ = mesh_.C()[cellI_];
134 subCellVolume_ = mesh_.V()[cellI_];
140 cutFaceCentres_.append(faceCentre_);
141 cutFaceAreas_.append(faceArea_);
152 VOF_ = subCellVolume_ / mesh_.V()[cellI_];
157 subCellCentre_ =
Zero;
164 subCellCentre_ = mesh_.C()[cellI_];
165 subCellVolume_ = mesh_.V()[cellI_];
175 if (facePoints_.empty())
178 calcIsoFacePointsFromEdges
195 cutFaceCentres_.clear();
196 cutFaceAreas_.clear();
197 isoFaceEdges_.clear();
201 subCellCentre_ =
Zero;
202 subCellVolume_ = -10;