Contents
It is a numerical integration technique. From MathWorld, Ueberhuber (1997, p. 71) and Krommer and Ueberhuber (1998, pp. 49 and 155-165) use the word "quadrature" to mean numerical computation of a univariate integral, and "cubature" to mean numerical computation of a multiple integral.
Please, see the module documentation here.
From the Nanostructures and Computation Wiki at MIT, Steven W. Johnson has written a simple C package for adaptive multidimensional integration (cubature) of vector-valued functions over hypercubes and this is a Python wrapper for the referred C package.
If you are changing _cubature.pyx
, you must have Cython installed in order
to create a new _cubature.c
file (the same is valid for
_test_integrands.pyx
). The setup.py
script will automatically try to
use the Cython compiler first.
To install in the site-packages
directory and make it importable from
anywhere:
python setup.py install
Windows:
set prefix=anydirectory set PYTHONPATH=%anydirectory%\Lib\site-packages;%PYTHONPATH% mkdir %anydirectory% python setup.py install --prefix="%anydirectory%"
Linux:
export prefix=anydirectory export PYTHONPATH=$anydirectory/Lib/site-packages;$PYTHONPATH mkdir $anydirectory python setup.py install --prefix=$anydirectory
It will create an .egg
file that will go to
$anydirectory\Lib\site-packages
. This file can be unzipped to obtain the
importable module, OR, the full path to this .egg
can be added to
$PYTHONPATH
(which can be also done inside a script through
sys.path.append().
If you want to build it locally (without installing in Python's
site-packages
) just type:
python setup.py build_ext --inplace clean
The Python wrapper has been proven using test integrands from the C package and some additional testing functions from Genz. The integrands were implemented in Cython and verified with Mathematica.
After building cubature, run the unit tests with the `pytest`
package in
the package directory. Be aware that this takes several minutes:
python -m pytest test_cubature.py
Some examples are given in ./examples.
You are welcome to fork this repository and modify it in whatever way you want. It will also be nice if you could send a pull request here in case you think your modifications are valuable for another person.
This wrapper follows the GNU-GPL license terms discribed in the C Package.