<font size=2 face="sans-serif">Hi,</font>
<br>
<br><font size=2 face="sans-serif">Who can I talk to about the implementation
of the TSTheta Implicit time stepper?</font>
<br>
<br><font size=2 face="sans-serif">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.</font>
<br>
<br><font size=2 face="sans-serif">Thanks,</font>
<br>
<br><font size=2 face="sans-serif">Mike Henderson<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>