From e6b29af8a314d3a6a538222a98383c215d410ddf Mon Sep 17 00:00:00 2001 From: Ulrich Drepper Date: Sun, 28 Oct 2007 01:14:40 +0000 Subject: [BZ #2549] 2006-04-14 H.J. Lu [BZ #2549] * math/libm-test.inc (check_float_internal): Support denormalized return. --- math/libm-test.inc | 23 ++++++++++++++++++----- 1 file changed, 18 insertions(+), 5 deletions(-) (limited to 'math') diff --git a/math/libm-test.inc b/math/libm-test.inc index da68d40bb4..e698923f4e 100644 --- a/math/libm-test.inc +++ b/math/libm-test.inc @@ -487,11 +487,24 @@ check_float_internal (const char *test_name, FLOAT computed, FLOAT expected, else { diff = FUNC(fabs) (computed - expected); - /* ilogb (0) isn't allowed. */ - if (expected == 0.0) - ulp = diff / FUNC(ldexp) (1.0, - MANT_DIG); - else - ulp = diff / FUNC(ldexp) (1.0, FUNC(ilogb) (expected) - MANT_DIG); + switch (fpclassify (expected)) + { + case FP_ZERO: + /* ilogb (0) isn't allowed. */ + ulp = diff / FUNC(ldexp) (1.0, - MANT_DIG); + break; + case FP_NORMAL: + ulp = diff / FUNC(ldexp) (1.0, FUNC(ilogb) (expected) - MANT_DIG); + break; + case FP_SUBNORMAL: + ulp = (FUNC(ldexp) (diff, MANT_DIG) + / FUNC(ldexp) (1.0, FUNC(ilogb) (expected))); + break; + default: + /* It should never happen. */ + abort (); + break; + } set_max_error (ulp, curr_max_error); print_diff = 1; if ((exceptions & IGNORE_ZERO_INF_SIGN) == 0 -- cgit 1.4.1