Use stableNorm in ComplexEigenSolver.

(cherry picked from commit 0f0c76dc2985fc6722406ccdf7488359f6f23400)
This commit is contained in:
Antonio Sánchez 2024-01-29 23:46:23 +00:00 committed by Antonio Sanchez
parent 42b04a08c4
commit 287c801780

View File

@ -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();
}
}