Implementing hybrid explicit/implicit ODE stepper

Rhys Ulerich
Wed Jun 3 14:32:00 GMT 2009

Hi all,

I'm looking into implementing a hybrid explicit, implicit ODE
integrator within the GSL ODE framework.
This one splits the right hand side of u_t = f(u) into an implicit
linear operator L(u) and an explicit nonlinear operator N(u).  The
scheme requires L(u) be available in a matrix or matrix-like form.

What's the "right way" to shoehorn this into the gsl_odeiv_system
struct (
L(u) as gsl_odeiv_system.function with N(u) passed within
gsl_odeiv_system.params ?  The reverse?  Neither?


P.S. It's Spalart, Moser, & Rogers 1991 appendix A for anyone interested.

More information about the Gsl-discuss mailing list