stabilise expansion calculation for if no point on arc
authormattijs <mattijs>
Wed, 22 Jul 2009 10:21:45 +0000 (22 11:21 +0100)
committermattijs <mattijs>
Wed, 22 Jul 2009 10:21:45 +0000 (22 11:21 +0100)
applications/utilities/mesh/generation/blockMesh/setEdge.C

index 44e5ca5..3b5c9ea 100644 (file)
@@ -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();