[petsc-dev] Fw: TSTheta for DAE's

[Michael E Henderson]

Hi,

Who can I talk to about the implementation of the TSTheta Implicit time 
stepper?

The implementation computes ts->vec_sol which solves the implicit system. 
But the monitor returns ts->vec_sol+dt*th->Xdot, which does not 
necessarily satisfy the equations. At least that's how I read the code. I 
think that the call to VecAXPY which changes ts->vec_sol needs to be moved 
to after the call to TSMonitor in src/ts/impls/implicit/theta/theta.c 
Purely a matter of what is passed to the monitor routine.

Thanks,

Mike Henderson
------------------------------------------------------------------------------------------------------------------------------------
Mathematical Sciences, TJ Watson Research Center
mhender at watson.ibm.com
http://www.research.ibm.com/people/h/henderson/
http://multifario.sourceforge.net/
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20101221/7032add1/attachment.html>