| /bionic/libm/src/ |
| D | e_acos.c | 64 double z,p,q,r,w,s,c,df; in __ieee754_acos() local
79 z = x*x; in __ieee754_acos()
80 p = z*(pS0+z*(pS1+z*(pS2+z*(pS3+z*(pS4+z*pS5))))); in __ieee754_acos()
81 q = one+z*(qS1+z*(qS2+z*(qS3+z*qS4))); in __ieee754_acos()
85 z = (one+x)*0.5; in __ieee754_acos()
86 p = z*(pS0+z*(pS1+z*(pS2+z*(pS3+z*(pS4+z*pS5))))); in __ieee754_acos()
87 q = one+z*(qS1+z*(qS2+z*(qS3+z*qS4))); in __ieee754_acos()
88 s = sqrt(z); in __ieee754_acos()
93 z = (one-x)*0.5; in __ieee754_acos()
94 s = sqrt(z); in __ieee754_acos()
[all …]
|
| D | e_acosf.c | 42 float z,p,q,r,w,s,c,df; in __ieee754_acosf() local
54 z = x*x; in __ieee754_acosf()
55 p = z*(pS0+z*(pS1+z*(pS2+z*(pS3+z*(pS4+z*pS5))))); in __ieee754_acosf()
56 q = one+z*(qS1+z*(qS2+z*(qS3+z*qS4))); in __ieee754_acosf()
60 z = (one+x)*(float)0.5; in __ieee754_acosf()
61 p = z*(pS0+z*(pS1+z*(pS2+z*(pS3+z*(pS4+z*pS5))))); in __ieee754_acosf()
62 q = one+z*(qS1+z*(qS2+z*(qS3+z*qS4))); in __ieee754_acosf()
63 s = __ieee754_sqrtf(z); in __ieee754_acosf()
69 z = (one-x)*(float)0.5; in __ieee754_acosf()
70 s = __ieee754_sqrtf(z); in __ieee754_acosf()
[all …]
|
| D | e_j1.c | 89 double z, s,c,ss,cc,r,u,v,y; in __ieee754_j1() local
102 z = cos(y+y); in __ieee754_j1()
103 if ((s*c)>zero) cc = z/ss; in __ieee754_j1()
104 else ss = z/cc; in __ieee754_j1()
110 if(ix>0x48000000) z = (invsqrtpi*cc)/sqrt(y); in __ieee754_j1()
113 z = invsqrtpi*(u*cc-v*ss)/sqrt(y); in __ieee754_j1()
115 if(hx<0) return -z; in __ieee754_j1()
116 else return z; in __ieee754_j1()
121 z = x*x; in __ieee754_j1()
122 r = z*(r00+z*(r01+z*(r02+z*r03))); in __ieee754_j1()
[all …]
|
| D | e_j1f.c | 47 float z, s,c,ss,cc,r,u,v,y; in __ieee754_j1f() local
60 z = cosf(y+y); in __ieee754_j1f()
61 if ((s*c)>zero) cc = z/ss; in __ieee754_j1f()
62 else ss = z/cc; in __ieee754_j1f()
68 if(((uint32_t)hx)>0x80000000) z = (invsqrtpi*cc)/sqrtf(y); in __ieee754_j1f()
71 z = invsqrtpi*(u*cc-v*ss)/sqrtf(y); in __ieee754_j1f()
73 if(hx<0) return -z; in __ieee754_j1f()
74 else return z; in __ieee754_j1f()
79 z = x*x; in __ieee754_j1f()
80 r = z*(r00+z*(r01+z*(r02+z*r03))); in __ieee754_j1f()
[all …]
|
| D | w_cabs.c | 17 cabs(z) in cabs() argument
18 double complex z; in cabs()
20 return hypot(creal(z), cimag(z));
24 z_abs(z) in z_abs() argument
25 double complex *z; in z_abs()
27 return hypot(creal(*z), cimag(*z));
|
| D | e_j0f.c | 46 float z, s,c,ss,cc,r,u,v; in __ieee754_j0f() local
59 z = -cosf(x+x); in __ieee754_j0f()
60 if ((s*c)<zero) cc = z/ss; in __ieee754_j0f()
61 else ss = z/cc; in __ieee754_j0f()
67 if(((uint32_t)hx)>0x80000000) z = (invsqrtpi*cc)/sqrtf(x); in __ieee754_j0f()
70 z = invsqrtpi*(u*cc-v*ss)/sqrtf(x); in __ieee754_j0f()
72 return z; in __ieee754_j0f()
80 z = x*x; in __ieee754_j0f()
81 r = z*(R02+z*(R03+z*(R04+z*R05))); in __ieee754_j0f()
82 s = one+z*(S01+z*(S02+z*(S03+z*S04))); in __ieee754_j0f()
[all …]
|
| D | e_j0.c | 88 double z, s,c,ss,cc,r,u,v; in __ieee754_j0() local
101 z = -cos(x+x); in __ieee754_j0()
102 if ((s*c)<zero) cc = z/ss; in __ieee754_j0()
103 else ss = z/cc; in __ieee754_j0()
109 if(ix>0x48000000) z = (invsqrtpi*cc)/sqrt(x); in __ieee754_j0()
112 z = invsqrtpi*(u*cc-v*ss)/sqrt(x); in __ieee754_j0()
114 return z; in __ieee754_j0()
122 z = x*x; in __ieee754_j0()
123 r = z*(R02+z*(R03+z*(R04+z*R05))); in __ieee754_j0()
124 s = one+z*(S01+z*(S02+z*(S03+z*S04))); in __ieee754_j0()
[all …]
|
| D | k_cos.c | 72 double hz,z,r,w; in __kernel_cos() local
74 z = x*x; in __kernel_cos()
75 r = z*(C1+z*(C2+z*(C3+z*(C4+z*(C5+z*C6))))); in __kernel_cos()
76 hz = (float)0.5*z; in __kernel_cos()
78 return w + (((one-w)-hz) + (z*r-x*y)); in __kernel_cos()
|
| D | k_sin.c | 63 double z,r,v; in __kernel_sin() local
65 z = x*x; in __kernel_sin()
66 v = z*x; in __kernel_sin()
67 r = S2+z*(S3+z*(S4+z*(S5+z*S6))); 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 | k_cosf.c | 40 double r, w, z; in __kernel_cosdf() local
43 z = x*x; in __kernel_cosdf()
44 w = z*z; in __kernel_cosdf()
45 r = C2+z*C3; in __kernel_cosdf()
46 return ((one+z*C0) + w*C1) + (w*z)*r; in __kernel_cosdf()
|
| D | k_sinf.c | 39 double r, s, w, z; in __kernel_sindf() local
42 z = x*x; in __kernel_sindf()
43 w = z*z; in __kernel_sindf()
44 r = S3+z*S4; in __kernel_sindf()
45 s = z*x; in __kernel_sindf()
46 return (x + s*(S1+z*S2)) + s*w*r; in __kernel_sindf()
|
| D | k_tanf.c | 42 double z,r,w,s,t,u; in __kernel_tandf() local
44 z = x*x; in __kernel_tandf()
59 r = T[4]+z*T[5]; in __kernel_tandf()
60 t = T[2]+z*T[3]; in __kernel_tandf()
61 w = z*z; in __kernel_tandf()
62 s = z*x; in __kernel_tandf()
63 u = T[0]+z*T[1]; in __kernel_tandf()
|
| D | w_cabsf.c | 18 cabsf(z) in cabsf() argument
19 float complex z; in cabsf()
22 return hypotf(crealf(z), cimagf(z));
|
| D | k_tan.c | 81 double z, r, v, w, s; in __kernel_tan() local
91 z = pio4 - x; in __kernel_tan()
93 x = z + w; in __kernel_tan()
96 z = x * x; in __kernel_tan()
97 w = z * z; in __kernel_tan()
105 v = z * (T[2] + w * (T[4] + w * (T[6] + w * (T[8] + w * (T[10] + in __kernel_tan()
107 s = z * x; in __kernel_tan()
108 r = y + z * (s * (r + v) + y); in __kernel_tan()
125 z = w; in __kernel_tan()
126 SET_LOW_WORD(z,0); in __kernel_tan()
[all …]
|
| D | e_rem_pio2.c | 78 double z,w,t,r,fn; in __ieee754_rem_pio2() local
89 z = x - pio2_1; in __ieee754_rem_pio2()
91 y[0] = z - pio2_1t; in __ieee754_rem_pio2()
92 y[1] = (z-y[0])-pio2_1t; in __ieee754_rem_pio2()
94 z -= pio2_2; in __ieee754_rem_pio2()
95 y[0] = z - pio2_2t; in __ieee754_rem_pio2()
96 y[1] = (z-y[0])-pio2_2t; in __ieee754_rem_pio2()
100 z = x + pio2_1; in __ieee754_rem_pio2()
102 y[0] = z + pio2_1t; in __ieee754_rem_pio2()
103 y[1] = (z-y[0])+pio2_1t; in __ieee754_rem_pio2()
[all …]
|
| D | s_erf.c | 194 double R,S,P,Q,s,y,z,r; in erf() local
208 z = x*x; in erf()
209 r = pp0+z*(pp1+z*(pp2+z*(pp3+z*pp4))); in erf()
210 s = one+z*(qq1+z*(qq2+z*(qq3+z*(qq4+z*qq5)))); in erf()
236 z = x; in erf()
237 SET_LOW_WORD(z,0); in erf()
238 r = __ieee754_exp(-z*z-0.5625)*__ieee754_exp((z-x)*(z+x)+R/S); in erf()
246 double R,S,P,Q,s,y,z,r; in erfc() local
257 z = x*x; in erfc()
258 r = pp0+z*(pp1+z*(pp2+z*(pp3+z*pp4))); in erfc()
[all …]
|
| D | e_lgammaf_r.c | 95 float y,z; in sin_pif() local
108 z = floorf(y); in sin_pif()
109 if(z!=y) { /* inexact anyway */ in sin_pif()
117 if(ix<0x4b000000) z = y+two23; /* exact */ in sin_pif()
118 GET_FLOAT_WORD(n,z); in sin_pif()
141 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 | e_atan2f.c | 34 float z; in __ieee754_atan2f() local
82 if(k > 60) z=pi_o_2+(float)0.5*pi_lo; /* |y/x| > 2**60 */ in __ieee754_atan2f()
83 else if(hx<0&&k<-60) z=0.0; /* |y|/x < -2**60 */ in __ieee754_atan2f()
84 else z=atanf(fabsf(y/x)); /* safe to do y/x */ in __ieee754_atan2f()
86 case 0: return z ; /* atan(+,+) */ in __ieee754_atan2f()
89 GET_FLOAT_WORD(zh,z); in __ieee754_atan2f()
90 SET_FLOAT_WORD(z,zh ^ 0x80000000); in __ieee754_atan2f()
92 return z ; /* atan(-,+) */ in __ieee754_atan2f()
93 case 2: return pi-(z-pi_lo);/* atan(+,-) */ in __ieee754_atan2f()
95 return (z-pi_lo)-pi;/* atan(-,-) */ in __ieee754_atan2f()
|
| D | s_erff.c | 102 float R,S,P,Q,s,y,z,r; in erff() local
117 z = x*x; in erff()
118 r = pp0+z*(pp1+z*(pp2+z*(pp3+z*pp4))); in erff()
119 s = one+z*(qq1+z*(qq2+z*(qq3+z*(qq4+z*qq5)))); in erff()
146 SET_FLOAT_WORD(z,ix&0xfffff000); in erff()
147 r = __ieee754_expf(-z*z-(float)0.5625)*__ieee754_expf((z-x)*(z+x)+R/S); in erff()
155 float R,S,P,Q,s,y,z,r; in erfcf() local
166 z = x*x; in erfcf()
167 r = pp0+z*(pp1+z*(pp2+z*(pp3+z*pp4))); in erfcf()
168 s = one+z*(qq1+z*(qq2+z*(qq3+z*(qq4+z*qq5)))); in erfcf()
[all …]
|
| D | s_fma.c | 53 fma(double x, double y, double z) in fma() argument
63 if (z == 0.0) in fma()
66 return (x * y + z); in fma()
69 if (!isfinite(x) || !isfinite(y) || !isfinite(z)) in fma()
70 return (x * y + z); in fma()
74 zs = frexp(z, &ez); in fma()
90 if (x > 0.0 ^ y < 0.0 ^ z < 0.0) in fma()
99 if (z > 0.0) in fma()
108 if (z < 0.0) in fma()
120 if (!isnormal(z)) in fma()
[all …]
|
| D | e_lgamma_r.c | 161 double y,z; in sin_pi() local
174 z = floor(y); in sin_pi()
175 if(z!=y) { /* inexact anyway */ in sin_pi()
183 if(ix<0x43300000) z = y+two52; /* exact */ in sin_pi()
184 GET_LOW_WORD(n,z); in sin_pi()
207 double t,y,z,nadj,p,p1,p2,p3,q,r,w; in __ieee754_lgamma_r() local
250 z = y*y; in __ieee754_lgamma_r()
251 p1 = a0+z*(a2+z*(a4+z*(a6+z*(a8+z*a10)))); in __ieee754_lgamma_r()
252 p2 = z*(a1+z*(a3+z*(a5+z*(a7+z*(a9+z*a11))))); in __ieee754_lgamma_r()
256 z = y*y; in __ieee754_lgamma_r()
[all …]
|
| D | s_conjf.c | 32 conjf(float complex z) in conjf() argument
34 return crealf(z) - I * cimagf(z); in conjf()
|
| D | s_conjl.c | 32 conjl(long double complex z) in conjl() argument
34 return creall(z) - I * cimagl(z); in conjl()
|
| D | s_conj.c | 32 conj(double complex z) in conj() argument
34 return creal(z) - I * cimag(z); in conj()
|
| /bionic/libm/bsdsrc/ |
| D | b_exp.c | 94 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))));
|