Home
last modified time | relevance | path

Searched refs:z (Results 1 – 25 of 135) sorted by relevance

123456

/bionic/libm/upstream-freebsd/lib/msun/src/
Dk_sincosl.h54 long double hz, r, v, w, z; in __kernel_sincosl() local
56 z = x * x; in __kernel_sincosl()
57 v = z * x; in __kernel_sincosl()
62 r = S2 + z * (S3 + z * (S4 + z * (S5 + z * (S6 + z * (S7 + z * S8))))); in __kernel_sincosl()
65 *sn = x + v * (S1 + z * r); in __kernel_sincosl()
67 *sn = x - ((z * (y / 2 - v * r) - y) - v * S1); in __kernel_sincosl()
69 hz = z / 2; in __kernel_sincosl()
71 r = z * (C1 + z * (C2 + z * (C3 + z * (C4 + z * (C5 + z * (C6 + in __kernel_sincosl()
72 z * C7)))))); in __kernel_sincosl()
73 *cs = w + (((1 - w) - hz) + (z * r - x * y)); in __kernel_sincosl()
[all …]
De_j1.c89 double z, s,c,ss,cc,r,u,v,y; in __ieee754_j1() local
101 z = cos(y+y); in __ieee754_j1()
102 if ((s*c)>zero) cc = z/ss; in __ieee754_j1()
103 else ss = z/cc; in __ieee754_j1()
109 if(ix>0x48000000) z = (invsqrtpi*cc)/sqrt(y); in __ieee754_j1()
112 z = invsqrtpi*(u*cc-v*ss)/sqrt(y); in __ieee754_j1()
114 if(hx<0) return -z; in __ieee754_j1()
115 else return z; in __ieee754_j1()
120 z = x*x; in __ieee754_j1()
121 r = z*(r00+z*(r01+z*(r02+z*r03))); in __ieee754_j1()
[all …]
De_j1f.c51 float z, s,c,ss,cc,r,u,v,y; in __ieee754_j1f() local
63 z = cosf(y+y); in __ieee754_j1f()
64 if ((s*c)>zero) cc = z/ss; in __ieee754_j1f()
65 else ss = z/cc; in __ieee754_j1f()
71 if(ix>0x58000000) z = (invsqrtpi*cc)/sqrtf(y); /* |x|>2**49 */ in __ieee754_j1f()
74 z = invsqrtpi*(u*cc-v*ss)/sqrtf(y); in __ieee754_j1f()
76 if(hx<0) return -z; in __ieee754_j1f()
77 else return z; in __ieee754_j1f()
82 z = x*x; in __ieee754_j1f()
83 r = z*(r00+z*(r01+z*(r02+z*r03))); in __ieee754_j1f()
[all …]
De_j0f.c50 float z, s,c,ss,cc,r,u,v; in __ieee754_j0f() local
62 z = -cosf(x+x); in __ieee754_j0f()
63 if ((s*c)<zero) cc = z/ss; in __ieee754_j0f()
64 else ss = z/cc; in __ieee754_j0f()
70 if(ix>0x58000000) z = (invsqrtpi*cc)/sqrtf(x); /* |x|>2**49 */ in __ieee754_j0f()
73 z = invsqrtpi*(u*cc-v*ss)/sqrtf(x); in __ieee754_j0f()
75 return z; in __ieee754_j0f()
83 z = x*x; in __ieee754_j0f()
84 r = z*(R02+z*(R03+z*(R04+z*R05))); in __ieee754_j0f()
85 s = one+z*(S01+z*(S02+z*(S03+z*S04))); in __ieee754_j0f()
[all …]
De_j0.c88 double z, s,c,ss,cc,r,u,v; in __ieee754_j0() local
100 z = -cos(x+x); in __ieee754_j0()
101 if ((s*c)<zero) cc = z/ss; in __ieee754_j0()
102 else ss = z/cc; in __ieee754_j0()
108 if(ix>0x48000000) z = (invsqrtpi*cc)/sqrt(x); in __ieee754_j0()
111 z = invsqrtpi*(u*cc-v*ss)/sqrt(x); in __ieee754_j0()
113 return z; in __ieee754_j0()
121 z = x*x; in __ieee754_j0()
122 r = z*(R02+z*(R03+z*(R04+z*R05))); in __ieee754_j0()
123 s = one+z*(S01+z*(S02+z*(S03+z*S04))); in __ieee754_j0()
[all …]
De_acos.c67 double z,p,q,r,w,s,c,df; in __ieee754_acos() local
82 z = x*x; in __ieee754_acos()
83 p = z*(pS0+z*(pS1+z*(pS2+z*(pS3+z*(pS4+z*pS5))))); in __ieee754_acos()
84 q = one+z*(qS1+z*(qS2+z*(qS3+z*qS4))); in __ieee754_acos()
88 z = (one+x)*0.5; in __ieee754_acos()
89 p = z*(pS0+z*(pS1+z*(pS2+z*(pS3+z*(pS4+z*pS5))))); in __ieee754_acos()
90 q = one+z*(qS1+z*(qS2+z*(qS3+z*qS4))); in __ieee754_acos()
91 s = sqrt(z); in __ieee754_acos()
96 z = (one-x)*0.5; in __ieee754_acos()
97 s = sqrt(z); in __ieee754_acos()
[all …]
De_rem_pio2.c57 double z,w,t,r,fn; in __ieee754_rem_pio2() local
73 z = x - pio2_1; /* one round good to 85 bits */ in __ieee754_rem_pio2()
74 y[0] = z - pio2_1t; in __ieee754_rem_pio2()
75 y[1] = (z-y[0])-pio2_1t; in __ieee754_rem_pio2()
78 z = x + pio2_1; in __ieee754_rem_pio2()
79 y[0] = z + pio2_1t; in __ieee754_rem_pio2()
80 y[1] = (z-y[0])+pio2_1t; in __ieee754_rem_pio2()
85 z = x - 2*pio2_1; in __ieee754_rem_pio2()
86 y[0] = z - 2*pio2_1t; in __ieee754_rem_pio2()
87 y[1] = (z-y[0])-2*pio2_1t; in __ieee754_rem_pio2()
[all …]
Dk_sincos.h36 double hz, r, v, w, z; in __kernel_sincos() local
38 z = x * x; in __kernel_sincos()
39 w = z * z; in __kernel_sincos()
40 r = S2 + z * (S3 + z * S4) + z * w * (S5 + z * S6); in __kernel_sincos()
41 v = z * x; in __kernel_sincos()
44 *sn = x + v * (S1 + z * r); in __kernel_sincos()
46 *sn = x - ((z * (y / 2 - v * r) - y) - v * S1); in __kernel_sincos()
48 r = z * (C1 + z * (C2 + z * C3)) + w * w * (C4 + z * (C5 + z * C6)); in __kernel_sincos()
49 hz = z / 2; in __kernel_sincos()
51 *cs = w + (((1 - w) - hz) + (z * r - x * y)); in __kernel_sincos()
Ds_erf.c200 double R,S,P,Q,s,y,z,r; in erf() local
214 z = x*x; in erf()
215 r = pp0+z*(pp1+z*(pp2+z*(pp3+z*pp4))); in erf()
216 s = one+z*(qq1+z*(qq2+z*(qq3+z*(qq4+z*qq5)))); in erf()
239 z = x; in erf()
240 SET_LOW_WORD(z,0); in erf()
241 r = __ieee754_exp(-z*z-0.5625)*__ieee754_exp((z-x)*(z+x)+R/S); in erf()
253 double R,S,P,Q,s,y,z,r; in erfc() local
264 z = x*x; in erfc()
265 r = pp0+z*(pp1+z*(pp2+z*(pp3+z*pp4))); in erfc()
[all …]
De_acosf.c37 float z,p,q,r,w,s,c,df; in __ieee754_acosf() local
50 z = x*x; in __ieee754_acosf()
51 p = z*(pS0+z*(pS1+z*pS2)); in __ieee754_acosf()
52 q = one+z*qS1; in __ieee754_acosf()
56 z = (one+x)*(float)0.5; in __ieee754_acosf()
57 p = z*(pS0+z*(pS1+z*pS2)); in __ieee754_acosf()
58 q = one+z*qS1; in __ieee754_acosf()
59 s = sqrtf(z); in __ieee754_acosf()
65 z = (one-x)*(float)0.5; in __ieee754_acosf()
66 s = sqrtf(z); in __ieee754_acosf()
[all …]
De_lgamma_r.c167 double y,z; in sin_pi() local
173 z = vz-0x1p52; /* rint(y) for the above range */ in sin_pi()
174 if (z == y) in sin_pi()
179 z = vz-0x1p50; /* y rounded to a multiple of 0.25 */ in sin_pi()
180 if (z > y) { in sin_pi()
181 z -= 0.25; /* adjust to round down */ in sin_pi()
185 y = y - z + n * 0.25; /* y mod 2 */ in sin_pi()
204 double nadj,p,p1,p2,p3,q,r,t,w,y,z; in __ieee754_lgamma_r() local
252 z = y*y; in __ieee754_lgamma_r()
253 p1 = a0+z*(a2+z*(a4+z*(a6+z*(a8+z*a10)))); in __ieee754_lgamma_r()
[all …]
Ds_erff.c84 float R,S,P,Q,s,y,z,r; in erff() local
98 z = x*x; in erff()
99 r = pp0+z*(pp1+z*pp2); in erff()
100 s = one+z*(qq1+z*(qq2+z*qq3)); in erff()
122 SET_FLOAT_WORD(z,hx&0xffffe000); in erff()
123 r = expf(-z*z-0.5625F)*expf((z-x)*(z+x)+R/S); in erff()
131 float R,S,P,Q,s,y,z,r; in erfcf() local
142 z = x*x; in erfcf()
143 r = pp0+z*(pp1+z*pp2); in erfcf()
144 s = one+z*(qq1+z*(qq2+z*qq3)); in erfcf()
[all …]
De_lgammaf_r.c88 float y,z; in sin_pif() local
94 z = vz-0x1p23F; /* rintf(y) for the above range */ in sin_pif()
95 if (z == y) in sin_pif()
100 z = vz-0x1p21F; /* y rounded to a multiple of 0.25 */ in sin_pif()
101 if (z > y) { in sin_pif()
102 z -= 0.25F; /* adjust to round down */ in sin_pif()
106 y = y - z + n * 0.25F; /* y mod 2 */ in sin_pif()
125 float nadj,p,p1,p2,q,r,t,w,y,z; in __ieee754_lgammaf_r() local
173 z = y*y; in __ieee754_lgammaf_r()
174 p1 = a0+z*(a2+z*a4); in __ieee754_lgammaf_r()
[all …]
De_powf.c61 float z,ax,z_h,z_l,p_h,p_l; in __ieee754_powf() local
117 z = ax; /*x is +-0,+-inf,+-1*/ in __ieee754_powf()
118 if(hy<0) z = one/z; /* z = (1/|x|) */ in __ieee754_powf()
121 z = (z-z)/(z-z); /* (-1)**non-int is NaN */ in __ieee754_powf()
123 z = -z; /* (x<0)**odd = -(|x|**odd) */ in __ieee754_powf()
125 return z; in __ieee754_powf()
210 z = p_l+p_h; in __ieee754_powf()
211 GET_FLOAT_WORD(j,z); in __ieee754_powf()
215 if(p_l+ovt>z-p_h) return sn*huge*huge; /* overflow */ in __ieee754_powf()
220 if(p_l<=z-p_h) return sn*tiny*tiny; /* underflow */ in __ieee754_powf()
[all …]
De_pow.c103 double z,ax,z_h,z_l,p_h,p_l; in __ieee754_pow() local
168 z = ax; /*x is +-0,+-inf,+-1*/ in __ieee754_pow()
169 if(hy<0) z = one/z; /* z = (1/|x|) */ in __ieee754_pow()
172 z = (z-z)/(z-z); /* (-1)**non-int is NaN */ in __ieee754_pow()
174 z = -z; /* (x<0)**odd = -(|x|**odd) */ in __ieee754_pow()
176 return z; in __ieee754_pow()
265 z = p_l+p_h; in __ieee754_pow()
266 EXTRACT_WORDS(j,i,z); in __ieee754_pow()
271 if(p_l+ovt>z-p_h) return s*huge*huge; /* overflow */ in __ieee754_pow()
277 if(p_l<=z-p_h) return s*tiny*tiny; /* underflow */ in __ieee754_pow()
[all …]
Dk_cos.c71 double hz,z,r,w; in __kernel_cos() local
73 z = x*x; in __kernel_cos()
74 w = z*z; in __kernel_cos()
75 r = z*(C1+z*(C2+z*C3)) + w*w*(C4+z*(C5+z*C6)); in __kernel_cos()
76 hz = 0.5*z; in __kernel_cos()
78 return w + (((one-w)-hz) + (z*r-x*y)); in __kernel_cos()
/bionic/libm/upstream-freebsd/lib/msun/ld128/
Dk_cosl.c53 long double hz,z,r,w; in __kernel_cosl() local
55 z = x*x; in __kernel_cosl()
56 r = z*(C1+z*(C2+z*(C3+z*(C4+z*(C5+z*(C6+z*(C7+ in __kernel_cosl()
57 z*(C8+z*(C9+z*(C10+z*C11)))))))))); in __kernel_cosl()
58 hz = 0.5*z; in __kernel_cosl()
60 return w + (((one-w)-hz) + (z*r-x*y)); in __kernel_cosl()
Ds_erfl.c196 long double ax,R,S,P,Q,s,y,z,r; in erfl() local
215 z = x*x; in erfl()
216 r = pp0+z*(pp1+z*(pp2+z*(pp3+z*(pp4+z*(pp5+z*(pp6+z*(pp7+ in erfl()
217 z*(pp8+z*pp9)))))))); in erfl()
218 s = one+z*(qq1+z*(qq2+z*(qq3+z*(qq4+z*(qq5+z*(qq6+z*(qq7+ in erfl()
219 z*(qq8+z*qq9)))))))); in erfl()
250 z = (float)ax; in erfl()
251 r = expl(-z*z-0.5625)*expl((z-ax)*(z+ax)+R/S); in erfl()
258 long double ax,R,S,P,Q,s,y,z,r; in erfcl() local
273 z = x*x; in erfcl()
[all …]
Dk_sinl.c51 long double z,r,v; in __kernel_sinl() local
53 z = x*x; in __kernel_sinl()
54 v = z*x; in __kernel_sinl()
55 r = S2+z*(S3+z*(S4+z*(S5+z*(S6+z*(S7+z*(S8+ in __kernel_sinl()
56 z*(S9+z*(S10+z*(S11+z*S12))))))))); in __kernel_sinl()
57 if(iy==0) return x+v*(S1+z*r); in __kernel_sinl()
58 else return x-((z*(half*y-v*r)-y)-v*S1); in __kernel_sinl()
De_lgammal_r.c190 long double y,z; in sin_pil() local
197 z = vz-0x1.p112; in sin_pil()
198 if (z == y) in sin_pil()
203 z = vz-0x1.p110; in sin_pil()
204 if (z > y) { in sin_pil()
205 z -= 0.25; in sin_pil()
209 y = y - z + n * 0.25; in sin_pil()
227 long double nadj,p,p1,p2,p3,q,r,t,w,y,z; in lgammal_r() local
276 z = y*y; in lgammal_r()
277 p1 = a0+z*(a2+z*(a4+z*(a6+z*(a8+z*(a10+z*(a12+z*(a14+z*(a16+ in lgammal_r()
[all …]
De_powl.c150 long double z, ax, z_h, z_l, p_h, p_l; in powl() local
201 z = 0.5 * y; in powl()
202 if (floorl (z) == z) in powl()
245 z = ax; /*x is +-0,+-inf,+-1 */ in powl()
247 z = one / z; /* z = (1/|x|) */ in powl()
252 z = (z - z) / (z - z); /* (-1)**non-int is NaN */ in powl()
255 z = -z; /* (x<0)**odd = -(|x|**odd) */ in powl()
257 return z; in powl()
378 z = p_l + p_h; in powl()
379 o.value = z; in powl()
[all …]
/bionic/libm/upstream-freebsd/lib/msun/bsdsrc/
Db_tgamma.c161 double z, p; local
164 z = one/(x*x);
165 p = Pa0+z*(Pa1+z*(Pa2+z*(Pa3+z*(Pa4+z*(Pa5+z*(Pa6+z*Pa7))))));
250 ratfun_gam(z, c) in ratfun_gam() argument
251 double z, c; in ratfun_gam()
256 q = Q0 +z*(Q1+z*(Q2+z*(Q3+z*(Q4+z*(Q5+z*(Q6+z*(Q7+z*Q8)))))));
257 p = P0 + z*(P1 + z*(P2 + z*(P3 + z*P4)));
261 t.a = z, TRUNC(t.a); /* t ~= z + c */
262 t.b = (z - t.a) + c;
263 t.b *= (t.a + z);
[all …]
Db_exp.c94 double z,hi,lo,c;
113 z=x*x;
114 c= x - z*(p1+z*(p2+z*(p3+z*(p4+z*p5))));
140 double z,hi,lo; local
149 z = invln2*x;
150 k = z + copysign(.5, x);
157 z=x*x;
158 c= x - z*(p1+z*(p2+z*(p3+z*(p4+z*p5))));
/bionic/libc/upstream-openbsd/lib/libc/gdtoa/
Dmisc.c181 ULong xi, z; local
199 z = (xi >> 16) * m + (y >> 16);
200 carry = z >> 16;
201 *x++ = (z << 16) + (y & 0xffff);
290 ULLong carry, z; local
292 ULong carry, z; local
326 z = *x++ * (ULLong)y + *xc + carry;
327 carry = z >> 32;
328 *xc++ = z & 0xffffffffUL;
342 z = (*x & 0xffff) * y + (*xc & 0xffff) + carry;
[all …]
/bionic/tests/
Dcomplex_test.cpp320 double complex z; in TEST() local
323 z = ctanh(nan("") + 0i); in TEST()
324 ASSERT_TRUE(isnan(creal(z))); in TEST()
325 ASSERT_EQ(0.0, cimag(z)); in TEST()
328 z = ctanh(nan("") + 2.0i); in TEST()
329 ASSERT_TRUE(isnan(creal(z))); in TEST()
330 ASSERT_TRUE(isnan(cimag(z))); in TEST()
333 z = ctanh(nan("") + nan("") * I); in TEST()
334 ASSERT_TRUE(isnan(creal(z))); in TEST()
335 ASSERT_TRUE(isnan(cimag(z))); in TEST()
[all …]

123456