diff options
| author |  Gael Guennebaud <g.gael@free.fr> | 2010-06-24 16:42:43 +0200 |
|---|---|---|
| committer | Gael Guennebaud <g.gael@free.fr> | 2010-06-24 16:42:43 +0200 |
| commit | 905beb0953c2ae5064f4b7b5bf534555a0010220 (patch) | |
| tree | d5375d5cafdff6c812bf31dad5e1e7a27622128a /Eigen/src/Core/products/SelfadjointMatrixMatrix.h | |
| parent | af38bccd3d5a5b2a67e3a0ee39a48ad1cf6e2084 (diff) | |
fix symm
Diffstat (limited to 'Eigen/src/Core/products/SelfadjointMatrixMatrix.h')
| -rw-r--r-- | Eigen/src/Core/products/SelfadjointMatrixMatrix.h | 2 |
1 files changed, 2 insertions, 0 deletions
diff --git a/Eigen/src/Core/products/SelfadjointMatrixMatrix.h b/Eigen/src/Core/products/SelfadjointMatrixMatrix.h index cc9333384..563feb094 100644 --- a/Eigen/src/Core/products/SelfadjointMatrixMatrix.h +++ b/Eigen/src/Core/products/SelfadjointMatrixMatrix.h @@ -262,6 +262,8 @@ struct ei_product_selfadjoint_matrix<Scalar,Index,LhsStorageOrder,true,Conjugate Index mc = rows; // cache block size along the M direction Index nc = cols; // cache block size along the N direction computeProductBlockingSizes<Scalar,Scalar>(kc, mc, nc); + // kc must smaller than mc + kc = std::min(kc,mc); Scalar* blockA = ei_aligned_stack_new(Scalar, kc*mc); std::size_t sizeB = kc*Blocking::PacketSize*Blocking::nr + kc*cols; |