From c8b2c603b0fe391d24ae0c95a1e20b157bfae33a 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! --- 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 7275db2cc..72e1740c1 100644 --- a/Eigen/src/SparseCholesky/SimplicialCholesky_impl.h +++ b/Eigen/src/SparseCholesky/SimplicialCholesky_impl.h @@ -131,7 +131,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]];