diff options
author | Gael Guennebaud <g.gael@free.fr> | 2017-02-20 13:44:37 +0100 |
---|---|---|
committer | Gael Guennebaud <g.gael@free.fr> | 2017-02-20 13:44:37 +0100 |
commit | 65728257036652fe1cb337a19ee68d8ec01462a3 (patch) | |
tree | d3396100754698d167842ec743905ad4b99b0d3a /Eigen/src/Householder | |
parent | a811a0469655f7ad24a01d219c0afe752214255e (diff) |
bug #1395: fix the use of compile-time vectors as inputs of JacobiSVD.
Diffstat (limited to 'Eigen/src/Householder')
-rw-r--r-- | Eigen/src/Householder/BlockHouseholder.h | 3 |
1 files changed, 2 insertions, 1 deletions
diff --git a/Eigen/src/Householder/BlockHouseholder.h b/Eigen/src/Householder/BlockHouseholder.h index 39bf8c83d..01a7ed188 100644 --- a/Eigen/src/Householder/BlockHouseholder.h +++ b/Eigen/src/Householder/BlockHouseholder.h @@ -87,7 +87,8 @@ void apply_block_householder_on_the_left(MatrixType& mat, const VectorsType& vec const TriangularView<const VectorsType, UnitLower> V(vectors); // A -= V T V^* A - Matrix<typename MatrixType::Scalar,VectorsType::ColsAtCompileTime,MatrixType::ColsAtCompileTime,0, + Matrix<typename MatrixType::Scalar,VectorsType::ColsAtCompileTime,MatrixType::ColsAtCompileTime, + (VectorsType::MaxColsAtCompileTime==1 && MatrixType::MaxColsAtCompileTime!=1)?RowMajor:ColMajor, VectorsType::MaxColsAtCompileTime,MatrixType::MaxColsAtCompileTime> tmp = V.adjoint() * mat; // FIXME add .noalias() once the triangular product can work inplace if(forward) tmp = T.template triangularView<Upper>() * tmp; |