devtools/eigen
1
0
Fork 0
mirror of https://gitlab.com/libeigen/eigen.git synced 2026-04-10 11:34:33 +08:00
Code Issues Packages Projects Releases Wiki Activity

* rewrite of the QR decomposition:

Browse Source 
- works for complex
  - allows direct access to the matrix R
* removed the scale by the matrix dimensions in MatrixBase::isMuchSmallerThan(scalar)
This commit is contained in:
Gael Guennebaud
2008-06-07 22:47:11 +00:00
parent eb7b7b2cfc
commit 4dd57b585d
3 changed files with 65 additions and 44 deletions
Download Patch File Download Diff File Expand all files Collapse all files

2
Eigen/src/QR/Tridiagonalization.h
View File

@@ -226,7 +226,7 @@ Tridiagonalization<MatrixType>::matrixQ(void) const
m_matrix.const_cast_derived().coeffRef(i+1,i) = 1;
matQ.corner(BottomRight,n-i-1,n-i-1) -=
((m_hCoeffs[i] * m_matrix.col(i).end(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();
m_matrix.const_cast_derived().coeffRef(i+1,i) = tmp;