commit af38bccd3d5a5b2a67e3a0ee39a48ad1cf6e2084
author Gael Guennebaud <g.gael@free.fr> Thu Jun 24 16:26:27 2010 +0200
committer Gael Guennebaud <g.gael@free.fr> Thu Jun 24 16:26:27 2010 +0200
tree 38768b36672241a08807fa8254034086ebf0a95b
parent e499646c74ec6853dd043043b5f6ffad22bb70e8

fix syrk

Eigen/src/Core/products/SelfadjointProduct.h

diff --git a/Eigen/src/Core/products/SelfadjointProduct.h b/Eigen/src/Core/products/SelfadjointProduct.h
index 8ce797c..9eb5de7 100644
--- a/Eigen/src/Core/products/SelfadjointProduct.h
+++ b/Eigen/src/Core/products/SelfadjointProduct.h
@@ -74,6 +74,8 @@
     Index mc = size;  // cache block size along the M direction
     Index nc = size;  // cache block size along the N direction
     computeProductBlockingSizes<Scalar,Scalar>(kc, mc, nc);
+    // !!! mc must be a multiple of nr:
+    mc = (mc/Blocking::nr)*Blocking::nr;
 
     Scalar* blockA = ei_aligned_stack_new(Scalar, kc*mc);
     std::size_t sizeB = kc*Blocking::PacketSize*Blocking::nr + kc*size;
@@ -142,7 +144,7 @@
 
   ei_selfadjoint_product<Scalar, Index,
     _ActualUType::Flags&RowMajorBit ? RowMajor : ColMajor,
-    ei_traits<MatrixType>::Flags&RowMajorBit ? RowMajor : ColMajor,
+    MatrixType::Flags&RowMajorBit ? RowMajor : ColMajor,
     !UBlasTraits::NeedToConjugate, UpLo>
     ::run(_expression().cols(), actualU.cols(), &actualU.coeff(0,0), actualU.outerStride(),
           const_cast<Scalar*>(_expression().data()), _expression().outerStride(), actualAlpha);