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

Fix eigensolver preconditioning #276

Open
aprokop opened this issue Aug 31, 2020 · 4 comments
Open

Fix eigensolver preconditioning #276

aprokop opened this issue Aug 31, 2020 · 4 comments

Comments

@aprokop
Copy link
Collaborator

aprokop commented Aug 31, 2020

From this comment:

Hey @aprokop , I've got everything working and tested except for a failing anasazi solve in the 4-processor eigen_handle tests (both C++ and Fortran wrapper):

p=0: *** Caught standard std::exception of type 'std::logic_error' :

 ../src/fortrilinos_hl/eigen_handle.cpp:140:

 Throw number = 1

 Throw test that evaluated to true: !(r == 0)

 Error!

The failing line is after

    Anasazi::ReturnType r = solver_->solve();

I tried updating the MueLu solver verbosity to see if that would help shed light on the error, but no luck there.

@aprokop
Copy link
Collaborator Author

aprokop commented Aug 31, 2020

The cause is completely unclear. Running the test without preconditioning works. However, using either Ifpack2 or MueLu does not converge.

@sethrj
Copy link
Collaborator

sethrj commented Sep 1, 2020

Is the fact that it passes on a single processor surprising?

@aprokop
Copy link
Collaborator Author

aprokop commented Sep 1, 2020

Is the fact that it passes on a single processor surprising?

Kinda.

@sethrj
Copy link
Collaborator

sethrj commented Nov 25, 2020

This is at least better-understood and worked around in #282. The check for convergence has been removed in favor of an accessor on the solved values.

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

No branches or pull requests

2 participants