* sysdeps/alpha/fpu/s_llround.c: New file.
* sysdeps/alpha/fpu/s_llroundf.c: New file. * sysdeps/alpha/fpu/s_lround.c: New file. * sysdeps/alpha/fpu/s_lroundf.c: New file. * sysdeps/alpha/fpu/s_round.c: New file. * sysdeps/alpha/fpu/s_roundf.c: New file. * sysdeps/alpha/fpu/s_trunc.c: New file. * sysdeps/alpha/fpu/s_truncf.c: New file. * sysdeps/alpha/fpu/s_llround.c: New file. * sysdeps/alpha/fpu/s_llroundf.c: New file. * sysdeps/alpha/fpu/s_lround.c: New file. * sysdeps/alpha/fpu/s_lroundf.c: New file. * sysdeps/alpha/fpu/s_round.c: New file. * sysdeps/alpha/fpu/s_roundf.c: New file. * sysdeps/alpha/fpu/s_trunc.c: New file. * sysdeps/alpha/fpu/s_truncf.c: New file.
This commit is contained in:
Richard Henderson
parent
aea1d19e2f
commit
b6df547c34
|
|
@ -1,5 +1,14 @@
|
|||
2007-03-14 Richard Henderson <rth@redhat.com>
|
||||
|
||||
* sysdeps/alpha/fpu/s_llround.c: New file.
|
||||
* sysdeps/alpha/fpu/s_llroundf.c: New file.
|
||||
* sysdeps/alpha/fpu/s_lround.c: New file.
|
||||
* sysdeps/alpha/fpu/s_lroundf.c: New file.
|
||||
* sysdeps/alpha/fpu/s_round.c: New file.
|
||||
* sysdeps/alpha/fpu/s_roundf.c: New file.
|
||||
* sysdeps/alpha/fpu/s_trunc.c: New file.
|
||||
* sysdeps/alpha/fpu/s_truncf.c: New file.
|
||||
|
||||
* sysdeps/alpha/fpu/s_ceil.c: Rewrite without branches.
|
||||
* sysdeps/alpha/fpu/s_ceilf.c: Likewise.
|
||||
* sysdeps/alpha/fpu/s_floor.c: Likewise.
|
||||
|
|
|
|||
|
|
@ -21,9 +21,7 @@
|
|||
#include <math_ldbl_opt.h>
|
||||
|
||||
|
||||
/* Use the -inf rounding mode conversion instructions to implement
|
||||
floor. We note when the exponent is large enough that the value
|
||||
must be integral, as this avoids unpleasant integer overflows. */
|
||||
/* Use the -inf rounding mode conversion instructions to implement floor. */
|
||||
|
||||
double
|
||||
__floor (double x)
|
||||
|
|
|
|||
|
|
@ -20,9 +20,7 @@
|
|||
#include <math.h>
|
||||
|
||||
|
||||
/* Use the -inf rounding mode conversion instructions to implement
|
||||
floor. We note when the exponent is large enough that the value
|
||||
must be integral, as this avoids unpleasant integer overflows. */
|
||||
/* Use the -inf rounding mode conversion instructions to implement floor. */
|
||||
|
||||
float
|
||||
__floorf (float x)
|
||||
|
|
|
|||
|
|
@ -0,0 +1 @@
|
|||
/* In s_lround.c. */
|
||||
|
|
@ -0,0 +1 @@
|
|||
/* In s_lroundf.c. */
|
||||
|
|
@ -0,0 +1,48 @@
|
|||
/* Copyright (C) 2007 Free Software Foundation, Inc.
|
||||
This file is part of the GNU C Library.
|
||||
|
||||
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, write to the Free
|
||||
Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA
|
||||
02111-1307 USA. */
|
||||
|
||||
#define __llround not___llround
|
||||
#define llround not_llround
|
||||
#include <math.h>
|
||||
#include <math_ldbl_opt.h>
|
||||
#undef __llround
|
||||
#undef llround
|
||||
|
||||
long int
|
||||
__lround (double x)
|
||||
{
|
||||
double adj;
|
||||
|
||||
adj = 0x1.fffffffffffffp-2; /* nextafter (0.5, 0.0) */
|
||||
adj = copysign (adj, x);
|
||||
return x + adj;
|
||||
}
|
||||
|
||||
strong_alias (__lround, __llround)
|
||||
weak_alias (__lround, lround)
|
||||
weak_alias (__llround, llround)
|
||||
#ifdef NO_LONG_DOUBLE
|
||||
strong_alias (__lround, __lroundl)
|
||||
strong_alias (__lround, __llroundl)
|
||||
weak_alias (__lroundl, lroundl)
|
||||
weak_alias (__llroundl, llroundl)
|
||||
#endif
|
||||
#if LONG_DOUBLE_COMPAT(libm, GLIBC_2_0)
|
||||
compat_symbol (libm, __lround, lroundl, GLIBC_2_0);
|
||||
compat_symbol (libm, __llround, llroundl, GLIBC_2_0);
|
||||
#endif
|
||||
|
|
@ -0,0 +1,38 @@
|
|||
/* Copyright (C) 2007 Free Software Foundation, Inc.
|
||||
This file is part of the GNU C Library.
|
||||
|
||||
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, write to the Free
|
||||
Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA
|
||||
02111-1307 USA. */
|
||||
|
||||
#define __llroundf not___llroundf
|
||||
#define llroundf not_llroundf
|
||||
#include <math.h>
|
||||
#undef __llroundf
|
||||
#undef llroundf
|
||||
|
||||
|
||||
long int
|
||||
__lroundf (float x)
|
||||
{
|
||||
float adj;
|
||||
|
||||
adj = 0x1.fffffep-2; /* nextafterf (0.5f, 0.0f) */
|
||||
adj = copysignf (adj, x);
|
||||
return x + adj;
|
||||
}
|
||||
|
||||
strong_alias (__lroundf, __llroundf)
|
||||
weak_alias (__lroundf, lroundf)
|
||||
weak_alias (__llroundf, llroundf)
|
||||
|
|
@ -0,0 +1,49 @@
|
|||
/* Copyright (C) 2007 Free Software Foundation, Inc.
|
||||
This file is part of the GNU C Library.
|
||||
|
||||
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, write to the Free
|
||||
Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA
|
||||
02111-1307 USA. */
|
||||
|
||||
#include <math.h>
|
||||
#include <math_ldbl_opt.h>
|
||||
|
||||
|
||||
double
|
||||
__round (double x)
|
||||
{
|
||||
const double almost_half = 0x1.fffffffffffffp-2;
|
||||
const double two52 = 0x1.0p52;
|
||||
double tmp, r;
|
||||
|
||||
__asm (
|
||||
#ifdef _IEEE_FP_INEXACT
|
||||
"addt/suic %2, %3, %1\n\tsubt/suic %1, %3, %0"
|
||||
#else
|
||||
"addt/suc %2, %3, %1\n\tsubt/suc %1, %3, %0"
|
||||
#endif
|
||||
: "=&f"(r), "=&f"(tmp)
|
||||
: "f"(fabs (x) + almost_half), "f"(two52));
|
||||
|
||||
return copysign (r, x);
|
||||
}
|
||||
|
||||
weak_alias (__round, round)
|
||||
#ifdef NO_LONG_DOUBLE
|
||||
strong_alias (__round, __roundl)
|
||||
weak_alias (__roundl, roundl)
|
||||
#endif
|
||||
#if LONG_DOUBLE_COMPAT(libm, GLIBC_2_0)
|
||||
compat_symbol (libm, __round, roundl, GLIBC_2_0);
|
||||
#endif
|
||||
|
|
@ -0,0 +1,44 @@
|
|||
/* Copyright (C) 2007 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, write to the Free
|
||||
Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA
|
||||
02111-1307 USA. */
|
||||
|
||||
#include <math.h>
|
||||
|
||||
|
||||
float
|
||||
__roundf (float x)
|
||||
{
|
||||
const float almost_half = 0x1.fffffep-2;
|
||||
const float two23 = 0x1.0p23;
|
||||
float r, tmp;
|
||||
|
||||
__asm (
|
||||
#ifdef _IEEE_FP_INEXACT
|
||||
"adds/suic %2, %3, %1\n\tsubs/suic %1, %3, %0"
|
||||
#else
|
||||
"adds/suc %2, %3, %1\n\tsubs/suc %1, %3, %0"
|
||||
#endif
|
||||
: "=&f"(r), "=&f"(tmp)
|
||||
: "f"(fabsf (x) + almost_half), "f"(two23));
|
||||
|
||||
/* round(-0) == -0, and in general we'll always have the same
|
||||
sign as our input. */
|
||||
return copysignf (r, x);
|
||||
}
|
||||
|
||||
weak_alias (__roundf, roundf)
|
||||
|
|
@ -0,0 +1,53 @@
|
|||
/* Copyright (C) 2007 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, write to the Free
|
||||
Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA
|
||||
02111-1307 USA. */
|
||||
|
||||
#include <math.h>
|
||||
#include <math_ldbl_opt.h>
|
||||
|
||||
|
||||
/* Use the chopped rounding mode conversion instructions to implement trunc. */
|
||||
|
||||
double
|
||||
__trunc (double x)
|
||||
{
|
||||
double two52 = copysign (0x1.0p52, x);
|
||||
double r, tmp;
|
||||
|
||||
__asm (
|
||||
#ifdef _IEEE_FP_INEXACT
|
||||
"addt/suic %2, %3, %1\n\tsubt/suic %1, %3, %0"
|
||||
#else
|
||||
"addt/suc %2, %3, %1\n\tsubt/suc %1, %3, %0"
|
||||
#endif
|
||||
: "=&f"(r), "=&f"(tmp)
|
||||
: "f"(x), "f"(two52));
|
||||
|
||||
/* trunc(-0) == -0, and in general we'll always have the same
|
||||
sign as our input. */
|
||||
return copysign (r, x);
|
||||
}
|
||||
|
||||
weak_alias (__trunc, trunc)
|
||||
#ifdef NO_LONG_DOUBLE
|
||||
strong_alias (__trunc, __truncl)
|
||||
weak_alias (__trunc, truncl)
|
||||
#endif
|
||||
#if LONG_DOUBLE_COMPAT(libm, GLIBC_2_0)
|
||||
compat_symbol (libm, __trunc, truncl, GLIBC_2_0);
|
||||
#endif
|
||||
|
|
@ -0,0 +1,45 @@
|
|||
/* Copyright (C) 2007 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, write to the Free
|
||||
Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA
|
||||
02111-1307 USA. */
|
||||
|
||||
#include <math.h>
|
||||
|
||||
|
||||
/* Use the chopped rounding mode conversion instructions to implement trunc. */
|
||||
|
||||
float
|
||||
__truncf (float x)
|
||||
{
|
||||
float two23 = copysignf (0x1.0p23, x);
|
||||
float r, tmp;
|
||||
|
||||
__asm (
|
||||
#ifdef _IEEE_FP_INEXACT
|
||||
"adds/suic %2, %3, %1\n\tsubs/suic %1, %3, %0"
|
||||
#else
|
||||
"adds/suc %2, %3, %1\n\tsubs/suc %1, %3, %0"
|
||||
#endif
|
||||
: "=&f"(r), "=&f"(tmp)
|
||||
: "f"(x), "f"(two23));
|
||||
|
||||
/* trunc(-0) == -0, and in general we'll always have the same
|
||||
sign as our input. */
|
||||
return copysignf (r, x);
|
||||
}
|
||||
|
||||
weak_alias (__truncf, truncf)
|
||||
Loading…
Reference in New Issue