diff options
author | 2010-04-02 21:33:34 +0100 | |
---|---|---|
committer | 2010-04-02 21:33:34 +0100 | |
commit | d88d1cfa626dcc10ffb0ec53284bdc9f3bb991c7 (patch) | |
tree | 41eb52b518b2b49601e31225c9e24a05a54c015e /Eigen/src/Eigenvalues/RealSchur.h | |
parent | 79e1ce609319e749d8d9a9aa2ffbcf0600ab1b93 (diff) | |
parent | 9d6afdeb22d1ccc17a2d97966163c6d8f7651047 (diff) |
Merge.
Diffstat (limited to 'Eigen/src/Eigenvalues/RealSchur.h')
-rw-r--r-- | Eigen/src/Eigenvalues/RealSchur.h | 7 |
1 files changed, 1 insertions, 6 deletions
diff --git a/Eigen/src/Eigenvalues/RealSchur.h b/Eigen/src/Eigenvalues/RealSchur.h index eb400c815..395b80089 100644 --- a/Eigen/src/Eigenvalues/RealSchur.h +++ b/Eigen/src/Eigenvalues/RealSchur.h @@ -131,17 +131,12 @@ void RealSchur<MatrixType>::hqr2() Scalar exshift = 0.0; Scalar p=0,q=0,r=0,s=0,z=0,w,x,y; - // Store roots isolated by balanc and compute matrix norm + // Compute matrix norm // FIXME to be efficient the following would requires a triangular reduxion code // Scalar norm = m_matT.upper().cwiseAbs().sum() + m_matT.corner(BottomLeft,n,n).diagonal().cwiseAbs().sum(); Scalar norm = 0.0; for (int j = 0; j < size; ++j) { - // FIXME what's the purpose of the following since the condition is always false - if ((j < low) || (j > high)) - { - m_eivalues.coeffRef(j) = ComplexScalar(m_matT.coeff(j,j), 0.0); - } norm += m_matT.row(j).segment(std::max(j-1,0), size-std::max(j-1,0)).cwiseAbs().sum(); } |