Options
Rigorous and portable standard functions
Publikationstyp
Journal Article
Date Issued
2001
Sprache
English
Author(s)
Institut
TORE-URI
Journal
Volume
41
Issue
3
Start Page
540
End Page
562
Citation
BIT Numerical Mathematics 41 (3): 540-562 (2001)
Publisher DOI
Scopus ID
Today's floating point implementations of elementary transcendental functions are usually very accurate. However, with few exceptions, the actual accuracy is not known. In the present paper we describe a rigorous, accurate, fast and portable implementation of the elementary standard functions based on some existing approximate standard functions. The scheme is outlined for IEEE 754, but not difficult to adapt to other floating point formats. A Matlab implementation is available on the net. Accuracy of the proposed algorithms can be rigorously estimated. As an example we prove that the relative accuracy of the exponential function is better than 2.07 eps in a slightly reduced argument range (eps denoting the relative rounding error unit). Otherwise, extensive computational tests suggest for all elementary functions and all suitable arguments an accuracy better than about 3 eps.
Subjects
Accuracy
Elementary standard functions
Floating point implementation
DDC Class
004: Informatik
510: Mathematik