On Sat, Feb 11, 2012 at 10:58, Konstantinos Kontzialis < [email protected]> 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); }