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/ldbl-96/e_hypotl.c | |
| parent | e4a110f59e6317a52dde765ceb3571763ca07a9a (diff) | |
| download | glibc-c88769dda4199d15267259fd3c9b986e6690a901.tar.gz | |
Fix hypot handling of subnormals (bug 16316, bug 16330).
Diffstat (limited to 'sysdeps/ieee754/ldbl-96/e_hypotl.c')
| -rw-r--r-- | sysdeps/ieee754/ldbl-96/e_hypotl.c | 11 |
1 files changed, 11 insertions, 0 deletions
diff --git a/sysdeps/ieee754/ldbl-96/e_hypotl.c b/sysdeps/ieee754/ldbl-96/e_hypotl.c index 7895488848..d3152f91e5 100644 --- a/sysdeps/ieee754/ldbl-96/e_hypotl.c +++ b/sysdeps/ieee754/ldbl-96/e_hypotl.c @@ -89,6 +89,17 @@ long double __ieee754_hypotl(long double x, long double y) b *= t1; a *= t1; k -= 16382; + GET_LDOUBLE_EXP (ea, a); + GET_LDOUBLE_EXP (eb, b); + if (eb > ea) + { + t1 = a; + a = b; + b = t1; + j = ea; + ea = eb; + eb = j; + } } else { /* scale a and b by 2^9600 */ ea += 0x2580; /* a *= 2^9600 */ eb += 0x2580; /* b *= 2^9600 */ |