Implementing hybrid explicit/implicit ODE stepper
Thu Jun 4 18:52:00 GMT 2009
On 06/04/2009 08:40 PM, Rhys Ulerich wrote:
> This integrator tends to get used in DNS of Navier-Stokes since the
> memory savings become significant for very large fields.
Ok, now I understand a bit more. I don't think the current GSL
ode-initval is directly suitable for that kind of specialized problem.
The framework does not support separation of f(u) into two components.
You would first have to make changes to gsl_odeiv_system to provide
that, and then modify all the rest accordingly. Might be doable, I don't
Would your application benefit from step size control, or evolution
functionality similar to those in ode-initval? If yes, then you might
try to start from ode-initval and modify it. Otherwise I suggest you
look for how others have designed these kinds of solvers and then create
a new framework suitable for this kind of a problem.
Of course, if you end up implementing only that one stepper, you can use
code snippets from ode-initval for it.
More information about the Gsl-discuss