diff options
author | Gael Guennebaud <g.gael@free.fr> | 2019-02-11 22:56:19 +0100 |
---|---|---|
committer | Gael Guennebaud <g.gael@free.fr> | 2019-02-11 22:56:19 +0100 |
commit | bdcb5f33043c559c7100f8fd5eb55fbbd0cdfc69 (patch) | |
tree | e0fe7726640e96a397a23888bdaaf0be3ccd6a48 | |
parent | 2edfc6807d84c3eff5abfeca809eca676c29ed7a (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...)
-rw-r--r-- | Eigen/src/LU/PartialPivLU.h | 4 |
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; } |