Accurate Padé Global Approximations for the Mittag-Leffler Function, Its Inverse, and Its Partial Derivatives to Efficiently Compute Convergent Power Series

Carson Ingo*, Thomas R. Barrick, Andrew G. Webb, Itamar Ronen

*Corresponding author for this work

Research output: Contribution to journalArticlepeer-review

12 Scopus citations

Abstract

The fractional derivative operator provides a mathematical means to concisely describe a heterogeneous and relatively complex system that exhibits non-local, power-law behavior. Discretization of a fractional partial differential equation is non-trivial and computationally intensive. Furthermore, the closed form solution, particularly in the case of the fractional time derivative, comes in the form of a convergent power series of the Mittag-Leffler type, which involves considerable computational burden for accurate and efficient calculation. Here, we extend the Padé global approximation technique in order to accurately compute the Mittag-Leffler function and its inverse operation. Furthermore, we introduce the use of the Padé technique for the partial derivatives of the Mittag-Leffler function, which holistically eliminates the need for the calculation of any convergent power series when fitting to data sets. The Padé approximations provide great utility, particularly in the field of diffusion-weighted magnetic resonance imaging, in which Mittag-Leffler function and its partial derivatives are computed on a voxel-wise basis, with each data set containing on the order of 10 7 voxels. Finally, the Padé approximation codes are provided freely for those interested in modeling and fitting data for processes encompassed by this convergent power series.

Original languageEnglish (US)
Pages (from-to)347-362
Number of pages16
JournalInternational Journal of Applied and Computational Mathematics
Volume3
Issue number2
DOIs
StatePublished - Jun 1 2017

Keywords

  • Anomalous diffusion
  • Fractional derivative
  • Magnetic resonance imaging
  • Mittag-Leffler function
  • Padé approximation
  • Power law dynamics

ASJC Scopus subject areas

  • Computational Mathematics
  • Applied Mathematics

Fingerprint

Dive into the research topics of 'Accurate Padé Global Approximations for the Mittag-Leffler Function, Its Inverse, and Its Partial Derivatives to Efficiently Compute Convergent Power Series'. Together they form a unique fingerprint.

Cite this