aboutsummaryrefslogtreecommitdiffhomepage
path: root/test/adjoint.cpp
diff options
context:
space:
mode:
authorGravatar Benoit Jacob <jacob.benoit.1@gmail.com>2008-05-28 04:38:16 +0000
committerGravatar Benoit Jacob <jacob.benoit.1@gmail.com>2008-05-28 04:38:16 +0000
commitaebecae510dd29f5573d3f86dfed526e6d8be9a8 (patch)
treeb9bb9d2d15011409fe5f66e1fbb56f1998f5cc91 /test/adjoint.cpp
parent559233c73e86474d67f5730f9b20e46ccea93b7f (diff)
* find the proper way of nesting the expression in Flagged:
finally that's more subtle than just using ei_nested, because when flagging with NestByValueBit we want to store the expression by value already, regardless of whether it already had the NestByValueBit set. * rename temporary() ----> nestByValue() * move the old Product.h to disabled/, replace by what was ProductWIP.h * tweak -O and -g flags for tests and examples * reorder the tests -- basic things go first * simplifications, e.g. in many methoeds return derived() and count on implicit casting to the actual return type. * strip some not-really-useful stuff from the heaviest tests
Diffstat (limited to 'test/adjoint.cpp')
-rw-r--r--test/adjoint.cpp14
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