qiurui
/
glibc
mirror of git://sourceware.org/git/glibc.git
1
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

alpha: Remove nearbyint and nearbyintf implementations 

The original implementation was written for EV5, which does not
record inexact in the status register for /SU (but no /I) insns.
But EV6 does record the inexact status; the lack of /I simply
means that the exception is suppressed.

Adding feholdexcept becomes the bulk of the overhead, so we might
as well use the default implementation.
This commit is contained in:
Richard Henderson 2014-07-03 10:37:59 -07:00
parent d5e16ade14
commit 69862e2511
3 changed files with 5 additions and 72 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
ChangeLog
View File

@ -1,5 +1,9 @@
2014-07-03 Richard Henderson <rth@redhat.com> 2014-07-03 Richard Henderson <rth@redhat.com>
* sysdeps/alpha/fpu/s_nearbyintf.c: Remove file.
* sysdeps/alpha/fpu/s_nearbyint.c (__nearbyint): Remove;
include sysdeps/ieee754/dbl-64/wordsize-64/s_nearbyint.c.
* sysdeps/alpha/fpu/math_private.h (math_opt_barrier): New. * sysdeps/alpha/fpu/math_private.h (math_opt_barrier): New.
(math_force_eval): New. (math_force_eval): New.

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

@ -16,33 +16,8 @@
License along with the GNU C Library. If not, see License along with the GNU C Library. If not, see
<http://www.gnu.org/licenses/>. */ <http://www.gnu.org/licenses/>. */
#include <math.h> #include <sysdeps/ieee754/dbl-64/wordsize-64/s_nearbyint.c>
#include <math_ldbl_opt.h>
double
__nearbyint (double x)
{
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. */
x = copysign(new_x, x);
}
return x;
}
weak_alias (__nearbyint, nearbyint)
#ifdef NO_LONG_DOUBLE
strong_alias (__nearbyint, __nearbyintl)
weak_alias (__nearbyint, nearbyintl)
#endif
#if LONG_DOUBLE_COMPAT(libm, GLIBC_2_1) #if LONG_DOUBLE_COMPAT(libm, GLIBC_2_1)
compat_symbol (libm, __nearbyint, nearbyintl, GLIBC_2_1); compat_symbol (libm, __nearbyint, nearbyintl, GLIBC_2_1);
#endif #endif

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

@ -1,46 +0,0 @@
/* Copyright (C) 2007-2014 Free Software Foundation, Inc.
This file is part of the GNU C Library.
Contributed by Richard Henderson.
The GNU C Library is free software; you can redistribute it and/or
modify it under the terms of the GNU Lesser General Public
License as published by the Free Software Foundation; either
version 2.1 of the License, or (at your option) any later version.
The GNU C Library is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
Lesser General Public License for more details.
You should have received a copy of the GNU Lesser General Public
License along with the GNU C Library. If not, see
<http://www.gnu.org/licenses/>. */
#include <math.h>
float
__nearbyintf (float x)
{
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. */
float tmp1, tmp2, new_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)