On Jan 26, 2010, at 1:30 PM, Jed Brown wrote:
> Note that TSGL is just like any other TS implementation, your code
> doesn't need to make any explicit reference to it, just set up a TS,
> call TSSolve, and then it can be used via -ts_type gl. Admittedly, many
> of the TS implementations are pretty naive, though -ts_type sundials is
> robust, -ts_type pseudo is useful for steady-state computations, and
> -ts_type ssp provides several optimal strong stability preserving
> schemes (not really useful to people doing continuous Galerkin, perhaps
> important with DG).
Exactly... shouldn't be hard to adapt the current SNES wrapper... only real difference is that the residual and jacobian calls pass t and u_t in addition to the stuff SNES normally passes.