Commit fd076755 authored by Kirill Terekhov's avatar Kirill Terekhov

Fix svd breakdown in bcgsl

parent e4076bbf
......@@ -206,6 +206,7 @@ namespace INMOST
flag = 0;
break;
}
if( l == 0 ) break;
if (fabs(w[nm]) + anorm == anorm)
break;
}
......@@ -698,8 +699,8 @@ namespace INMOST
//Recompute rho1
#else
reason = "r_tilde perturbation algorithm is disabled";
halt = true;
break;
halt = true;
break;
#endif
}
#if defined(USE_OMP)
......
Markdown is supported
0%
or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment