Emil,

I think that TSInterpolate_RK is broken.

>From the manual page of TSRKRegister

http://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/TS/TSRKRegister.html

binterp is supposed to be a matrix but we are passing in the b vector in
TSRKRegisterAll. Also, the NULL case for binterp is not taken into account.

Could you please provide a fix in maint?

Thanks,
-- 
Stefano

Reply via email to