[petsc-users] Question about how to solve the DAE step by step using PETSC4PY.

[Xiong, Jing]

Good afternoon,

Thanks for all your help.
I got a question about how to solve the DAE step by step using PETSC4PY.
I got two DAE systems, let's say f1(x1, u1) and f2(x2, u2). The simplified algorithm I need to implement is as follows:
Step 1: solve f1 for 1 step.
Step 2: use part of x1 as input for f2, which means u2 = x1[:n].
Step 3: solve f2 for one step with u2 = x1[:n].
Step 4: use part of x2 as input for f1, which means u1 = x2[:m].

I'm not able to find any examples of how to use PETSC4PY in such scenarios. If using the "scikits.odes.dae" package, it is like:
daesolver.init_step(timeInit, XInit, XpInit)
daesolver.step(time)

Please let me know if there are any examples I can refer to. Thank you.

Best,
Jing
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20220325/986b3754/attachment.html>