diff --git a/src/OpenFOAM/meshes/primitiveMesh/primitiveMeshFaceCentresAndAreas.C b/src/OpenFOAM/meshes/primitiveMesh/primitiveMeshFaceCentresAndAreas.C index 8de9076603..589ff7a98b 100644 --- a/src/OpenFOAM/meshes/primitiveMesh/primitiveMeshFaceCentresAndAreas.C +++ b/src/OpenFOAM/meshes/primitiveMesh/primitiveMeshFaceCentresAndAreas.C @@ -118,13 +118,12 @@ void Foam::primitiveMesh::makeFaceCentresAndAreas sumAc += a*c; } - + // This is to deal with zero-area faces. Mark very small faces + // to be detected in e.g., processorPolyPatch. if (sumA < ROOTVSMALL) { - // Sum of area too small. No chance of reliably calculating - // centroid so fallback to average. fCtrs[facei] = fCentre; - fAreas[facei] = 0.5*sumN; + fAreas[facei] = vector::zero; } else {