From 1e2a6b5eada5ec9c3aeafe2fe0804d12204d4dfa Mon Sep 17 00:00:00 2001
From: mattijs <mattijs>
Date: Thu, 18 Apr 2013 10:07:52 +0100
Subject: [PATCH] BUG: polyMeshTools: volume ratio caclulation
---
src/OpenFOAM/meshes/polyMesh/polyMeshCheck/polyMeshTools.C | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/src/OpenFOAM/meshes/polyMesh/polyMeshCheck/polyMeshTools.C b/src/OpenFOAM/meshes/polyMesh/polyMeshCheck/polyMeshTools.C
index 1884485a697..59947b0fbe2 100644
--- a/src/OpenFOAM/meshes/polyMesh/polyMeshCheck/polyMeshTools.C
+++ b/src/OpenFOAM/meshes/polyMesh/polyMeshCheck/polyMeshTools.C
@@ -264,7 +264,7 @@ Foam::tmp<Foam::scalarField> Foam::polyMeshTools::volRatio
scalar volOwn = vol[own[faceI]];
scalar volNei = vol[nei[faceI]];
- ratio[faceI] = min(volOwn,volNei)/(volOwn+volNei+VSMALL);
+ ratio[faceI] = min(volOwn,volNei)/(max(volOwn, volNei)+VSMALL);
}
@@ -286,7 +286,7 @@ Foam::tmp<Foam::scalarField> Foam::polyMeshTools::volRatio
scalar volOwn = vol[own[faceI]];
scalar volNei = neiVol[bFaceI];
- ratio[faceI] = min(volOwn,volNei)/(volOwn+volNei+VSMALL);
+ ratio[faceI] = min(volOwn,volNei)/(max(volOwn, volNei)+VSMALL);
}
}
}
--
GitLab