-
Notifications
You must be signed in to change notification settings - Fork 29
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
eigs() doesn't work in my computer, but works in another one. (Arpack 0.5.4) #147
Comments
Generalized EVP with symmetric positive (semi-)definite sparse matrices. Arpack 0.5.4 fails to converge, 0.5.3 converges. Called as
The error is
|
@AlvaroPY : Could you try 0.5.3 please? |
Thank you!! With Arpack 0.5.3, at every first run I have another type of error:
But at the following runs everything goes well with Arpack 0.5.3, in both of my computers. |
I've been having the same error but for the :LM option. My sparse matrix will trigger the automatic shift-and-invert (on both 5.3 and 5.4 versions), and I think that could be related to the commit on libarpack shown by @PetrKryslUCSD. The omission of the dot when inverting on the composed dmap function may be causing issues? (Not sure if at that spot it deals with a matrix or a scalar) |
@AlvaroPY Could you rename present issue to reflect that 0.5.4 version is broken. And create a new issue for 0.5.3? |
This repository appears to have lost its maintainers. I would suggest to switch to ArnoldiMethod.jl. |
I will do it now! And thanks for the suggestion. Since my last post here I decided to continue on with ArnoldiMethod.jl (KrylovKit.jl having some limitations with non-Hermitian generalized EVP) rather than downgrading to Arpack 0.5.3, despite some efficiency advantage with the latter. |
I had some free time and tried some tests. Seems that (at least for my case) the culprit was after performing the shift and invert, and switching from LM to SM (to preserve the originally intended eigenvalue search), a new conditional (the second capture from @PetrKryslUCSD ) would also change the I've created a fork reverting that conditional, and for my case now it works. @AlvaroPY, could you try if the fork works for your case? |
Hi. Is 0.5.0 Version around? I need it :( |
You should certainly be able to install version v0.5.0 as |
I take it back. ArnoldiMethod apparently does not converge for symmetric, spd pencil. The eigenvectors are bad. JuliaLinearAlgebra/ArnoldiMethod.jl#114 |
Probably this is not related with the package itself, but I would greatly appreciate your help. I have the same code running in two different computers with the same version of Julia (1.7.2) and Arpack (0.5.4). In one computer, everything works well (I have two instances of EVPs, one generalized with dense complex non-Hermitian matrices and the other, generalized with sparse real symmetric positive definite matrices). I've compared to eigen() from LinearAlgebra and ArnoldiMethod and the solutions are ok. But in the other computer, I have constantly the following error:
Since eigs() from Arpack is faster than ArnoldiMethod, I would like to use it in my computer.
The text was updated successfully, but these errors were encountered: