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.
diff --git a/test/nomalloc.cpp b/test/nomalloc.cpp
index 060276a..d85e9e5 100644
--- a/test/nomalloc.cpp
+++ b/test/nomalloc.cpp
@@ -81,11 +81,12 @@
m2.template selfadjointView<Lower>().rankUpdate(m1.row(0),-1);
// The following fancy matrix-matrix products are not safe yet regarding static allocation
-// m1 += m1.template triangularView<Upper>() * m2.col(;
-// m1.template selfadjointView<Lower>().rankUpdate(m2);
-// m1 += m1.template triangularView<Upper>() * m2;
+// m1.col(1) += m1.template triangularView<Upper>() * m2.col(0);
+ m2.template selfadjointView<Lower>().rankUpdate(m1);
+ m2 += m2.template triangularView<Upper>() * m1;
+ m2.template triangularView<Upper>() = m2 * m2;
// m1 += m1.template selfadjointView<Lower>() * m2;
-// VERIFY_IS_APPROX(m1,m1);
+ VERIFY_IS_APPROX(m2,m2);
}
template<typename Scalar>