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

(grafted from c8b2c603b0fe391d24ae0c95a1e20b157bfae33a
)
This commit is contained in:
Gael Guennebaud 2019-09-30 16:14:34 +02:00
parent 52d159c19f
commit 165db26dc0

View File

@ -156,7 +156,7 @@ void SimplicialCholeskyBase<Derived>::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]];