pygslodeiv2 package¶

Python binding for odeiv2 in GNU Scientific Library (GSL).

pygslodeiv2.integrate_adaptive(rhs, jac, y0, x0, xend, dx0, atol, rtol, dx_min=0.0, dx_max=0.0, nsteps=500, nderiv=0, check_callable=False, check_indexing=False, **kwargs)[source]¶

Integrates a system of ordinary differential equations.

Parameters:

rhs: callable

Function with signature f(t, y, fout) which modifies fout inplace.

jac: callable

Function with signature j(t, y, jmat_out, dfdx_out) which modifies jmat_out and dfdx_out inplace.

y0: array_like

initial values of the dependent variables

x0: float

initial value of the independent variable

xend: float

stopping value for the independent variable

dx0: float

initial step-size

atol: float

absolute tolerance

rtol: float

relative tolerance

dx_min: float

minimum step (default: 0.0)

dx_max: float

maximum step (default: 0.0)

nsteps: int

maximum number of steps (default: 500)

nderiv: int

number of derivatives (0 or 1) (default: 0)

check_callable: bool (default: False)

perform signature sanity checks on rhs and jac

check_indexing: bool (default: False)

perform item setting sanity checks on rhs and jac.

**kwargs:

‘method’: str

One of: ‘rk2’, ‘rk4’, ‘rkf45’, ‘rkck’, ‘rk8pd’, ‘rk1imp’,

‘rk2imp’, ‘rk4imp’, ‘bsimp’, ‘msadams’, ‘msbdf’

Returns:

(xout, yout, info):

xout: 1-dimensional array of values for the independent variable yout: 2-dimensional array of the dependent variables (axis 1) for

values corresponding to xout (axis 0)

info: dictionary with information about the integration

pygslodeiv2.integrate_predefined(rhs, jac, y0, xout, dx0, atol, rtol, dx_min=0.0, dx_max=0.0, nsteps=500, nderiv=0, check_callable=False, check_indexing=False, **kwargs)[source]¶