We have two "splitting" methods in PETSc, additive Runge-Kutta IMEX (a nonlinearly implicit method) and Rosenbrock-W (a more convenient linearly implicit method). Unlike almost all ad-hoc splitting methods, these have high order accuracy and provably good stability properties. Read the man pages and use the IFunction/RHSFunction/IJacobian interface, see src/ts/examples/tutorials/ex22.c and other examples.<br>
<br><div class="gmail_quote">On Mon, Jun 18, 2012 at 1:54 PM, Gaetano Esposito <span dir="ltr"><<a href="mailto:gaetano@email.virginia.edu" target="_blank">gaetano@email.virginia.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
I was wondering whether there is a suggested strategy to implement<br>
something like a timestep splitting method for systems of PDE's (NS)<br>
in Petsc.<br>
By timestep splitting I mean this: advance independently in time each<br>
term of the equations, then add all the solutions to get the total<br>
effects.<br>
>From my basic understanding of the TS object, I could basically define<br>
only one RHS function.<br>
Is creating several TS objects a recommended approach? What are the<br>
downsides of this?<br>
<br>
Many Thanks,<br>
<br>
--g<br>
</blockquote></div><br>