mirror of
https://sourceware.org/git/glibc.git
synced 2024-12-22 10:50:07 +00:00
* sysdeps/alpha/fpu/fraiseexcpt.c: Use get/set_fp_control instead of arithmetic instructions.
* sysdeps/alpha/fpu/s_ceil.c: Use round to -inf instead of playing
with the fpcr. Protect from INV exception.
* sysdeps/alpha/fpu/s_ceilf.c: Likewise.
* sysdeps/alpha/fpu/s_floor.c: Protect from INV exception.
* sysdeps/alpha/fpu/s_floorf.c: Likewise.
* sysdeps/alpha/fpu/s_copysign.c: New.
* sysdeps/alpha/fpu/s_copysignf.c: New.
* sysdeps/alpha/fpu/s_fabs.c: New.
* sysdeps/alpha/fpu/s_fabsf.c: New.
* sysdeps/alpha/fpu/s_rint.c: New.
* sysdeps/alpha/fpu/s_rintf.c: New.
* sysdeps/alpha/fpu/fraiseexcpt.c: Use get/set_fp_control instead
of arithmetic instructions.
* sysdeps/alpha/fpu/s_ceil.c: Use round to -inf instead of playing
with the fpcr. Protect from INV exception.
* sysdeps/alpha/fpu/s_ceilf.c: Likewise.
* sysdeps/alpha/fpu/s_floor.c: Protect from INV exception.
* sysdeps/alpha/fpu/s_floorf.c: Likewise.
* sysdeps/alpha/fpu/s_copysign.c: New.
* sysdeps/alpha/fpu/s_copysignf.c: New.
* sysdeps/alpha/fpu/s_fabs.c: New.
* sysdeps/alpha/fpu/s_fabsf.c: New.
* sysdeps/alpha/fpu/s_rint.c: New.
* sysdeps/alpha/fpu/s_rintf.c: New.
This commit is contained in:
Richard Henderson
parent
45159f0082
commit
b7434f8243
16
ChangeLog
16
ChangeLog
@ -1,5 +1,21 @@
|
||||
2000-03-20 Richard Henderson <rth@cygnus.com>
|
||||
|
||||
* sysdeps/alpha/fpu/fraiseexcpt.c: Use get/set_fp_control instead
|
||||
of arithmetic instructions.
|
||||
|
||||
* sysdeps/alpha/fpu/s_ceil.c: Use round to -inf instead of playing
|
||||
with the fpcr. Protect from INV exception.
|
||||
* sysdeps/alpha/fpu/s_ceilf.c: Likewise.
|
||||
* sysdeps/alpha/fpu/s_floor.c: Protect from INV exception.
|
||||
* sysdeps/alpha/fpu/s_floorf.c: Likewise.
|
||||
|
||||
* sysdeps/alpha/fpu/s_copysign.c: New.
|
||||
* sysdeps/alpha/fpu/s_copysignf.c: New.
|
||||
* sysdeps/alpha/fpu/s_fabs.c: New.
|
||||
* sysdeps/alpha/fpu/s_fabsf.c: New.
|
||||
* sysdeps/alpha/fpu/s_rint.c: New.
|
||||
* sysdeps/alpha/fpu/s_rintf.c: New.
|
||||
|
||||
* sysdeps/unix/sysv/linux/alpha/bits/mman.h: Add MADV constants.
|
||||
|
||||
* sysdeps/unix/sysv/linux/alpha/Makefile (sysdep_routines):
|
||||
|
||||
@ -24,43 +24,16 @@
|
||||
int
|
||||
__feraiseexcept (int excepts)
|
||||
{
|
||||
double tmp;
|
||||
double dummy;
|
||||
unsigned long int tmp;
|
||||
|
||||
/* Raise exceptions represented by EXPECTS. But we must raise only
|
||||
one signal at a time. It is important the if the overflow/underflow
|
||||
exception and the inexact exception are given at the same time,
|
||||
the overflow/underflow exception precedes the inexact exception. */
|
||||
/* Get the current exception state. */
|
||||
tmp = __ieee_get_fp_control ();
|
||||
|
||||
/* We do these bits in assembly to be certain GCC doesn't optimize
|
||||
away something important. */
|
||||
/* Set all the bits that were called for. */
|
||||
tmp |= (excepts & FE_ALL_EXCEPT);
|
||||
|
||||
/* First: invalid exception. */
|
||||
if (FE_INVALID & excepts)
|
||||
/* One example of a invalid operation is 0 * Infinity. */
|
||||
__asm__ __volatile__("mult/sui $f31,%1,%0; trapb"
|
||||
: "=&f" (tmp) : "f" (HUGE_VAL));
|
||||
|
||||
/* Next: division by zero. */
|
||||
if (FE_DIVBYZERO & excepts)
|
||||
__asm__ __volatile__("cmpteq $f31,$f31,%1; divt/sui %1,$f31,%0; trapb"
|
||||
: "=&f" (tmp), "=f" (dummy));
|
||||
|
||||
/* Next: overflow. */
|
||||
if (FE_OVERFLOW & excepts)
|
||||
__asm__ __volatile__("mult/sui %1,%1,%0; trapb"
|
||||
: "=&f" (tmp) : "f" (DBL_MAX));
|
||||
|
||||
/* Next: underflow. */
|
||||
if (FE_UNDERFLOW & excepts)
|
||||
__asm__ __volatile__("divt/sui %1,%2,%0; trapb"
|
||||
: "=&f" (tmp) : "f" (DBL_MIN),
|
||||
"f" ((double) (1UL << 60)));
|
||||
|
||||
/* Last: inexact. */
|
||||
if (FE_INEXACT & excepts)
|
||||
__asm__ __volatile__("divt/sui %1,%2,%0; trapb"
|
||||
: "=&f" (tmp) : "f" (1.0), "f" (M_PI));
|
||||
/* And store it back. */
|
||||
__ieee_set_fp_control (tmp);
|
||||
|
||||
/* Success. */
|
||||
return 0;
|
||||
|
||||
@ -1,4 +1,4 @@
|
||||
/* Copyright (C) 1998 Free Software Foundation, Inc.
|
||||
/* Copyright (C) 1998, 2000 Free Software Foundation, Inc.
|
||||
This file is part of the GNU C Library.
|
||||
Contributed by Richard Henderson.
|
||||
|
||||
@ -19,34 +19,30 @@
|
||||
|
||||
#include <math.h>
|
||||
|
||||
/* Use the -inf rounding mode conversion instructions to implement
|
||||
ceil, via something akin to -floor(-x). This is much faster than
|
||||
playing with the fpcr to achieve +inf rounding mode. */
|
||||
|
||||
double
|
||||
__ceil (double x)
|
||||
{
|
||||
if (x != 0 && fabs (x) < 9007199254740992.0) /* 1 << DBL_MANT_DIG */
|
||||
if (isless (fabs (x), 9007199254740992.0)) /* 1 << DBL_MANT_DIG */
|
||||
{
|
||||
double tmp1;
|
||||
unsigned long fpcr0, fpcr1;
|
||||
unsigned long pinf = 3UL << 58;
|
||||
double tmp1, new_x;
|
||||
|
||||
/* Set round to +inf. */
|
||||
__asm __volatile("excb; mf_fpcr %0" : "=f"(fpcr0));
|
||||
__asm __volatile("mt_fpcr %0; excb" : : "f"(fpcr0 | pinf));
|
||||
|
||||
/* Calculate! */
|
||||
new_x = -x;
|
||||
__asm (
|
||||
#ifdef _IEEE_FP_INEXACT
|
||||
__asm("cvttq/svid %2,%1\n\tcvtqt/suid %1,%0"
|
||||
: "=f"(x), "=&f"(tmp1)
|
||||
: "f"(x));
|
||||
"cvttq/svim %2,%1\n\t"
|
||||
#else
|
||||
__asm("cvttq/svd %2,%1\n\tcvtqt/d %1,%0"
|
||||
: "=f"(x), "=&f"(tmp1)
|
||||
: "f"(x));
|
||||
"cvttq/svm %2,%1\n\t"
|
||||
#endif
|
||||
"cvtqt/m %1,%0\n\t"
|
||||
: "=f"(new_x), "=&f"(tmp1)
|
||||
: "f"(new_x));
|
||||
|
||||
/* 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));
|
||||
/* Fix up the negation we did above, as well as handling -0 properly. */
|
||||
x = copysign(new_x, x);
|
||||
}
|
||||
return x;
|
||||
}
|
||||
|
||||
@ -1,4 +1,4 @@
|
||||
/* Copyright (C) 1998 Free Software Foundation, Inc.
|
||||
/* Copyright (C) 1998, 2000 Free Software Foundation, Inc.
|
||||
This file is part of the GNU C Library.
|
||||
Contributed by Richard Henderson.
|
||||
|
||||
@ -19,39 +19,35 @@
|
||||
|
||||
#include <math.h>
|
||||
|
||||
/* Use the -inf rounding mode conversion instructions to implement
|
||||
ceil, via something akin to -floor(-x). This is much faster than
|
||||
playing with the fpcr to achieve +inf rounding mode. */
|
||||
|
||||
float
|
||||
__ceilf (float x)
|
||||
{
|
||||
if (x != 0 && fabsf (x) < 16777216.0f) /* 1 << FLT_MANT_DIG */
|
||||
if (isless (fabsf (x), 16777216.0f)) /* 1 << FLT_MANT_DIG */
|
||||
{
|
||||
float tmp1, tmp2;
|
||||
unsigned long fpcr0, fpcr1;
|
||||
unsigned long pinf = 3UL << 58;
|
||||
/* 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. */
|
||||
|
||||
/* 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. */
|
||||
float tmp1, tmp2, new_x;
|
||||
|
||||
new_x = -x;
|
||||
__asm ("cvtst/s %3,%2\n\t"
|
||||
#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));
|
||||
"cvttq/svim %2,%1\n\t"
|
||||
#else
|
||||
__asm("cvtst/s %3,%2\n\tcvttq/svd %2,%1\n\tcvtqt/d %1,%0"
|
||||
: "=f"(x), "=&f"(tmp1), "=&f"(tmp2)
|
||||
: "f"(x));
|
||||
"cvttq/svm %2,%1\n\t"
|
||||
#endif
|
||||
"cvtqt/m %1,%0\n\t"
|
||||
: "=f"(new_x), "=&f"(tmp1), "=&f"(tmp2)
|
||||
: "f"(new_x));
|
||||
|
||||
/* 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));
|
||||
/* Fix up the negation we did above, as well as handling -0 properly. */
|
||||
x = copysignf(new_x, x);
|
||||
}
|
||||
return x;
|
||||
}
|
||||
|
||||
33
sysdeps/alpha/fpu/s_copysign.c
Normal file
33
sysdeps/alpha/fpu/s_copysign.c
Normal file
@ -0,0 +1,33 @@
|
||||
/* Copyright (C) 2000 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
|
||||
__copysign (double x, double y)
|
||||
{
|
||||
__asm ("cpys %1, %2, %0" : "=f" (x) : "f" (y), "f" (x));
|
||||
return x;
|
||||
}
|
||||
|
||||
weak_alias (__copysign, copysign)
|
||||
#ifdef NO_LONG_DOUBLE
|
||||
strong_alias (__copysign, __copysignl)
|
||||
weak_alias (__copysign, copysignl)
|
||||
#endif
|
||||
29
sysdeps/alpha/fpu/s_copysignf.c
Normal file
29
sysdeps/alpha/fpu/s_copysignf.c
Normal file
@ -0,0 +1,29 @@
|
||||
/* Copyright (C) 2000 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
|
||||
__copysignf (float x, float y)
|
||||
{
|
||||
__asm ("cpys %1, %2, %0" : "=f" (x) : "f" (y), "f" (x));
|
||||
return x;
|
||||
}
|
||||
|
||||
weak_alias (__copysignf, copysignf)
|
||||
37
sysdeps/alpha/fpu/s_fabs.c
Normal file
37
sysdeps/alpha/fpu/s_fabs.c
Normal file
@ -0,0 +1,37 @@
|
||||
/* Copyright (C) 2000 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
|
||||
__fabs (double x)
|
||||
{
|
||||
#if __GNUC_PREREQ (2, 8)
|
||||
return __builtin_fabs (x);
|
||||
#else
|
||||
__asm ("cpys $f31, %1, %0" : "=f" (x) : "f" (x));
|
||||
return x;
|
||||
#endif
|
||||
}
|
||||
|
||||
weak_alias (__fabs, fabs)
|
||||
#ifdef NO_LONG_DOUBLE
|
||||
strong_alias (__fabs, __fabsl)
|
||||
weak_alias (__fabs, fabsl)
|
||||
#endif
|
||||
33
sysdeps/alpha/fpu/s_fabsf.c
Normal file
33
sysdeps/alpha/fpu/s_fabsf.c
Normal file
@ -0,0 +1,33 @@
|
||||
/* Copyright (C) 2000 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
|
||||
__fabsf (float x)
|
||||
{
|
||||
#if __GNUC_PREREQ (2, 8)
|
||||
return __builtin_fabsf (x);
|
||||
#else
|
||||
__asm ("cpys $f31, %1, %0" : "=f" (x) : "f" (x));
|
||||
return x;
|
||||
#endif
|
||||
}
|
||||
|
||||
weak_alias (__fabsf, fabsf)
|
||||
@ -1,4 +1,4 @@
|
||||
/* Copyright (C) 1998, 1999 Free Software Foundation, Inc.
|
||||
/* Copyright (C) 1998, 1999, 2000 Free Software Foundation, Inc.
|
||||
This file is part of the GNU C Library.
|
||||
Contributed by Richard Henderson.
|
||||
|
||||
@ -27,10 +27,10 @@
|
||||
double
|
||||
__floor (double x)
|
||||
{
|
||||
/* Check not zero since floor(-0) == -0. */
|
||||
if (x != 0 && fabs (x) < 9007199254740992.0) /* 1 << DBL_MANT_DIG */
|
||||
if (isless (fabs (x), 9007199254740992.0)) /* 1 << DBL_MANT_DIG */
|
||||
{
|
||||
double __tmp1;
|
||||
double tmp1, new_x;
|
||||
|
||||
__asm (
|
||||
#ifdef _IEEE_FP_INEXACT
|
||||
"cvttq/svim %2,%1\n\t"
|
||||
@ -38,8 +38,12 @@ __floor (double x)
|
||||
"cvttq/svm %2,%1\n\t"
|
||||
#endif
|
||||
"cvtqt/m %1,%0\n\t"
|
||||
: "=f"(x), "=&f"(__tmp1)
|
||||
: "=f"(new_x), "=&f"(tmp1)
|
||||
: "f"(x));
|
||||
|
||||
/* floor(-0) == -0, and in general we'll always have the same
|
||||
sign as our input. */
|
||||
x = copysign(new_x, x);
|
||||
}
|
||||
return x;
|
||||
}
|
||||
|
||||
@ -1,4 +1,4 @@
|
||||
/* Copyright (C) 1998, 1999 Free Software Foundation, Inc.
|
||||
/* Copyright (C) 1998, 1999, 2000 Free Software Foundation, Inc.
|
||||
This file is part of the GNU C Library.
|
||||
Contributed by Richard Henderson.
|
||||
|
||||
@ -27,15 +27,14 @@
|
||||
float
|
||||
__floorf (float x)
|
||||
{
|
||||
/* Check not zero since floor(-0) == -0. */
|
||||
if (x != 0 && fabsf (x) < 16777216.0f) /* 1 << FLT_MANT_DIG */
|
||||
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;
|
||||
float tmp1, tmp2, new_x;
|
||||
|
||||
__asm ("cvtst/s %3,%2\n\t"
|
||||
#ifdef _IEEE_FP_INEXACT
|
||||
@ -44,8 +43,12 @@ __floorf (float x)
|
||||
"cvttq/svm %2,%1\n\t"
|
||||
#endif
|
||||
"cvtqt/m %1,%0\n\t"
|
||||
: "=f"(x), "=&f"(tmp1), "=&f"(tmp2)
|
||||
: "=f"(new_x), "=&f"(tmp1), "=&f"(tmp2)
|
||||
: "f"(x));
|
||||
|
||||
/* floor(-0) == -0, and in general we'll always have the same
|
||||
sign as our input. */
|
||||
x = copysignf(new_x, x);
|
||||
}
|
||||
return x;
|
||||
}
|
||||
|
||||
50
sysdeps/alpha/fpu/s_rint.c
Normal file
50
sysdeps/alpha/fpu/s_rint.c
Normal file
@ -0,0 +1,50 @@
|
||||
/* Copyright (C) 2000 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
|
||||
__rint (double x)
|
||||
{
|
||||
if (isless (fabs (x), 9007199254740992.0)) /* 1 << DBL_MANT_DIG */
|
||||
{
|
||||
double tmp1, new_x;
|
||||
__asm (
|
||||
#ifdef _IEEE_FP_INEXACT
|
||||
"cvttq/svid %2,%1\n\t"
|
||||
#else
|
||||
"cvttq/svd %2,%1\n\t"
|
||||
#endif
|
||||
"cvtqt/d %1,%0\n\t"
|
||||
: "=f"(new_x), "=&f"(tmp1)
|
||||
: "f"(x));
|
||||
|
||||
/* rint(-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 (__rint, rint)
|
||||
#ifdef NO_LONG_DOUBLE
|
||||
strong_alias (__rint, __rintl)
|
||||
weak_alias (__rint, rintl)
|
||||
#endif
|
||||
52
sysdeps/alpha/fpu/s_rintf.c
Normal file
52
sysdeps/alpha/fpu/s_rintf.c
Normal file
@ -0,0 +1,52 @@
|
||||
/* Copyright (C) 2000 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
|
||||
__rintf (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"
|
||||
#ifdef _IEEE_FP_INEXACT
|
||||
"cvttq/svid %2,%1\n\t"
|
||||
#else
|
||||
"cvttq/svd %2,%1\n\t"
|
||||
#endif
|
||||
"cvtqt/d %1,%0\n\t"
|
||||
: "=f"(new_x), "=&f"(tmp1), "=&f"(tmp2)
|
||||
: "f"(x));
|
||||
|
||||
/* rint(-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 (__rintf, rintf)
|
||||
Loading…
Reference in New Issue
Block a user