devtools/eigen
1
0
Fork 0
mirror of https://gitlab.com/libeigen/eigen.git synced 2026-04-10 11:34:33 +08:00
Code Issues Packages Projects Releases Wiki Activity

don't try to finish early

Browse Source
This commit is contained in:
Benoit Jacob
2009-10-15 16:09:17 -04:00
parent 924b55e9a9
commit 41e942d3fb
1 changed files with 2 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

4
Eigen/src/LU/LU.h
View File

@@ -361,8 +361,8 @@ LU<MatrixType>& LU<MatrixType>::compute(const MatrixType& matrix)
col_of_biggest_in_corner += k; col_of_biggest_in_corner += k;
if(k==0) biggest = biggest_in_corner; if(k==0) biggest = biggest_in_corner;
// if the corner is negligible, then we have less than full rank, and we can finish early // if the corner is exactly zero, terminate to avoid generating NaN values
if(ei_isMuchSmallerThan(biggest_in_corner, biggest, m_precision)) if(biggest_in_corner == RealScalar(0))
{ {
m_rank = k; m_rank = k;
for(int i = k; i < size; i++) for(int i = k; i < size; i++)