diff options
| author |  Gael Guennebaud <g.gael@free.fr> | 2015-03-06 16:11:22 +0100 |
|---|---|---|
| committer | Gael Guennebaud <g.gael@free.fr> | 2015-03-06 16:11:22 +0100 |
| commit | 87681e508f1da5064f5b2699f17049306ac9a2ca (patch) | |
| tree | 2cf24d11a02f75283af62dbf663fabf831decddb /Eigen/src/Core/products/GeneralMatrixMatrix.h | |
| parent | 4c8eeeaed6e5fe3f62c06e743971cefb5eaff359 (diff) | |
bug #978: early return for vanishing products
Diffstat (limited to 'Eigen/src/Core/products/GeneralMatrixMatrix.h')
| -rw-r--r-- | Eigen/src/Core/products/GeneralMatrixMatrix.h | 2 |
1 files changed, 2 insertions, 0 deletions
diff --git a/Eigen/src/Core/products/GeneralMatrixMatrix.h b/Eigen/src/Core/products/GeneralMatrixMatrix.h index c76f48154..7fd707ed7 100644 --- a/Eigen/src/Core/products/GeneralMatrixMatrix.h +++ b/Eigen/src/Core/products/GeneralMatrixMatrix.h @@ -457,6 +457,8 @@ struct generic_product_impl<Lhs,Rhs,DenseShape,DenseShape,GemmProduct> static void scaleAndAddTo(Dest& dst, const Lhs& a_lhs, const Rhs& a_rhs, const Scalar& alpha) { eigen_assert(dst.rows()==a_lhs.rows() && dst.cols()==a_rhs.cols()); + if(a_lhs.cols()==0 || a_lhs.rows()==0 || a_rhs.cols()==0) + return; typename internal::add_const_on_value_type<ActualLhsType>::type lhs = LhsBlasTraits::extract(a_lhs); typename internal::add_const_on_value_type<ActualRhsType>::type rhs = RhsBlasTraits::extract(a_rhs); |