From 7054a8be58e0a9ff38c20d2c26cf5a56000b08d3 Mon Sep 17 00:00:00 2001
From: mattijs <mattijs>
Date: Mon, 27 Jan 2014 12:54:48 +0000
Subject: [PATCH] ENH: primitiveMeshTools: handle zero/neg volume cells

---
 .../primitiveMesh/primitiveMeshCheck/primitiveMeshTools.C   | 6 ++++--
 1 file changed, 4 insertions(+), 2 deletions(-)

diff --git a/src/OpenFOAM/meshes/primitiveMesh/primitiveMeshCheck/primitiveMeshTools.C b/src/OpenFOAM/meshes/primitiveMesh/primitiveMeshCheck/primitiveMeshTools.C
index 241f4eef8e..fa83d3ae9d 100644
--- a/src/OpenFOAM/meshes/primitiveMesh/primitiveMeshCheck/primitiveMeshTools.C
+++ b/src/OpenFOAM/meshes/primitiveMesh/primitiveMeshCheck/primitiveMeshTools.C
@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2012-2013 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2012-2014 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -315,10 +315,12 @@ void Foam::primitiveMeshTools::cellClosedness
         scalar aspectRatio = maxCmpt/(minCmpt + VSMALL);
         if (nDims == 3)
         {
+            scalar v = max(VSMALL, vols[cellI]);
+
             aspectRatio = max
             (
                 aspectRatio,
-                1.0/6.0*cmptSum(sumMagClosed[cellI])/pow(vols[cellI], 2.0/3.0)
+                1.0/6.0*cmptSum(sumMagClosed[cellI])/pow(v, 2.0/3.0)
             );
         }