Computation of Elementary Functions on the IBM RISC System/6000 Processors

The additional speed and precision of the IBM RISC System/6000 floating-point unit have motivated reexamination of algorithms to perform division, square root, and the elementary functions. New results are obtained which avoid the necessity of doing special testing to get the last bit rounded correctly in accordance with all of the IEEE rounding modes in the case of division and square root. For the elementary function library, a technique is described for always getting the last bit rounded correctly in the selected rounding mode.