Skip to main content

Collection of algorithms for numerically calculating fractional derivatives.

Project description

differint

This package is used for numerically calculating fractional derivatives and integrals (differintegrals). Options for varying definitions of the differintegral are available, including the Grunwald-Letnikov (GL), the 'improved' Grunwald-Letnikov (GLI), the Riemann-Liouville (RL), and the Caputo (L1, L2, and L2C). Through the API, you can compute differintegrals at a point or over an array of function values.

See below for an example of how to use this package, or check out the wiki for references, signatures, and examples for each function.

Motivation

There is little in the way of readily available, easy-to-use code for numerical fractional calculus. What is currently available are functions that are generally either smart parts of a much larger package, or only offer one numerical algorithm. The differint package offers a variety of algorithms for computing differintegrals and several auxiliary functions relating to generalized binomial coefficients.

Installation

This project requires Python 3+ and NumPy to run.

Installation from the Python Packaging index (https://pypi.python.org/pypi) is simple using pip.

pip install differint

Included Files

Core File Description
differint/differint.py Contains algorithms for fractional differentiation and integration.
tests/test.py Testing suite containing all unit tests.

Both of the above files have corresponding __init__.py files.

Setup File Description
.gitignore List of files to ignore during git push/pull requests.
CONTRIBUTING.md Instructions for potential contributors to the differint project.
LICENSE MIT license agreement.
MANIFEST.in Selects the README file for uploading to PyPi.
README.md This README file.
README.rst This README file in ReStructuredText format.
init.py __init__ file for overall package.
changelog.txt List of updates to package.
setup.py Script for downloading package from pip.

Example Usage

Taking a fractional derivative is easy with the differint package. Let's take the 1/2 derivative of the square root function on the interval [0,1], using the Riemann-Liouville definition of the fractional derivative.

import numpy as np
import differint.differint as df

def f(x):
   return x**0.5

DF = df.RL(0.5, f)
print(DF)

You can also specify the endpoints of the domain and the number of points used as follows.

DF = df.RL(0.5, f, 0, 1, 128)

For a description of all functions, their signatures, and more usage examples, see the project's wiki.

Tests

All tests can be run with nose from the command line. Setup will automatically install nose if it is not present on your machine.

python setup.py tests

Alternatively, you can run the test script directly.

cd <file_path>/differint/tests/
python test.py

API Reference

In this section we cover the usage of the various functions within the differint package.

Main Function Usage
GLpoint Computes the GL differintegral at a point
GL Computes the GL differintegral over an entire array of function values using the Fast Fourier Transform
GLI Computes the improved GL differintegral over an entire array of function values
CRONE Calculates the GL derivative approximation using the CRONE operator.
RLpoint Computes the RL differintegral at a point
RL Computes the RL differintegral over an entire array of function values using matrix methods
CaputoL1point Computes the Caputo differintegral at a point using the L1 algorithm
CaputoL2point Computes the Caputo differintegral at a point using the L2 algorithm
CaputoL2Cpoint Computes the Caputo differintegral at a point using the L2C algorithm
PCsolver Solves IVPs for fractional ODEs of the form ${}^CD^\alpha[y(x)]=f(x,y(x))$ using the predictor-corrector method
Auxiliary Function Usage
isInteger Determine if a number is an integer
isPositiveInteger Determine if a number is an integer, and if it is greater than 0
checkValues Used to check for valid algorithm input types
GLIinterpolat Define interpolating coefficients for the improved GL algorithm
functionCheck Determines if algorithm function input is callable or an array of numbers
poch Computes the Pochhammer symbol
Gamma Computes the gamma function, an extension of the factorial to complex numbers
Beta Computes the beta function, a function related to the binomial coefficient
MittagLeffler Computes the two parameter Mittag-Leffler function, which is important in the solution of fractional ODEs
GLcoeffs Determines the convolution filter composed of generalized binomial coefficients used in the GL algorithm
RLcoeffs Calculates the coefficients used in the RLpoint and RL algorithms
RLmatrix Determines the matrix used in the RL algorithm
PCcoeffs Determines the coefficients used in the PC algorithm

Contribute

To contribute to this project, see the contributing guidelines.

Credits

Baleanu, D., Diethelm, K., Scalas, E., & Trujillo, J.J. (2012). Fractional Calculus: Models and Numerical Methods. World Scientific.

Oldham, K.B. & Spanier, J. (1974). The Fractional Calculus: Theory and Applications of Differentiation and Integration to Arbitrary Order. Academic Press Inc.

Karniadakis, G.E.. (2019). Handbook of Fractional Calculus with Applications Volume 3: Numerical Methods. De Gruyter.

License

MIT © Matthew Adams

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

differint-1.0.0.tar.gz (16.9 kB view details)

Uploaded Source

Built Distribution

differint-1.0.0-py3-none-any.whl (11.8 kB view details)

Uploaded Python 3

File details

Details for the file differint-1.0.0.tar.gz.

File metadata

  • Download URL: differint-1.0.0.tar.gz
  • Upload date:
  • Size: 16.9 kB
  • Tags: Source
  • Uploaded using Trusted Publishing? No
  • Uploaded via: twine/4.0.2 CPython/3.9.16

File hashes

Hashes for differint-1.0.0.tar.gz
Algorithm Hash digest
SHA256 dc4e00bc54bdc0eeea864da496e5a3c08bf9246da1882f885846a29e70d5c449
MD5 680f1ea52283d9852689861b1e90c714
BLAKE2b-256 483595cfa80aba1966dc365377e841584420ac10f8cd49d29e0257ca7255ff80

See more details on using hashes here.

File details

Details for the file differint-1.0.0-py3-none-any.whl.

File metadata

  • Download URL: differint-1.0.0-py3-none-any.whl
  • Upload date:
  • Size: 11.8 kB
  • Tags: Python 3
  • Uploaded using Trusted Publishing? No
  • Uploaded via: twine/4.0.2 CPython/3.9.16

File hashes

Hashes for differint-1.0.0-py3-none-any.whl
Algorithm Hash digest
SHA256 f15361a394a4b10098914b705680094597d498e2ddffc4f429539a0006da7c11
MD5 e3bbf36efdb07760d97671de61778dde
BLAKE2b-256 71ae8a79f26573f9868d6275ea975634e709c6a683522176841d98a4ccfc154b

See more details on using hashes here.

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