Abstract
Special functions are pervasive in all fields of science. The most well-known application areas are in physics, engineering, chemistry, computer science and statistics. Because of their importance, several books and a large collection of papers have been devoted to the numerical computation of these functions. But up to this date, even environments such as Maple, Mathematica, MATLAB and libraries such as IMSL, CERN and NAG offer no routines for the reliable evaluation of special functions. Here the notion reliable indicates that, together with the function evaluation a guaranteed upper bound on the total error or, equivalently, an enclosure for the exact result, is computed.
We point out how limit-periodic continued fraction representations of these functions can be helpful in this respect. The newly developed (and implemented) scalable precision technique is mainly based on the use of sharpened a priori truncation error and round-off error upper bounds for real continued fraction representations of special functions of a real variable. The implementation is reliable in the sense that it returns a sharp interval enclosure for the requested function evaluation, at the same cost as the evaluation.
Access this chapter
Tax calculation will be finalised at checkout
Purchases are for personal use only
Preview
Unable to display preview. Download preview PDF.
References
Colman, M., Cuyt, A.: Gauss and confluent hypergeometric functions accurate to the last digit. ACM TOMS (in preparation, 2006)
Cuyt, A., Becuwe, S.: Reliable software for the evaluation of several special functions. ACM TOMS (in preparation, 2006)
Cuyt, A., Verdonk, B.: Computer arithmetic: basic theory. SIAM, Philadelphia (in preparation, 2006)
Cuyt, A., Verdonk, B., Waadeland, H.: Efficient and reliable multiprecision implementation of elementary and special functions. SIAM Journal on Scientific Computing (to appear, 2006)
Gautschi, W.: Computational aspects of three-term recurrence relations. SIAM Rev. 9, 24–82 (1987)
Jones, W.B., Thron, W.J.: Numerical stability in evaluating continued fractions. Math. Comp. 28, 795–810 (1974)
Lorentzen, L., Waadeland, H.: Continued fractions with applications. North-Holland Publishing Company, Amsterdam (1992)
Thron, W.J., Waadeland, H.: Accelerating convergence of limit periodic continued fractions K(a n /1). Numer. Math. 34, 155–170 (1980)
Author information
Authors and Affiliations
Editor information
Editors and Affiliations
Rights and permissions
Copyright information
© 2006 Springer-Verlag Berlin Heidelberg
About this paper
Cite this paper
Cuyt, A., Becuwe, S. (2006). Towards Reliable Software for the Evaluation of a Class of Special Functions. In: Iglesias, A., Takayama, N. (eds) Mathematical Software - ICMS 2006. ICMS 2006. Lecture Notes in Computer Science, vol 4151. Springer, Berlin, Heidelberg. https://doi.org/10.1007/11832225_31
Download citation
DOI: https://doi.org/10.1007/11832225_31
Publisher Name: Springer, Berlin, Heidelberg
Print ISBN: 978-3-540-38084-9
Online ISBN: 978-3-540-38086-3
eBook Packages: Computer ScienceComputer Science (R0)