Skip to main content

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

Project description

Build status PyPI version License

pygslodeiv2 provides a Python binding to the Ordinary Differential Equation integration routines exposed by the odeiv2 interface of GSL - GNU Scientific Library. The odeiv2 interface allows a user to numerically integrate (systems of) differential equations.

The following stepping functions are available:

  • rk2

  • rk4

  • rkf45

  • rkck

  • rk8pd

  • rk1imp

  • rk2imp

  • rk4imp

  • bsimp

  • msadams

  • msbdf

Note that all implicit steppers (those ending with “imp”) and msbdf require a user supplied callback for calculating the jacobian.

Documentation

Autogenerated API documentation is found here: https://bjodah.github.com/pygslodeiv2

Installation

Binary distribution is available here: https://anaconda.org/bjodah/pygslodeiv2

Source distribution is available here: https://pypi.python.org/pypi/pygslodeiv2

Example

The classic van der Pol oscillator (see examples/van_der_pol.py)

>>> import numpy as np
>>> from pygslodeiv2 import integrate_predefined  # also: integrate_adaptive
>>> mu = 1.0
>>> def f(t, y, dydt):
...     dydt[0] = y[1]
...     dydt[1] = -y[0] + mu*y[1]*(1 - y[0]**2)
...
>>> def j(t, y, Jmat, dfdt):
...     Jmat[0, 0] = 0
...     Jmat[0, 1] = 1
...     Jmat[1, 0] = -1 -mu*2*y[1]*y[0]
...     Jmat[1, 1] = mu*(1 - y[0]**2)
...     dfdt[0] = 0
...     dfdt[1] = 0
...
>>> y0 = [1, 0]; dt0=1e-8; t0=0.0; atol=1e-8; rtol=1e-8
>>> tout = np.linspace(0, 10.0, 200)
>>> yout = integrate_predefined(f, j, y0, tout, atol, rtol, dt0,
...                             method='bsimp')  # Implicit Bulirsch-Stoer
>>> import matplotlib.pyplot as plt
>>> plt.plot(tout, yout)
https://raw.githubusercontent.com/bjodah/pygslodeiv2/master/examples/van_der_pol.png

License

The source code is Open Source and is released under GNU GPL v3. See LICENSE for further details. Contributors are welcome to suggest improvements at https://github.com/bjodah/pygslodeiv2

Author

Björn I. Dahlgren, contact:

  • gmail address: bjodah

  • kth.se address: bda

Project details


Download files

Download the file for your platform. If you're not sure which to choose, learn more about installing packages.

Source Distribution

pygslodeiv2-0.3.1.tar.gz (66.8 kB view hashes)

Uploaded source

Supported by

AWS AWS Cloud computing and Security Sponsor Datadog Datadog Monitoring Fastly Fastly CDN Google Google Download Analytics Microsoft Microsoft PSF Sponsor Pingdom Pingdom Monitoring Sentry Sentry Error logging StatusPage StatusPage Status page