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

ARPACK feature simplification #129

Open
jarlebring opened this issue May 9, 2021 · 0 comments
Open

ARPACK feature simplification #129

jarlebring opened this issue May 9, 2021 · 0 comments

Comments

@jarlebring
Copy link
Contributor

jarlebring commented May 9, 2021

Based on discussion in #120 we want to make a cleaner interface with fortran arpack, in order to lessen the number of corner cases which are making the code difficult to maintain.

We can reduce the number of target choices, e.g., LM, LR, SR, LI, SI (and possible general eigsorter like in KrylovKit). In particular, in order to get functionality SM or shift-and-invert one would instead need to explicitly provide a transformed problem, e.g., with LinearMap:

   Afact=factorize(A-s*I);
   Atransformed=LinearMap(x->Afact\x,size(A,1))
   (d,v)=eigs(Atransformed,which=:LM)
   d = 1 ./(d .+ s)

The advantage is that we would be able to remove features associated with sigma and explicittransform.

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

1 participant