<div class="gmail_quote">On Sat, Feb 11, 2012 at 10:58, Konstantinos Kontzialis <span dir="ltr">&lt;<a href="mailto:ckontzialis@lycos.com">ckontzialis@lycos.com</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Dear Jed,<br>
<br>
 I use TSSolve() and version 3.2-p5.<br></blockquote><div><br></div><div>The loop in TSSolve() is simple. Maybe you can use a debugger to find out what happened to post-step function pointer.</div><div><br></div><div><div>
    while (!ts-&gt;reason) {</div><div>      ierr = TSPreStep(ts);CHKERRQ(ierr);</div><div>      ierr = TSStep(ts);CHKERRQ(ierr);</div><div>      ierr = TSPostStep(ts);CHKERRQ(ierr);</div><div>      ierr = TSMonitor(ts,ts-&gt;steps,ts-&gt;ptime,ts-&gt;vec_sol);CHKERRQ(ierr);</div>
<div>    }</div></div></div>