| /sysdeps/ieee754/flt-32/ |
| A D | k_tanf.c | 45 float z,r,v,w,s; in __kernel_tanf() local
64 w = pio4lo-y; in __kernel_tanf()
65 x = z+w; y = 0.0; in __kernel_tanf()
70 w = z*z; in __kernel_tanf()
75 r = T[1]+w*(T[3]+w*(T[5]+w*(T[7]+w*(T[9]+w*T[11])))); in __kernel_tanf()
76 v = z*(T[2]+w*(T[4]+w*(T[6]+w*(T[8]+w*(T[10]+w*T[12]))))); in __kernel_tanf()
80 w = x+r; in __kernel_tanf()
83 return (float)(1-((hx>>30)&2))*(v-(float)2.0*(x-(w*w/(w+v)-r))); in __kernel_tanf()
85 if(iy==1) return w; in __kernel_tanf()
91 z = w; in __kernel_tanf()
[all …]
|
| A D | s_atanf.c | 59 float w,s1,s2,z; in __atanf() local
92 w = z*z; in __atanf()
94 s1 = z*(aT[0]+w*(aT[2]+w*(aT[4]+w*(aT[6]+w*(aT[8]+w*aT[10]))))); in __atanf()
95 s2 = w*(aT[1]+w*(aT[3]+w*(aT[5]+w*(aT[7]+w*aT[9])))); in __atanf()
|
| A D | e_asinf.c | 66 float t,w,p,q,c,r,s; in __ieee754_asinf() local
81 w = t * (p0 + t * (p1 + t * (p2 + t * (p3 + t * p4)))); in __ieee754_asinf()
82 return x+x*w; in __ieee754_asinf()
86 w = one-fabsf(x); in __ieee754_asinf()
87 t = w*0.5f; in __ieee754_asinf()
94 w = s; in __ieee754_asinf()
95 GET_FLOAT_WORD(iw,w); in __ieee754_asinf()
96 SET_FLOAT_WORD(w,iw&0xfffff000); in __ieee754_asinf()
97 c = (t-w*w)/(s+w); in __ieee754_asinf()
100 q = pio4_hi-2.0f*w; in __ieee754_asinf()
|
| A D | e_coshf.c | 26 float t,w; in __ieee754_coshf() local
38 w = one+t; in __ieee754_coshf()
39 return one+(t*t)/(w+w); in __ieee754_coshf()
52 w = __ieee754_expf(half*fabsf(x)); in __ieee754_coshf()
53 t = half*w; in __ieee754_coshf()
54 return t*w; in __ieee754_coshf()
|
| A D | s_rintf.c | 34 float w, t; in __rintf()
42 w = TWO23[sx] + x; in __rintf()
43 t = w - TWO23[sx]; in __rintf()
56 w = TWO23[sx] + x; in __rintf()
57 return w - TWO23[sx]; in __rintf()
|
| A D | s_asinhf.c | 29 float w; in __asinhf() local
39 w = __ieee754_logf(fabsf(x))+ln2; in __asinhf()
43 w = __ieee754_logf(2.0f*xa+one/(sqrtf(xa*xa+one)+xa)); in __asinhf()
46 w =__log1pf(xa+t/(one+sqrtf(one+t))); in __asinhf()
49 return copysignf(w, x); in __asinhf()
|
| /sysdeps/ieee754/ldbl-128ibm/ |
| A D | e_hypotl.c | 71 w = a+b; /* for sNaN */ in __ieee754_hypotl()
73 return w; in __ieee754_hypotl()
75 w = a; in __ieee754_hypotl()
77 w = b; in __ieee754_hypotl()
78 return w; in __ieee754_hypotl()
101 w = a-b; in __ieee754_hypotl()
102 if (w>b) { in __ieee754_hypotl()
129 w = sqrtl(a1*b1-(w*(-w)-(a1*b2+a2*b))); in __ieee754_hypotl()
133 w *= kld; in __ieee754_hypotl()
135 return w; in __ieee754_hypotl()
[all …]
|
| A D | k_tanl.c | 91 long double z, r, v, w, s; in __kernel_tanl() local
124 w = pio4lo - y; in __kernel_tanl()
125 x = z + w; in __kernel_tanl()
136 w = x + r; in __kernel_tanl()
140 w = (v - 2.0 * (x - (w * w / (w + v) - r))); in __kernel_tanl()
148 w = -w; in __kernel_tanl()
150 return w; in __kernel_tanl()
153 return w; in __kernel_tanl()
160 u1 = ldbl_high (w); in __kernel_tanl()
162 z = -1.0 / w; in __kernel_tanl()
|
| A D | e_coshl.c | 43 long double t,w; in __ieee754_coshl() local
59 w = one+t; in __ieee754_coshl()
60 return one+(t*t)/(w+w); in __ieee754_coshl()
74 w = __ieee754_expl(half*fabsl(x)); in __ieee754_coshl()
75 t = half*w; in __ieee754_coshl()
76 return t*w; in __ieee754_coshl()
|
| A D | s_asinhl.c | 41 long double t,w; in __asinhl() local
54 w = __ieee754_logl(fabsl(x))+ln2; in __asinhl()
57 w = __ieee754_logl(2.0*t+one/(sqrtl(x*x+one)+t)); in __asinhl()
60 w =__log1pl(fabsl(x)+t/(one+sqrtl(one+t))); in __asinhl()
62 if(hx>0) return w; else return -w; in __asinhl()
|
| A D | e_sinhl.c | 42 long double t,w,h; in __ieee754_sinhl() local
64 w = t/(t+one); in __ieee754_sinhl()
65 return h*(t+w); in __ieee754_sinhl()
73 w = __ieee754_expl(0.5*fabsl(x)); in __ieee754_sinhl()
74 t = h*w; in __ieee754_sinhl()
75 return t*w; in __ieee754_sinhl()
|
| A D | e_asinl.c | 137 long double a, t, w, p, q, c, r, s; in __ieee754_asinl() local
199 w = one - a; in __ieee754_asinl()
200 t = w * 0.5; in __ieee754_asinl()
227 w = p / q; in __ieee754_asinl()
228 return x + x * w; in __ieee754_asinl()
234 w = p / q; in __ieee754_asinl()
235 t = pio2_hi - (2.0 * (s + s * w) - pio2_lo); in __ieee754_asinl()
239 w = ldbl_high (s); in __ieee754_asinl()
240 c = (t - w * w) / (s + w); in __ieee754_asinl()
243 q = pio4_hi - 2.0 * w; in __ieee754_asinl()
|
| /sysdeps/ieee754/ldbl-96/ |
| A D | k_tanl.c | 91 long double z, r, v, w, s; in __kernel_tanl() local
121 w = pio4lo - y; in __kernel_tanl()
122 x = z + w; in __kernel_tanl()
133 w = x + r; in __kernel_tanl()
137 w = (v - 2.0 * (x - (w * w / (w + v) - r))); in __kernel_tanl()
145 w = -w; in __kernel_tanl()
147 return w; in __kernel_tanl()
150 return w; in __kernel_tanl()
|
| A D | e_asinl.c | 99 long double t, w, p, q, c, r, s; in __ieee754_asinl() local
128 w = p / q; in __ieee754_asinl()
129 return x + x * w; in __ieee754_asinl()
133 w = one - fabsl (x); in __ieee754_asinl()
134 t = w * 0.5; in __ieee754_asinl()
140 w = p / q; in __ieee754_asinl()
141 t = pio2_hi - (2.0 * (s + s * w) - pio2_lo); in __ieee754_asinl()
147 SET_LDOUBLE_WORDS (w,k,i0,i1); in __ieee754_asinl()
148 c = (t - w * w) / (s + w); in __ieee754_asinl()
151 q = pio4_hi - 2.0 * w; in __ieee754_asinl()
|
| A D | e_coshl.c | 46 long double t,w; in __ieee754_coshl() local
60 w = one+t; in __ieee754_coshl()
61 return one+(t*t)/(w+w); in __ieee754_coshl()
77 w = __ieee754_expl(half*fabsl(x)); in __ieee754_coshl()
78 t = half*w; in __ieee754_coshl()
79 return t*w; in __ieee754_coshl()
|
| A D | s_fmal.c | 39 w.d = z; in __fmal()
80 if (w.ieee.exponent >= 3) in __fmal()
91 || (w.ieee.exponent == 1 in __fmal()
165 if (w.ieee.exponent) in __fmal()
168 w.d *= 0x1p130L; in __fmal()
176 z = w.d; in __fmal()
278 w.d = a1 + u.d; in __fmal()
280 return w.d * 0x1p-130L; in __fmal()
285 w.d = 0.0L; in __fmal()
290 w.d *= 0x1p-2L; in __fmal()
[all …]
|
| A D | s_asinhl.c | 43 long double t,w; in __asinhl() local
53 w = __ieee754_logl(fabsl(x))+ln2; in __asinhl()
57 w = __ieee754_logl(2.0*xa+one/(sqrtl(xa*xa+one)+xa)); in __asinhl()
60 w =__log1pl(xa+t/(one+sqrtl(one+t))); in __asinhl()
63 return copysignl(w, x); in __asinhl()
|
| /sysdeps/ieee754/ldbl-128/ |
| A D | k_tanl.c | 91 _Float128 z, r, v, w, s; in __kernel_tanl() local
124 w = pio4lo - y; in __kernel_tanl()
125 x = z + w; in __kernel_tanl()
136 w = x + r; in __kernel_tanl()
140 w = (v - 2.0 * (x - (w * w / (w + v) - r))); in __kernel_tanl()
148 w = -w; in __kernel_tanl()
150 return w; in __kernel_tanl()
153 return w; in __kernel_tanl()
158 u1.value = w; in __kernel_tanl()
162 z = -1.0 / w; in __kernel_tanl()
|
| A D | e_coshl.c | 64 _Float128 t, w; in __ieee754_coshl() local
84 w = one + t; in __ieee754_coshl()
86 return one + (t * t) / (w + w); in __ieee754_coshl()
103 w = __ieee754_expl (half * u.value); in __ieee754_coshl()
104 t = half * w; in __ieee754_coshl()
105 return t * w; in __ieee754_coshl()
|
| A D | s_fmal.c | 47 w.d = z; in __fmal()
88 if (w.ieee.exponent >= 3) in __fmal()
99 || (w.ieee.exponent == 1 in __fmal()
175 if (w.ieee.exponent) in __fmal()
178 w.d *= L(0x1p228); in __fmal()
186 z = w.d; in __fmal()
288 w.d = a1 + u.d; in __fmal()
290 return w.d * L(0x1p-228); in __fmal()
295 w.d = 0; in __fmal()
300 w.d *= L(0x1p-2); in __fmal()
[all …]
|
| A D | s_asinhl.c | 44 _Float128 t, w; in __asinhl() local
62 w = __ieee754_logl (u.value) + ln2; in __asinhl()
67 w = __ieee754_logl (2.0 * t + one / (sqrtl (x * x + one) + t)); in __asinhl()
72 w = __log1pl (u.value + t / (one + sqrtl (one + t))); in __asinhl()
75 return -w; in __asinhl()
77 return w; in __asinhl()
|
| A D | s_rintl.c | 49 _Float128 w, t; in __rintl()
57 w = TWO112[sx] + x; in __rintl()
58 t = w - TWO112[sx]; in __rintl()
71 w = TWO112[sx] + x; in __rintl()
72 return w - TWO112[sx]; in __rintl()
|
| /sysdeps/ieee754/dbl-64/ |
| A D | e_cosh.c | 43 double t,w; in __ieee754_cosh() local
57 w = one+t; in __ieee754_cosh()
58 return one+(t*t)/(w+w); in __ieee754_cosh()
74 w = __ieee754_exp(half*fabs(x)); in __ieee754_cosh()
75 t = half*w; in __ieee754_cosh()
76 return t*w; in __ieee754_cosh()
|
| A D | s_fma.c | 50 w.d = z; in __fma()
89 if (w.ieee.exponent >= 3) in __fma()
100 || (w.ieee.exponent == 1 in __fma()
174 if (w.ieee.exponent) in __fma()
177 w.d *= 0x1p108; in __fma()
185 z = w.d; in __fma()
287 w.d = a1 + u.d; in __fma()
289 return w.d * 0x1p-108; in __fma()
294 w.d = 0.0; in __fma()
299 w.d *= 0x1p-2; in __fma()
[all …]
|
| A D | s_asinh.c | 38 double w; in __asinh() local
52 w = __ieee754_log (fabs (x)) + ln2; in __asinh()
59 w = __ieee754_log (2.0 * xa + one / (sqrt (xa * xa + one) + in __asinh()
65 w = __log1p (xa + t / (one + sqrt (one + t))); in __asinh()
68 return copysign (w, x); in __asinh()
|