From ad1dc50b5779f4e6f87eb9e5865ce6fcfcc5a828 Mon Sep 17 00:00:00 2001 From: Gael Guennebaud Date: Sun, 3 Nov 2013 13:19:55 +0100 Subject: [PATCH] Check for minimal norm solutions --- test/jacobisvd.cpp | 8 +------- 1 file changed, 1 insertion(+), 7 deletions(-) diff --git a/test/jacobisvd.cpp b/test/jacobisvd.cpp index e378de477..a0d8562c1 100644 --- a/test/jacobisvd.cpp +++ b/test/jacobisvd.cpp @@ -125,13 +125,7 @@ void jacobisvd_solve(const MatrixType& m, unsigned int computationOptions) RhsType3 rhs3 = C * rhs2; JacobiSVD svd3(m3, computationOptions); SolutionType x3 = svd3.solve(rhs3); - if(svd3.rank()!=rank) { - std::cout << m3 << "\n\n"; - std::cout << svd3.singularValues().transpose() << "\n"; - std::cout << svd3.rank() << " == " << rank << "\n"; - std::cout << x21.norm() << " == " << x3.norm() << "\n"; - } -// VERIFY_IS_APPROX(m3*x3, rhs3); + VERIFY_IS_APPROX(m3*x3, rhs3); VERIFY_IS_APPROX(m3*x21, rhs3); VERIFY_IS_APPROX(m2*x3, rhs2);