diff options
| author |  Gael Guennebaud <g.gael@free.fr> | 2017-11-16 17:55:24 +0100 |
|---|---|---|
| committer | Gael Guennebaud <g.gael@free.fr> | 2017-11-16 17:55:24 +0100 |
| commit | 672bdc126b0923e6228a024ce62d1f18b05840ea (patch) | |
| tree | cb5788075f7db8d4ecaceea9819a03e0ffa332f8 /Eigen/src/Cholesky | |
| parent | 95ecb2b5d670f80133ffc3fb34215ab7e33288ac (diff) | |
bug #1479: fix failure detection in LDLT
Diffstat (limited to 'Eigen/src/Cholesky')
| -rw-r--r-- | Eigen/src/Cholesky/LDLT.h | 2 |
1 files changed, 2 insertions, 0 deletions
diff --git a/Eigen/src/Cholesky/LDLT.h b/Eigen/src/Cholesky/LDLT.h index 968427b3a..13a8f6d14 100644 --- a/Eigen/src/Cholesky/LDLT.h +++ b/Eigen/src/Cholesky/LDLT.h @@ -375,6 +375,8 @@ template<> struct ldlt_inplace<Lower> if((rs>0) && pivot_is_valid) A21 /= realAkk; + else if(rs>0) + ret = ret && (A21.array()==Scalar(0)).all(); if(found_zero_pivot && pivot_is_valid) ret = false; // factorization failed else if(!pivot_is_valid) found_zero_pivot = true; |