diff options
| author | Szabolcs Nagy <szabolcs.nagy@arm.com> | 2022-10-13 11:10:10 +0100 |
|---|---|---|
| committer | Szabolcs Nagy <szabolcs.nagy@arm.com> | 2022-10-13 11:43:54 +0100 |
| commit | 4ad4ea848b335d803bcf4754a1cea78f697eaec5 (patch) | |
| tree | 415fff2c47389cdbe5f40fe3c1820ac4530c0f9f | |
| parent | f660eb9cc7ddaed3900d0cb3b2468816790855ff (diff) | |
| download | glibc-4ad4ea848b335d803bcf4754a1cea78f697eaec5.tar glibc-4ad4ea848b335d803bcf4754a1cea78f697eaec5.tar.gz glibc-4ad4ea848b335d803bcf4754a1cea78f697eaec5.tar.bz2 glibc-4ad4ea848b335d803bcf4754a1cea78f697eaec5.zip | |
math: Fix asin and acos invalid exception
This works around a gcc issue where it const folds inf/inf into nan,
preventing the invalid exception signal to be raised.
(x-x)/(x-x) is more robust against optimizations and works for x==nan
too.
The issue should be fixed in gcc-11.3.0 and gcc-12, but glibc supports
older compilers.
| -rw-r--r-- | sysdeps/ieee754/dbl-64/e_asin.c | 18 |
1 files changed, 2 insertions, 16 deletions
diff --git a/sysdeps/ieee754/dbl-64/e_asin.c b/sysdeps/ieee754/dbl-64/e_asin.c index e7ea0cbe8f..6b7c971e76 100644 --- a/sysdeps/ieee754/dbl-64/e_asin.c +++ b/sysdeps/ieee754/dbl-64/e_asin.c @@ -165,14 +165,7 @@ __ieee754_asin(double x){ /*---------------------------- |x|>=1 -------------------------------*/ else if (k==0x3ff00000 && u.i[LOW_HALF]==0) return (m>0)?hp0.x:-hp0.x; else - if (k>0x7ff00000 || (k == 0x7ff00000 && u.i[LOW_HALF] != 0)) return x + x; - else { - u.i[HIGH_HALF]=0x7ff00000; - v.i[HIGH_HALF]=0x7ff00000; - u.i[LOW_HALF]=0; - v.i[LOW_HALF]=0; - return u.x/v.x; /* NaN */ - } + return (x - x) / (x - x); } #ifndef __ieee754_asin libm_alias_finite (__ieee754_asin, __asin) @@ -334,14 +327,7 @@ __ieee754_acos(double x) else if (k==0x3ff00000 && u.i[LOW_HALF]==0) return (m>0)?0:2.0*hp0.x; else - if (k>0x7ff00000 || (k == 0x7ff00000 && u.i[LOW_HALF] != 0)) return x + x; - else { - u.i[HIGH_HALF]=0x7ff00000; - v.i[HIGH_HALF]=0x7ff00000; - u.i[LOW_HALF]=0; - v.i[LOW_HALF]=0; - return u.x/v.x; - } + return (x - x) / (x - x); } #ifndef __ieee754_acos libm_alias_finite (__ieee754_acos, __acos) |