zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

ENH: if the face area is less than ROOTVSMALL then return vector::zero

Browse Source
This commit is contained in:
laurence
2012-12-11 17:05:03 +00:00
parent 3ad25a12d4
commit 3d49ababd4
1 changed files with 3 additions and 4 deletions
Download Patch File Download Diff File Expand all files Collapse all files

7
src/OpenFOAM/meshes/primitiveMesh/primitiveMeshFaceCentresAndAreas.C
View File

@ -118,13 +118,12 @@ void Foam::primitiveMesh::makeFaceCentresAndAreas
sumAc += a*c; 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) if (sumA < ROOTVSMALL)
{ {
// Sum of area too small. No chance of reliably calculating
// centroid so fallback to average.
fCtrs[facei] = fCentre; fCtrs[facei] = fCentre;
fAreas[facei] = 0.5*sumN; fAreas[facei] = vector::zero;
} }
else else
{ {