aboutsummaryrefslogtreecommitdiffhomepage
path: root/Eigen/src/LU/PartialPivLU.h
diff options
context:
space:
mode:
authorGravatar Gael Guennebaud <g.gael@free.fr>2019-02-11 22:56:19 +0100
committerGravatar Gael Guennebaud <g.gael@free.fr>2019-02-11 22:56:19 +0100
commitbdcb5f33043c559c7100f8fd5eb55fbbd0cdfc69 (patch)
treee0fe7726640e96a397a23888bdaaf0be3ccd6a48 /Eigen/src/LU/PartialPivLU.h
parent2edfc6807d84c3eff5abfeca809eca676c29ed7a (diff)
Let's properly use Score instead of std::abs, and remove deprecated FIXME ( a /= b does a/b and not a * (1/b) as it was a long time ago...)
Diffstat (limited to 'Eigen/src/LU/PartialPivLU.h')
-rw-r--r--Eigen/src/LU/PartialPivLU.h4
1 files changed, 1 insertions, 3 deletions
diff --git a/Eigen/src/LU/PartialPivLU.h b/Eigen/src/LU/PartialPivLU.h
index 12e72880d..cba7a0fcf 100644
--- a/Eigen/src/LU/PartialPivLU.h
+++ b/Eigen/src/LU/PartialPivLU.h
@@ -387,8 +387,6 @@ struct partial_lu_impl
++nb_transpositions;
}
- // FIXME shall we introduce a safe quotient expression in cas 1/lu.coeff(k,k)
- // overflow but not the actual quotient?
lu.col(k).tail(fix<RRows>(rrows)) /= lu.coeff(k,k);
}
else if(first_zero_pivot==-1)
@@ -407,7 +405,7 @@ struct partial_lu_impl
{
Index k = endk;
row_transpositions[k] = PivIndex(k);
- if (std::abs(lu(k, k)) == 0 && first_zero_pivot == -1)
+ if (Scoring()(lu(k, k)) == Score(0) && first_zero_pivot == -1)
first_zero_pivot = k;
}