aboutsummaryrefslogtreecommitdiff
path: root/lib/msun/ld80
diff options
context:
space:
mode:
Diffstat (limited to 'lib/msun/ld80')
-rw-r--r--lib/msun/ld80/e_lgammal_r.c2
-rw-r--r--lib/msun/ld80/e_powl.c26
-rw-r--r--lib/msun/ld80/e_rem_pio2l.h4
-rw-r--r--lib/msun/ld80/invtrig.c1
-rw-r--r--lib/msun/ld80/k_cosl.c2
-rw-r--r--lib/msun/ld80/k_expl.h1
-rw-r--r--lib/msun/ld80/k_sinl.c2
-rw-r--r--lib/msun/ld80/k_tanl.c3
-rw-r--r--lib/msun/ld80/s_cexpl.c1
-rw-r--r--lib/msun/ld80/s_erfl.c2
-rw-r--r--lib/msun/ld80/s_exp2l.c1
-rw-r--r--lib/msun/ld80/s_expl.c1
-rw-r--r--lib/msun/ld80/s_logl.c3
13 files changed, 14 insertions, 35 deletions
diff --git a/lib/msun/ld80/e_lgammal_r.c b/lib/msun/ld80/e_lgammal_r.c
index cfd73b408b87..9f8b8f0f98bb 100644
--- a/lib/msun/ld80/e_lgammal_r.c
+++ b/lib/msun/ld80/e_lgammal_r.c
@@ -1,4 +1,3 @@
-/* @(#)e_lgamma_r.c 1.3 95/01/18 */
/*
* ====================================================
* Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
@@ -10,7 +9,6 @@
* ====================================================
*/
-#include <sys/cdefs.h>
/*
* See e_lgamma_r.c for complete comments.
*
diff --git a/lib/msun/ld80/e_powl.c b/lib/msun/ld80/e_powl.c
index 2daf26aab93b..d04d77a5ee98 100644
--- a/lib/msun/ld80/e_powl.c
+++ b/lib/msun/ld80/e_powl.c
@@ -14,7 +14,6 @@
* OR IN CONNECTION WITH THE USE OR PERFORMANCE OF THIS SOFTWARE.
*/
-#include <sys/cdefs.h>
#include <math.h>
#include "math_private.h"
@@ -24,10 +23,10 @@
* P[0] x^n + P[1] x^(n-1) + ... + P[n]
*/
static inline long double
-__polevll(long double x, long double *PP, int n)
+__polevll(long double x, const long double *PP, int n)
{
long double y;
- long double *P;
+ const long double *P;
P = PP;
y = *P++;
@@ -43,10 +42,10 @@ __polevll(long double x, long double *PP, int n)
* x^n + P[0] x^(n-1) + P[1] x^(n-2) + ... + P[n]
*/
static inline long double
-__p1evll(long double x, long double *PP, int n)
+__p1evll(long double x, const long double *PP, int n)
{
long double y;
- long double *P;
+ const long double *P;
P = PP;
n -= 1;
@@ -114,7 +113,6 @@ __p1evll(long double x, long double *PP, int n)
*
*/
-#include <sys/cdefs.h>
#include <float.h>
#include <math.h>
@@ -128,13 +126,13 @@ __p1evll(long double x, long double *PP, int n)
/* log(1+x) = x - .5x^2 + x^3 * P(z)/Q(z)
* on the domain 2^(-1/32) - 1 <= x <= 2^(1/32) - 1
*/
-static long double P[] = {
+static const long double P[] = {
8.3319510773868690346226E-4L,
4.9000050881978028599627E-1L,
1.7500123722550302671919E0L,
1.4000100839971580279335E0L,
};
-static long double Q[] = {
+static const long double Q[] = {
/* 1.0000000000000000000000E0L,*/
5.2500282295834889175431E0L,
8.4000598057587009834666E0L,
@@ -143,7 +141,7 @@ static long double Q[] = {
/* A[i] = 2^(-i/32), rounded to IEEE long double precision.
* If i is even, A[i] + B[i/2] gives additional accuracy.
*/
-static long double A[33] = {
+static const long double A[33] = {
1.0000000000000000000000E0L,
9.7857206208770013448287E-1L,
9.5760328069857364691013E-1L,
@@ -178,7 +176,7 @@ static long double A[33] = {
5.1094857432705833910408E-1L,
5.0000000000000000000000E-1L,
};
-static long double B[17] = {
+static const long double B[17] = {
0.0000000000000000000000E0L,
2.6176170809902549338711E-20L,
-1.0126791927256478897086E-20L,
@@ -201,7 +199,7 @@ static long double B[17] = {
/* 2^x = 1 + x P(x),
* on the interval -1/32 <= x <= 0
*/
-static long double R[] = {
+static const long double R[] = {
1.5089970579127659901157E-5L,
1.5402715328927013076125E-4L,
1.3333556028915671091390E-3L,
@@ -232,13 +230,13 @@ static long double R[] = {
static const long double MAXLOGL = 1.1356523406294143949492E4L;
static const long double MINLOGL = -1.13994985314888605586758E4L;
static const long double LOGE2L = 6.9314718055994530941723E-1L;
-static volatile long double z;
-static long double w, W, Wa, Wb, ya, yb, u;
+static _Thread_local volatile long double z;
+static _Thread_local long double w, W, Wa, Wb, ya, yb, u;
static const long double huge = 0x1p10000L;
#if 0 /* XXX Prevent gcc from erroneously constant folding this. */
static const long double twom10000 = 0x1p-10000L;
#else
-static volatile long double twom10000 = 0x1p-10000L;
+static _Thread_local volatile long double twom10000 = 0x1p-10000L;
#endif
static long double reducl( long double );
diff --git a/lib/msun/ld80/e_rem_pio2l.h b/lib/msun/ld80/e_rem_pio2l.h
index d9e4d3a0941f..b39e63b72afc 100644
--- a/lib/msun/ld80/e_rem_pio2l.h
+++ b/lib/msun/ld80/e_rem_pio2l.h
@@ -1,4 +1,3 @@
-/* From: @(#)e_rem_pio2.c 1.4 95/01/18 */
/*
* ====================================================
* Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
@@ -13,7 +12,6 @@
* Optimized by Bruce D. Evans.
*/
-#include <sys/cdefs.h>
/* ld80 version of __ieee754_rem_pio2l(x,y)
*
* return the remainder of x rem pi/2 in y[0]+y[1]
@@ -68,7 +66,7 @@ pio2_2t = 6.36831716351095013979e-25L, /* 0xc51701b839a25205.0p-144 */
pio2_3t = -2.75299651904407171810e-37L; /* -0xbb5bf6c7ddd660ce.0p-185 */
#endif
-static inline __always_inline int
+static __always_inline int
__ieee754_rem_pio2l(long double x, long double *y)
{
union IEEEl2bits u,u1;
diff --git a/lib/msun/ld80/invtrig.c b/lib/msun/ld80/invtrig.c
index 0e2d49546391..76eb681454d3 100644
--- a/lib/msun/ld80/invtrig.c
+++ b/lib/msun/ld80/invtrig.c
@@ -26,7 +26,6 @@
* SUCH DAMAGE.
*/
-#include <sys/cdefs.h>
#include "invtrig.h"
/*
diff --git a/lib/msun/ld80/k_cosl.c b/lib/msun/ld80/k_cosl.c
index e76d429b90ef..371c360aa2ce 100644
--- a/lib/msun/ld80/k_cosl.c
+++ b/lib/msun/ld80/k_cosl.c
@@ -1,4 +1,3 @@
-/* From: @(#)k_cos.c 1.3 95/01/18 */
/*
* ====================================================
* Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
@@ -11,7 +10,6 @@
* ====================================================
*/
-#include <sys/cdefs.h>
/*
* ld80 version of k_cos.c. See ../src/k_cos.c for most comments.
*/
diff --git a/lib/msun/ld80/k_expl.h b/lib/msun/ld80/k_expl.h
index 8c3617e4b177..fabf5a792205 100644
--- a/lib/msun/ld80/k_expl.h
+++ b/lib/msun/ld80/k_expl.h
@@ -30,7 +30,6 @@
* Optimized by Bruce D. Evans.
*/
-#include <sys/cdefs.h>
/*
* See s_expl.c for more comments about __k_expl().
*
diff --git a/lib/msun/ld80/k_sinl.c b/lib/msun/ld80/k_sinl.c
index 21c1946b5e1b..ffc00be23f1c 100644
--- a/lib/msun/ld80/k_sinl.c
+++ b/lib/msun/ld80/k_sinl.c
@@ -1,4 +1,3 @@
-/* From: @(#)k_sin.c 1.3 95/01/18 */
/*
* ====================================================
* Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
@@ -11,7 +10,6 @@
* ====================================================
*/
-#include <sys/cdefs.h>
/*
* ld80 version of k_sin.c. See ../src/k_sin.c for most comments.
*/
diff --git a/lib/msun/ld80/k_tanl.c b/lib/msun/ld80/k_tanl.c
index 4965c669979f..b3a4d9e05255 100644
--- a/lib/msun/ld80/k_tanl.c
+++ b/lib/msun/ld80/k_tanl.c
@@ -1,5 +1,3 @@
-/* From: @(#)k_tan.c 1.5 04/04/22 SMI */
-
/*
* ====================================================
* Copyright 2004 Sun Microsystems, Inc. All Rights Reserved.
@@ -11,7 +9,6 @@
* ====================================================
*/
-#include <sys/cdefs.h>
/*
* ld80 version of k_tan.c. See ../src/k_tan.c for most comments.
*/
diff --git a/lib/msun/ld80/s_cexpl.c b/lib/msun/ld80/s_cexpl.c
index 81a4d727561b..cbe9dcd64099 100644
--- a/lib/msun/ld80/s_cexpl.c
+++ b/lib/msun/ld80/s_cexpl.c
@@ -28,7 +28,6 @@
* src/s_cexp.c converted to long double complex by Steven G. Kargl
*/
-#include <sys/cdefs.h>
#include <complex.h>
#include <float.h>
#ifdef __i386__
diff --git a/lib/msun/ld80/s_erfl.c b/lib/msun/ld80/s_erfl.c
index 362b4f07d6bc..047077ca4470 100644
--- a/lib/msun/ld80/s_erfl.c
+++ b/lib/msun/ld80/s_erfl.c
@@ -1,4 +1,3 @@
-/* @(#)s_erf.c 5.1 93/09/24 */
/*
* ====================================================
* Copyright (C) 1993 by Sun Microsystems, Inc. All rights reserved.
@@ -10,7 +9,6 @@
* ====================================================
*/
-#include <sys/cdefs.h>
/*
* See s_erf.c for complete comments.
*
diff --git a/lib/msun/ld80/s_exp2l.c b/lib/msun/ld80/s_exp2l.c
index 94aecf1843cd..9d979994ece8 100644
--- a/lib/msun/ld80/s_exp2l.c
+++ b/lib/msun/ld80/s_exp2l.c
@@ -26,7 +26,6 @@
* SUCH DAMAGE.
*/
-#include <sys/cdefs.h>
#include <float.h>
#include <stdint.h>
diff --git a/lib/msun/ld80/s_expl.c b/lib/msun/ld80/s_expl.c
index 065e02fae1b0..03d7b366151a 100644
--- a/lib/msun/ld80/s_expl.c
+++ b/lib/msun/ld80/s_expl.c
@@ -28,7 +28,6 @@
* Optimized by Bruce D. Evans.
*/
-#include <sys/cdefs.h>
/**
* Compute the exponential of x for Intel 80-bit format. This is based on:
*
diff --git a/lib/msun/ld80/s_logl.c b/lib/msun/ld80/s_logl.c
index 459374d7d164..7141abe086d5 100644
--- a/lib/msun/ld80/s_logl.c
+++ b/lib/msun/ld80/s_logl.c
@@ -26,7 +26,6 @@
* THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
*/
-#include <sys/cdefs.h>
/**
* Implementation of the natural logarithm of x for Intel 80-bit format.
*
@@ -445,7 +444,7 @@ struct ld {
#endif
#ifdef STRUCT_RETURN
-static inline __always_inline void
+static __always_inline void
k_logl(long double x, struct ld *rp)
#else
long double