From 165db26dc07d9219524e10b69aa937536172e21c Mon Sep 17 00:00:00 2001 From: Gael Guennebaud Date: Mon, 30 Sep 2019 16:14:34 +0200 Subject: [PATCH] Fix speed issue with SimplicialLDLT for complexes: the diagonal is real! (grafted from c8b2c603b0fe391d24ae0c95a1e20b157bfae33a ) --- Eigen/src/SparseCholesky/SimplicialCholesky_impl.h | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/Eigen/src/SparseCholesky/SimplicialCholesky_impl.h b/Eigen/src/SparseCholesky/SimplicialCholesky_impl.h index 31e06995b..7b6183d08 100644 --- a/Eigen/src/SparseCholesky/SimplicialCholesky_impl.h +++ b/Eigen/src/SparseCholesky/SimplicialCholesky_impl.h @@ -156,7 +156,7 @@ void SimplicialCholeskyBase::factorize_preordered(const CholMatrixType& /* the nonzero entry L(k,i) */ Scalar l_ki; if(DoLDLT) - l_ki = yi / m_diag[i]; + l_ki = yi / numext::real(m_diag[i]); else yi = l_ki = yi / Lx[Lp[i]];