80 lines
2.8 KiB
Diff
80 lines
2.8 KiB
Diff
|
diff --git a/sysdeps/ieee754/ldbl-96/Makefile b/sysdeps/ieee754/ldbl-96/Makefile
|
||
|
|
index 995e90d6da..318628aed6 100644
|
||
|
|
--- a/sysdeps/ieee754/ldbl-96/Makefile
|
||
|
|
+++ b/sysdeps/ieee754/ldbl-96/Makefile
|
||
|
|
@@ -17,5 +17,6 @@
|
||
|
|
# <https://www.gnu.org/licenses/>.
|
||
|
|
|
||
|
|
ifeq ($(subdir),math)
|
||
|
|
-tests += test-canonical-ldbl-96 test-totalorderl-ldbl-96
|
||
|
|
+tests += test-canonical-ldbl-96 test-totalorderl-ldbl-96 test-sinl-pseudo
|
||
|
|
+CFLAGS-test-sinl-pseudo.c += -fstack-protector-all
|
||
|
|
endif
|
||
|
|
diff --git a/sysdeps/ieee754/ldbl-96/e_rem_pio2l.c b/sysdeps/ieee754/ldbl-96/e_rem_pio2l.c
|
||
|
|
index 5f742321ae..bcdf20179f 100644
|
||
|
|
--- a/sysdeps/ieee754/ldbl-96/e_rem_pio2l.c
|
||
|
|
+++ b/sysdeps/ieee754/ldbl-96/e_rem_pio2l.c
|
||
|
|
@@ -210,6 +210,18 @@ __ieee754_rem_pio2l (long double x, long double *y)
|
||
|
|
return 0;
|
||
|
|
}
|
||
|
|
|
||
|
|
+ if ((i0 & 0x80000000) == 0)
|
||
|
|
+ {
|
||
|
|
+ /* Pseudo-zero and unnormal representations are not valid
|
||
|
|
+ representations of long double. We need to avoid stack
|
||
|
|
+ corruption in __kernel_rem_pio2, which expects input in a
|
||
|
|
+ particular normal form, but those representations do not need
|
||
|
|
+ to be consistently handled like any particular floating-point
|
||
|
|
+ value. */
|
||
|
|
+ y[1] = y[0] = __builtin_nanl ("");
|
||
|
|
+ return 0;
|
||
|
|
+ }
|
||
|
|
+
|
||
|
|
/* Split the 64 bits of the mantissa into three 24-bit integers
|
||
|
|
stored in a double array. */
|
||
|
|
exp = j0 - 23;
|
||
|
|
--- /dev/null
|
||
|
|
+++ b/sysdeps/ieee754/ldbl-96/test-sinl-pseudo.c
|
||
|
|
@@ -0,0 +1,41 @@
|
||
|
|
+/* Test sinl for pseudo-zeros and unnormals for ldbl-96 (bug 25487).
|
||
|
|
+ Copyright (C) 2020 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, see
|
||
|
|
+ <https://www.gnu.org/licenses/>. */
|
||
|
|
+
|
||
|
|
+#include <math.h>
|
||
|
|
+#include <math_ldbl.h>
|
||
|
|
+#include <stdint.h>
|
||
|
|
+
|
||
|
|
+static int
|
||
|
|
+do_test (void)
|
||
|
|
+{
|
||
|
|
+ for (int i = 0; i < 64; i++)
|
||
|
|
+ {
|
||
|
|
+ uint64_t sig = i == 63 ? 0 : 1ULL << i;
|
||
|
|
+ long double ld;
|
||
|
|
+ SET_LDOUBLE_WORDS (ld, 0x4141,
|
||
|
|
+ sig >> 32, sig & 0xffffffffULL);
|
||
|
|
+ /* The requirement is that no stack overflow occurs when the
|
||
|
|
+ pseudo-zero or unnormal goes through range reduction. */
|
||
|
|
+ volatile long double ldr;
|
||
|
|
+ ldr = sinl (ld);
|
||
|
|
+ (void) ldr;
|
||
|
|
+ }
|
||
|
|
+ return 0;
|
||
|
|
+}
|
||
|
|
+
|
||
|
|
+#include <support/test-driver.c>
|