/bionic/libm/upstream-freebsd/lib/msun/src/ |
D | e_acos.c | 67 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 …]
|
D | e_j0f.c | 44 float z, s,c,ss,cc,r,u,v; in __ieee754_j0f() local 57 z = -cosf(x+x); in __ieee754_j0f() 58 if ((s*c)<zero) cc = z/ss; in __ieee754_j0f() 59 else ss = z/cc; in __ieee754_j0f() 65 if(ix>0x80000000) z = (invsqrtpi*cc)/sqrtf(x); in __ieee754_j0f() 68 z = invsqrtpi*(u*cc-v*ss)/sqrtf(x); in __ieee754_j0f() 70 return z; in __ieee754_j0f() 78 z = x*x; in __ieee754_j0f() 79 r = z*(R02+z*(R03+z*(R04+z*R05))); in __ieee754_j0f() 80 s = one+z*(S01+z*(S02+z*(S03+z*S04))); in __ieee754_j0f() [all …]
|
D | e_j1f.c | 45 float z, s,c,ss,cc,r,u,v,y; in __ieee754_j1f() local 58 z = cosf(y+y); in __ieee754_j1f() 59 if ((s*c)>zero) cc = z/ss; in __ieee754_j1f() 60 else ss = z/cc; in __ieee754_j1f() 66 if(ix>0x80000000) z = (invsqrtpi*cc)/sqrtf(y); in __ieee754_j1f() 69 z = invsqrtpi*(u*cc-v*ss)/sqrtf(y); in __ieee754_j1f() 71 if(hx<0) return -z; in __ieee754_j1f() 72 else return z; in __ieee754_j1f() 77 z = x*x; in __ieee754_j1f() 78 r = z*(r00+z*(r01+z*(r02+z*r03))); in __ieee754_j1f() [all …]
|
D | e_j1.c | 88 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 …]
|
D | e_j0.c | 87 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 …]
|
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_acosf.c | 37 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 …]
|
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 | k_cosf.c | 39 double r, w, z; in __kernel_cosdf() local 42 z = x*x; in __kernel_cosdf() 43 w = z*z; in __kernel_cosdf() 44 r = C2+z*C3; in __kernel_cosdf() 45 return ((one+z*C0) + w*C1) + (w*z)*r; in __kernel_cosdf()
|
D | k_tanf.c | 41 double z,r,w,s,t,u; in __kernel_tandf() local 43 z = x*x; in __kernel_tandf() 58 r = T[4]+z*T[5]; in __kernel_tandf() 59 t = T[2]+z*T[3]; in __kernel_tandf() 60 w = z*z; in __kernel_tandf() 61 s = z*x; in __kernel_tandf() 62 u = T[0]+z*T[1]; in __kernel_tandf()
|
D | k_sinf.c | 38 double r, s, w, z; in __kernel_sindf() local 41 z = x*x; in __kernel_sindf() 42 w = z*z; in __kernel_sindf() 43 r = S3+z*S4; in __kernel_sindf() 44 s = z*x; in __kernel_sindf() 45 return (x + s*(S1+z*S2)) + s*w*r; in __kernel_sindf()
|
D | s_cprojf.c | 36 cprojf(float complex z) in cprojf() argument 39 if (!isinf(crealf(z)) && !isinf(cimagf(z))) in cprojf() 40 return (z); in cprojf() 42 return (cpackf(INFINITY, copysignf(0.0, cimagf(z)))); in cprojf()
|
D | s_cprojl.c | 36 cprojl(long double complex z) in cprojl() argument 39 if (!isinf(creall(z)) && !isinf(cimagl(z))) in cprojl() 40 return (z); in cprojl() 42 return (cpackl(INFINITY, copysignl(0.0, cimagl(z)))); in cprojl()
|
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 …]
|
D | s_cproj.c | 36 cproj(double complex z) in cproj() argument 39 if (!isinf(creal(z)) && !isinf(cimag(z))) in cproj() 40 return (z); in cproj() 42 return (cpack(INFINITY, copysign(0.0, cimag(z)))); in cproj()
|
D | e_acosl.c | 47 long double z,p,q,r,w,s,c,df; in acosl() local 61 z = x*x; in acosl() 62 p = P(z); in acosl() 63 q = Q(z); in acosl() 67 z = (one+x)*0.5; in acosl() 68 p = P(z); in acosl() 69 q = Q(z); in acosl() 70 s = sqrtl(z); in acosl() 75 z = (one-x)*0.5; in acosl() 76 s = sqrtl(z); in acosl() [all …]
|
D | s_tanl.c | 54 union IEEEl2bits z; in tanl() local 59 z.e = x; in tanl() 60 s = z.bits.sign; in tanl() 61 z.bits.sign = 0; in tanl() 64 if (z.bits.exp == 0) in tanl() 68 if (z.bits.exp == 32767) in tanl() 74 if (z.e < M_PI_4) { in tanl() 75 hi = __kernel_tanl(z.e, 0, 0); in tanl()
|
D | s_erf.c | 193 double R,S,P,Q,s,y,z,r; in erf() local 207 z = x*x; in erf() 208 r = pp0+z*(pp1+z*(pp2+z*(pp3+z*pp4))); in erf() 209 s = one+z*(qq1+z*(qq2+z*(qq3+z*(qq4+z*qq5)))); in erf() 235 z = x; in erf() 236 SET_LOW_WORD(z,0); in erf() 237 r = __ieee754_exp(-z*z-0.5625)*__ieee754_exp((z-x)*(z+x)+R/S); in erf() 245 double R,S,P,Q,s,y,z,r; in erfc() local 256 z = x*x; in erfc() 257 r = pp0+z*(pp1+z*(pp2+z*(pp3+z*pp4))); in erfc() [all …]
|
D | w_cabsf.c | 18 cabsf(z) in cabsf() argument 19 float complex z; in cabsf() 22 return hypotf(crealf(z), cimagf(z));
|
D | e_lgammaf_r.c | 94 float y,z; in sin_pif() local 107 z = floorf(y); in sin_pif() 108 if(z!=y) { /* inexact anyway */ in sin_pif() 116 if(ix<0x4b000000) z = y+two23; /* exact */ in sin_pif() 117 GET_FLOAT_WORD(n,z); in sin_pif() 140 float t,y,z,nadj,p,p1,p2,p3,q,r,w; in __ieee754_lgammaf_r() local 184 z = y*y; in __ieee754_lgammaf_r() 185 p1 = a0+z*(a2+z*(a4+z*(a6+z*(a8+z*a10)))); in __ieee754_lgammaf_r() 186 p2 = z*(a1+z*(a3+z*(a5+z*(a7+z*(a9+z*a11))))); in __ieee754_lgammaf_r() 190 z = y*y; in __ieee754_lgammaf_r() [all …]
|
D | s_erff.c | 101 float R,S,P,Q,s,y,z,r; in erff() local 116 z = x*x; in erff() 117 r = pp0+z*(pp1+z*(pp2+z*(pp3+z*pp4))); in erff() 118 s = one+z*(qq1+z*(qq2+z*(qq3+z*(qq4+z*qq5)))); in erff() 145 SET_FLOAT_WORD(z,ix&0xfffff000); in erff() 146 r = __ieee754_expf(-z*z-(float)0.5625)*__ieee754_expf((z-x)*(z+x)+R/S); in erff() 154 float R,S,P,Q,s,y,z,r; in erfcf() local 165 z = x*x; in erfcf() 166 r = pp0+z*(pp1+z*(pp2+z*(pp3+z*pp4))); in erfcf() 167 s = one+z*(qq1+z*(qq2+z*(qq3+z*(qq4+z*qq5)))); in erfcf() [all …]
|
D | s_sinl.c | 48 union IEEEl2bits z; in sinl() local 53 z.e = x; in sinl() 54 s = z.bits.sign; in sinl() 55 z.bits.sign = 0; in sinl() 58 if (z.bits.exp == 0) in sinl() 62 if (z.bits.exp == 32767) in sinl() 68 if (z.e < M_PI_4) { in sinl() 69 hi = __kernel_sinl(z.e, 0, 0); in sinl()
|
/bionic/libm/upstream-freebsd/lib/msun/bsdsrc/ |
D | b_exp.c | 92 double z,hi,lo,c; 111 z=x*x; 112 c= x - z*(p1+z*(p2+z*(p3+z*(p4+z*p5)))); 138 double z,hi,lo; local 147 z = invln2*x; 148 k = z + copysign(.5, x); 155 z=x*x; 156 c= x - z*(p1+z*(p2+z*(p3+z*(p4+z*p5))));
|
D | b_tgamma.c | 159 double z, p; local 162 z = one/(x*x); 163 p = Pa0+z*(Pa1+z*(Pa2+z*(Pa3+z*(Pa4+z*(Pa5+z*(Pa6+z*Pa7)))))); 248 ratfun_gam(z, c) in ratfun_gam() argument 249 double z, c; in ratfun_gam() 254 q = Q0 +z*(Q1+z*(Q2+z*(Q3+z*(Q4+z*(Q5+z*(Q6+z*(Q7+z*Q8))))))); 255 p = P0 + z*(P1 + z*(P2 + z*(P3 + z*P4))); 259 t.a = z, TRUNC(t.a); /* t ~= z + c */ 260 t.b = (z - t.a) + c; 261 t.b *= (t.a + z); [all …]
|