Home
last modified time | relevance | path

Searched refs:k (Results 1 – 25 of 60) sorted by relevance

123

/sysdeps/ieee754/flt-32/
A Ds_scalbnf.c27 int32_t k,ix; in __scalbnf() local
29 k = (ix&0x7f800000)>>23; /* extract exponent */ in __scalbnf()
30 if (__builtin_expect(k==0, 0)) { /* 0 or subnormal x */ in __scalbnf()
34 k = ((ix&0x7f800000)>>23) - 25; in __scalbnf()
39 if (__builtin_expect(n> 50000 || k+n > 0xfe, 0)) in __scalbnf()
43 k = k+n; in __scalbnf()
44 if (__builtin_expect(k > 0, 1)) /* normal result */ in __scalbnf()
45 {SET_FLOAT_WORD(x,(ix&0x807fffff)|(k<<23)); return x;} in __scalbnf()
46 if (k <= -25) in __scalbnf()
48 k += 25; /* subnormal result */ in __scalbnf()
[all …]
A Ds_scalblnf.c27 int32_t k,ix; in __scalblnf() local
29 k = (ix&0x7f800000)>>23; /* extract exponent */ in __scalblnf()
30 if (__builtin_expect(k==0, 0)) { /* 0 or subnormal x */ in __scalblnf()
34 k = ((ix&0x7f800000)>>23) - 25; in __scalblnf()
39 if (__builtin_expect(n> 50000 || k+n > 0xfe, 0)) in __scalblnf()
43 k = k+n; in __scalblnf()
44 if (__builtin_expect(k > 0, 1)) /* normal result */ in __scalblnf()
45 {SET_FLOAT_WORD(x,(ix&0x807fffff)|(k<<23)); return x;} in __scalblnf()
46 if (k <= -25) in __scalblnf()
48 k += 25; /* subnormal result */ in __scalblnf()
[all …]
A Ds_log1pf.c40 int32_t k,hx,hu,ax; in __log1pf() local
45 k = 1; in __log1pf()
65 if(k!=0) { in __log1pf()
69 k = (hu>>23)-127; in __log1pf()
76 k = (hu>>23)-127; in __log1pf()
83 k += 1; in __log1pf()
92 if(k==0) return zero; in __log1pf()
93 else {c += k*ln2_lo; return k*ln2_hi+c;} in __log1pf()
96 if(k==0) return f-R; else in __log1pf()
97 return k*ln2_hi-((R-(k*ln2_lo+c))-f); in __log1pf()
[all …]
A Ds_expm1f.c43 int32_t k,xsb; in __expm1f() local
73 {hi = x - ln2_hi; lo = ln2_lo; k = 1;} in __expm1f()
75 {hi = x + ln2_hi; lo = -ln2_lo; k = -1;} in __expm1f()
77 k = invln2*x+((xsb==0)?(float)0.5:(float)-0.5); in __expm1f()
78 t = k; in __expm1f()
90 else k = 0; in __expm1f()
98 if(k==0) return x - (x*e-hxs); /* c is 0 */ in __expm1f()
103 if(k==1) { in __expm1f()
107 if (k <= -2 || k>56) { /* suffice to return exp(x)-1 */ in __expm1f()
115 if(k<23) { in __expm1f()
[all …]
A De_log10f.c30 int32_t i,k,hx; in __ieee754_log10f() local
34 k=0; in __ieee754_log10f()
40 k -= 25; x *= two25; /* subnormal number, scale up x */ in __ieee754_log10f()
44 k += (hx>>23)-127; in __ieee754_log10f()
45 i = ((uint32_t)k&0x80000000)>>31; in __ieee754_log10f()
47 y = (float)(k+i); in __ieee754_log10f()
/sysdeps/ieee754/ldbl-96/
A Ds_scalblnl.c34 int32_t k,es,hx,lx; in __scalblnl() local
36 k = es&0x7fff; /* extract exponent */ in __scalblnl()
37 if (__builtin_expect(k==0, 0)) { /* 0 or subnormal x */ in __scalblnl()
41 k = (es&0x7fff) - 63; in __scalblnl()
46 if (__builtin_expect(n> 50000 || k+n > 0x7ffe, 0)) in __scalblnl()
50 k = k+n; in __scalblnl()
51 if (__builtin_expect(k > 0, 1)) /* normal result */ in __scalblnl()
52 {SET_LDOUBLE_EXP(x,(es&0x8000)|k); return x;} in __scalblnl()
53 if (k <= -64) in __scalblnl()
55 k += 64; /* subnormal result */ in __scalblnl()
[all …]
/sysdeps/ieee754/dbl-64/
A Ds_scalbln.c32 int64_t k; in __scalbln() local
34 k = (ix >> 52) & 0x7ff; /* extract exponent */ in __scalbln()
35 if (__builtin_expect(k==0, 0)) { /* 0 or subnormal x */ in __scalbln()
39 k = ((ix >> 52) & 0x7ff) - 54; in __scalbln()
41 if (__builtin_expect(k==0x7ff, 0)) return x+x; /* NaN or Inf */ in __scalbln()
44 if (__builtin_expect(n> 50000 || k+n > 0x7fe, 0)) in __scalbln()
48 k = k+n; in __scalbln()
49 if (__builtin_expect(k > 0, 1)) /* normal result */ in __scalbln()
52 if (k <= -54) in __scalbln()
54 k += 54; /* subnormal result */ in __scalbln()
[all …]
A Ds_scalbn.c32 int64_t k; in __scalbn() local
34 k = (ix >> 52) & 0x7ff; /* extract exponent */ in __scalbn()
35 if (__builtin_expect(k==0, 0)) { /* 0 or subnormal x */ in __scalbn()
39 k = ((ix >> 52) & 0x7ff) - 54; in __scalbn()
41 if (__builtin_expect(k==0x7ff, 0)) return x+x; /* NaN or Inf */ in __scalbn()
44 if (__builtin_expect(n> 50000 || k+n > 0x7fe, 0)) in __scalbn()
48 k = k+n; in __scalbn()
49 if (__builtin_expect(k > 0, 1)) /* normal result */ in __scalbn()
52 if (k <= -54) in __scalbn()
54 k += 54; /* subnormal result */ in __scalbn()
[all …]
A Ds_log1p.c106 int32_t k, hx, hu, ax; in __log1p() local
111 k = 1; in __log1p()
134 k = 0; f = x; hu = 1; in __log1p()
139 if (k != 0) in __log1p()
163 k += 1; in __log1p()
174 if (k == 0) in __log1p()
178 c += k * ln2_lo; return k * ln2_hi + c; in __log1p()
182 if (k == 0) in __log1p()
185 return k * ln2_hi - ((R - (k * ln2_lo + c)) - f); in __log1p()
194 if (k == 0) in __log1p()
[all …]
A Ds_expm1.c137 int32_t k, xsb; in __expm1() local
182 hi = x - ln2_hi; lo = ln2_lo; k = 1; in __expm1()
186 hi = x + ln2_hi; lo = -ln2_lo; k = -1; in __expm1()
191 k = invln2 * x + ((xsb == 0) ? 0.5 : -0.5); in __expm1()
192 t = k; in __expm1()
206 k = 0; in __expm1()
217 if (k == 0) in __expm1()
223 if (k == -1) in __expm1()
225 if (k == 1) in __expm1()
232 if (k <= -2 || k > 56) /* suffice to return exp(x)-1 */ in __expm1()
[all …]
A Ds_sincos.c34 int k; in __sincos() local
39 k = u.i[HIGH_HALF] & 0x7fffffff; in __sincos()
41 if (k < 0x400368fd) in __sincos()
45 if (k < 0x3e400000) in __sincos()
47 if (k < 0x3e500000) in __sincos()
54 else if (k < 0x3feb6000) in __sincos()
70 if (k < 0x7ff00000) in __sincos()
76 n = k < 0x419921FB ? reduce_sincos (x, &a, &da) : __branred (x, &a, &da); in __sincos()
A De_log10.c62 int32_t k; in __ieee754_log10() local
66 k = 0; in __ieee754_log10()
73 k -= 54; in __ieee754_log10()
80 k += (hx >> 52) - 1023; in __ieee754_log10()
81 i = ((uint64_t) k & UINT64_C(0x8000000000000000)) >> 63; in __ieee754_log10()
83 y = (double) (k + i); in __ieee754_log10()
A Ds_sin.c204 int4 k, m, n; in __sin() local
218 else if (k < 0x3feb6000) in __sin()
225 else if (k < 0x400368fd) in __sin()
233 else if (k < 0x419921FB) in __sin()
240 else if (k < 0x7ff00000) in __sin()
268 int4 k, m, n; in __cos() local
276 k = 0x7fffffff & m; in __cos()
279 if (k < 0x3e400000) in __cos()
282 else if (k < 0x3feb6000) in __cos()
288 else if (k < 0x400368fd) in __cos()
[all …]
A De_asin.c50 int4 k,m,n; in __ieee754_asin() local
56 if (k < 0x3e500000) in __ieee754_asin()
63 if (k < 0x3fc00000) { in __ieee754_asin()
72 if (k<0x3fd00000) n = 11*((k&0x000fffff)>>15); in __ieee754_asin()
86 if (k < 0x3fe80000) { in __ieee754_asin()
150 k=v.i[HIGH_HALF]; in __ieee754_asin()
151 t=inroot[(k&0x001fffff)>>14]*powtwo[511-(k>>21)]; in __ieee754_asin()
194 int4 k,m,n; in libm_alias_finite() local
197 k = 0x7fffffff&m; in libm_alias_finite()
215 if (k<0x3fd00000) n = 11*((k&0x000fffff)>>15); in libm_alias_finite()
[all …]
/sysdeps/ieee754/ldbl-128/
A Ds_scalblnl.c38 int64_t k,hx,lx; in __scalblnl() local
40 k = (hx>>48)&0x7fff; /* extract exponent */ in __scalblnl()
41 if (k==0) { /* 0 or subnormal x */ in __scalblnl()
45 k = ((hx>>48)&0x7fff) - 114; in __scalblnl()
47 if (k==0x7fff) return x+x; /* NaN or Inf */ in __scalblnl()
49 if (n> 50000 || k+n > 0x7ffe) in __scalblnl()
53 k = k+n; in __scalblnl()
54 if (k > 0) /* normal result */ in __scalblnl()
56 if (k <= -114) in __scalblnl()
58 k += 114; /* subnormal result */ in __scalblnl()
[all …]
A Ds_scalbnl.c38 int64_t k,hx,lx; in __scalbnl() local
40 k = (hx>>48)&0x7fff; /* extract exponent */ in __scalbnl()
41 if (k==0) { /* 0 or subnormal x */ in __scalbnl()
45 k = ((hx>>48)&0x7fff) - 114; in __scalbnl()
47 if (k==0x7fff) return x+x; /* NaN or Inf */ in __scalbnl()
49 if (n> 50000 || k+n > 0x7ffe) in __scalbnl()
53 k = k+n; in __scalbnl()
54 if (k > 0) /* normal result */ in __scalbnl()
56 if (k <= -114) in __scalbnl()
58 k += 114; /* subnormal result */ in __scalbnl()
[all …]
A Dldbl2mpn.c103 int j, k, l; in __mpn_extract_long_double() local
118 for (k = N - 1; k >= l; k--) in __mpn_extract_long_double()
119 res_ptr[k] = res_ptr[k-l]; in __mpn_extract_long_double()
122 for (k = N - 1; k > l; k--) in __mpn_extract_long_double()
123 res_ptr[k] = res_ptr[k-l] << cnt in __mpn_extract_long_double()
124 | res_ptr[k-l-1] >> (BITS_PER_MP_LIMB - cnt); in __mpn_extract_long_double()
125 res_ptr[k--] = res_ptr[0] << cnt; in __mpn_extract_long_double()
128 for (; k >= 0; k--) in __mpn_extract_long_double()
129 res_ptr[k] = 0; in __mpn_extract_long_double()
A Ds_atanl.c179 int k, sign; in __atanl() local
184 k = s.parts32.w0; in __atanl()
185 if (k & 0x80000000) in __atanl()
191 k &= 0x7fffffff; in __atanl()
192 if (k >= 0x7fff0000) in __atanl()
205 if (k <= 0x3fc50000) /* |x| < 2**-58 */ in __atanl()
225 if (k >= 0x40024800) /* 10.25 */ in __atanl()
227 k = 83; in __atanl()
235 k = 8.0 * x + 0.25; in __atanl()
236 u = L(0.125) * k; in __atanl()
[all …]
/sysdeps/ieee754/ldbl-128ibm/
A Ds_atanl.c178 int32_t k, sign, lx; in __atanl() local
183 EXTRACT_WORDS (k, lx, xhi); in __atanl()
184 sign = k & 0x80000000; in __atanl()
187 k &= 0x7fffffff; in __atanl()
188 if (k >= 0x7ff00000) in __atanl()
191 if (((k - 0x7ff00000) | lx) != 0) in __atanl()
221 if (k >= 0x40248000) /* 10.25 */ in __atanl()
223 k = 83; in __atanl()
231 k = 8.0 * x + 0.25; in __atanl()
232 u = 0.125 * k; in __atanl()
[all …]
A Ds_scalblnl.c41 int64_t k,l,hx,lx; in __scalblnl() local
48 k = (hx>>52)&0x7ff; /* extract exponent */ in __scalblnl()
50 if (k==0) { /* 0 or subnormal x */ in __scalblnl()
55 k = ((hx>>52)&0x7ff) - 54; in __scalblnl()
59 if (n> 50000 || k+n > 0x7fe) in __scalblnl()
63 k = k+n; in __scalblnl()
64 if (k > 0) { /* normal result */ in __scalblnl()
65 hx = (hx&0x800fffffffffffffULL)|(k<<52); in __scalblnl()
94 if (k <= -54) in __scalblnl()
96 k += 54; /* subnormal result */ in __scalblnl()
[all …]
A Ds_scalbnl.c41 int64_t k,l,hx,lx; in __scalbnl() local
48 k = (hx>>52)&0x7ff; /* extract exponent */ in __scalbnl()
50 if (k==0) { /* 0 or subnormal x */ in __scalbnl()
55 k = ((hx>>52)&0x7ff) - 54; in __scalbnl()
59 if (n> 50000 || k+n > 0x7fe) in __scalbnl()
63 k = k+n; in __scalbnl()
64 if (k > 0) { /* normal result */ in __scalbnl()
65 hx = (hx&0x800fffffffffffffULL)|(k<<52); in __scalbnl()
94 if (k <= -54) in __scalbnl()
96 k += 54; /* subnormal result */ in __scalbnl()
[all …]
A De_logl.c195 int k, e; in __ieee754_logl() local
204 k = m & 0x7fffffff; in __ieee754_logl()
206 if ((k | lx) == 0) in __ieee754_logl()
216 if (k >= 0x7ff00000) in __ieee754_logl()
227 k = 64; in __ieee754_logl()
250 k = (m - 0xff000) >> 13; in __ieee754_logl()
253 INSERT_WORDS (xhi, 0x3ff00000 + (k << 13), 0); in __ieee754_logl()
257 k += 64; in __ieee754_logl()
261 k = (m - 0xfe000) >> 14; in __ieee754_logl()
262 INSERT_WORDS (xhi, 0x3fe00000 + (k << 14), 0); in __ieee754_logl()
[all …]
A Dldbl2mpn.c160 int j, k, l; in __mpn_extract_long_double() local
175 for (k = N - 1; k >= l; k--) in __mpn_extract_long_double()
176 res_ptr[k] = res_ptr[k-l]; in __mpn_extract_long_double()
179 for (k = N - 1; k > l; k--) in __mpn_extract_long_double()
180 res_ptr[k] = res_ptr[k-l] << cnt in __mpn_extract_long_double()
181 | res_ptr[k-l-1] >> (BITS_PER_MP_LIMB - cnt); in __mpn_extract_long_double()
182 res_ptr[k--] = res_ptr[0] << cnt; in __mpn_extract_long_double()
185 for (; k >= 0; k--) in __mpn_extract_long_double()
186 res_ptr[k] = 0; in __mpn_extract_long_double()
A De_sqrtl.c56 uint64_t k, l; in __ieee754_sqrtl() local
61 k=a.i[0] & INT64_C(0x7fffffffffffffff); in __ieee754_sqrtl()
63 if (k>INT64_C(0x000fffff00000000) && k<INT64_C(0x7ff0000000000000)) { in __ieee754_sqrtl()
65 l = (k&INT64_C(0x001fffffffffffff))|INT64_C(0x3fe0000000000000); in __ieee754_sqrtl()
67 n = (int64_t) ((l - k) * 2) >> 53; in __ieee754_sqrtl()
87 c.i[0] = INT64_C(0x2000000000000000)+((k&INT64_C(0x7fe0000000000000))>>1); in __ieee754_sqrtl()
95 if (k>=INT64_C(0x7ff0000000000000)) in __ieee754_sqrtl()
/sysdeps/i386/
A Dtst-auditmod3a.c24 audit1_test (int i, int j, int k) in audit1_test() argument
26 if (i != 1 || j != 2 || k != 3) in audit1_test()
33 audit2_test (int i, int j, int k) in audit2_test() argument
35 if (i != 1 || j != 2 || k != 3) in audit2_test()

Completed in 28 milliseconds

123