diff options
| author | Antonio Sánchez <cantonios@google.com> | 2024-01-29 23:46:23 +0000 |
|---|---|---|
| committer | Antonio Sanchez <cantonios@google.com> | 2024-01-30 08:39:35 -0800 |
| commit | 287c8017808fb4c7a651ba4c02363a773e2f0c46 (patch) | |
| tree | 985f9d786707dd07f9347ddf876bfb10edb70448 | |
| parent | 42b04a08c40b223197bbced29101f2edb6925034 (diff) | |
| download | eigen-287c8017808fb4c7a651ba4c02363a773e2f0c46.tar.gz | |
Use stableNorm in ComplexEigenSolver.
(cherry picked from commit 0f0c76dc2985fc6722406ccdf7488359f6f23400)
| -rw-r--r-- | Eigen/src/Eigenvalues/ComplexEigenSolver.h | 5 |
1 files changed, 2 insertions, 3 deletions
diff --git a/Eigen/src/Eigenvalues/ComplexEigenSolver.h b/Eigen/src/Eigenvalues/ComplexEigenSolver.h index 081e918f1..6a99f4c5e 100644 --- a/Eigen/src/Eigenvalues/ComplexEigenSolver.h +++ b/Eigen/src/Eigenvalues/ComplexEigenSolver.h @@ -316,9 +316,8 @@ void ComplexEigenSolver<MatrixType>::doComputeEigenvectors(RealScalar matrixnorm // Compute V as V = U X; now A = U T U^* = U X D X^(-1) U^* = V D V^(-1) m_eivec.noalias() = m_schur.matrixU() * m_matX; // .. and normalize the eigenvectors - for(Index k=0 ; k<n ; k++) - { - m_eivec.col(k).normalize(); + for (Index k = 0; k < n; k++) { + m_eivec.col(k).stableNormalize(); } } |