diff options
author | Christoph Hertzberg <chtz@informatik.uni-bremen.de> | 2018-12-03 16:18:15 +0100 |
---|---|---|
committer | Christoph Hertzberg <chtz@informatik.uni-bremen.de> | 2018-12-03 16:18:15 +0100 |
commit | 919414b9fe2ad7fdcb0f2b2cbdf6b5322d0f2034 (patch) | |
tree | b5269508221dd5d95b38386423619254d780d937 /Eigen/src/Cholesky/LDLT.h | |
parent | 0ea7ae72130cac7334823ec442f0a8a6772c9ab8 (diff) |
bug #785: Make Cholesky decomposition work for empty matrices
Diffstat (limited to 'Eigen/src/Cholesky/LDLT.h')
-rw-r--r-- | Eigen/src/Cholesky/LDLT.h | 3 |
1 files changed, 2 insertions, 1 deletions
diff --git a/Eigen/src/Cholesky/LDLT.h b/Eigen/src/Cholesky/LDLT.h index 2dfeac333..6831eab3d 100644 --- a/Eigen/src/Cholesky/LDLT.h +++ b/Eigen/src/Cholesky/LDLT.h @@ -304,7 +304,8 @@ template<> struct ldlt_inplace<Lower> if (size <= 1) { transpositions.setIdentity(); - if (numext::real(mat.coeff(0,0)) > static_cast<RealScalar>(0) ) sign = PositiveSemiDef; + if(size==0) sign = ZeroSign; + else if (numext::real(mat.coeff(0,0)) > static_cast<RealScalar>(0) ) sign = PositiveSemiDef; else if (numext::real(mat.coeff(0,0)) < static_cast<RealScalar>(0)) sign = NegativeSemiDef; else sign = ZeroSign; return true; |