Skip to content
Snippets Groups Projects
Commit 2ea070db authored by Christophe Geuzaine's avatar Christophe Geuzaine
Browse files

lowered tolerance in sys3x3_with_tol -> check benchmarks/bug/bug_ruled.geo

parent ef585516
No related branches found
No related tags found
No related merge requests found
// $Id: Numeric.cpp,v 1.16 2001-07-26 21:26:34 geuzaine Exp $
// $Id: Numeric.cpp,v 1.17 2001-08-01 21:18:21 geuzaine Exp $
#include "Gmsh.h"
#include "Const.h"
......@@ -130,7 +130,7 @@ int sys3x3_with_tol (double mat[3][3], double b[3], double res[3], double *det){
DSQR (mat[2][0]) + DSQR (mat[2][1]) + DSQR (mat[2][2]) ;
// TOLERANCE ! WARNING WARNING
if (norm == 0.0 || fabs (*det) / norm < 1.e-07){
if (norm == 0.0 || fabs (*det) / norm < 1.e-12){
Msg(DEBUG, "Assuming 3x3 matrix is singular (det/norm == %g)", fabs(*det)/norm);
res[0] = res[1] = res[2] = 0.0 ;
return 0;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment