From de3a9a2f77891b7073498759f03390386b2a2d5a Mon Sep 17 00:00:00 2001
From: mattijs <mattijs>
Date: Wed, 22 Jul 2009 11:21:45 +0100
Subject: [PATCH] stabilise expansion calculation for if no point on arc

---
 applications/utilities/mesh/generation/blockMesh/setEdge.C | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/applications/utilities/mesh/generation/blockMesh/setEdge.C b/applications/utilities/mesh/generation/blockMesh/setEdge.C
index 44e5ca5..3b5c9ea 100644
--- a/applications/utilities/mesh/generation/blockMesh/setEdge.C
+++ b/applications/utilities/mesh/generation/blockMesh/setEdge.C
@@ -100,7 +100,7 @@ void blockDescriptor::setEdge(label edgeI, label start, label end, label dim)
                 scalar gExp = calcGexp(expand_[edgeI], dim);
 
                 // divide the line
-                lineDivide divEdge(curvedEdges_[nCEI], dim, 1.0/gExp);
+                lineDivide divEdge(curvedEdges_[nCEI], dim, 1.0/(gExp+SMALL));
 
                 pointField p = divEdge.points();
                 scalarList d = divEdge.lambdaDivisions();
-- 
2.11.4.GIT