diff options
| author | Adhemerval Zanella <azanella@linux.vnet.ibm.com> | 2012-06-01 10:10:18 -0300 |
|---|---|---|
| committer | Adhemerval Zanella <azanella@linux.vnet.ibm.com> | 2012-06-01 10:10:18 -0300 |
| commit | 73a68f94d650ac485db6b95242e1cb6591b193d3 (patch) | |
| tree | fb07542aa7411ef1ddbf8d213521d72e208b370c /sysdeps/powerpc/fpu/k_cosf.c | |
| parent | ebc64a18c07742ebb85d7ddeca376995c0e67593 (diff) | |
| download | glibc-73a68f94d650ac485db6b95242e1cb6591b193d3.tar.gz | |
PowerPC: Fix for POWER7 sinf/cosf
This patch fixes some sinf/cosf calculations that generated unexpected
underflows exceptions.
Diffstat (limited to 'sysdeps/powerpc/fpu/k_cosf.c')
| -rw-r--r-- | sysdeps/powerpc/fpu/k_cosf.c | 5 |
1 files changed, 3 insertions, 2 deletions
diff --git a/sysdeps/powerpc/fpu/k_cosf.c b/sysdeps/powerpc/fpu/k_cosf.c index e2571d6802..d5668e3127 100644 --- a/sysdeps/powerpc/fpu/k_cosf.c +++ b/sysdeps/powerpc/fpu/k_cosf.c @@ -18,6 +18,7 @@ not, see <http://www.gnu.org/licenses/>. */ #include <math.h> +#include <fenv.h> #include <math_private.h> static const float twom27 = 7.4505806e-09; @@ -40,8 +41,8 @@ __kernel_cosf (float x, float y) ix = __builtin_fabsf (x); if (ix < twom27) { /* |x| < 2**-27 */ - if (x == 0.0) - return one; + __feraiseexcept (FE_INEXACT); + return one; } z = x * x; r = z * (C1 + z * (C2 + z * (C3 + z * (C4 + z * (C5 + z * C6))))); |