Fix rank condition check, using default MATLAB tolerance (thanks to Jesper Linde and Malin Adolfson)
parent
bc8f546cbb
commit
61d6e49dfd
|
@ -53,7 +53,7 @@ end
|
||||||
JJ1 = JJ(:,ind1);
|
JJ1 = JJ(:,ind1);
|
||||||
[eu,ee2,ee1] = svd( JJ1, 0 );
|
[eu,ee2,ee1] = svd( JJ1, 0 );
|
||||||
condJ= cond(JJ1);
|
condJ= cond(JJ1);
|
||||||
rankJ = rank(JJ./norm(JJ),1.e-10);
|
rankJ = rank(JJ);
|
||||||
rankJJ = rankJ;
|
rankJJ = rankJ;
|
||||||
% if hess_flag==0,
|
% if hess_flag==0,
|
||||||
% rankJJ = rank(JJ'*JJ);
|
% rankJJ = rank(JJ'*JJ);
|
||||||
|
|
Loading…
Reference in New Issue