diff options
author | 2009-02-03 17:50:35 +0000 | |
---|---|---|
committer | 2009-02-03 17:50:35 +0000 | |
commit | b9a82be7271f21f78c1bce858bace71407c070c1 (patch) | |
tree | 5c2170fe648883a57864019a7ca3433620a9dfe3 /test | |
parent | 6c5868cc997c50ca3f61c81ea551e5b234f45cd7 (diff) |
Add full pivoting to LDLT decomposition.
Diffstat (limited to 'test')
-rw-r--r-- | test/cholesky.cpp | 3 | ||||
-rw-r--r-- | test/main.h | 2 |
2 files changed, 3 insertions, 2 deletions
diff --git a/test/cholesky.cpp b/test/cholesky.cpp index 2e3353d21..b3e0df438 100644 --- a/test/cholesky.cpp +++ b/test/cholesky.cpp @@ -83,7 +83,8 @@ template<typename MatrixType> void cholesky(const MatrixType& m) { LDLT<SquareMatrixType> ldlt(symm); VERIFY(ldlt.isPositiveDefinite()); - VERIFY_IS_APPROX(symm, ldlt.matrixL() * ldlt.vectorD().asDiagonal() * ldlt.matrixL().adjoint()); + // TODO(keir): This doesn't make sense now that LDLT pivots. + //VERIFY_IS_APPROX(symm, ldlt.matrixL() * ldlt.vectorD().asDiagonal() * ldlt.matrixL().adjoint()); ldlt.solve(vecB, &vecX); VERIFY_IS_APPROX(symm * vecX, vecB); ldlt.solve(matB, &matX); diff --git a/test/main.h b/test/main.h index 0356f60bf..489f78bff 100644 --- a/test/main.h +++ b/test/main.h @@ -44,7 +44,7 @@ namespace Eigen #define EI_PP_MAKE_STRING2(S) #S #define EI_PP_MAKE_STRING(S) EI_PP_MAKE_STRING2(S) - +#define EIGEN_DEFAULT_IO_FORMAT IOFormat(4, AlignCols, " ", "\n", "", "", "", "") #ifndef EIGEN_NO_ASSERTION_CHECKING |