| /glibc-2.36/sysdeps/ieee754/dbl-64/ |
| D | e_asin.c | 47 __ieee754_asin(double x){ in __ieee754_asin() argument
52 u.x = x; in __ieee754_asin()
58 math_check_force_underflow (x); in __ieee754_asin()
59 return x; /* for x->0 => sin(x)=x */ in __ieee754_asin()
64 x2 = x*x; in __ieee754_asin()
65 t = (((((f6*x2 + f5)*x2 + f4)*x2 + f3)*x2 + f2)*x2 + f1)*(x2*x); in __ieee754_asin()
66 res = x+t; /* res=arcsin(x) according to Taylor series */ in __ieee754_asin()
74 if (m>0) xx = x - asncs.x[n]; in __ieee754_asin()
75 else xx = -x - asncs.x[n]; in __ieee754_asin()
76 t = asncs.x[n+1]*xx; in __ieee754_asin()
[all …]
|
| D | e_remainder.c | 45 __ieee754_remainder (double x, double y) in __ieee754_remainder() argument
50 u.x = x; in __ieee754_remainder()
51 t.x = y; in __ieee754_remainder()
60 return x; in __ieee754_remainder()
63 z = x / t.x; in __ieee754_remainder()
65 d = (z + big.x) - big.x; in __ieee754_remainder()
66 xx = (x - d * v.x) - d * (t.x - v.x); in __ieee754_remainder()
68 return (xx != 0) ? xx : ((x > 0) ? ZERO.x : nZERO.x); in __ieee754_remainder()
71 if (fabs (xx) > 0.5 * t.x) in __ieee754_remainder()
72 return (z > d) ? xx - t.x : xx + t.x; in __ieee754_remainder()
[all …]
|
| D | s_sin.c | 60 #define TAYLOR_SIN(xx, x, dx) \ argument
62 double t = ((POLYNOMIAL (xx) * (x) - 0.5 * (dx)) * (xx) + (dx)); \
63 double res = (x) + t; \
70 sn = __sincostab.x[k]; \
71 ssn = __sincostab.x[k + 1]; \
72 cs = __sincostab.x[k + 2]; \
73 ccs = __sincostab.x[k + 3]; \
83 double x[440]; member
93 int __branred (double x, double *a, double *aa);
100 do_cos (double x, double dx) in do_cos() argument
[all …]
|
| D | dla.h | 47 #define EADD(x,y,z,zz) \ argument
48 z=(x)+(y); zz=(fabs(x)>fabs(y)) ? (((x)-(z))+(y)) : (((y)-(z))+(x));
55 #define ESUB(x,y,z,zz) \ argument
56 z=(x)-(y); zz=(fabs(x)>fabs(y)) ? (((x)-(z))-(y)) : ((x)-((y)+(z)));
60 # define DLA_FMS(x, y, z) __builtin_fma (x, y, -(z)) argument
69 # define EMULV(x, y, z, zz) \ argument
70 z = x * y; zz = DLA_FMS (x, y, z);
72 # define EMULV(x, y, z, zz) \ argument
73 ({ __typeof__ (x) __p, hx, tx, hy, ty; \
74 __p = CN * (x); hx = ((x) - __p) + __p; tx = (x) - hx; \
[all …]
|
| /glibc-2.36/math/ |
| D | test-tgmath.c | 184 #define x dx macro
193 #define x fx macro
203 #define x lx macro
234 a = cos (cos (x)); in F()
236 a = sin (sin (x)); in F()
238 a = tan (tan (x)); in F()
240 c = atan2 (atan2 (a, c), atan2 (b, x)); in F()
241 a = cosh (cosh (x)); in F()
243 a = sinh (sinh (x)); in F()
245 a = tanh (tanh (x)); in F()
[all …]
|
| D | s_csqrt_template.c | 26 M_DECL_FUNC (__csqrt) (CFLOAT x) in M_DECL_FUNC()
29 int rcls = fpclassify (__real__ x); in M_DECL_FUNC()
30 int icls = fpclassify (__imag__ x); in M_DECL_FUNC()
37 __imag__ res = __imag__ x; in M_DECL_FUNC()
41 if (__real__ x < 0) in M_DECL_FUNC()
44 __imag__ res = M_COPYSIGN (M_HUGE_VAL, __imag__ x); in M_DECL_FUNC()
48 __real__ res = __real__ x; in M_DECL_FUNC()
50 ? M_NAN : M_COPYSIGN (0, __imag__ x)); in M_DECL_FUNC()
63 if (__real__ x < 0) in M_DECL_FUNC()
66 __imag__ res = M_COPYSIGN (M_SQRT (-__real__ x), __imag__ x); in M_DECL_FUNC()
[all …]
|
| /glibc-2.36/string/ |
| D | endian.h | 39 # define htobe16(x) __bswap_16 (x) argument
40 # define htole16(x) __uint16_identity (x) argument
41 # define be16toh(x) __bswap_16 (x) argument
42 # define le16toh(x) __uint16_identity (x) argument
44 # define htobe32(x) __bswap_32 (x) argument
45 # define htole32(x) __uint32_identity (x) argument
46 # define be32toh(x) __bswap_32 (x) argument
47 # define le32toh(x) __uint32_identity (x) argument
49 # define htobe64(x) __bswap_64 (x) argument
50 # define htole64(x) __uint64_identity (x) argument
[all …]
|
| /glibc-2.36/sysdeps/ieee754/ldbl-128ibm/ |
| D | s_log1pl.c | 126 long double x, y, z, r, s; in __log1pl() local
142 x = xm1; in __log1pl()
144 x = xm1 + 1.0L; in __log1pl()
147 if (x <= 0.0L) in __log1pl()
149 if (x == 0.0L) in __log1pl()
152 return (zero / (x - x)); in __log1pl()
158 x = __frexpl (x, &e); in __log1pl()
164 if (x < sqrth) in __log1pl()
167 z = x - 0.5L; in __log1pl()
172 z = x - 0.5L; in __log1pl()
[all …]
|
| D | s_cbrtl.c | 34 __cbrtl (long double x) in __cbrtl() argument
39 if (!isfinite (x)) in __cbrtl()
40 return x + x; in __cbrtl()
42 if (x == 0) in __cbrtl()
43 return (x); in __cbrtl()
45 if (x > 0) in __cbrtl()
50 x = -x; in __cbrtl()
53 z = x; in __cbrtl()
55 x = __frexpl (x, &e); in __cbrtl()
59 x = ((((1.3584464340920900529734e-1L * x in __cbrtl()
[all …]
|
| D | e_log2l.c | 140 neval (long double x, const long double *p, int n) in neval() argument
148 y = y * x + *p--; in neval()
158 deval (long double x, const long double *p, int n) in deval() argument
163 y = x + *p--; in deval()
166 y = y * x + *p--; in deval()
175 __ieee754_log2l (long double x) in __ieee754_log2l() argument
184 xhi = ldbl_high (x); in __ieee754_log2l()
187 return (-1.0L / fabsl (x)); /* log2l(+-0)=-inf */ in __ieee754_log2l()
189 return (x - x) / (x - x); in __ieee754_log2l()
191 return (x + x); in __ieee754_log2l()
[all …]
|
| D | e_lgammal_r.c | 675 neval (long double x, const long double *p, int n) in neval() argument
683 y = y * x + *p--; in neval()
693 deval (long double x, const long double *p, int n) in deval() argument
698 y = x + *p--; in deval()
701 y = y * x + *p--; in deval()
709 __ieee754_lgammal_r (long double x, int *signgamp) in __ieee754_lgammal_r() argument
716 if (! isfinite (x)) in __ieee754_lgammal_r()
717 return x * x; in __ieee754_lgammal_r()
719 if (x == 0) in __ieee754_lgammal_r()
721 if (signbit (x)) in __ieee754_lgammal_r()
[all …]
|
| D | e_log10l.c | 146 neval (long double x, const long double *p, int n) in neval() argument
154 y = y * x + *p--; in neval()
164 deval (long double x, const long double *p, int n) in deval() argument
169 y = x + *p--; in deval()
172 y = y * x + *p--; in deval()
181 __ieee754_log10l (long double x) in __ieee754_log10l() argument
190 xhi = ldbl_high (x); in __ieee754_log10l()
193 return (-1.0L / fabsl (x)); /* log10l(+-0)=-inf */ in __ieee754_log10l()
195 return (x - x) / (x - x); in __ieee754_log10l()
197 return (x + x); in __ieee754_log10l()
[all …]
|
| /glibc-2.36/bits/ |
| D | floatn-common.h | 83 # define __f16(x) ((_Float16) x##f) argument
85 # define __f16(x) x##f16 argument
91 # define __f32(x) x##f argument
93 # define __f32(x) x##f32 argument
100 # define __f64(x) x##l argument
102 # define __f64(x) x argument
105 # define __f64(x) x##f64 argument
111 # define __f32x(x) x argument
113 # define __f32x(x) x##f32x argument
120 # define __f64x(x) x##l argument
[all …]
|
| /glibc-2.36/nptl/ |
| D | tst-audit-threads.h | 26 #define NUM(x, y) CONCAT (x, y) argument
28 #define FUNC10(x) \ argument
29 FUNC (NUM (x, 0)); \
30 FUNC (NUM (x, 1)); \
31 FUNC (NUM (x, 2)); \
32 FUNC (NUM (x, 3)); \
33 FUNC (NUM (x, 4)); \
34 FUNC (NUM (x, 5)); \
35 FUNC (NUM (x, 6)); \
36 FUNC (NUM (x, 7)); \
[all …]
|
| /glibc-2.36/sysdeps/generic/ |
| D | a.out.h | 81 #define N_BADMAG(x) (N_MAGIC(x) != OMAGIC && N_MAGIC(x) != ZMAGIC) argument
82 #define N_ADDRADJ(x) (N_MAGIC(x) == ZMAGIC ? SEGMENT_SIZE : 0) argument
83 #define N_DATOFF(x) (N_TXTOFF(x) + (x).a_text - N_ADDRADJ(x)) argument
84 #define N_TRELOFF(x) (N_DATOFF(x) + (x).a_data + (x).a_shdata) argument
85 #define N_SYMOFF(x) (N_DRELOFF(x) + (x).a_drsize + (x).a_shdrsize) argument
86 #define N_TXTADDR(x) SEGMENT_SIZE argument
128 #define N_BADMAG(x) \ argument
129 (N_MAGIC(x) != OMAGIC && N_MAGIC(x) != NMAGIC \
130 && N_MAGIC(x) != ZMAGIC)
133 #define _N_BADMAG(x) \ argument
[all …]
|
| /glibc-2.36/sysdeps/ieee754/ldbl-128/ |
| D | s_cbrtl.c | 68 __cbrtl (_Float128 x) in __cbrtl() argument
73 if (!isfinite (x)) in __cbrtl()
74 return x + x; in __cbrtl()
76 if (x == 0) in __cbrtl()
77 return (x); in __cbrtl()
79 if (x > 0) in __cbrtl()
84 x = -x; in __cbrtl()
87 z = x; in __cbrtl()
89 x = __frexpl (x, &e); in __cbrtl()
93 x = ((((L(1.3584464340920900529734e-1) * x in __cbrtl()
[all …]
|
| D | s_log1pl.c | 126 _Float128 x, y, z, r, s; in __log1pl() local
150 x = xm1; in __log1pl()
152 x = xm1 + 1; in __log1pl()
155 if (x <= 0) in __log1pl()
157 if (x == 0) in __log1pl()
160 return (zero / (x - x)); in __log1pl()
166 x = __frexpl (x, &e); in __log1pl()
172 if (x < sqrth) in __log1pl()
175 z = x - L(0.5); in __log1pl()
180 z = x - L(0.5); in __log1pl()
[all …]
|
| D | e_log2l.c | 140 neval (_Float128 x, const _Float128 *p, int n) in neval() argument
148 y = y * x + *p--; in neval()
158 deval (_Float128 x, const _Float128 *p, int n) in deval() argument
163 y = x + *p--; in deval()
166 y = y * x + *p--; in deval()
175 __ieee754_log2l (_Float128 x) in __ieee754_log2l() argument
183 GET_LDOUBLE_WORDS64 (hx, lx, x); in __ieee754_log2l()
185 return (-1 / fabsl (x)); /* log2l(+-0)=-inf */ in __ieee754_log2l()
187 return (x - x) / (x - x); in __ieee754_log2l()
189 return (x + x); in __ieee754_log2l()
[all …]
|
| D | e_log10l.c | 146 neval (_Float128 x, const _Float128 *p, int n) in neval() argument
154 y = y * x + *p--; in neval()
164 deval (_Float128 x, const _Float128 *p, int n) in deval() argument
169 y = x + *p--; in deval()
172 y = y * x + *p--; in deval()
181 __ieee754_log10l (_Float128 x) in __ieee754_log10l() argument
189 GET_LDOUBLE_WORDS64 (hx, lx, x); in __ieee754_log10l()
191 return (-1 / fabsl (x)); /* log10l(+-0)=-inf */ in __ieee754_log10l()
193 return (x - x) / (x - x); in __ieee754_log10l()
195 return (x + x); in __ieee754_log10l()
[all …]
|
| /glibc-2.36/sysdeps/ieee754/flt-32/ |
| D | e_j0f.c | 196 j0f_asympt (float x) in j0f_asympt() argument
200 if (x == 0x1.4665d2p+24f) in j0f_asympt()
202 if (x == 0x1.a9afdep+7f) in j0f_asympt()
204 double y = 1.0 / (double) x; in j0f_asympt()
210 h = reduce_aux (x, &n, alpha0); in j0f_asympt()
215 float t = cst / sqrtf (x) * (float) beta0; in j0f_asympt()
231 j0f_near_root (float x, float z) in j0f_near_root() argument
236 index_f = roundf ((x - FIRST_ZERO_J0) / M_PIf); in j0f_near_root()
240 return j0f_asympt (x); in j0f_near_root()
246 if (! (x0 <= x && x <= x1)) in j0f_near_root()
[all …]
|
| D | e_j1f.c | 200 j1f_asympt (float x) in j1f_asympt() argument
203 if (x < 0) in j1f_asympt()
205 x = -x; in j1f_asympt()
208 double y = 1.0 / (double) x; in j1f_asympt()
215 h = reduce_aux (x, &n, alpha1); in j1f_asympt()
220 float t = cst / sqrtf (x) * (float) beta1; in j1f_asympt()
236 j1f_near_root (float x, float z) in j1f_near_root() argument
241 if (x < 0) in j1f_near_root()
243 x = -x; in j1f_near_root()
246 index_f = roundf ((x - FIRST_ZERO_J1) / M_PIf); in j1f_near_root()
[all …]
|
| D | s_remquof.c | 29 __remquof (float x, float y, int *quo) in __remquof() argument
35 GET_FLOAT_WORD (hx, x); in __remquof()
44 return (x * y) / (x * y); /* y = 0 */ in __remquof()
47 return (x * y) / (x * y); in __remquof()
50 x = __ieee754_fmodf (x, 8 * y); /* now x < 8y */ in __remquof()
55 return zero * x; in __remquof()
58 x = fabsf (x); in __remquof()
62 if (hy <= 0x7e7fffff && x >= 4 * y) in __remquof()
64 x -= 4 * y; in __remquof()
67 if (hy <= 0x7effffff && x >= 2 * y) in __remquof()
[all …]
|
| /glibc-2.36/sysdeps/unix/sysv/linux/sparc/ |
| D | a.out.h | 65 #define N_BADMAG(x) \ argument
66 (N_MAGIC(x) != OMAGIC && N_MAGIC(x) != NMAGIC \
67 && N_MAGIC(x) != ZMAGIC && N_MAGIC(x) != QMAGIC)
68 #define _N_HDROFF(x) (1024 - sizeof (struct exec)) argument
69 #define N_TXTOFF(x) \ argument
70 (N_MAGIC(x) == ZMAGIC ? 0 : sizeof (struct exec))
71 #define N_DATOFF(x) (N_TXTOFF(x) + (x).a_text) argument
72 #define N_TRELOFF(x) (N_DATOFF(x) + (x).a_data) argument
73 #define N_DRELOFF(x) (N_TRELOFF(x) + N_TRSIZE(x)) argument
74 #define N_SYMOFF(x) \ argument
[all …]
|
| /glibc-2.36/sysdeps/unix/sysv/linux/ |
| D | a.out.h | 62 #define N_BADMAG(x) \ argument
63 (N_MAGIC(x) != OMAGIC && N_MAGIC(x) != NMAGIC \
64 && N_MAGIC(x) != ZMAGIC && N_MAGIC(x) != QMAGIC)
65 #define _N_HDROFF(x) (1024 - sizeof (struct exec)) argument
66 #define N_TXTOFF(x) \ argument
67 (N_MAGIC(x) == ZMAGIC ? _N_HDROFF((x)) + sizeof (struct exec) \
68 : (N_MAGIC(x) == QMAGIC ? 0 : sizeof (struct exec)))
69 #define N_DATOFF(x) (N_TXTOFF(x) + (x).a_text) argument
70 #define N_TRELOFF(x) (N_DATOFF(x) + (x).a_data) argument
71 #define N_DRELOFF(x) (N_TRELOFF(x) + N_TRSIZE(x)) argument
[all …]
|
| /glibc-2.36/sysdeps/ia64/fpu/ |
| D | libm_support.h | 100 #define DOUBLE_HEX(hi,lo) 0x##lo,0x##hi /*LITTLE_ENDIAN*/
257 # define I64CONST(bits) 0x##bits##i64
258 # define U64CONST(bits) 0x##bits##ui64
260 # define I64CONST(bits) 0x##bits##L
261 # define U64CONST(bits) 0x##bits##uL
263 # define I64CONST(bits) 0x##bits##LL
264 # define U64CONST(bits) 0x##bits##uLL
284 #define __LIBM_ADDL1_K80( rhi,rlo,x,y, t1 ) \ argument
285 rhi = x + y; \
286 rlo = rhi - x; \
[all …]
|