diff options
author | Gael Guennebaud <g.gael@free.fr> | 2016-06-11 14:41:53 +0200 |
---|---|---|
committer | Gael Guennebaud <g.gael@free.fr> | 2016-06-11 14:41:53 +0200 |
commit | a3a4714abac02ba48a683c3c3967cebee2833188 (patch) | |
tree | 52e896504c71bc7d3ab7f960bd9835a79a9c558e /test/real_qz.cpp | |
parent | 83904a21c11ffdb88f3ad8a65ded7bf46c1a068a (diff) |
Add debug output.
Diffstat (limited to 'test/real_qz.cpp')
-rw-r--r-- | test/real_qz.cpp | 9 |
1 files changed, 9 insertions, 0 deletions
diff --git a/test/real_qz.cpp b/test/real_qz.cpp index a1766c6d9..45ae8d763 100644 --- a/test/real_qz.cpp +++ b/test/real_qz.cpp @@ -49,11 +49,20 @@ template<typename MatrixType> void real_qz(const MatrixType& m) for (Index i=0; i<A.cols(); i++) for (Index j=0; j<i; j++) { if (abs(qz.matrixT()(i,j))!=Scalar(0.0)) + { + std::cerr << "Error: T(" << i << "," << j << ") = " << qz.matrixT()(i,j) << std::endl; all_zeros = false; + } if (j<i-1 && abs(qz.matrixS()(i,j))!=Scalar(0.0)) + { + std::cerr << "Error: S(" << i << "," << j << ") = " << qz.matrixS()(i,j) << std::endl; all_zeros = false; + } if (j==i-1 && j>0 && abs(qz.matrixS()(i,j))!=Scalar(0.0) && abs(qz.matrixS()(i-1,j-1))!=Scalar(0.0)) + { + std::cerr << "Error: S(" << i << "," << j << ") = " << qz.matrixS()(i,j) << " && S(" << i-1 << "," << j-1 << ") = " << qz.matrixS()(i-1,j-1) << std::endl; all_zeros = false; + } } VERIFY_IS_EQUAL(all_zeros, true); VERIFY_IS_APPROX(qz.matrixQ()*qz.matrixS()*qz.matrixZ(), A); |