diff options
| author | Joseph Myers <joseph@codesourcery.com> | 2013-05-24 20:52:55 +0000 |
|---|---|---|
| committer | Joseph Myers <joseph@codesourcery.com> | 2013-05-24 20:52:55 +0000 |
| commit | 0323d08657f111267efa47bd448fbf6cd76befe8 (patch) | |
| tree | 7ebd164c7460e42866c7fe5e31ab9c4b556f23ab | |
| parent | dd4259b9f76fa8fd3690a91608d2e3a94e2d6713 (diff) | |
| download | glibc-0323d08657f111267efa47bd448fbf6cd76befe8.tar.gz | |
Fix ldbl-96 hypotl of subnormals (bug 15529).
| -rw-r--r-- | ChangeLog | 5 | ||||
| -rw-r--r-- | NEWS | 3 | ||||
| -rw-r--r-- | math/libm-test.inc | 13 | ||||
| -rw-r--r-- | sysdeps/ieee754/ldbl-96/e_hypotl.c | 2 |
4 files changed, 21 insertions, 2 deletions
@@ -1,5 +1,10 @@ 2013-05-24 Joseph Myers <joseph@codesourcery.com> + [BZ #15529] + * sysdeps/ieee754/ldbl-96/e_hypotl.c (__ieee754_hypotl): Set high + bit of mantissa of 2^16382. + * math/libm-test.inc (hypot_test_data): Add more tests. + * math/libm-test.inc: Add drem and pow10 to list of tested functions. (pow10_test): New function. @@ -18,7 +18,8 @@ Version 2.18 15287, 15304, 15305, 15307, 15309, 15327, 15330, 15335, 15336, 15337, 15339, 15342, 15346, 15359, 15361, 15366, 15380, 15381, 15394, 15395, 15405, 15406, 15409, 15416, 15418, 15419, 15423, 15424, 15426, 15429, - 15441, 15442, 15448, 15480, 15485, 15488, 15490, 15493, 15497, 15506. + 15441, 15442, 15448, 15480, 15485, 15488, 15490, 15493, 15497, 15506, + 15529. * CVE-2013-0242 Buffer overrun in regexp matcher has been fixed (Bugzilla #15078). diff --git a/math/libm-test.inc b/math/libm-test.inc index 91af668e94..7a6bf09a66 100644 --- a/math/libm-test.inc +++ b/math/libm-test.inc @@ -9387,6 +9387,19 @@ static const struct test_ff_f_data hypot_test_data[] = #if defined TEST_LDOUBLE && LDBL_MAX_EXP >= 16384 && !defined TEST_INLINE TEST_ff_f (hypot, 0x3p16381L, 0x4p16381L, 0x5p16381L), #endif + + TEST_ff_f (hypot, 0x1p-149L, 0x1p-149L, 1.9817352931807469938024533350782879785095e-45L, UNDERFLOW_EXCEPTION_FLOAT), + +#ifndef TEST_FLOAT + TEST_ff_f (hypot, 0x1p-1074L, 0x1p-1074L, 6.9871433705131320800651344656990806305791e-324L, UNDERFLOW_EXCEPTION_DOUBLE), +#endif + +#if defined TEST_LDOUBLE && LDBL_MIN_EXP <= -16381 && !defined TEST_INLINE + TEST_ff_f (hypot, 0x1p-16445L, 0x1p-16445L, 5.1550906155442528702558159159596215039925e-4951L, UNDERFLOW_EXCEPTION), +# if LDBL_MANT_DIG >= 113 + TEST_ff_f (hypot, 0x1p-16494L, 0x1p-16494L, 9.1572804726500807075521065242888978445857e-4966L, UNDERFLOW_EXCEPTION), +# endif +#endif }; static void diff --git a/sysdeps/ieee754/ldbl-96/e_hypotl.c b/sysdeps/ieee754/ldbl-96/e_hypotl.c index 306f92924c..7895488848 100644 --- a/sysdeps/ieee754/ldbl-96/e_hypotl.c +++ b/sysdeps/ieee754/ldbl-96/e_hypotl.c @@ -85,7 +85,7 @@ long double __ieee754_hypotl(long double x, long double y) u_int32_t high,low; GET_LDOUBLE_WORDS(exp,high,low,b); if((high|low)==0) return a; - SET_LDOUBLE_WORDS(t1, 0x7ffd, 0, 0); /* t1=2^16382 */ + SET_LDOUBLE_WORDS(t1, 0x7ffd, 0x80000000, 0); /* t1=2^16382 */ b *= t1; a *= t1; k -= 16382; |