mirror of
https://sourceware.org/git/glibc.git
synced 2024-12-18 00:41:06 +00:00
* sysdeps/alpha/fpu/bits/mathinline.h (__floorf, __floor):
Early out for -0. Optimize for !_IEEE_FP_INEXACT.
* sysdeps/alpha/fpu/s_floor.c: New.
* sysdeps/alpha/fpu/s_floorf.c: New.
* sysdeps/alpha/fpu/s_ceil.c: New.
* sysdeps/alpha/fpu/s_ceilf.c: New.
This commit is contained in:
Richard Henderson
parent
9b1370b857
commit
d7bb0de178
@ -75,7 +75,7 @@ __inline_copysign(copysign, double)
|
||||
#undef __MATH_INLINE_copysign
|
||||
|
||||
|
||||
#if defined __GNUC__ && (__GNUC__ > 2 || __GNUC__ == 2 && __GNUC_MINOR__ >= 8)
|
||||
#if defined __GNUC__ && (__GNUC__ > 2 || (__GNUC__ == 2 && __GNUC_MINOR__ >= 8))
|
||||
__MATH_INLINE float __fabsf (float __x) { return __builtin_fabsf (__x); }
|
||||
__MATH_INLINE float fabsf (float __x) { return __builtin_fabsf (__x); }
|
||||
__MATH_INLINE double __fabs (double __x) { return __builtin_fabs (__x); }
|
||||
@ -106,7 +106,8 @@ __inline_fabs(fabs, double)
|
||||
__MATH_INLINE float
|
||||
__floorf (float __x)
|
||||
{
|
||||
if (fabsf (__x) < 16777216.0f) /* 1 << FLT_MANT_DIG */
|
||||
/* Check not zero since floor(-0) == -0. */
|
||||
if (__x != 0 && 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
|
||||
@ -116,10 +117,13 @@ __floorf (float __x)
|
||||
float __tmp1, __tmp2;
|
||||
|
||||
__asm ("cvtst/s %3,%2\n\t"
|
||||
#ifdef _IEEE_FP_INEXACT
|
||||
"cvttq/svim %2,%1\n\t"
|
||||
"cvtqt/suim %1,%0\n\t"
|
||||
"trapb"
|
||||
: "=&f"(__x), "=&f"(__tmp1), "=&f"(__tmp2)
|
||||
#else
|
||||
"cvttq/svm %2,%1\n\t"
|
||||
#endif
|
||||
"cvtqt/m %1,%0\n\t"
|
||||
: "=f"(__x), "=&f"(__tmp1), "=&f"(__tmp2)
|
||||
: "f"(__x));
|
||||
}
|
||||
return __x;
|
||||
@ -128,13 +132,17 @@ __floorf (float __x)
|
||||
__MATH_INLINE double
|
||||
__floor (double __x)
|
||||
{
|
||||
if (fabs (__x) < 9007199254740992.0) /* 1 << DBL_MANT_DIG */
|
||||
if (__x != 0 && fabs (__x) < 9007199254740992.0) /* 1 << DBL_MANT_DIG */
|
||||
{
|
||||
double __tmp1;
|
||||
__asm ("cvttq/svim %2,%1\n\t"
|
||||
"cvtqt/suim %1,%0\n\t"
|
||||
"trapb"
|
||||
: "=&f"(__x), "=&f"(__tmp1)
|
||||
__asm (
|
||||
#ifdef _IEEE_FP_INEXACT
|
||||
"cvttq/svim %2,%1\n\t"
|
||||
#else
|
||||
"cvttq/svm %2,%1\n\t"
|
||||
#endif
|
||||
"cvtqt/m %1,%0\n\t"
|
||||
: "=f"(__x), "=&f"(__tmp1)
|
||||
: "f"(__x));
|
||||
}
|
||||
return __x;
|
||||
|
||||
58
sysdeps/alpha/fpu/s_ceil.c
Normal file
58
sysdeps/alpha/fpu/s_ceil.c
Normal file
@ -0,0 +1,58 @@
|
||||
/* Copyright (C) 1998 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 Library General Public License as
|
||||
published by the Free Software Foundation; either version 2 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
|
||||
Library General Public License for more details.
|
||||
|
||||
You should have received a copy of the GNU Library General Public
|
||||
License along with the GNU C Library; see the file COPYING.LIB. If not,
|
||||
write to the Free Software Foundation, Inc., 59 Temple Place - Suite 330,
|
||||
Boston, MA 02111-1307, USA. */
|
||||
|
||||
#include <math.h>
|
||||
|
||||
double
|
||||
__ceil (double x)
|
||||
{
|
||||
if (x != 0 && fabs (x) < 9007199254740992.0) /* 1 << DBL_MANT_DIG */
|
||||
{
|
||||
double tmp1;
|
||||
unsigned long fpcr0, fpcr1;
|
||||
unsigned long pinf = 3UL << 58;
|
||||
|
||||
/* Set round to +inf. */
|
||||
__asm __volatile("excb; mf_fpcr %0" : "=f"(fpcr0));
|
||||
__asm __volatile("mt_fpcr %0; excb" : : "f"(fpcr0 | pinf));
|
||||
|
||||
/* Calculate! */
|
||||
#ifdef _IEEE_FP_INEXACT
|
||||
__asm("cvttq/svid %2,%1\n\tcvtqt/suid %1,%0"
|
||||
: "=f"(x), "=&f"(tmp1)
|
||||
: "f"(x));
|
||||
#else
|
||||
__asm("cvttq/svd %2,%1\n\tcvtqt/d %1,%0"
|
||||
: "=f"(x), "=&f"(tmp1)
|
||||
: "f"(x));
|
||||
#endif
|
||||
|
||||
/* Reset rounding mode, while retaining new exception bits. */
|
||||
__asm __volatile("excb; mf_fpcr %0" : "=f"(fpcr1));
|
||||
fpcr0 = (fpcr0 & pinf) | (fpcr1 & ~pinf);
|
||||
__asm __volatile("mt_fpcr %0; excb" : : "f"(fpcr0));
|
||||
}
|
||||
return x;
|
||||
}
|
||||
|
||||
weak_alias (__ceil, ceil)
|
||||
#ifdef NO_LONG_DOUBLE
|
||||
strong_alias (__ceil, __ceill)
|
||||
weak_alias (__ceil, ceill)
|
||||
#endif
|
||||
59
sysdeps/alpha/fpu/s_ceilf.c
Normal file
59
sysdeps/alpha/fpu/s_ceilf.c
Normal file
@ -0,0 +1,59 @@
|
||||
/* Copyright (C) 1998 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 Library General Public License as
|
||||
published by the Free Software Foundation; either version 2 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
|
||||
Library General Public License for more details.
|
||||
|
||||
You should have received a copy of the GNU Library General Public
|
||||
License along with the GNU C Library; see the file COPYING.LIB. If not,
|
||||
write to the Free Software Foundation, Inc., 59 Temple Place - Suite 330,
|
||||
Boston, MA 02111-1307, USA. */
|
||||
|
||||
#include <math.h>
|
||||
|
||||
float
|
||||
__ceilf (float x)
|
||||
{
|
||||
if (x != 0 && fabsf (x) < 16777216.0f) /* 1 << FLT_MANT_DIG */
|
||||
{
|
||||
float tmp1, tmp2;
|
||||
unsigned long fpcr0, fpcr1;
|
||||
unsigned long pinf = 3UL << 58;
|
||||
|
||||
/* Set round to +inf. */
|
||||
__asm __volatile("excb; mf_fpcr %0" : "=f"(fpcr0));
|
||||
__asm __volatile("mt_fpcr %0; excb" : : "f"(fpcr0 | pinf));
|
||||
|
||||
/* Calculate!
|
||||
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. */
|
||||
|
||||
#ifdef _IEEE_FP_INEXACT
|
||||
__asm("cvtst/s %3,%2\n\tcvttq/svid %2,%1\n\tcvtqt/suid %1,%0"
|
||||
: "=f"(x), "=&f"(tmp1), "=&f"(tmp2)
|
||||
: "f"(x));
|
||||
#else
|
||||
__asm("cvtst/s %3,%2\n\tcvttq/svd %2,%1\n\tcvtqt/d %1,%0"
|
||||
: "=f"(x), "=&f"(tmp1), "=&f"(tmp2)
|
||||
: "f"(x));
|
||||
#endif
|
||||
|
||||
/* Reset rounding mode, while retaining new exception bits. */
|
||||
__asm __volatile("excb; mf_fpcr %0" : "=f"(fpcr1));
|
||||
fpcr0 = (fpcr0 & pinf) | (fpcr1 & ~pinf);
|
||||
__asm __volatile("mt_fpcr %0; excb" : : "f"(fpcr0));
|
||||
}
|
||||
return x;
|
||||
}
|
||||
|
||||
weak_alias (__ceilf, ceilf)
|
||||
39
sysdeps/alpha/fpu/s_floor.c
Normal file
39
sysdeps/alpha/fpu/s_floor.c
Normal file
@ -0,0 +1,39 @@
|
||||
/* Copyright (C) 1998 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 Library General Public License as
|
||||
published by the Free Software Foundation; either version 2 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
|
||||
Library General Public License for more details.
|
||||
|
||||
You should have received a copy of the GNU Library General Public
|
||||
License along with the GNU C Library; see the file COPYING.LIB. If not,
|
||||
write to the Free Software Foundation, Inc., 59 Temple Place - Suite 330,
|
||||
Boston, MA 02111-1307, USA. */
|
||||
|
||||
#ifndef __USE_EXTERN_INLINES
|
||||
#define __USE_EXTERN_INLINES
|
||||
#endif
|
||||
#define __floor __i_floor
|
||||
|
||||
#include <math.h>
|
||||
|
||||
#undef __floor
|
||||
|
||||
double
|
||||
__floor (double x)
|
||||
{
|
||||
return __i_floor(x);
|
||||
}
|
||||
|
||||
weak_alias (__floor, floor)
|
||||
#ifdef NO_LONG_DOUBLE
|
||||
strong_alias (__floor, __floorl)
|
||||
weak_alias (__floor, floorl)
|
||||
#endif
|
||||
35
sysdeps/alpha/fpu/s_floorf.c
Normal file
35
sysdeps/alpha/fpu/s_floorf.c
Normal file
@ -0,0 +1,35 @@
|
||||
/* Copyright (C) 1998 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 Library General Public License as
|
||||
published by the Free Software Foundation; either version 2 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
|
||||
Library General Public License for more details.
|
||||
|
||||
You should have received a copy of the GNU Library General Public
|
||||
License along with the GNU C Library; see the file COPYING.LIB. If not,
|
||||
write to the Free Software Foundation, Inc., 59 Temple Place - Suite 330,
|
||||
Boston, MA 02111-1307, USA. */
|
||||
|
||||
#ifndef __USE_EXTERN_INLINES
|
||||
#define __USE_EXTERN_INLINES
|
||||
#endif
|
||||
#define __floorf __i_floorf
|
||||
|
||||
#include <math.h>
|
||||
|
||||
#undef __floorf
|
||||
|
||||
float
|
||||
__floorf (float x)
|
||||
{
|
||||
return __i_floorf(x);
|
||||
}
|
||||
|
||||
weak_alias (__floorf, floorf)
|
||||
Loading…
Reference in New Issue
Block a user