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

Michael E Henderson mhender at us.ibm.com
Tue Dec 21 11:48:13 CST 2010


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

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.


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

More information about the petsc-dev mailing list