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>
