diff options
Diffstat (limited to 'Eigen/src/Householder/BlockHouseholder.h')
-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 c17ca10e0..60da512e1 100644 --- a/Eigen/src/Householder/BlockHouseholder.h +++ b/Eigen/src/Householder/BlockHouseholder.h @@ -66,7 +66,8 @@ void ei_apply_block_householder_on_the_left(MatrixType& mat, const VectorsType& const TriangularView<VectorsType, UnitLower>& V(vectors); // A -= V T V^* A - Matrix<typename MatrixType::Scalar,Dynamic,Dynamic> tmp = V.adjoint() * mat; + Matrix<typename MatrixType::Scalar,VectorsType::ColsAtCompileTime,MatrixType::ColsAtCompileTime,0, + VectorsType::MaxColsAtCompileTime,MatrixType::MaxColsAtCompileTime> tmp = V.adjoint() * mat; // FIXME add .noalias() once the triangular product can work inplace tmp = T.template triangularView<Upper>().adjoint() * tmp; mat.noalias() -= V * tmp; |