Publications Meetings The Profession Membership Programs Math Samplings Policy & Advocacy In the News About the AMS
   
Mobile Device Pairing
Green Open Access
Mathematics of Computation
Mathematics of Computation
ISSN 1088-6842(online) ISSN 0025-5718(print)

 

Sharp ULP rounding error bound
for the hypotenuse function


Author: Abraham Ziv
Journal: Math. Comp. 68 (1999), 1143-1148
MSC (1991): Primary 65G05; Secondary 65D20
Published electronically: February 13, 1999
MathSciNet review: 1648423
Full-text PDF Free Access

Abstract | References | Similar Articles | Additional Information

Abstract: The hypotenuse function, $z=\sqrt{x^2+y^2}$, is sometimes included in math library packages. Assuming that it is being computed by a straightforward algorithm, in a binary floating point environment, with round to nearest rounding mode, a sharp roundoff error bound is derived, for arbitrary precision. For IEEE single precision, or higher, the bound implies that $|\overline z-z|<1.222\, ulp(z)$ and $|\overline z-z|<1.222\, ulp(\overline z)$. Numerical experiments indicate that this bound is sharp and cannot be improved.


References [Enhancements On Off] (What's this?)


Similar Articles

Retrieve articles in Mathematics of Computation of the American Mathematical Society with MSC (1991): 65G05, 65D20

Retrieve articles in all journals with MSC (1991): 65G05, 65D20


Additional Information

Abraham Ziv
Affiliation: IBM Israel, Science and Technology, Matam–Advanced Technology Center, Haifa 31905, Israel
Email: ziv@haifasc3.vnet.ibm.com

DOI: http://dx.doi.org/10.1090/S0025-5718-99-01103-5
PII: S 0025-5718(99)01103-5
Keywords: Rounding error, error analysis, relative error, error bound, floating point, ULP, hypotenuse function, math library
Received by editor(s): December 1, 1997
Published electronically: February 13, 1999
Article copyright: © Copyright 1999 American Mathematical Society