From 7fdf2189516e2d4548c04ed4e18b7a28c28eb77c Mon Sep 17 00:00:00 2001
From: Gael Guennebaud <g.gael@free.fr>
Date: Thu, 17 Jun 2010 10:17:22 +0200
Subject: makes trmv works with the triangular matrix on the right

---
 test/product_trmv.cpp | 6 ++++++
 1 file changed, 6 insertions(+)

(limited to 'test/product_trmv.cpp')

diff --git a/test/product_trmv.cpp b/test/product_trmv.cpp
index f0962557a..2f5743187 100644
--- a/test/product_trmv.cpp
+++ b/test/product_trmv.cpp
@@ -76,6 +76,12 @@ template<typename MatrixType> void trmv(const MatrixType& m)
   VERIFY((m3.adjoint() * (s1*v1.conjugate())).isApprox(m1.adjoint().template triangularView<Eigen::Upper>() * (s1*v1.conjugate()), largerEps));
   m3 = m1.template triangularView<Eigen::UnitUpper>();
 
+  // check transposed cases:
+  m3 = m1.template triangularView<Eigen::Lower>();
+  VERIFY((v1.transpose() * m3).isApprox(v1.transpose() * m1.template triangularView<Eigen::Lower>(), largerEps));
+  VERIFY((v1.adjoint() * m3).isApprox(v1.adjoint() * m1.template triangularView<Eigen::Lower>(), largerEps));
+  VERIFY((v1.adjoint() * m3.adjoint()).isApprox(v1.adjoint() * m1.template triangularView<Eigen::Lower>().adjoint(), largerEps));
+
   // TODO check with sub-matrices
 }
 
-- 
cgit v1.2.3