Table-driven implementation of the Expm1 function in IEEE floating-point arithmetic

Abstract

Algorithms and implementation details for the function e x - 1 in both single and double precision of IEEE 754 arithmetic are presented here. With a table of moderate size, the implementations need only working-precision arithmetic and are provably accurate to within 0.58 ulp.