AuroraMiddleware/glibc
1
0
Fork 0
mirror of https://sourceware.org/git/glibc.git synced 2024-11-10 15:20:10 +00:00
Code Issues Packages Projects Releases Wiki Activity

[BZ #13848] alpha: Fix s_nearbyint implementation.

Browse Source
This commit is contained in:
Richard Henderson 2012-06-15 12:17:11 -07:00
parent 45c8de68d5
commit f56ed78d4a
2 changed files with 34 additions and 26 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

29
sysdeps/alpha/fpu/s_nearbyint.c
View File

@ -1,4 +1,4 @@
/* Copyright (C) 2007 Free Software Foundation, Inc.
/* Copyright (C) 2000-2012 Free Software Foundation, Inc.
This file is part of the GNU C Library.
Contributed by Richard Henderson.
@ -19,22 +19,23 @@
#include <math.h>
#include <math_ldbl_opt.h>
#ifdef _IEEE_FP_INEXACT
#error "Don't compile with -mieee-with-inexact"
#endif
double
__nearbyint (double x)
{
double two52 = copysign (0x1.0p52, x);
double r;
r = x + two52;
r = r - two52;
if (isless (fabs (x), 9007199254740992.0)) /* 1 << DBL_MANT_DIG */
{
double tmp1, new_x;
__asm ("cvttq/svd %2,%1\n\t"
"cvtqt/d %1,%0\n\t"
: "=f"(new_x), "=&f"(tmp1)
: "f"(x));
/* nearbyint(-0.1) == -0, and in general we'll always have the same sign
as our input. */
return copysign (r, x);
/* nearbyint(-0.1) == -0, and in general we'll always have the same
sign as our input. */
x = copysign(new_x, x);
}
return x;
}
weak_alias (__nearbyint, nearbyint)
@ -42,6 +43,6 @@ weak_alias (__nearbyint, nearbyint)
strong_alias (__nearbyint, __nearbyintl)
weak_alias (__nearbyint, nearbyintl)
#endif
#if LONG_DOUBLE_COMPAT(libm, GLIBC_2_1)
compat_symbol (libm, __nearbyint, nearbyintl, GLIBC_2_1);
#if LONG_DOUBLE_COMPAT(libm, GLIBC_2_0)
compat_symbol (libm, __nearbyint, nearbyintl, GLIBC_2_0);
#endif

31
sysdeps/alpha/fpu/s_nearbyintf.c
View File

@ -1,4 +1,4 @@
/* Copyright (C) 2007 Free Software Foundation, Inc.
/* Copyright (C) 2007-2012 Free Software Foundation, Inc.
This file is part of the GNU C Library.
Contributed by Richard Henderson.
@ -18,22 +18,29 @@
#include <math.h>
#ifdef _IEEE_FP_INEXACT
#error "Don't compile with -mieee-with-inexact"
#endif
float
__nearbyintf (float x)
{
float two23 = copysignf (0x1.0p23, x);
float r;
if (isless (fabsf (x), 16777216.0f)) /* 1 << FLT_MANT_DIG */
{
/* Note that Alpha S_Floating is stored in registers in a
restricted T_Floating format, so we don't even need to
convert back to S_Floating in the end. The initial
conversion to T_Floating is needed to handle denormals. */
r = x + two23;
r = r - two23;
float tmp1, tmp2, new_x;
/* nearbyint(-0.1) == -0, and in general we'll always have the same sign
as our input. */
return copysign (r, x);
__asm ("cvtst/s %3,%2\n\t"
"cvttq/svd %2,%1\n\t"
"cvtqt/d %1,%0\n\t"
: "=f"(new_x), "=&f"(tmp1), "=&f"(tmp2)
: "f"(x));
/* nearbyintf(-0.1) == -0, and in general we'll always have the same
sign as our input. */
x = copysignf(new_x, x);
}
return x;
}
weak_alias (__nearbyintf, nearbyintf)