sblauth / cashocs

computational adjoint-based shape optimization and optimal control software for python
https://cashocs.readthedocs.io/
GNU General Public License v3.0
50 stars 10 forks source link

Pseudo-transient continuation solver for globalizing Newton's method #229

Open sblauth opened 1 year ago

sblauth commented 1 year ago

Pseudo-Transient Continuation seems to be a robust way of globalizing Newton's method. It should be possible to implement this for general problems in cashocs.

One note: Users should have the option (optional!) to supply their own time derivative. This can be very useful when using mixed problems, where only one variable is truly time-dependent. This is analogous to the DAE formulation of Pseudo-Transient-Continuation. Another possible implementation would be for the user to indicate which parts of the (mixed) solution should receive time-derivatives.

sblauth commented 3 months ago

Just a comment: This is already implemented with the PETSc TS module - maybe we can use this? This is in line with the discussion of #441