From 1578421ed14c23fa5c7ab3c818a069f1c1cefb8a Mon Sep 17 00:00:00 2001 From: Gael Guennebaud Date: Wed, 15 Jul 2009 14:20:45 +0200 Subject: fix issue #25 : the problem was that we assumed Dynamic was a multiple of a packet size (also disable the test of blueNorm) --- Eigen/src/QR/Tridiagonalization.h | 1 + 1 file changed, 1 insertion(+) (limited to 'Eigen/src/QR/Tridiagonalization.h') diff --git a/Eigen/src/QR/Tridiagonalization.h b/Eigen/src/QR/Tridiagonalization.h index bd8ff4fe3..59bffc12a 100644 --- a/Eigen/src/QR/Tridiagonalization.h +++ b/Eigen/src/QR/Tridiagonalization.h @@ -279,6 +279,7 @@ Tridiagonalization::matrixQ(void) const Scalar tmp = m_matrix.coeff(i+1,i); m_matrix.const_cast_derived().coeffRef(i+1,i) = 1; + // TODO this product could be optimized by processing the submatrix per panel of at least 4 columns matQ.corner(BottomRight,n-i-1,n-i-1) -= ((m_hCoeffs.coeff(i) * m_matrix.col(i).end(n-i-1)) * (m_matrix.col(i).end(n-i-1).adjoint() * matQ.corner(BottomRight,n-i-1,n-i-1)).lazy()).lazy(); -- cgit v1.2.3