mirror of
https://github.com/crosstool-ng/crosstool-ng.git
synced 2024-12-22 22:22:23 +00:00
4b173045e9
Woo... It seems the glibc guys finally decided that tarballs were not deprecated, in fact. The patchset was vampirised from Gentoo (kudos, guys!), and applies to glibc+ports, so that's why it's been added as a patchset against ports, not against glibc.
271 lines
9.0 KiB
Diff
271 lines
9.0 KiB
Diff
http://sources.redhat.com/bugzilla/show_bug.cgi?id=5350
|
|
https://bugs.gentoo.org/264335
|
|
|
|
diff -durN glibc-2.10.1.orig/glibc-ports-2.10.1/sysdeps/alpha/fpu/s_ceil.c glibc-2.10.1/glibc-ports-2.10.1/sysdeps/alpha/fpu/s_ceil.c
|
|
--- glibc-2.10.1.orig/glibc-ports-2.10.1/sysdeps/alpha/fpu/s_ceil.c 2009-05-16 10:36:20.000000000 +0200
|
|
+++ glibc-2.10.1/glibc-ports-2.10.1/sysdeps/alpha/fpu/s_ceil.c 2009-11-13 00:50:59.000000000 +0100
|
|
@@ -27,20 +27,25 @@
|
|
double
|
|
__ceil (double x)
|
|
{
|
|
- double two52 = copysign (0x1.0p52, x);
|
|
- double r, tmp;
|
|
-
|
|
- __asm (
|
|
+ if (isless (fabs (x), 9007199254740992.0)) /* 1 << DBL_MANT_DIG */
|
|
+ {
|
|
+ double tmp1, new_x;
|
|
+
|
|
+ new_x = -x;
|
|
+ __asm (
|
|
#ifdef _IEEE_FP_INEXACT
|
|
- "addt/suim %2, %3, %1\n\tsubt/suim %1, %3, %0"
|
|
+ "cvttq/svim %2,%1\n\t"
|
|
#else
|
|
- "addt/sum %2, %3, %1\n\tsubt/sum %1, %3, %0"
|
|
+ "cvttq/svm %2,%1\n\t"
|
|
#endif
|
|
- : "=&f"(r), "=&f"(tmp)
|
|
- : "f"(-x), "f"(-two52));
|
|
+ "cvtqt/m %1,%0\n\t"
|
|
+ : "=f"(new_x), "=&f"(tmp1)
|
|
+ : "f"(new_x));
|
|
|
|
- /* Fix up the negation we did above, as well as handling -0 properly. */
|
|
- return copysign (r, x);
|
|
+ /* Fix up the negation we did above, as well as handling -0 properly. */
|
|
+ x = copysign(new_x, x);
|
|
+ }
|
|
+ return x;
|
|
}
|
|
|
|
weak_alias (__ceil, ceil)
|
|
diff -durN glibc-2.10.1.orig/glibc-ports-2.10.1/sysdeps/alpha/fpu/s_ceilf.c glibc-2.10.1/glibc-ports-2.10.1/sysdeps/alpha/fpu/s_ceilf.c
|
|
--- glibc-2.10.1.orig/glibc-ports-2.10.1/sysdeps/alpha/fpu/s_ceilf.c 2009-05-16 10:36:20.000000000 +0200
|
|
+++ glibc-2.10.1/glibc-ports-2.10.1/sysdeps/alpha/fpu/s_ceilf.c 2009-11-13 00:50:59.000000000 +0100
|
|
@@ -26,20 +26,30 @@
|
|
float
|
|
__ceilf (float x)
|
|
{
|
|
- float two23 = copysignf (0x1.0p23, x);
|
|
- float r, tmp;
|
|
-
|
|
- __asm (
|
|
+ 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;
|
|
+
|
|
+ new_x = -x;
|
|
+ __asm ("cvtst/s %3,%2\n\t"
|
|
#ifdef _IEEE_FP_INEXACT
|
|
- "adds/suim %2, %3, %1\n\tsubs/suim %1, %3, %0"
|
|
+ "cvttq/svim %2,%1\n\t"
|
|
#else
|
|
- "adds/sum %2, %3, %1\n\tsubs/sum %1, %3, %0"
|
|
+ "cvttq/svm %2,%1\n\t"
|
|
#endif
|
|
- : "=&f"(r), "=&f"(tmp)
|
|
- : "f"(-x), "f"(-two23));
|
|
+ "cvtqt/m %1,%0\n\t"
|
|
+ : "=f"(new_x), "=&f"(tmp1), "=&f"(tmp2)
|
|
+ : "f"(new_x));
|
|
|
|
- /* Fix up the negation we did above, as well as handling -0 properly. */
|
|
- return copysignf (r, x);
|
|
+ /* Fix up the negation we did above, as well as handling -0 properly. */
|
|
+ x = copysignf(new_x, x);
|
|
+ }
|
|
+ return x;
|
|
}
|
|
|
|
weak_alias (__ceilf, ceilf)
|
|
diff -durN glibc-2.10.1.orig/glibc-ports-2.10.1/sysdeps/alpha/fpu/s_floor.c glibc-2.10.1/glibc-ports-2.10.1/sysdeps/alpha/fpu/s_floor.c
|
|
--- glibc-2.10.1.orig/glibc-ports-2.10.1/sysdeps/alpha/fpu/s_floor.c 2009-05-16 10:36:20.000000000 +0200
|
|
+++ glibc-2.10.1/glibc-ports-2.10.1/sysdeps/alpha/fpu/s_floor.c 2009-11-13 00:50:59.000000000 +0100
|
|
@@ -21,26 +21,32 @@
|
|
#include <math_ldbl_opt.h>
|
|
|
|
|
|
-/* Use the -inf rounding mode conversion instructions to implement floor. */
|
|
+/* 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. */
|
|
|
|
double
|
|
__floor (double x)
|
|
{
|
|
- double two52 = copysign (0x1.0p52, x);
|
|
- double r, tmp;
|
|
-
|
|
- __asm (
|
|
+ if (isless (fabs (x), 9007199254740992.0)) /* 1 << DBL_MANT_DIG */
|
|
+ {
|
|
+ double tmp1, new_x;
|
|
+
|
|
+ __asm (
|
|
#ifdef _IEEE_FP_INEXACT
|
|
- "addt/suim %2, %3, %1\n\tsubt/suim %1, %3, %0"
|
|
+ "cvttq/svim %2,%1\n\t"
|
|
#else
|
|
- "addt/sum %2, %3, %1\n\tsubt/sum %1, %3, %0"
|
|
+ "cvttq/svm %2,%1\n\t"
|
|
#endif
|
|
- : "=&f"(r), "=&f"(tmp)
|
|
- : "f"(x), "f"(two52));
|
|
+ "cvtqt/m %1,%0\n\t"
|
|
+ : "=f"(new_x), "=&f"(tmp1)
|
|
+ : "f"(x));
|
|
|
|
- /* floor(-0) == -0, and in general we'll always have the same
|
|
- sign as our input. */
|
|
- return copysign (r, x);
|
|
+ /* floor(-0) == -0, and in general we'll always have the same
|
|
+ sign as our input. */
|
|
+ x = copysign(new_x, x);
|
|
+ }
|
|
+ return x;
|
|
}
|
|
|
|
weak_alias (__floor, floor)
|
|
diff -durN glibc-2.10.1.orig/glibc-ports-2.10.1/sysdeps/alpha/fpu/s_floorf.c glibc-2.10.1/glibc-ports-2.10.1/sysdeps/alpha/fpu/s_floorf.c
|
|
--- glibc-2.10.1.orig/glibc-ports-2.10.1/sysdeps/alpha/fpu/s_floorf.c 2009-05-16 10:36:20.000000000 +0200
|
|
+++ glibc-2.10.1/glibc-ports-2.10.1/sysdeps/alpha/fpu/s_floorf.c 2009-11-13 00:50:59.000000000 +0100
|
|
@@ -20,26 +20,37 @@
|
|
#include <math.h>
|
|
|
|
|
|
-/* Use the -inf rounding mode conversion instructions to implement floor. */
|
|
+/* 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. */
|
|
|
|
float
|
|
__floorf (float x)
|
|
{
|
|
- float two23 = copysignf (0x1.0p23, x);
|
|
- float r, tmp;
|
|
-
|
|
- __asm (
|
|
+ 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
|
|
- "adds/suim %2, %3, %1\n\tsubs/suim %1, %3, %0"
|
|
+ "cvttq/svim %2,%1\n\t"
|
|
#else
|
|
- "adds/sum %2, %3, %1\n\tsubs/sum %1, %3, %0"
|
|
+ "cvttq/svm %2,%1\n\t"
|
|
#endif
|
|
- : "=&f"(r), "=&f"(tmp)
|
|
- : "f"(x), "f"(two23));
|
|
+ "cvtqt/m %1,%0\n\t"
|
|
+ : "=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. */
|
|
- return copysignf (r, x);
|
|
+ /* floor(-0) == -0, and in general we'll always have the same
|
|
+ sign as our input. */
|
|
+ x = copysignf(new_x, x);
|
|
+ }
|
|
+ return x;
|
|
}
|
|
|
|
weak_alias (__floorf, floorf)
|
|
diff -durN glibc-2.10.1.orig/glibc-ports-2.10.1/sysdeps/alpha/fpu/s_rint.c glibc-2.10.1/glibc-ports-2.10.1/sysdeps/alpha/fpu/s_rint.c
|
|
--- glibc-2.10.1.orig/glibc-ports-2.10.1/sysdeps/alpha/fpu/s_rint.c 2009-05-16 10:36:20.000000000 +0200
|
|
+++ glibc-2.10.1/glibc-ports-2.10.1/sysdeps/alpha/fpu/s_rint.c 2009-11-13 00:50:59.000000000 +0100
|
|
@@ -24,15 +24,24 @@
|
|
double
|
|
__rint (double x)
|
|
{
|
|
- double two52 = copysign (0x1.0p52, x);
|
|
- double r;
|
|
-
|
|
- r = x + two52;
|
|
- r = r - two52;
|
|
+ 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. */
|
|
- return copysign (r, 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)
|
|
diff -durN glibc-2.10.1.orig/glibc-ports-2.10.1/sysdeps/alpha/fpu/s_rintf.c glibc-2.10.1/glibc-ports-2.10.1/sysdeps/alpha/fpu/s_rintf.c
|
|
--- glibc-2.10.1.orig/glibc-ports-2.10.1/sysdeps/alpha/fpu/s_rintf.c 2009-05-16 10:36:20.000000000 +0200
|
|
+++ glibc-2.10.1/glibc-ports-2.10.1/sysdeps/alpha/fpu/s_rintf.c 2009-11-13 00:50:59.000000000 +0100
|
|
@@ -23,15 +23,30 @@
|
|
float
|
|
__rintf (float x)
|
|
{
|
|
- float two23 = copysignf (0x1.0p23, x);
|
|
- float r;
|
|
+ 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. */
|
|
|
|
- r = x + two23;
|
|
- r = r - two23;
|
|
+ float tmp1, tmp2, new_x;
|
|
|
|
- /* rint(-0.1) == -0, and in general we'll always have the same sign
|
|
- as our input. */
|
|
- return copysign (r, 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)
|
|
diff -durN glibc-2.10.1.orig/ports/sysdeps/alpha/fpu/s_ceil.c glibc-2.10.1/ports/sysdeps/alpha/fpu/s_ceil.c
|
|
diff -durN glibc-2.10.1.orig/ports/sysdeps/alpha/fpu/s_ceilf.c glibc-2.10.1/ports/sysdeps/alpha/fpu/s_ceilf.c
|
|
diff -durN glibc-2.10.1.orig/ports/sysdeps/alpha/fpu/s_floor.c glibc-2.10.1/ports/sysdeps/alpha/fpu/s_floor.c
|
|
diff -durN glibc-2.10.1.orig/ports/sysdeps/alpha/fpu/s_floorf.c glibc-2.10.1/ports/sysdeps/alpha/fpu/s_floorf.c
|
|
diff -durN glibc-2.10.1.orig/ports/sysdeps/alpha/fpu/s_rint.c glibc-2.10.1/ports/sysdeps/alpha/fpu/s_rint.c
|
|
diff -durN glibc-2.10.1.orig/ports/sysdeps/alpha/fpu/s_rintf.c glibc-2.10.1/ports/sysdeps/alpha/fpu/s_rintf.c
|