diff options
| author |  Gael Guennebaud <g.gael@free.fr> | 2013-08-02 00:36:06 +0200 |
|---|---|---|
| committer | Gael Guennebaud <g.gael@free.fr> | 2013-08-02 00:36:06 +0200 |
| commit | e90229a429c74b302368e788a36d8bbef4d1e617 (patch) | |
| tree | 8bea0d275189a64fbfda461d686a4ab80f94c1da /test | |
| parent | cf884a9815c51c27f9172465afb32f3d1d899953 (diff) | |
reduce cancellation probablity
Diffstat (limited to 'test')
| -rw-r--r-- | test/product_trmm.cpp | 1 |
1 files changed, 1 insertions, 0 deletions
diff --git a/test/product_trmm.cpp b/test/product_trmm.cpp index 506a1aeb9..d715b9a36 100644 --- a/test/product_trmm.cpp +++ b/test/product_trmm.cpp @@ -51,6 +51,7 @@ void trmm(int rows=internal::random<int>(1,EIGEN_TEST_MAX_SIZE), ge_xs_save = ge_xs; VERIFY_IS_APPROX( (ge_xs_save + s1*triTr.conjugate() * (s2*ge_left.adjoint())).eval(), ge_xs.noalias() += (s1*mat.adjoint()).template triangularView<Mode>() * (s2*ge_left.adjoint()) ); + ge_sx.setRandom(); ge_sx_save = ge_sx; VERIFY_IS_APPROX( ge_sx_save - (ge_right.adjoint() * (-s1 * triTr).conjugate()).eval(), ge_sx.noalias() -= (ge_right.adjoint() * (-s1 * mat).adjoint().template triangularView<Mode>()).eval()); |