Optimization: enable cache-efficient application of HouseholderSequence.

This commit is contained in:
Gael Guennebaud
2014-09-04 09:15:59 +02:00
parent 80993b95d3
commit 8846aa6d1b
4 changed files with 89 additions and 28 deletions

View File

@@ -215,10 +215,10 @@ void upperbidiagonalization_blocked_helper(MatrixType& A,
if(k) u_k -= U_k1.adjoint() * X.row(k).head(k).adjoint();
}
// 5 - construct right Householder transform in-placecols
// 5 - construct right Householder transform in-place
u_k.makeHouseholderInPlace(tau_u, upper_diagonal[k]);
// this eases the application of Householder transforAions
// this eases the application of Householder transformations
// A(k,k+1) will store tau_u later
A(k,k+1) = Scalar(1);