From a9eabed42140c8704cfac1a7b02c1612d562b623 Mon Sep 17 00:00:00 2001 From: Benoit Jacob Date: Tue, 2 Feb 2010 07:06:15 -0500 Subject: [PATCH] Patch by 'Wolf' from the issue tracker: Fix bug #90, missing type cast in LU, allow to use LU with MPFR. --- Eigen/src/LU/LU.h | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/Eigen/src/LU/LU.h b/Eigen/src/LU/LU.h index 60b62a8e8..ef1487c6a 100644 --- a/Eigen/src/LU/LU.h +++ b/Eigen/src/LU/LU.h @@ -509,7 +509,7 @@ bool LU::solve( if(!isSurjective()) { // is c is in the image of U ? - RealScalar biggest_in_c = m_rank>0 ? c.corner(TopLeft, m_rank, c.cols()).cwise().abs().maxCoeff() : 0; + RealScalar biggest_in_c = m_rank>0 ? c.corner(TopLeft, m_rank, c.cols()).cwise().abs().maxCoeff() : RealScalar(0); for(int col = 0; col < c.cols(); ++col) for(int row = m_rank; row < c.rows(); ++row) if(!ei_isMuchSmallerThan(c.coeff(row,col), biggest_in_c, m_precision))