diff options
Diffstat (limited to 'test/adjoint.cpp')
| -rw-r--r-- | test/adjoint.cpp | 14 |
1 files changed, 3 insertions, 11 deletions
diff --git a/test/adjoint.cpp b/test/adjoint.cpp index ae1b002e5..9a15c4986 100644 --- a/test/adjoint.cpp +++ b/test/adjoint.cpp @@ -51,23 +51,15 @@ template<typename MatrixType> void adjoint(const MatrixType& m) Scalar s1 = ei_random<Scalar>(), s2 = ei_random<Scalar>(); - // check involutivity of adjoint, transpose, conjugate - VERIFY_IS_APPROX(m1.transpose().transpose(), m1); - VERIFY_IS_APPROX(m1.conjugate().conjugate(), m1); - VERIFY_IS_APPROX(m1.adjoint().adjoint(), m1); - // check basic compatibility of adjoint, transpose, conjugate VERIFY_IS_APPROX(m1.transpose().conjugate().adjoint(), m1); VERIFY_IS_APPROX(m1.adjoint().conjugate().transpose(), m1); - if(!NumTraits<Scalar>::IsComplex) - VERIFY_IS_APPROX(m1.adjoint().transpose(), m1); // check multiplicative behavior - VERIFY_IS_APPROX((m1.transpose() * m2).transpose(), m2.transpose() * m1); + std::cout << (m1.adjoint() * m2).adjoint() << std::endl; + std::cout << "------------------------------" << std::endl; + std::cout << m2.adjoint() * m1 << std::endl; VERIFY_IS_APPROX((m1.adjoint() * m2).adjoint(), m2.adjoint() * m1); - VERIFY_IS_APPROX((m1.transpose() * m2).conjugate(), m1.adjoint() * m2.conjugate()); - VERIFY_IS_APPROX((s1 * m1).transpose(), s1 * m1.transpose()); - VERIFY_IS_APPROX((s1 * m1).conjugate(), ei_conj(s1) * m1.conjugate()); VERIFY_IS_APPROX((s1 * m1).adjoint(), ei_conj(s1) * m1.adjoint()); // check basic properties of dot, norm, norm2 |