diff options
author | Joseph Myers <joseph@codesourcery.com> | 2013-12-17 13:42:13 +0000 |
---|---|---|
committer | Joseph Myers <joseph@codesourcery.com> | 2013-12-17 13:42:13 +0000 |
commit | c88769dda4199d15267259fd3c9b986e6690a901 (patch) | |
tree | 743ef394beae012b285cc655faa375293445361f /sysdeps/ieee754/dbl-64/e_hypot.c | |
parent | e4a110f59e6317a52dde765ceb3571763ca07a9a (diff) | |
download | glibc-c88769dda4199d15267259fd3c9b986e6690a901.tar glibc-c88769dda4199d15267259fd3c9b986e6690a901.tar.gz glibc-c88769dda4199d15267259fd3c9b986e6690a901.tar.bz2 glibc-c88769dda4199d15267259fd3c9b986e6690a901.zip |
Fix hypot handling of subnormals (bug 16316, bug 16330).
Diffstat (limited to 'sysdeps/ieee754/dbl-64/e_hypot.c')
-rw-r--r-- | sysdeps/ieee754/dbl-64/e_hypot.c | 11 |
1 files changed, 11 insertions, 0 deletions
diff --git a/sysdeps/ieee754/dbl-64/e_hypot.c b/sysdeps/ieee754/dbl-64/e_hypot.c index 9f4321efe1..500658d076 100644 --- a/sysdeps/ieee754/dbl-64/e_hypot.c +++ b/sysdeps/ieee754/dbl-64/e_hypot.c @@ -102,6 +102,17 @@ __ieee754_hypot (double x, double y) b *= t1; a *= t1; k -= 1022; + GET_HIGH_WORD (ha, a); + GET_HIGH_WORD (hb, b); + if (hb > ha) + { + t1 = a; + a = b; + b = t1; + j = ha; + ha = hb; + hb = j; + } } else /* scale a and b by 2^600 */ { |