Fix speed issue with SimplicialLDLT for complexes: the diagonal is real!

This commit is contained in:
Gael Guennebaud 2019-09-30 16:14:34 +02:00
parent 13ef08e5ac
commit c8b2c603b0

View File

@ -131,7 +131,7 @@ void SimplicialCholeskyBase<Derived>::factorize_preordered(const CholMatrixType&
/* the nonzero entry L(k,i) */ /* the nonzero entry L(k,i) */
Scalar l_ki; Scalar l_ki;
if(DoLDLT) if(DoLDLT)
l_ki = yi / m_diag[i]; l_ki = yi / numext::real(m_diag[i]);
else else
yi = l_ki = yi / Lx[Lp[i]]; yi = l_ki = yi / Lx[Lp[i]];