GitHub-Proxy/eigen
1
0
Fork 0
mirror of https://gitlab.com/libeigen/eigen.git synced 2025-08-13 20:26:03 +08:00
Code Issues Packages Projects Releases Wiki Activity

Fix bug #540: SelfAdjointEigenSolver improperly used the upper triangular part to extract the scaling factor.

Browse Source
This commit is contained in:
Sergey Popov 2013-01-12 12:07:49 +01:00
parent 7262cf783c
commit 761b3bbb69
1 changed files with 3 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
Eigen/src/Eigenvalues/SelfAdjointEigenSolver.h
View File

@ -409,9 +409,10 @@ SelfAdjointEigenSolver<MatrixType>& SelfAdjointEigenSolver<MatrixType>
MatrixType& mat = m_eivec; MatrixType& mat = m_eivec;
// map the matrix coefficients to [-1:1] to avoid over- and underflow. // map the matrix coefficients to [-1:1] to avoid over- and underflow.
RealScalar scale = matrix.cwiseAbs().maxCoeff(); mat = matrix.template triangularView<Lower>();
RealScalar scale = mat.cwiseAbs().maxCoeff();
if(scale==RealScalar(0)) scale = RealScalar(1); if(scale==RealScalar(0)) scale = RealScalar(1);
mat = matrix / scale; mat.template triangularView<Lower>() /= scale;
m_subdiag.resize(n-1); m_subdiag.resize(n-1);
internal::tridiagonalization_inplace(mat, diag, m_subdiag, computeEigenvectors); internal::tridiagonalization_inplace(mat, diag, m_subdiag, computeEigenvectors);