diff options
| author |  Gael Guennebaud <g.gael@free.fr> | 2017-06-09 14:38:04 +0200 |
|---|---|---|
| committer | Gael Guennebaud <g.gael@free.fr> | 2017-06-09 14:38:04 +0200 |
| commit | ba5cab576a7615af12389ff159c6ed57b5312d5e (patch) | |
| tree | e4200fcea398221324069d837cdd68a2863dbd96 /Eigen/src/Core/products/GeneralMatrixMatrixTriangular_BLAS.h | |
| parent | 90168c003d80236e935589d56f9fbd8e9c7aaa75 (diff) | |
bug #1405: enable StrictlyLower/StrictlyUpper triangularView as the destination of matrix*matrix products.
Diffstat (limited to 'Eigen/src/Core/products/GeneralMatrixMatrixTriangular_BLAS.h')
| -rw-r--r-- | Eigen/src/Core/products/GeneralMatrixMatrixTriangular_BLAS.h | 2 |
1 files changed, 1 insertions, 1 deletions
diff --git a/Eigen/src/Core/products/GeneralMatrixMatrixTriangular_BLAS.h b/Eigen/src/Core/products/GeneralMatrixMatrixTriangular_BLAS.h index 5b7c15cca..41e18ff07 100644 --- a/Eigen/src/Core/products/GeneralMatrixMatrixTriangular_BLAS.h +++ b/Eigen/src/Core/products/GeneralMatrixMatrixTriangular_BLAS.h @@ -52,7 +52,7 @@ struct general_matrix_matrix_triangular_product<Index,Scalar,LhsStorageOrder,Con static EIGEN_STRONG_INLINE void run(Index size, Index depth,const Scalar* lhs, Index lhsStride, \ const Scalar* rhs, Index rhsStride, Scalar* res, Index resStride, Scalar alpha, level3_blocking<Scalar, Scalar>& blocking) \ { \ - if (lhs==rhs) { \ + if ( lhs==rhs && ((UpLo&(Lower|Upper)==UpLo)) ) { \ general_matrix_matrix_rankupdate<Index,Scalar,LhsStorageOrder,ConjugateLhs,ColMajor,UpLo> \ ::run(size,depth,lhs,lhsStride,rhs,rhsStride,res,resStride,alpha,blocking); \ } else { \ |