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

Do not create unnecessary failing assembler for compressible Newton models with full A block preconditioner #4465

Merged

Conversation

gassmoeller
Copy link
Member

This was a bug in debug mode I noticed in one of @elodie-kendall's models. If we set up a compressible model with a Newton solver and use the full A block of the Stokes matrix for the preconditioner, then the NewtonHandler still adds an assembler for the compressible terms in the preconditioner, although the preconditioner matrix is not used (the full matrix is used instead). Additionally the added assembler will throw an assertion. This change unifies the code paths between the Newton solver and the regular solver (which does not add the assembler for these models).

@gassmoeller gassmoeller changed the title Do not create unnecessary failing assembler. Do not create unnecessary failing assembler for compressible Newton models with full A block preconditioner Jan 27, 2022
@gassmoeller gassmoeller force-pushed the avoid_creating_unnecessary_assembler branch from 43a9323 to 709ffe8 Compare January 28, 2022 16:53
Copy link
Contributor

@jdannberg jdannberg left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

This looks good to me. Thanks for the fix!
Since it's apparently a bugfix, do you want to add a changelog entry?

@gassmoeller
Copy link
Member Author

Changelog added.

@jdannberg jdannberg merged commit 748a11a into geodynamics:main Feb 3, 2022
@gassmoeller gassmoeller deleted the avoid_creating_unnecessary_assembler branch February 3, 2022 17:18
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

Successfully merging this pull request may close these issues.

None yet

2 participants