summaryrefslogtreecommitdiff
path: root/sysdeps/alpha/fpu/bits/mathinline.h
diff options
context:
space:
mode:
authorRichard Henderson <rth@redhat.com>1998-08-23 04:10:13 +0000
committerRichard Henderson <rth@redhat.com>1998-08-23 04:10:13 +0000
commit7d84a06715a3d378eb3ee0ee5d770a5ea4566a87 (patch)
treee7af13d76db60df9a34bac32508ff6d5aca81ab2 /sysdeps/alpha/fpu/bits/mathinline.h
parentba3904fdf8e5ae26c1a931e09e3c7a3896593d92 (diff)
¸
1998-08-23 Ricahrd Henderson <rth@cygnus.com> * sysdeps/alpha/elf/crtbegin.S: Fix .prologue; no pv used. * sysdeps/alpha/elf/crtend.S: Likewise. * sysdeps/alpha/elf/start.S: Likewise. * sysdeps/unix/sysv/linux/alpha/rt_sigaction.S: Fix .prologue; non-standard pv usage. * sysdeps/unix/sysv/linux/alpha/brk.S: Use jmp macro for relaxation. * sysdeps/unix/sysv/linux/alpha/getitimer.S: Likewise. * sysdeps/unix/sysv/linux/alpha/getrusage.S: Likewise. * sysdeps/unix/sysv/linux/alpha/gettimeofday.S: Likewise. * sysdeps/unix/sysv/linux/alpha/ieee_get_fp_control.S: Likewise. * sysdeps/unix/sysv/linux/alpha/ieee_set_fp_control.S: Likewise. * sysdeps/unix/sysv/linux/alpha/select.S: Likewise. * sysdeps/unix/sysv/linux/alpha/setitimer.S: Likewise. * sysdeps/unix/sysv/linux/alpha/settimeofday.S: Likewise. * sysdeps/unix/sysv/linux/alpha/utimes.S: Likewise. * sysdeps/unix/sysv/linux/alpha/wait4.S: Likewise. * sysdeps/alpha/fpu/e_sqrt.c: Use the asm version when the input is a finite non-denormal, deferring to the full IEEE version otherwise. * 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.
Diffstat (limited to 'sysdeps/alpha/fpu/bits/mathinline.h')
-rw-r--r--sysdeps/alpha/fpu/bits/mathinline.h28
1 files changed, 18 insertions, 10 deletions
diff --git a/sysdeps/alpha/fpu/bits/mathinline.h b/sysdeps/alpha/fpu/bits/mathinline.h
index 492d9f18a7..681ea70cf9 100644
--- a/sysdeps/alpha/fpu/bits/mathinline.h
+++ b/sysdeps/alpha/fpu/bits/mathinline.h
@@ -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;