<br><font size=2 face="sans-serif">Barry, Jed,</font>
<br>
<br><font size=2 face="sans-serif">Thanks. That's good to hear. I'll give
it a look.</font>
<br>
<br><font size=2 face="sans-serif">Mike<br>
------------------------------------------------------------------------------------------------------------------------------------<br>
Mathematical Sciences, TJ Watson Research Center<br>
mhender@watson.ibm.com<br>
</font><a href=http://www.research.ibm.com/people/h/henderson/><font size=2 face="sans-serif">http://www.research.ibm.com/people/h/henderson/</font></a><font size=2 face="sans-serif"><br>
</font><a href=http://multifario.sourceforge.net/><font size=2 face="sans-serif">http://multifario.sourceforge.net/</font></a><font size=2 face="sans-serif"><br>
</font>
<br>
<br>
<br>
<table width=100%>
<tr valign=top>
<td><font size=1 color=#5f5f5f face="sans-serif">From:</font>
<td><font size=1 face="sans-serif">Jed Brown <jed@59A2.org></font>
<tr valign=top>
<td><font size=1 color=#5f5f5f face="sans-serif">To:</font>
<td><font size=1 face="sans-serif">Barry Smith <bsmith@mcs.anl.gov>,
petsc-maint@mcs.anl.gov, Michael E Henderson/Watson/IBM@IBMUS</font>
<tr>
<td valign=top><font size=1 color=#5f5f5f face="sans-serif">Cc:</font>
<td><font size=1 face="sans-serif">For users of the development version
of PETSc <petsc-dev@mcs.anl.gov></font>
<tr valign=top>
<td><font size=1 color=#5f5f5f face="sans-serif">Date:</font>
<td><font size=1 face="sans-serif">03/15/10 03:31 PM</font>
<tr valign=top>
<td><font size=1 color=#5f5f5f face="sans-serif">Subject:</font>
<td><font size=1 face="sans-serif">Re: [petsc-dev] [petsc-maint #43201]
Questions about TS functionality.</font>
<tr valign=top>
<td><font size=1 color=#5f5f5f face="sans-serif">Sent by:</font>
<td><font size=1 face="sans-serif">Jed Brown <five9a2@gmail.com></font></table>
<br>
<hr noshade>
<br>
<br>
<br><tt><font size=2>Hi Michael, to elaborate on Barry's comments:<br>
<br>
PETSc-dev (and the forthcoming 3.1 release) has support for DAEs using<br>
the interface described here<br>
<br>
  </font></tt><a href="http://www.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-dev/docs/manualpages/TS/TSSetIFunction.html"><tt><font size=2>http://www.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-dev/docs/manualpages/TS/TSSetIFunction.html</font></tt></a><tt><font size=2><br>
<br>
Three of the TS implementations currently support this form:<br>
<br>
TSGL: very new general linear methods (A- and L-stable, order 1-5 with
full stage order, adaptive step and order)<br>
<br>
TSTHETA: classic Theta method, Theta=0.5 is Trapezoid rule (Crank-Nicolson),
Theta=1 is Backward Euler.<br>
<br>
TSPSEUDO: pseudo-transient continuation<br>
<br>
<br>
I wrote a section in the User's manual on solving DAE (6.1.2):<br>
<br>
  </font></tt><a href="http://www.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-dev/docs/manual.pdf"><tt><font size=2>http://www.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-dev/docs/manual.pdf</font></tt></a><tt><font size=2><br>
<br>
<br>
The controllers for TSGL are pretty naive so it will probably take a<br>
while before they are especially robust.<br>
<br>
As Barry says, we plan to make a number of improvements to TS (mostly in<br>
implementation and extensibility) after the release.<br>
<br>
Jed<br>
</font></tt>
<br>
<br>