X-Git-Url: http://git.salome-platform.org/gitweb/?p=modules%2Fsmesh.git;a=blobdiff_plain;f=src%2FStdMeshers%2FStdMeshers_ViscousLayers.cxx;h=6a6d871c43d65d25128c00a3f466aa903662e998;hp=5fb3736b05536d083c82379ef827d62b571e757f;hb=e2a638a02647937174d628352feeef9640c14ec0;hpb=423f2b873169ff0745e554b17c97d40de844b0fa diff --git a/src/StdMeshers/StdMeshers_ViscousLayers.cxx b/src/StdMeshers/StdMeshers_ViscousLayers.cxx index 5fb3736b0..6a6d871c4 100644 --- a/src/StdMeshers/StdMeshers_ViscousLayers.cxx +++ b/src/StdMeshers/StdMeshers_ViscousLayers.cxx @@ -801,10 +801,10 @@ namespace double u1 = intervals( i ); double u2 = intervals( i+1 ); curve.D2( 0.5*( u1+u2 ), p, drv1, drv2 ); - double cross = drv2 ^ drv1; + double cross = drv2 * drv1; //drv2 ^ drv1; if ( E.Orientation() == TopAbs_REVERSED ) cross = -cross; - isConvex = ( cross < 1e-9 ); + isConvex = ( cross > -1e-9 ); } // check if concavity is strong enough to care about it //const double maxAngle = 5 * Standard_PI180; @@ -865,9 +865,11 @@ namespace const char* fname = "/tmp/viscous.py"; cout << "execfile('"< 0 ) return error(SMESH_Comment("Can't shrink 2D mesh on face ") << f2sd->first ); - if ( shriStep > 20 ) + if ( shriStep > 200 ) return error(SMESH_Comment("Infinite loop at shrinking 2D mesh on face ") << f2sd->first ); }