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.

Sign up for GitHub

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

Jump to bottom

Julia and eigs in NEP-PACK #35

Closed
jarlebring opened this issue Jun 1, 2018 · 2 comments
Closed

Julia and eigs in NEP-PACK #35

jarlebring opened this issue Jun 1, 2018 · 2 comments

Comments

@jarlebring
Copy link
Member

jarlebring commented Jun 1, 2018 •
edited
Loading

As we have observed several times e.g. JuliaLang/julia#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
Copy link
Member Author

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

@jarlebring
Copy link
Member Author

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

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Assignees
No one assigned
Labels
None yet
Projects
None yet
Milestone
No milestone
Development

No branches or pull requests
1 participant
@jarlebring