diff options
| author |  Christoph Hertzberg <chtz@informatik.uni-bremen.de> | 2018-12-03 16:18:15 +0100 |
|---|---|---|
| committer | Christoph Hertzberg <chtz@informatik.uni-bremen.de> | 2018-12-03 16:18:15 +0100 |
| commit | 919414b9fe2ad7fdcb0f2b2cbdf6b5322d0f2034 (patch) | |
| tree | b5269508221dd5d95b38386423619254d780d937 /Eigen/src/Core/ConditionEstimator.h | |
| parent | 0ea7ae72130cac7334823ec442f0a8a6772c9ab8 (diff) | |
bug #785: Make Cholesky decomposition work for empty matrices
Diffstat (limited to 'Eigen/src/Core/ConditionEstimator.h')
| -rw-r--r-- | Eigen/src/Core/ConditionEstimator.h | 2 |
1 files changed, 1 insertions, 1 deletions
diff --git a/Eigen/src/Core/ConditionEstimator.h b/Eigen/src/Core/ConditionEstimator.h index aa7efdc76..8c1a89129 100644 --- a/Eigen/src/Core/ConditionEstimator.h +++ b/Eigen/src/Core/ConditionEstimator.h @@ -160,7 +160,7 @@ rcond_estimate_helper(typename Decomposition::RealScalar matrix_norm, const Deco { typedef typename Decomposition::RealScalar RealScalar; eigen_assert(dec.rows() == dec.cols()); - if (dec.rows() == 0) return RealScalar(1); + if (dec.rows() == 0) return RealScalar(1)/RealScalar(0); if (matrix_norm == RealScalar(0)) return RealScalar(0); if (dec.rows() == 1) return RealScalar(1); const RealScalar inverse_matrix_norm = rcond_invmatrix_L1_norm_estimate(dec); |