diff options
author | Benoit Jacob <jacob.benoit.1@gmail.com> | 2008-08-04 04:45:59 +0000 |
---|---|---|
committer | Benoit Jacob <jacob.benoit.1@gmail.com> | 2008-08-04 04:45:59 +0000 |
commit | c2f8ecf46683adcab0db05199ee2ebe15e6ada4a (patch) | |
tree | 34e909defb2401c96332663c127225be6d08700c /Eigen/src/LU/Determinant.h | |
parent | f81dfcf00b8fb26bd21495023799118fa444870a (diff) |
* LU decomposition, supporting all rectangular matrices, with full
pivoting for better numerical stability. For now the only application is
determinant.
* New determinant unit-test.
* Disable most of Swap.h for now as it makes LU fail (mysterious).
Anyway Swap needs a big overhaul as proposed on IRC.
* Remnants of old class Inverse removed.
* Some warnings fixed.
Diffstat (limited to 'Eigen/src/LU/Determinant.h')
-rw-r--r-- | Eigen/src/LU/Determinant.h | 4 |
1 files changed, 2 insertions, 2 deletions
diff --git a/Eigen/src/LU/Determinant.h b/Eigen/src/LU/Determinant.h index fdb1673d8..96e400564 100644 --- a/Eigen/src/LU/Determinant.h +++ b/Eigen/src/LU/Determinant.h @@ -63,7 +63,7 @@ const typename Derived::Scalar ei_bruteforce_det(const MatrixBase<Derived>& m) - ei_bruteforce_det4_helper(m,1,3,0,2) + ei_bruteforce_det4_helper(m,2,3,0,1); default: - assert(false); + ei_internal_assert(false); } } @@ -85,7 +85,7 @@ typename ei_traits<Derived>::Scalar MatrixBase<Derived>::determinant() const return derived().diagonal().redux(ei_scalar_product_op<Scalar>()); } else if(rows() <= 4) return ei_bruteforce_det(derived()); - else assert(false); // unimplemented for now + else return lu().determinant(); } #endif // EIGEN_DETERMINANT_H |