-
Notifications
You must be signed in to change notification settings - Fork 36
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
AdditiveSchwarz preconditioner not working in structure module #655
Comments
again, this is more a reminder for myself to fix this. But any hints are appreciated. |
It seems like the update is not done correctly as you guessed. What the error message says is that we try to factorize a matrix which is already factorized. For an update we need to first set the matrices and then factorize them. |
Can you check if https://github.com/bergbauer/exadg/tree/reinit_matrix fixes the problem? |
exactly the fix I have in my branch rn |
The problem is the |
I stumbled upon the issue regarding the
AdditiveSchwarz
preconditioner when looking at #652Issue #652 (steady linear elasticity case) is fixed via #654, but the problem with the AS prec is that after one initial use it becomes invalid.
See the table in
#652 (comment)
and the stacktrace in
#652 (comment)
Probably the update is not done correctly as well here, or the internal state variables are not updated correctly.
@bergbauer could you comment on that guess of mine, I think you implemented that, right?
... I can search for the bug if you agree that something like this could be the reason.
I further assume that in the Poisson module this did not come up, since we do a single solve (similar to the working steady linear case, as seen in #652 (comment)).
The text was updated successfully, but these errors were encountered: