Rump, Siegfried M.Siegfried M.Rump2021-02-182021-02-182001BIT Numerical Mathematics 41 (3): 540-562 (2001)http://hdl.handle.net/11420/8846Today'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.en1572-9125BIT20013540562AccuracyElementary standard functionsFloating point implementationInformatikMathematikRigorous and portable standard functionsJournal Article10.1023/A:1021971313412Other