GitHub-Proxy/eigen
1
0
Fork 0
mirror of https://gitlab.com/libeigen/eigen.git synced 2025-05-02 00:34:14 +08:00
Code Issues Packages Projects Releases Wiki Activity

Widen short-circuiting ReciprocalConditionNumberEstimate so we don't call InverseMatrixL1NormEstimate for dec.rows() <= 1.

Browse Source
This commit is contained in:
Rasmus Munk Larsen 2016-04-07 16:45:40 -07:00
parent d51803a728
commit 283c51cd5e
1 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
Eigen/src/Core/ConditionEstimator.h
View File

@ -66,7 +66,7 @@ typename Decomposition::RealScalar ReciprocalConditionNumberEstimate(
eigen_assert(matrix.rows() == dec.rows()); eigen_assert(matrix.rows() == dec.rows());
eigen_assert(matrix.cols() == dec.cols()); eigen_assert(matrix.cols() == dec.cols());
eigen_assert(matrix.rows() == matrix.cols()); eigen_assert(matrix.rows() == matrix.cols());
if (dec.rows() == 0) { if (dec.rows() <= 1) {
return static_cast<typename Decomposition::RealScalar>(1); return static_cast<typename Decomposition::RealScalar>(1);
} }
return ReciprocalConditionNumberEstimate(MatrixL1Norm(matrix), dec); return ReciprocalConditionNumberEstimate(MatrixL1Norm(matrix), dec);
@ -93,7 +93,7 @@ typename Decomposition::RealScalar ReciprocalConditionNumberEstimate(
typename Decomposition::RealScalar matrix_norm, const Decomposition& dec) { typename Decomposition::RealScalar matrix_norm, const Decomposition& dec) {
typedef typename Decomposition::RealScalar RealScalar; typedef typename Decomposition::RealScalar RealScalar;
eigen_assert(dec.rows() == dec.cols()); eigen_assert(dec.rows() == dec.cols());
if (dec.rows() == 0) { if (dec.rows() <= 1) {
return static_cast<RealScalar>(1); return static_cast<RealScalar>(1);
} }
if (matrix_norm == static_cast<RealScalar>(0)) { if (matrix_norm == static_cast<RealScalar>(0)) {