From e58827d2ed32cee5362e4d7d007da06a2bdc7309 Mon Sep 17 00:00:00 2001 From: Gael Guennebaud Date: Mon, 25 Jan 2016 17:16:33 +0100 Subject: bug #51: make general_matrix_matrix_triangular_product use L3-blocking helper so that general symmetric rank-updates and general-matrix-to-triangular products do not trigger dynamic memory allocation for fixed size matrices. --- test/nomalloc.cpp | 9 +++++---- 1 file changed, 5 insertions(+), 4 deletions(-) (limited to 'test/nomalloc.cpp') diff --git a/test/nomalloc.cpp b/test/nomalloc.cpp index 060276a20..d85e9e5bc 100644 --- a/test/nomalloc.cpp +++ b/test/nomalloc.cpp @@ -81,11 +81,12 @@ template void nomalloc(const MatrixType& m) m2.template selfadjointView().rankUpdate(m1.row(0),-1); // The following fancy matrix-matrix products are not safe yet regarding static allocation -// m1 += m1.template triangularView() * m2.col(; -// m1.template selfadjointView().rankUpdate(m2); -// m1 += m1.template triangularView() * m2; +// m1.col(1) += m1.template triangularView() * m2.col(0); + m2.template selfadjointView().rankUpdate(m1); + m2 += m2.template triangularView() * m1; + m2.template triangularView() = m2 * m2; // m1 += m1.template selfadjointView() * m2; -// VERIFY_IS_APPROX(m1,m1); + VERIFY_IS_APPROX(m2,m2); } template -- cgit v1.2.3