| /picolibc-3.7.0-3.6.0/newlib/libm/ld/ld128/ |
| D | e_lgammal_r.c | 759 long double p, q, w, z, nx; in lgammal_r() local
786 z = q - p; in lgammal_r()
787 if (z > 0.5L) in lgammal_r()
790 z = p - q; in lgammal_r()
792 z = q * sinl (PIL * z); in lgammal_r()
793 if (z == 0.0L) in lgammal_r()
796 z = logl (PIL / z) - w; in lgammal_r()
797 return (z); in lgammal_r()
815 z = x - 0.25L; in lgammal_r()
816 p = z * neval (z, RN1r25, NRN1r25) / deval (z, RD1r25, NRD1r25); in lgammal_r()
[all …]
|
| D | s_log1pl.c | 124 long double x, y, z, r, s; in log1pl() local
163 z = x - 0.5L; in log1pl()
164 y = 0.5L * z + 0.5L; in log1pl()
168 z = x - 0.5L; in log1pl()
169 z -= 0.5L; in log1pl()
172 x = z / y; in log1pl()
173 z = x * x; in log1pl()
174 r = ((((R5 * z in log1pl()
175 + R4) * z in log1pl()
176 + R3) * z in log1pl()
[all …]
|
| D | s_erfl.c | 749 long double a, y, z; in erfl() local
778 z = x * x; in erfl()
787 y = a + a * neval (z, TN1, NTN1) / deval (z, TD1, NTD1); in erfl()
803 long double y, z, p, r; in erfcl() local
832 z = x - 0.25L; in erfcl()
833 y = C13b + z * neval (z, RNr13, NRNr13) / deval (z, RDr13, NRDr13); in erfcl()
837 z = x - 0.375L; in erfcl()
838 y = C14b + z * neval (z, RNr14, NRNr14) / deval (z, RDr14, NRDr14); in erfcl()
842 z = x - 0.5L; in erfcl()
843 y = C15b + z * neval (z, RNr15, NRNr15) / deval (z, RDr15, NRDr15); in erfcl()
[all …]
|
| D | k_cosl.c | 51 long double hz,z,r,w; in __kernel_cosl() local
53 z = x*x; in __kernel_cosl()
54 r = z*(C1+z*(C2+z*(C3+z*(C4+z*(C5+z*(C6+z*(C7+ in __kernel_cosl()
55 z*((long double)C8+z*((long double)C9+z*((long double)C10+z*(long double)C11)))))))))); in __kernel_cosl()
56 hz = 0.5L*z; in __kernel_cosl()
58 return w + ((((long double)one-w)-hz) + (z*r-x*y)); in __kernel_cosl()
|
| D | k_sinl.c | 49 long double z,r,v; in __kernel_sinl() local
51 z = x*x; in __kernel_sinl()
52 v = z*x; in __kernel_sinl()
53 r = S2+z*(S3+z*(S4+z*(S5+z*(S6+z*(S7+z*(S8+ in __kernel_sinl()
54 z*((long double)S9+z*((long double)S10+z*((long double)S11+z*(long double)S12))))))))); in __kernel_sinl()
55 if(iy==0) return x+v*(S1+z*r); in __kernel_sinl()
56 else return x-((z*((long double)half*y-v*r)-y)-v*S1); in __kernel_sinl()
|
| D | e_log10l.c | 180 long double z; in log10l() local
210 z = x - 0.5L; in log10l()
211 y = 0.5L * z + 0.5L; in log10l()
215 z = x - 0.5L; in log10l()
216 z -= 0.5L; in log10l()
219 x = z / y; in log10l()
220 z = x * x; in log10l()
221 y = x * (z * neval (z, R, 5) / deval (z, S, 5)); in log10l()
237 z = x * x; in log10l()
238 y = x * (z * neval (x, P, 12) / deval (x, Q, 11)); in log10l()
[all …]
|
| D | e_log2l.c | 174 long double z; in log2l() local
204 z = x - 0.5L; in log2l()
205 y = 0.5L * z + 0.5L; in log2l()
209 z = x - 0.5L; in log2l()
210 z -= 0.5L; in log2l()
213 x = z / y; in log2l()
214 z = x * x; in log2l()
215 y = x * (z * neval (z, R, 5) / deval (z, S, 5)); in log2l()
231 z = x * x; in log2l()
232 y = x * (z * neval (x, P, 12) / deval (x, Q, 11)); in log2l()
[all …]
|
| D | e_logl.c | 190 long double z, y, w; in logl() local
259 z = x - 1.0L; in logl()
270 z = (u.value - t.value) / t.value; in logl()
273 w = z * z; in logl()
274 y = ((((((((((((l15 * z in logl()
275 + l14) * z in logl()
276 + l13) * z in logl()
277 + l12) * z in logl()
278 + l11) * z in logl()
279 + l10) * z in logl()
[all …]
|
| /picolibc-3.7.0-3.6.0/newlib/libm/math/ |
| D | sf_acos.c | 36 float z, p, q, r, w, s, c, df; in acosf() local
51 z = x * x; in acosf()
52 p = z * (pS0 + z * (pS1 + z * (pS2 + z * (pS3 + z * (pS4 + z * pS5))))); in acosf()
53 q = one + z * (qS1 + z * (qS2 + z * (qS3 + z * qS4))); in acosf()
57 z = (one + x) * (float)0.5; in acosf()
58 p = z * (pS0 + z * (pS1 + z * (pS2 + z * (pS3 + z * (pS4 + z * pS5))))); in acosf()
59 q = one + z * (qS1 + z * (qS2 + z * (qS3 + z * qS4))); in acosf()
60 s = sqrtf(z); in acosf()
66 z = (one - x) * (float)0.5; in acosf()
67 s = sqrtf(z); in acosf()
[all …]
|
| D | s_acos.c | 61 __float64 z, p, q, r, w, s, c, df; in acos64() local
79 z = x * x; in acos64()
80 p = z * (pS0 + z * (pS1 + z * (pS2 + z * (pS3 + z * (pS4 + z * pS5))))); in acos64()
81 q = one + z * (qS1 + z * (qS2 + z * (qS3 + z * qS4))); in acos64()
85 z = (one + x) * _F_64(0.5); in acos64()
86 p = z * (pS0 + z * (pS1 + z * (pS2 + z * (pS3 + z * (pS4 + z * pS5))))); in acos64()
87 q = one + z * (qS1 + z * (qS2 + z * (qS3 + z * qS4))); in acos64()
88 s = sqrt64(z); in acos64()
93 z = (one - x) * _F_64(0.5); in acos64()
94 s = sqrt64(z); in acos64()
[all …]
|
| D | sf_j1.c | 39 float z, s, c, ss, cc, r, u, v, y; in j1f() local
58 z = cosf(y + y); in j1f()
60 cc = z / ss; in j1f()
62 ss = z / cc; in j1f()
69 z = (invsqrtpi * cc) / sqrtf(y); in j1f()
73 z = invsqrtpi * (u * cc - v * ss) / sqrtf(y); in j1f()
76 return -z; in j1f()
78 return z; in j1f()
85 z = x * x; in j1f()
86 r = z * (r00 + z * (r01 + z * (r02 + z * r03))); in j1f()
[all …]
|
| D | sf_j0.c | 38 float z, s, c, ss, cc, r, u, v; in j0f() local
56 z = -cosf(x + x); in j0f()
58 cc = z / ss; in j0f()
60 ss = z / cc; in j0f()
67 z = (invsqrtpi * cc) / sqrtf(x); in j0f()
71 z = invsqrtpi * (u * cc - v * ss) / sqrtf(x); in j0f()
73 return z; in j0f()
83 z = x * x; in j0f()
84 r = z * (R02 + z * (R03 + z * (R04 + z * R05))); in j0f()
85 s = one + z * (S01 + z * (S02 + z * (S03 + z * S04))); in j0f()
[all …]
|
| D | s_j0.c | 85 __float64 z, s, c, ss, cc, r, u, v; in j064() local
103 z = -cos64(x + x); in j064()
105 cc = z / ss; in j064()
107 ss = z / cc; in j064()
114 z = (invsqrtpi * cc) / sqrt64(x); in j064()
118 z = invsqrtpi * (u * cc - v * ss) / sqrt64(x); in j064()
120 return z; in j064()
130 z = x * x; in j064()
131 r = z * (R02 + z * (R03 + z * (R04 + z * R05))); in j064()
132 s = one + z * (S01 + z * (S02 + z * (S03 + z * S04))); in j064()
[all …]
|
| D | s_j1.c | 87 __float64 z, s, c, ss, cc, r, u, v, y; in j164() local
105 z = cos64(y + y); in j164()
107 cc = z / ss; in j164()
109 ss = z / cc; in j164()
116 z = (invsqrtpi * cc) / sqrt64(y); in j164()
120 z = invsqrtpi * (u * cc - v * ss) / sqrt64(y); in j164()
123 return -z; in j164()
125 return z; in j164()
133 z = x * x; in j164()
134 r = z * (r00 + z * (r01 + z * (r02 + z * r03))); in j164()
[all …]
|
| D | kf_cos.c | 29 float a, hz, z, r, qx; in __kernel_cosf() local
35 z = x * x; in __kernel_cosf()
36 r = z * (C1 + z * (C2 + z * (C3 + z * (C4 + z * (C5 + z * C6))))); in __kernel_cosf()
38 return one - ((float)0.5 * z - (z * r - x * y)); in __kernel_cosf()
45 hz = (float)0.5 * z - qx; in __kernel_cosf()
47 return a - (hz - (z * r - x * y)); in __kernel_cosf()
|
| D | k_cos.c | 65 __float64 a, hz, z, r, qx; in __kernel_cos() local
71 z = x * x; in __kernel_cos()
72 r = z * (C1 + z * (C2 + z * (C3 + z * (C4 + z * (C5 + z * C6))))); in __kernel_cos()
74 return one - (_F_64(0.5) * z - (z * r - x * y)); in __kernel_cos()
81 hz = _F_64(0.5) * z - qx; in __kernel_cos()
83 return a - (hz - (z * r - x * y)); in __kernel_cos()
|
| D | sf_erf.c | 100 float R, S, P, Q, s, y, z, r; in erff() local
115 z = x * x; in erff()
116 r = pp0 + z * (pp1 + z * (pp2 + z * (pp3 + z * pp4))); in erff()
117 s = one + z * (qq1 + z * (qq2 + z * (qq3 + z * (qq4 + z * qq5)))); in erff()
160 SET_FLOAT_WORD(z, ix & 0xfffff000); in erff()
161 r = expf(-z * z - (float)0.5625) * expf((z - x) * (z + x) + R / S); in erff()
172 float R, S, P, Q, s, y, z, r; in erfcf() local
183 z = x * x; in erfcf()
184 r = pp0 + z * (pp1 + z * (pp2 + z * (pp3 + z * pp4))); in erfcf()
185 s = one + z * (qq1 + z * (qq2 + z * (qq3 + z * (qq4 + z * qq5)))); in erfcf()
[all …]
|
| /picolibc-3.7.0-3.6.0/newlib/libm/ld/ld80/ |
| D | e_logl.c | 117 long double y, z; in logl() local
148 z = x - 0.5L; in logl()
149 y = 0.5L * z + 0.5L; in logl()
153 z = x - 0.5L; in logl()
154 z -= 0.5L; in logl()
157 x = z / y; in logl()
158 z = x*x; in logl()
159 z = x * ( z * __polevll( z, R, 3 ) / __p1evll( z, S, 3 ) ); in logl()
160 z = z + e * C2; in logl()
161 z = z + x; in logl()
[all …]
|
| D | s_log1pl.c | 114 long double x, y, z; in log1pl() local
146 z = x - 0.5L; in log1pl()
147 y = 0.5L * z + 0.5L; in log1pl()
151 z = x - 0.5L; in log1pl()
152 z -= 0.5L; in log1pl()
155 x = z / y; in log1pl()
156 z = x*x; in log1pl()
157 z = x * ( z * __polevll( z, R, 3 ) / __p1evll( z, S, 3 ) ); in log1pl()
158 z = z + e * C2; in log1pl()
159 z = z + x; in log1pl()
[all …]
|
| D | e_log2l.c | 117 volatile long double z; in log2l() local
151 z = x - 0.5L; in log2l()
152 y = 0.5L * z + 0.5L; in log2l()
156 z = x - 0.5L; in log2l()
157 z -= 0.5L; in log2l()
160 x = z / y; in log2l()
161 z = x*x; in log2l()
162 y = x * ( z * __polevll( z, R, 3 ) / __p1evll( z, S, 3 ) ); in log2l()
178 z = x*x; in log2l()
179 y = x * ( z * __polevll( x, P, 6 ) / __p1evll( x, Q, 7 ) ); in log2l()
[all …]
|
| D | e_log10l.c | 123 volatile long double z; in log10l() local
154 z = x - 0.5L; in log10l()
155 y = 0.5L * z + 0.5L; in log10l()
159 z = x - 0.5L; in log10l()
160 z -= 0.5L; in log10l()
163 x = z / y; in log10l()
164 z = x*x; in log10l()
165 y = x * ( z * __polevll( z, R, 3 ) / __p1evll( z, S, 3 ) ); in log10l()
181 z = x*x; in log10l()
182 y = x * ( z * __polevll( x, P, 6 ) / __p1evll( x, Q, 7 ) ); in log10l()
[all …]
|
| D | e_tgammal.c | 217 long double p, q, z; in tgammal() local
246 z = q - p; in tgammal()
247 if( z > 0.5L ) in tgammal()
250 z = q - p; in tgammal()
252 z = q * sinl( PIL * z ); in tgammal()
253 z = fabsl(z) * stirf(q); in tgammal()
254 if( z <= PIL/LDBL_MAX ) in tgammal()
258 z = PIL/z; in tgammal()
262 z = stirf(x); in tgammal()
264 return( sign * z ); in tgammal()
[all …]
|
| D | s_erfl.c | 249 long double R, S, P, Q, s, y, z, r; in erfl() local
271 z = x * x; in erfl()
272 r = pp[0] + z * (pp[1] in erfl()
273 + z * (pp[2] + z * (pp[3] + z * (pp[4] + z * pp[5])))); in erfl()
274 s = qq[0] + z * (qq[1] in erfl()
275 + z * (qq[2] + z * (qq[3] + z * (qq[4] + z * (qq[5] + z))))); in erfl()
314 z = x; in erfl()
315 GET_LDOUBLE_WORDS (i, i0, i1, z); in erfl()
317 SET_LDOUBLE_WORDS (z, i, i0, i1); in erfl()
319 expl (-z * z - 0.5625L) * expl ((z - x) * (z + x) + R / S); in erfl()
[all …]
|
| D | k_cosl.c | 67 long double hz,z,r,w; in __kernel_cosl() local
69 z = x*x; in __kernel_cosl()
70 …r = z*(C1+z*((long double) C2+z*((long double) C3+z*((long double) C4+z*((long double) C5+z*((lon… in __kernel_cosl()
71 hz = 0.5l*z; in __kernel_cosl()
73 return w + (((1.0l-w)-hz) + (z*r-x*y)); in __kernel_cosl()
|
| D | k_sinl.c | 50 long double z,r,v; in __kernel_sinl() local
52 z = x*x; in __kernel_sinl()
53 v = z*x; in __kernel_sinl()
54 …r = (long double) S2+z*((long double) S3+z*((long double) S4+z*((long double) S5+z*((long double)… in __kernel_sinl()
55 if(iy==0) return x+v*(S1+z*r); in __kernel_sinl()
56 else return x-((z*(0.5l*y-v*r)-y)-v*S1); in __kernel_sinl()
|
