diff --git a/gtsam/base/tests/testCholesky.cpp b/gtsam/base/tests/testCholesky.cpp index dd26d2004..3fd42fa58 100644 --- a/gtsam/base/tests/testCholesky.cpp +++ b/gtsam/base/tests/testCholesky.cpp @@ -98,33 +98,6 @@ TEST(cholesky, choleskyPartial) { EXPECT(assert_equal(ublas::symmetric_adaptor(ABC), actual, 1e-9)); } -/* ************************************************************************* */ -TEST(cholesky, choleskyCareful) { - Matrix A = Matrix_(7,7, - 5184.9, -4381.3, 1137.3, -4448.6, 5127.2, 3264.2, 2.4894, - -4381.3, 1.1345e+05, -28036., -13136., -1.1277e+05, -76360., 19.118, - 1137.3, -28036., 1.3193e+05, 3192.3, 27877., -1.0613e+05, -38.623, - -4448.6, -13136., 3192.3, 6417.7, 12294., 8530., -5.4028, - 5127.2, -1.1277e+05, 27877., 12294., 1.1222e+05, 75952., -18.507, - 3264.2, -76360., -1.0613e+05, 8530., 75952., 1.7642e+05, 21.269, - 2.4894, 19.118, -38.623, -5.4028, -18.507, 21.269, 0.014511); - - Matrix expected = Matrix_(7,7, - 72.006, -60.846, 15.795, -61.781, 71.206, 45.332, 0.034572, - 0.0, 331.28, -81.728, -50.999, -327.33, -222.17, 0.064059, - 0.0, 0.0, 353.55, 3.6466e-06, 0.00014052, -353.55, -0.09598, - 0.0, 0.0, 0.0, 0.024427, -2.3538, 0.33228, -0.0004039, - 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, - 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, - 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0); - - MatrixColMajor actual = A; - choleskyCareful(actual); - actual = ublas::triangular_adaptor(actual); - - EXPECT(assert_equal(expected, Matrix(actual), 1e-3)); -} - /* ************************************************************************* */ int main() { TestResult tr;