[petsc-dev] [petsc-maint #43201] Questions about TS functionality.

[Jed Brown]

Hi Michael, to elaborate on Barry's comments:

PETSc-dev (and the forthcoming 3.1 release) has support for DAEs using
the interface described here

  http://www.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-dev/docs/manualpages/TS/TSSetIFunction.html

Three of the TS implementations currently support this form:

TSGL: very new general linear methods (A- and L-stable, order 1-5 with full stage order, adaptive step and order)

TSTHETA: classic Theta method, Theta=0.5 is Trapezoid rule (Crank-Nicolson), Theta=1 is Backward Euler.

TSPSEUDO: pseudo-transient continuation


I wrote a section in the User's manual on solving DAE (6.1.2):

  http://www.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-dev/docs/manual.pdf


The controllers for TSGL are pretty naive so it will probably take a
while before they are especially robust.

As Barry says, we plan to make a number of improvements to TS (mostly in
implementation and extensibility) after the release.

Jed