diff options
Diffstat (limited to 'test/cholesky.cpp')
| -rw-r--r-- | test/cholesky.cpp | 3 |
1 files changed, 2 insertions, 1 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); |