/bionic/libm/upstream-freebsd/lib/msun/src/ |
D | k_sincosl.h | 54 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 …]
|
D | e_acos.c | 67 double z,p,q,r,w,s,c,df; in acos() local 82 z = x*x; in acos() 83 p = z*(pS0+z*(pS1+z*(pS2+z*(pS3+z*(pS4+z*pS5))))); in acos() 84 q = one+z*(qS1+z*(qS2+z*(qS3+z*qS4))); in acos() 88 z = (one+x)*0.5; in acos() 89 p = z*(pS0+z*(pS1+z*(pS2+z*(pS3+z*(pS4+z*pS5))))); in acos() 90 q = one+z*(qS1+z*(qS2+z*(qS3+z*qS4))); in acos() 91 s = sqrt(z); in acos() 96 z = (one-x)*0.5; in acos() 97 s = sqrt(z); in acos() [all …]
|
D | e_j1.c | 89 double z, s,c,ss,cc,r,u,v,y; in j1() local 101 z = cos(y+y); in j1() 102 if ((s*c)>zero) cc = z/ss; in j1() 103 else ss = z/cc; in j1() 109 if(ix>0x48000000) z = (invsqrtpi*cc)/sqrt(y); in j1() 112 z = invsqrtpi*(u*cc-v*ss)/sqrt(y); in j1() 114 if(hx<0) return -z; in j1() 115 else return z; in j1() 120 z = x*x; in j1() 121 r = z*(r00+z*(r01+z*(r02+z*r03))); in j1() [all …]
|
D | e_j1f.c | 51 float z, s,c,ss,cc,r,u,v,y; in j1f() local 63 z = cosf(y+y); in j1f() 64 if ((s*c)>zero) cc = z/ss; in j1f() 65 else ss = z/cc; in j1f() 71 if(ix>0x58000000) z = (invsqrtpi*cc)/sqrtf(y); /* |x|>2**49 */ in j1f() 74 z = invsqrtpi*(u*cc-v*ss)/sqrtf(y); in j1f() 76 if(hx<0) return -z; in j1f() 77 else return z; in j1f() 82 z = x*x; in j1f() 83 r = z*(r00+z*(r01+z*(r02+z*r03))); in j1f() [all …]
|
D | e_j0f.c | 50 float z, s,c,ss,cc,r,u,v; in j0f() local 62 z = -cosf(x+x); in j0f() 63 if ((s*c)<zero) cc = z/ss; in j0f() 64 else ss = z/cc; in j0f() 70 if(ix>0x58000000) z = (invsqrtpi*cc)/sqrtf(x); /* |x|>2**49 */ in j0f() 73 z = invsqrtpi*(u*cc-v*ss)/sqrtf(x); in j0f() 75 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 | e_j0.c | 88 double z, s,c,ss,cc,r,u,v; in j0() local 100 z = -cos(x+x); in j0() 101 if ((s*c)<zero) cc = z/ss; in j0() 102 else ss = z/cc; in j0() 108 if(ix>0x48000000) z = (invsqrtpi*cc)/sqrt(x); in j0() 111 z = invsqrtpi*(u*cc-v*ss)/sqrt(x); in j0() 113 return z; in j0() 121 z = x*x; in j0() 122 r = z*(R02+z*(R03+z*(R04+z*R05))); in j0() 123 s = one+z*(S01+z*(S02+z*(S03+z*S04))); in j0() [all …]
|
D | k_sincos.h | 36 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()
|
D | e_rem_pio2.c | 57 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 …]
|
D | e_acosf.c | 37 float z,p,q,r,w,s,c,df; in acosf() local 50 z = x*x; in acosf() 51 p = z*(pS0+z*(pS1+z*pS2)); in acosf() 52 q = one+z*qS1; in acosf() 56 z = (one+x)*(float)0.5; in acosf() 57 p = z*(pS0+z*(pS1+z*pS2)); in acosf() 58 q = one+z*qS1; in acosf() 59 s = sqrtf(z); in acosf() 65 z = (one-x)*(float)0.5; in acosf() 66 s = sqrtf(z); in acosf() [all …]
|
D | s_erf.c | 200 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 = exp(-z*z-0.5625)*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 …]
|
D | e_lgamma_r.c | 167 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 lgamma_r() local 252 z = y*y; in lgamma_r() 253 p1 = a0+z*(a2+z*(a4+z*(a6+z*(a8+z*a10)))); in lgamma_r() [all …]
|
D | s_erff.c | 84 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 …]
|
D | e_lgammaf_r.c | 88 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 lgammaf_r() local 173 z = y*y; in lgammaf_r() 174 p1 = a0+z*(a2+z*a4); in lgammaf_r() [all …]
|
D | k_cos.c | 71 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()
|
D | k_sin.c | 62 double z,r,v,w; in __kernel_sin() local 64 z = x*x; in __kernel_sin() 65 w = z*z; in __kernel_sin() 66 r = S2+z*(S3+z*S4) + z*w*(S5+z*S6); in __kernel_sin() 67 v = z*x; in __kernel_sin() 68 if(iy==0) return x+v*(S1+z*r); in __kernel_sin() 69 else return x-((z*(half*y-v*r)-y)-v*S1); in __kernel_sin()
|
D | e_powf.c | 61 float z,ax,z_h,z_l,p_h,p_l; in powf() local 117 z = ax; /*x is +-0,+-inf,+-1*/ in powf() 118 if(hy<0) z = one/z; /* z = (1/|x|) */ in powf() 121 z = (z-z)/(z-z); /* (-1)**non-int is NaN */ in powf() 123 z = -z; /* (x<0)**odd = -(|x|**odd) */ in powf() 125 return z; in powf() 210 z = p_l+p_h; in powf() 211 GET_FLOAT_WORD(j,z); in powf() 215 if(p_l+ovt>z-p_h) return sn*huge*huge; /* overflow */ in powf() 220 if(p_l<=z-p_h) return sn*tiny*tiny; /* underflow */ in powf() [all …]
|
D | k_sincosf.h | 34 double r, s, w, z; in __kernel_sincosdf() local 36 z = x * x; in __kernel_sincosdf() 37 w = z * z; in __kernel_sincosdf() 38 r = S3 + z * S4; in __kernel_sincosdf() 39 s = z * x; in __kernel_sincosdf() 40 *sn = (x + s * (S1 + z * S2)) + s * w * r; in __kernel_sincosdf() 41 r = C2 + z * C3; in __kernel_sincosdf() 42 *cs = ((1 + z * C0) + w * C1) + (w * z) * r; in __kernel_sincosdf()
|
D | k_tan.c | 80 double z, r, v, w, s; in __kernel_tan() local 90 z = pio4 - x; in __kernel_tan() 92 x = z + w; in __kernel_tan() 95 z = x * x; in __kernel_tan() 96 w = z * z; in __kernel_tan() 104 v = z * (T[2] + w * (T[4] + w * (T[6] + w * (T[8] + w * (T[10] + in __kernel_tan() 106 s = z * x; in __kernel_tan() 107 r = y + z * (s * (r + v) + y); in __kernel_tan() 124 z = w; in __kernel_tan() 125 SET_LOW_WORD(z,0); in __kernel_tan() [all …]
|
/bionic/libm/upstream-freebsd/lib/msun/ld128/ |
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*(C8+z*(C9+z*(C10+z*(C11+z*C12))))))))))); in __kernel_cosl() 56 hz = 0.5*z; in __kernel_cosl() 58 return w + (((one-w)-hz) + (z*r-x*y)); in __kernel_cosl()
|
D | k_sinl.c | 51 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()
|
D | s_erfl.c | 196 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 …]
|
D | e_lgammal_r.c | 190 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 …]
|
/bionic/libm/upstream-freebsd/lib/msun/bsdsrc/ |
D | b_tgamma.c | 124 double p, z, thi, tlo, xhi, xlo; in large_gam() local 127 z = 1 / (x * x); in large_gam() 128 p = Pa0 + z * (Pa1 + z * (Pa2 + z * (Pa3 + z * (Pa4 + z * (Pa5 + in large_gam() 129 z * (Pa6 + z * Pa7)))))); in large_gam() 184 ratfun_gam(double z, double c) in ratfun_gam() argument 189 q = Q0 + z * (Q1 + z * (Q2 + z * (Q3 + z * (Q4 + z * (Q5 + in ratfun_gam() 190 z * (Q6 + z * (Q7 + z * Q8))))))); in ratfun_gam() 191 p = P0 + z * (P1 + z * (P2 + z * (P3 + z * P4))); in ratfun_gam() 195 thi = (float)z; in ratfun_gam() 196 tlo = (z - thi) + c; in ratfun_gam() [all …]
|
/bionic/tests/ |
D | complex_test.cpp | 315 double complex z; in TEST() local 318 z = ctanh(nan("") + 0i); in TEST() 319 ASSERT_TRUE(isnan(creal(z))); in TEST() 320 ASSERT_EQ(0.0, cimag(z)); in TEST() 323 z = ctanh(nan("") + 2.0i); in TEST() 324 ASSERT_TRUE(isnan(creal(z))); in TEST() 325 ASSERT_TRUE(isnan(cimag(z))); in TEST() 328 z = ctanh(nan("") + nan("") * I); in TEST() 329 ASSERT_TRUE(isnan(creal(z))); in TEST() 330 ASSERT_TRUE(isnan(cimag(z))); in TEST() [all …]
|
/bionic/libc/upstream-openbsd/lib/libc/gdtoa/ |
D | misc.c | 181 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 …]
|