Precise and fast computation of inverse Fermi–Dirac integral of order 1/2 by minimax rational function approximation

Toshio Fukushima

Applied Mathematics and Computation, 2015, vol. 259, issue C, 698-707

Abstract: The single and double precision procedures are developed for the inverse Fermi–Dirac integral of order 1/2 by the minimax rational function approximation. The maximum error of the new approximations is one and 7 machine epsilons in the single and double precision computations, respectively. Meanwhile, the CPU time of the new approximations is so small as to be comparable to that of elementary functions. As a result, the new double precision approximation achieves the 15 digit accuracy and runs 30–84% faster than Antia’s 28 bit precision approximation (Antia, 1993). Also, the new single precision approximation is of the 24 bit accuracy and runs 10–86% faster than Antia’s 15 bit precision approximation.

Keywords: Fermi–Dirac integral; Function approximation; Inverse Fermi–Dirac integral; Minimax approximation; Rational function approximation (search for similar items in EconPapers)
Date: 2015
References: View complete reference list from CitEc
Citations:

Downloads: (external link)
http://www.sciencedirect.com/science/article/pii/S0096300315003094
Full text for ScienceDirect subscribers only

Related works:
This item may be available elsewhere in EconPapers: Search for items with the same title.

Export reference: BibTeX RIS (EndNote, ProCite, RefMan) HTML/Text

Persistent link: https://EconPapers.repec.org/RePEc:eee:apmaco:v:259:y:2015:i:c:p:698-707

DOI: 10.1016/j.amc.2015.03.015

Access Statistics for this article

Applied Mathematics and Computation is currently edited by Theodore Simos

More articles in Applied Mathematics and Computation from Elsevier
Bibliographic data for series maintained by Catherine Liu ().