diff options
author | Gael Guennebaud <g.gael@free.fr> | 2016-06-11 14:41:36 +0200 |
---|---|---|
committer | Gael Guennebaud <g.gael@free.fr> | 2016-06-11 14:41:36 +0200 |
commit | 83904a21c11ffdb88f3ad8a65ded7bf46c1a068a (patch) | |
tree | a594bab75f5acd051d1e0aae7e00a1b411c2f25f /Eigen/src/Eigenvalues/RealQZ.h | |
parent | 65d33e5898ee7642c5a2226e1080ee1d95e3a045 (diff) |
Make sure T(i+1,i)==0 when diagonalizing T(i:i+1,i:i+1)
Diffstat (limited to 'Eigen/src/Eigenvalues/RealQZ.h')
-rw-r--r-- | Eigen/src/Eigenvalues/RealQZ.h | 6 |
1 files changed, 3 insertions, 3 deletions
diff --git a/Eigen/src/Eigenvalues/RealQZ.h b/Eigen/src/Eigenvalues/RealQZ.h index c4715b954..b3a910dd9 100644 --- a/Eigen/src/Eigenvalues/RealQZ.h +++ b/Eigen/src/Eigenvalues/RealQZ.h @@ -630,11 +630,11 @@ namespace Eigen { internal::real_2x2_jacobi_svd(m_T, i, i+1, &j_left, &j_right); // Apply resulting Jacobi rotations - m_T.applyOnTheLeft(i,i+1,j_left); - m_T.applyOnTheRight(i,i+1,j_right); m_S.applyOnTheLeft(i,i+1,j_left); m_S.applyOnTheRight(i,i+1,j_right); - m_T(i,i+1) = Scalar(0); + m_T.applyOnTheLeft(i,i+1,j_left); + m_T.applyOnTheRight(i,i+1,j_right); + m_T(i+1,i) = m_T(i,i+1) = Scalar(0); if(m_computeQZ) { m_Q.applyOnTheRight(i,i+1,j_left.transpose()); |