JuliaLinearAlgebra / ArnoldiMethod.jl

The Arnoldi Method with Krylov-Schur restart, natively in Julia.
https://julialinearalgebra.github.io/ArnoldiMethod.jl/dev
MIT License
96 stars 18 forks source link

Deal with new ritz values in 1:nlock range #117

Closed haampie closed 7 months ago

haampie commented 7 months ago

ArnoldiMethod.jl did not deal with repeated eigenvalues well.

The typical convergence behavior is that the second converges slower / after the first, so that when sorted by, a not yet converged eigenvalue may pop up in the range 1:nlock. I guess previously they all converged nicely in order.