<commit_before>========
pyodeint
========

.. image:: http://hera.physchem.kth.se:8080/github.com/bjodah/pyodeint/status.svg?branch=master
   :target: http://hera.physchem.kth.se:8080/github.com/bjodah/pyodeint
   :alt: Build status

``pyodeint`` provides a `Python <http://www.python.org>`_ binding to
`odeint <http://www.odint.com>`_. Currently, only the the 4th order
Rosenbrock stepper has been wrapped and exposed to python (using the
adpative integration routine).

Example
=======
The classic van der Pol oscillator (see `examlpes/van_der_pol.py <examlpes/van_der_pol.py>`_)

::

   >>> from pyodeint import 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
   ...
   >>> y0 = [1; 0]; tend=10.0; dt0=1e-8; t0=0.0; atol=1e-8; rtol=1e-8
   >>> tout, yout = integrate_adaptive(f, j, ny, atol, rtol, y0, t0, tend, dt0)


License
=======
The source code is Open Source and is released under the very permissive
"simplified (2-clause) BSD license". See ``LICENSE.txt`` for further details.<commit_msg>Fix typo in link in README<commit_after>========
pyodeint
========

.. image:: http://hera.physchem.kth.se:8080/github.com/bjodah/pyodeint/status.svg?branch=master
   :target: http://hera.physchem.kth.se:8080/github.com/bjodah/pyodeint
   :alt: Build status

``pyodeint`` provides a `Python <http://www.python.org>`_ binding to
`odeint <http://www.odint.com>`_. Currently, only the the 4th order
Rosenbrock stepper has been wrapped and exposed to python (using the
adpative integration routine).

Example
=======
The classic van der Pol oscillator (see `examples/van_der_pol.py <examples/van_der_pol.py>`_)

::

   >>> from pyodeint import 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
   ...
   >>> y0 = [1; 0]; tend=10.0; dt0=1e-8; t0=0.0; atol=1e-8; rtol=1e-8
   >>> tout, yout = integrate_adaptive(f, j, ny, atol, rtol, y0, t0, tend, dt0)


License
=======
The source code is Open Source and is released under the very permissive
"simplified (2-clause) BSD license". See ``LICENSE.txt`` for further details.
{'fraction_non_alphanumeric': 0.20609803164801235, 'fraction_numerical': 0.04399845619451949, 'mean_word_length': 2.963302752293578, 'pattern_counts': {'":': 0, '<': 9, '<?xml version=': 0, '>': 45, 'https://': 0, 'lorem ipsum': 0, 'www.': 4, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
{'dir': 'git-commits-cleaned', 'id': '2669867', 'n_tokens_mistral': 1117, 'n_tokens_neox': 1038, 'n_words': 258}