diff options
| author | Szabolcs Nagy <szabolcs.nagy@arm.com> | 2018-12-10 11:08:36 +0000 |
|---|---|---|
| committer | Szabolcs Nagy <szabolcs.nagy@arm.com> | 2018-12-11 10:01:43 +0000 |
| commit | 505b5b292293a5d6bd4046a6bc7f8c2381a33da4 (patch) | |
| tree | feeacf415d6fdbb306f2911deede485a6d92fda7 /math/auto-libm-test-in | |
| parent | 304c61a24f909168c16793ccf7c686237e53d003 (diff) | |
| download | glibc-505b5b292293a5d6bd4046a6bc7f8c2381a33da4.tar glibc-505b5b292293a5d6bd4046a6bc7f8c2381a33da4.tar.gz glibc-505b5b292293a5d6bd4046a6bc7f8c2381a33da4.tar.bz2 glibc-505b5b292293a5d6bd4046a6bc7f8c2381a33da4.zip | |
Fix powf overflow handling in non-nearest rounding mode [BZ #23961]
The threshold value at which powf overflows depends on the rounding mode
and the current check did not take this into account. So when the result
was rounded away from zero it could become infinity without setting
errno to ERANGE.
Example: pow(0x1.7ac7cp+5, 23) is 0x1.fffffep+127 + 0.1633ulp
If the result goes above 0x1.fffffep+127 + 0.5ulp then errno is set,
which is fine in nearest rounding mode, but
powf(0x1.7ac7cp+5, 23) is inf in upward rounding mode
powf(-0x1.7ac7cp+5, 23) is -inf in downward rounding mode
and the previous implementation did not set errno in these cases.
The fix tries to avoid affecting the common code path or calling a
function that may introduce a stack frame, so float arithmetics is used
to check the rounding mode and the threshold is selected accordingly.
[BZ #23961]
* math/auto-libm-test-in: Add new test case.
* math/auto-libm-test-out-pow: Regenerated.
* sysdeps/ieee754/flt-32/e_powf.c (__powf): Fix overflow check.
Diffstat (limited to 'math/auto-libm-test-in')
| -rw-r--r-- | math/auto-libm-test-in | 2 |
1 files changed, 2 insertions, 0 deletions
diff --git a/math/auto-libm-test-in b/math/auto-libm-test-in index 785328562f..c488e7225c 100644 --- a/math/auto-libm-test-in +++ b/math/auto-libm-test-in @@ -7154,6 +7154,8 @@ pow 0x1.430d4cp+0 0x5.0e462p+4 pow 0x9.8b82ap-4 -0x1.99907ap+12 pow 0xd.73035p-4 -0x1.47bb8p+8 pow 0x1.059c76p+0 0x1.ff80bep+11 +pow 0x1.7ac7cp+5 23 +pow -0x1.7ac7cp+5 23 sin 0 sin -0 |