From: mattijs Date: Wed, 22 Jul 2009 10:21:45 +0000 (+0100) Subject: stabilise expansion calculation for if no point on arc X-Git-Url: https://repo.or.cz/w/OpenFOAM-1.5.x.git/commitdiff_plain/de3a9a2f77891b7073498759f03390386b2a2d5a stabilise expansion calculation for if no point on arc --- 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();