Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Missing eigenvalues #90

Closed
ttx002000 opened this issue Jun 16, 2024 · 2 comments
Closed

Missing eigenvalues #90

ttx002000 opened this issue Jun 16, 2024 · 2 comments

Comments

@ttx002000
Copy link

ttx002000 commented Jun 16, 2024

I have the some big sparse matrix to be diagonalized. When compared with the Arpack results, it seems that there are multiple eigenvalues missing in the degenerate subspace using Krylovkit.

The only change I make in my code is changing from

MB_spectrum,ζ,info=eigsolve(sparse(matrix_index1,matrix_index2,matrix_value,length(state_can),length(state_can)),25,:SR)

to

MB_spectrum,ζ=eigs(sparse(matrix_index1,matrix_index2,matrix_value,length(state_can),length(state_can)), nev=25,which=:SR)

The second one use Arpack method, and here are the eigenvalues it obtained (the first 15 of it)
image

The first one use KrylovKit method, and here are the eigenvalues it obtained (the first 15 of it)
image

It seems that Krylovkit missed 2 of the 6 eigenvalues at 3.212527.

@lkdvos
Copy link
Collaborator

lkdvos commented Jun 17, 2024

Hi ttx,
Note that really, Krylov methods are not super well-suited for degenerate subspaces. (See for example the note in the docstring)
This being said, often the Krylov methods still work, but you typically have to play around a bit with the parameters. If you make the tol more strict, or the krylovdim a bit larger, it might then find these other eigenvectors as well.

@lkdvos
Copy link
Collaborator

lkdvos commented Jun 24, 2024

Did you manage to find the other eigenvalues with the higher tolerances? Can I close this issue?

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants