Julia and eigs in NEP-PACK

nep-pack / NonlinearEigenproblems.jl

Nonlinear eigenvalue problems in Julia: Iterative methods and benchmarks

https://nep-pack.github.io/NonlinearEigenproblems.jl

MIT License

92 stars 15 forks source link

Julia and eigs in NEP-PACK #35

Closed jarlebring closed 6 years ago

jarlebring commented 6 years ago

As we have observed several times e.g. #1, julias eigs is broken for a large class of generalized eigenvalue problems. As pointed out here https://github.com/JuliaLang/julia/issues/ 24668 this seems to be restricted to indefinite B-matrices, without a clear solution in sight. Our current workaround with MATLAB is not so nice for many reasons. We should try to come up with better work-arounds.

jarlebring commented 6 years ago

I have tried to do a work-around in b3ac014c605f31fa7db162d4eddd25275bee4693 by creating a LinearMap corresponding to a shift-and-inverted matrix.

jarlebring commented 6 years ago

Since the work around with a eigs+LinearMap has not caused any problems so far, I will close this issue.