TPSA restructure (3/3): Newton defaults for linear elasticity - #7359
Draft
svenn-t wants to merge 3 commits into
Draft
TPSA restructure (3/3): Newton defaults for linear elasticity#7359svenn-t wants to merge 3 commits into
svenn-t wants to merge 3 commits into
Conversation
- Old 7x7 block structure is broken up to several sub-matrices, purpose-made for preconditioners like Hypre. - New structure made to interface 7x7 structure in linearizer without changes.
- One preconditioner for whole TPSA system; three sub-block preconditioner can be set by user - Allows for Hypre, which is default if linked; else ilu0
- newton procedure now: check linearization, if error larger than tolerance, solve linear system - lowered linear solve tolerance - more clear fixed-stress conv. check - bug fix: num. dofs on owned cells only - moved calc. of info to mech. output to separate function, which is called at endTimeStep() - option to scale linear system for better conv. check with dune solvers (preconditioners are invariant)
Contributor
Author
|
jenkins build this please |
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
Add this suggestion to a batch that can be applied as a single commit.This suggestion is invalid because no changes were made to the code.Suggestions cannot be applied while the pull request is closed.Suggestions cannot be applied while viewing a subset of changes.Only one suggestion per line can be applied in a batch.Add this suggestion to a batch that can be applied as a single commit.Applying suggestions on deleted lines is not supported.You must change the existing code in this line in order to create a valid suggestion.Outdated suggestions cannot be applied.This suggestion has been applied or marked resolved.Suggestions cannot be applied from pending reviews.Suggestions cannot be applied on multi-line comments.Suggestions cannot be applied while the pull request is queued to merge.Suggestion cannot be applied right now. Please check back later.
Change Newton to default to 1 iteration, to speed up linear elasticity solves with TPSA. Now we only do one linearization, check the error, and proceed with one linear solve if linearization error is too high. Increased linear solver tolerance to compensate for only doing 1 Newton iteration.
Also includes some minor changes and fixes:
Includes both PRs #7357 and #7358, which should be merged first. Draft mode until then.