On Sat, Feb 11, 2012 at 10:58, Konstantinos Kontzialis <ckontzialis@lycos.com> wrote:

Dear Jed,

I use TSSolve() and version 3.2-p5.

The loop in TSSolve() is simple. Maybe you can use a debugger to find out what happened to post-step function pointer.

while (!ts->reason) {

ierr = TSPreStep(ts);CHKERRQ(ierr);

ierr = TSStep(ts);CHKERRQ(ierr);

ierr = TSPostStep(ts);CHKERRQ(ierr);

ierr = TSMonitor(ts,ts->steps,ts->ptime,ts->vec_sol);CHKERRQ(ierr);

}