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

[Development] Should we add a minimum number of iterations for the inexact Newton's method #290

Closed
sblauth opened this issue Aug 21, 2023 · 2 comments
Assignees

Comments

@sblauth
Copy link
Owner

sblauth commented Aug 21, 2023

For large scale problems, one issue with the inexact Newton's method is that the assembly of the Jacobian matrix takes lots of time. In the very early stages (first 2-5 iterations) the inexact Newton's method has a very low tolerance, so that the linear solver converges in about 1-10 iterations (this does not take a long time).

Hence, it might be an idea to let the inexact Newton's method use a minimum number of iterations (say 10-20) regardless of the tolerance, so that the information of the Jacobian matrix can be used more efficiently instead of re-assembling it (with very minor changes) so often, which is expensive.

@sblauth sblauth self-assigned this Aug 21, 2023
@sblauth
Copy link
Owner Author

sblauth commented Oct 26, 2023

This could be implemented, but as additional parameter (which then can be set via the configuration files for optimization problems)

@sblauth
Copy link
Owner Author

sblauth commented Feb 2, 2024

This is available starting with PETSc >= 3.20 - the user only has to specify "ksp_min_it" in the ksp_options for the solve. Hence, it is not necessary to provide a cashocs-based interface for this.

Closing this issue and the corresponding PR #343

@sblauth sblauth closed this as completed Feb 2, 2024
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
1 participant