diff options
| author |  Desire NUENTSA <desire.nuentsa_wakam@inria.fr> | 2013-06-28 22:27:45 +0200 |
|---|---|---|
| committer | Desire NUENTSA <desire.nuentsa_wakam@inria.fr> | 2013-06-28 22:27:45 +0200 |
| commit | 9f035c876a9d84f32b39bacba179e9b2a542b57c (patch) | |
| tree | c3815bcd69867236e3e0a3dd3e5c3ccc34c64fe8 /test/sparse_product.cpp | |
| parent | fc27cbd914b8afcf03b6120e3b23206bdc57b0a5 (diff) | |
Fiw bug #553: add support for sparse matrix time sparse self-adjoint view products
Diffstat (limited to 'test/sparse_product.cpp')
| -rw-r--r-- | test/sparse_product.cpp | 9 |
1 files changed, 9 insertions, 0 deletions
diff --git a/test/sparse_product.cpp b/test/sparse_product.cpp index 338c3d378..51eed428b 100644 --- a/test/sparse_product.cpp +++ b/test/sparse_product.cpp @@ -201,7 +201,16 @@ template<typename SparseMatrixType> void sparse_product() VERIFY_IS_APPROX(x=mUp.template selfadjointView<Upper>()*b, refX=refS*b); VERIFY_IS_APPROX(x=mLo.template selfadjointView<Lower>()*b, refX=refS*b); VERIFY_IS_APPROX(x=mS.template selfadjointView<Upper|Lower>()*b, refX=refS*b); + + // sparse selfadjointView * sparse + SparseMatrixType mSres(rows,rows); + VERIFY_IS_APPROX(mSres = mLo.template selfadjointView<Lower>()*mS, + refX = refLo.template selfadjointView<Lower>()*refS); + // sparse * sparse selfadjointview + VERIFY_IS_APPROX(mSres = mS * mLo.template selfadjointView<Lower>(), + refX = refS * refLo.template selfadjointView<Lower>()); } + } // New test for Bug in SparseTimeDenseProduct |