Validated evaluation of special mathematical functions

Because of the importance of special functions, several books and a large collection of papers have been devoted to their use and computation, the most well-known being the Abramowitz and Stegun handbook (Abramowitz and Stegun, 1964) 1] and its successor (Olver et al. 0000) 2]. However, until now no environment offers routines for the provable correct multiprecision and radix-independent evaluation of these special functions. We point out how we make good use of series and limit-periodic continued fraction representations in a package that is being developed at the University of Antwerp. Our scalable precision technique is mainly based on the use of sharpened a priori truncation and round-off error upper bounds for real arguments. The implementation is validated in the sense that it returns a sharp interval enclosure for the requested function evaluation, at the same cost as the evaluation.

[1]  David H. Bailey,et al.  A Fortran 90-based multiprecision system , 1995, TOMS.

[2]  Annie A. M. Cuyt,et al.  Validated computation of certain hypergeometric functions , 2012, TOMS.

[3]  Walter Gautschi,et al.  NUMERICAL EVALUATION OF SPECIAL FUNCTIONS , 2001 .

[4]  Jean-Michel Muller,et al.  Elementary Functions: Algorithms and Implementation , 1997 .

[5]  Richard P. Brent,et al.  Recent technical reports , 1977, SIGA.

[6]  Ansi Ieee,et al.  IEEE Standard for Binary Floating Point Arithmetic , 1985 .

[7]  Annie A. M. Cuyt,et al.  Handbook of Continued Fractions for Special Functions , 2008 .

[8]  Daniel W. Lozier,et al.  NIST Digital Library of Mathematical Functions , 2003, Annals of Mathematics and Artificial Intelligence.

[9]  S. Moshier,et al.  Computer approximations , 1986 .

[10]  Annie A. M. Cuyt,et al.  Efficient and Reliable Multiprecision Implementation of Elementary and Special Functions , 2006, SIAM J. Sci. Comput..

[11]  Haakon Waadeland,et al.  Continued fractions with applications , 1994 .

[12]  Guido D. Salvucci,et al.  Ieee standard for binary floating-point arithmetic , 1985 .

[13]  M. Abramowitz,et al.  Handbook of Mathematical Functions With Formulas, Graphs and Mathematical Tables (National Bureau of Standards Applied Mathematics Series No. 55) , 1965 .

[14]  William B. Jones,et al.  Numerical stability in evaluating continued fractions , 1974 .

[15]  John Gill Infinite compositions of Möbius transformations , 1973 .

[16]  D. Owen Handbook of Mathematical Functions with Formulas , 1965 .

[17]  David M. Smith,et al.  Algorithm 786: multiple-precision complex arithmetic and functions , 1998, TOMS.

[18]  Lisa Lorentzen A Priori Truncation Error Bounds for Continued Fractions , 2003 .

[19]  Ronald F. Boisvert,et al.  NIST Handbook of Mathematical Functions , 2010 .

[20]  Jackson B. Lackey,et al.  Errata: Handbook of mathematical functions with formulas, graphs, and mathematical tables (Superintendent of Documents, U. S. Government Printing Office, Washington, D. C., 1964) by Milton Abramowitz and Irene A. Stegun , 1977 .

[21]  Mei Han An,et al.  accuracy and stability of numerical algorithms , 1991 .

[22]  Vincent Lefèvre,et al.  MPFR: A multiple-precision binary floating-point library with correct rounding , 2007, TOMS.