diff options
| author | Szabolcs Nagy <szabolcs.nagy@arm.com> | 2017-12-05 18:09:22 +0000 |
|---|---|---|
| committer | Szabolcs Nagy <szabolcs.nagy@arm.com> | 2017-12-05 18:09:22 +0000 |
| commit | 00d54af7c86cdcabc18c00619c229f49f4978610 (patch) | |
| tree | 6bc09f21632f62fc82f9f5b057e4fa9493adf3ca | |
| parent | 1f70eae4cf0ca826ddbc2b61720e8aeae1e13de4 (diff) | |
| download | glibc-00d54af7c86cdcabc18c00619c229f49f4978610.zip glibc-00d54af7c86cdcabc18c00619c229f49f4978610.tar.gz glibc-00d54af7c86cdcabc18c00619c229f49f4978610.tar.bz2 | |
[PATCH] fix sinf(NAN)
sinf(NAN) should not signal invalid fp exception
so use isless instead of < where NAN is compared.
this makes the sinf tests pass on aarch64.
* sysdeps/ieee754/flt-32/s_sinf.c (sinf): Use isless.
| -rw-r--r-- | ChangeLog | 4 | ||||
| -rw-r--r-- | sysdeps/ieee754/flt-32/s_sinf.c | 6 |
2 files changed, 7 insertions, 3 deletions
@@ -1,3 +1,7 @@ +2017-12-05 Szabolcs Nagy <szabolcs.nagy@arm.com> + + * sysdeps/ieee754/flt-32/s_sinf.c (sinf): Use isless. + 2017-12-05 Joseph Myers <joseph@codesourcery.com> * stdlib/strfromd.c: Include <bits/floatn.h>. diff --git a/sysdeps/ieee754/flt-32/s_sinf.c b/sysdeps/ieee754/flt-32/s_sinf.c index 0b76477..8b98573 100644 --- a/sysdeps/ieee754/flt-32/s_sinf.c +++ b/sysdeps/ieee754/flt-32/s_sinf.c @@ -128,7 +128,7 @@ SINF_FUNC (float x) double theta = x; double abstheta = fabs (theta); /* If |x|< Pi/4. */ - if (abstheta < M_PI_4) + if (isless (abstheta, M_PI_4)) { if (abstheta >= 0x1p-5) /* |x| >= 2^-5. */ { @@ -162,8 +162,8 @@ SINF_FUNC (float x) } else /* |x| >= Pi/4. */ { - unsigned int signbit = (x < 0); - if (abstheta < 9 * M_PI_4) /* |x| < 9*Pi/4. */ + unsigned int signbit = isless (x, 0); + if (isless (abstheta, 9 * M_PI_4)) /* |x| < 9*Pi/4. */ { /* There are cases where FE_UPWARD rounding mode can produce a result of abstheta * inv_PI_4 == 9, |