Cast to diagonalSize to RealScalar instead Scalar.

This commit is contained in:
Rasmus Munk Larsen 2018-08-09 14:46:17 -07:00
parent 8603d80029
commit bfc5091dd5

View File

@ -320,7 +320,7 @@ template<typename _MatrixType> class FullPivLU
return m_usePrescribedThreshold ? m_prescribedThreshold
// this formula comes from experimenting (see "LU precision tuning" thread on the list)
// and turns out to be identical to Higham's formula used already in LDLt.
: NumTraits<Scalar>::epsilon() * Scalar(m_lu.diagonalSize());
: NumTraits<Scalar>::epsilon() * RealScalar(m_lu.diagonalSize());
}
/** \returns the rank of the matrix of which *this is the LU decomposition.