diff options
author | Gael Guennebaud <g.gael@free.fr> | 2010-07-15 20:39:20 +0200 |
---|---|---|
committer | Gael Guennebaud <g.gael@free.fr> | 2010-07-15 20:39:20 +0200 |
commit | 84fdbded4d8a7fae0cbd6fb5836f2aa0e442f356 (patch) | |
tree | 8b1f266641a2fd8c2c2f7424709dc82d09b52716 /test/product_trmm.cpp | |
parent | cf9edd995810de7575685d1164590f06edbcab3d (diff) |
add support for strictly triangular matrix in trmm though it is not really useful
Diffstat (limited to 'test/product_trmm.cpp')
-rw-r--r-- | test/product_trmm.cpp | 9 |
1 files changed, 8 insertions, 1 deletions
diff --git a/test/product_trmm.cpp b/test/product_trmm.cpp index e20b408c4..e4790fb66 100644 --- a/test/product_trmm.cpp +++ b/test/product_trmm.cpp @@ -35,7 +35,7 @@ template<typename Scalar> void trmm(int size,int /*othersize*/) DenseIndex cols = ei_random<DenseIndex>(1,size); MatrixColMaj triV(rows,cols), triH(cols,rows), upTri(cols,rows), loTri(rows,cols), - unitUpTri(cols,rows), unitLoTri(rows,cols); + unitUpTri(cols,rows), unitLoTri(rows,cols), strictlyUpTri(cols,rows), strictlyLoTri(rows,cols); MatrixColMaj ge1(rows,cols), ge2(cols,rows), ge3; MatrixRowMaj rge3; @@ -48,6 +48,8 @@ template<typename Scalar> void trmm(int size,int /*othersize*/) upTri = triH.template triangularView<Upper>(); unitLoTri = triV.template triangularView<UnitLower>(); unitUpTri = triH.template triangularView<UnitUpper>(); + strictlyLoTri = triV.template triangularView<StrictlyLower>(); + strictlyUpTri = triH.template triangularView<StrictlyUpper>(); ge1.setRandom(); ge2.setRandom(); @@ -72,6 +74,11 @@ template<typename Scalar> void trmm(int size,int /*othersize*/) VERIFY_IS_APPROX( rge3.noalias() = ge2 * triV.template triangularView<UnitLower>(), ge2 * unitLoTri); VERIFY_IS_APPROX( ge3 = ge2 * triV.template triangularView<UnitLower>(), ge2 * unitLoTri); VERIFY_IS_APPROX( ge3 = (s1*triV).adjoint().template triangularView<UnitUpper>() * ge2.adjoint(), ei_conj(s1) * unitLoTri.adjoint() * ge2.adjoint()); + + VERIFY_IS_APPROX( ge3 = triV.template triangularView<StrictlyLower>() * ge2, strictlyLoTri * ge2); + VERIFY_IS_APPROX( rge3.noalias() = ge2 * triV.template triangularView<StrictlyLower>(), ge2 * strictlyLoTri); + VERIFY_IS_APPROX( ge3 = ge2 * triV.template triangularView<StrictlyLower>(), ge2 * strictlyLoTri); + VERIFY_IS_APPROX( ge3 = (s1*triV).adjoint().template triangularView<StrictlyUpper>() * ge2.adjoint(), ei_conj(s1) * strictlyLoTri.adjoint() * ge2.adjoint()); } void test_product_trmm() |