Home
last modified time | relevance | path

Searched refs:ix (Results 1 – 25 of 151) sorted by relevance

1234567

/glibc-2.36/math/
Dk_casinh_template.c35 FLOAT rx, ix; in M_DECL_FUNC() local
40 ix = M_FABS (__imag__ x); in M_DECL_FUNC()
42 if (rx >= 1 / M_EPSILON || ix >= 1 / M_EPSILON) in M_DECL_FUNC()
49 __imag__ y = ix; in M_DECL_FUNC()
61 else if (rx >= M_LIT (0.5) && ix < M_EPSILON / 8) in M_DECL_FUNC()
69 __imag__ res = M_ATAN2 (ix, s); in M_DECL_FUNC()
71 else if (rx < M_EPSILON / 8 && ix >= M_LIT (1.5)) in M_DECL_FUNC()
73 FLOAT s = M_SQRT ((ix + 1) * (ix - 1)); in M_DECL_FUNC()
75 __real__ res = M_LOG (ix + s); in M_DECL_FUNC()
81 else if (ix > 1 && ix < M_LIT (1.5) && rx < M_LIT (0.5)) in M_DECL_FUNC()
[all …]
/glibc-2.36/sysdeps/ieee754/flt-32/
De_sqrtf.c29 int32_t ix,s,q,m,t,i; in __ieee754_sqrtf()
32 GET_FLOAT_WORD(ix,x); in __ieee754_sqrtf()
35 if((ix&0x7f800000)==0x7f800000) { in __ieee754_sqrtf()
40 if(ix<=0) { in __ieee754_sqrtf()
41 if((ix&(~sign))==0) return x;/* sqrt(+-0) = +-0 */ in __ieee754_sqrtf()
42 else if(ix<0) in __ieee754_sqrtf()
46 m = (ix>>23); in __ieee754_sqrtf()
48 for(i=0;(ix&0x00800000)==0;i++) ix<<=1; in __ieee754_sqrtf()
52 ix = (ix&0x007fffff)|0x00800000; in __ieee754_sqrtf()
54 ix += ix; in __ieee754_sqrtf()
[all …]
De_powf.c45 log2_inline (uint32_t ix) in log2_inline() argument
55 tmp = ix - OFF; in log2_inline()
58 iz = ix - top; in log2_inline()
141 zeroinfnan (uint32_t ix) in zeroinfnan() argument
143 return 2 * ix - 1 >= 2u * 0x7f800000 - 1; in zeroinfnan()
150 uint32_t ix, iy; in __powf() local
152 ix = asuint (x); in __powf()
154 if (__glibc_unlikely (ix - 0x00800000 >= 0x7f800000 - 0x00800000 in __powf()
162 if (ix == 0x3f800000) in __powf()
164 if (2 * ix > 2u * 0x7f800000 || 2 * iy > 2u * 0x7f800000) in __powf()
[all …]
De_j0f.c257 int32_t hx,ix; in __ieee754_j0f() local
260 ix = hx&0x7fffffff; in __ieee754_j0f()
261 if(ix>=0x7f800000) return one/(x*x); in __ieee754_j0f()
263 if(ix >= 0x40000000) { /* |x| >= 2.0 */ in __ieee754_j0f()
268 if (ix >= 0x7f000000) in __ieee754_j0f()
279 if (ix <= 0x5c000000) in __ieee754_j0f()
294 if(ix<0x39000000) { /* |x| < 2**-13 */ in __ieee754_j0f()
296 if(ix<0x32000000) return one; /* |x|<2**-27 */ in __ieee754_j0f()
302 if(ix < 0x3F800000) { /* |x| < 1.00 */ in __ieee754_j0f()
541 int32_t hx,ix; in __ieee754_y0f() local
[all …]
De_lgammaf_r.c95 int n,ix; in sin_pif() local
97 GET_FLOAT_WORD(ix,x); in sin_pif()
98 ix &= 0x7fffffff; in sin_pif()
100 if(ix<0x3e800000) return __sinf (pi*x); in sin_pif()
113 if(ix>=0x4b800000) { in sin_pif()
116 if(ix<0x4b000000) z = y+two23; /* exact */ in sin_pif()
141 int i,hx,ix; in __ieee754_lgammaf_r() local
147 ix = hx&0x7fffffff; in __ieee754_lgammaf_r()
148 if(__builtin_expect(ix>=0x7f800000, 0)) return x*x; in __ieee754_lgammaf_r()
149 if(__builtin_expect(ix==0, 0)) in __ieee754_lgammaf_r()
[all …]
De_j1f.c266 int32_t hx,ix; in __ieee754_j1f() local
269 ix = hx&0x7fffffff; in __ieee754_j1f()
270 if(__builtin_expect(ix>=0x7f800000, 0)) return one/x; in __ieee754_j1f()
272 if(ix >= 0x40000000) { /* |x| >= 2.0 */ in __ieee754_j1f()
277 if (ix >= 0x7f000000) in __ieee754_j1f()
288 if (ix <= 0x5c000000) in __ieee754_j1f()
305 if(__builtin_expect(ix<0x32000000, 0)) { /* |x|<2**-27 */ in __ieee754_j1f()
553 int32_t hx,ix; in __ieee754_y1f() local
556 ix = 0x7fffffff&hx; in __ieee754_y1f()
558 if(__builtin_expect(ix>=0x7f800000, 0)) return one/(x+x*x); in __ieee754_y1f()
[all …]
Ds_setpayloadf_main.c33 uint32_t ix; in FUNC() local
34 GET_FLOAT_WORD (ix, payload); in FUNC()
35 int exponent = ix >> EXPLICIT_MANT_DIG; in FUNC()
39 || (exponent < BIAS && !(SET_HIGH_BIT && ix == 0)) in FUNC()
40 || (ix & ((1U << (BIAS + EXPLICIT_MANT_DIG - exponent)) - 1)) != 0) in FUNC()
45 if (ix != 0) in FUNC()
47 ix &= (1U << EXPLICIT_MANT_DIG) - 1; in FUNC()
48 ix |= 1U << EXPLICIT_MANT_DIG; in FUNC()
49 ix >>= BIAS + EXPLICIT_MANT_DIG - exponent; in FUNC()
51 ix |= 0x7f800000 | (SET_HIGH_BIT ? 0x400000 : 0); in FUNC()
[all …]
Ds_fromfpf_main.c43 uint32_t ix; in FUNC() local
44 GET_FLOAT_WORD (ix, x); in FUNC()
45 bool negative = (ix & 0x80000000) != 0; in FUNC()
48 ix &= 0x7fffffff; in FUNC()
49 if (ix == 0) in FUNC()
51 int exponent = ix >> (MANT_DIG - 1); in FUNC()
57 ix &= ((1U << (MANT_DIG - 1)) - 1); in FUNC()
58 ix |= 1U << (MANT_DIG - 1); in FUNC()
63 uret = ix; in FUNC()
71 half_bit = (ix & h) != 0; in FUNC()
[all …]
Ds_erff.c104 int32_t hx,ix,i; in __erff() local
107 ix = hx&0x7fffffff; in __erff()
108 if(ix>=0x7f800000) { /* erf(nan)=nan */ in __erff()
113 if(ix < 0x3f580000) { /* |x|<0.84375 */ in __erff()
114 if(ix < 0x31800000) { /* |x|<2**-28 */ in __erff()
115 if (ix < 0x04000000) in __erff()
130 if(ix < 0x3fa00000) { /* 0.84375 <= |x| < 1.25 */ in __erff()
136 if (ix >= 0x40c00000) { /* inf>|x|>=6 */ in __erff()
141 if(ix< 0x4036DB6E) { /* |x| < 1/0.35 */ in __erff()
152 GET_FLOAT_WORD(ix,x); in __erff()
[all …]
De_log2f.c43 uint32_t ix, iz, top, tmp; in __log2f() local
46 ix = asuint (x); in __log2f()
49 if (__glibc_unlikely (ix == 0x3f800000)) in __log2f()
52 if (__glibc_unlikely (ix - 0x00800000 >= 0x7f800000 - 0x00800000)) in __log2f()
55 if (ix * 2 == 0) in __log2f()
57 if (ix == 0x7f800000) /* log2(inf) == inf. */ in __log2f()
59 if ((ix & 0x80000000) || ix * 2 >= 0xff000000) in __log2f()
62 ix = asuint (x * 0x1p23f); in __log2f()
63 ix -= 23 << 23; in __log2f()
69 tmp = ix - OFF; in __log2f()
[all …]
De_logf.c44 uint32_t ix, iz, tmp; in __logf() local
47 ix = asuint (x); in __logf()
50 if (__glibc_unlikely (ix == 0x3f800000)) in __logf()
53 if (__glibc_unlikely (ix - 0x00800000 >= 0x7f800000 - 0x00800000)) in __logf()
56 if (ix * 2 == 0) in __logf()
58 if (ix == 0x7f800000) /* log(inf) == inf. */ in __logf()
60 if ((ix & 0x80000000) || ix * 2 >= 0xff000000) in __logf()
63 ix = asuint (x * 0x1p23f); in __logf()
64 ix -= 23 << 23; in __logf()
70 tmp = ix - OFF; in __logf()
[all …]
Ds_roundevenf.c37 uint32_t ix, ux; in __roundevenf()
38 GET_FLOAT_WORD (ix, x); in __roundevenf()
39 ux = ix & 0x7fffffff; in __roundevenf()
60 if ((ix & (int_bit | (half_bit - 1))) != 0) in __roundevenf()
63 ix += half_bit; in __roundevenf()
64 ix &= ~(int_bit - 1); in __roundevenf()
68 ix = (ix & 0x80000000) | 0x3f800000; in __roundevenf()
71 ix &= 0x80000000; in __roundevenf()
72 SET_FLOAT_WORD (x, ix); in __roundevenf()
/glibc-2.36/sysdeps/ieee754/dbl-64/
De_j0.c89 int32_t hx, ix; in __ieee754_j0() local
92 ix = hx & 0x7fffffff; in __ieee754_j0()
93 if (ix >= 0x7ff00000) in __ieee754_j0()
96 if (ix >= 0x40000000) /* |x| >= 2.0 */ in __ieee754_j0()
101 if (ix < 0x7fe00000) /* make sure x+x not overflow */ in __ieee754_j0()
113 if (ix > 0x48000000) in __ieee754_j0()
122 if (ix < 0x3f200000) /* |x| < 2**-13 */ in __ieee754_j0()
125 if (ix < 0x3e400000) in __ieee754_j0()
137 if (ix < 0x3FF00000) /* |x| < 1.00 */ in __ieee754_j0()
166 int32_t hx, ix, lx; in __ieee754_y0() local
[all …]
De_j1.c93 int32_t hx, ix; in __ieee754_j1() local
96 ix = hx & 0x7fffffff; in __ieee754_j1()
97 if (__glibc_unlikely (ix >= 0x7ff00000)) in __ieee754_j1()
100 if (ix >= 0x40000000) /* |x| >= 2.0 */ in __ieee754_j1()
105 if (ix < 0x7fe00000) /* make sure y+y not overflow */ in __ieee754_j1()
117 if (ix > 0x48000000) in __ieee754_j1()
129 if (__glibc_unlikely (ix < 0x3e400000)) /* |x|<2**-27 */ in __ieee754_j1()
172 int32_t hx, ix, lx; in __ieee754_y1() local
175 ix = 0x7fffffff & hx; in __ieee754_y1()
177 if (__glibc_unlikely (ix >= 0x7ff00000)) in __ieee754_y1()
[all …]
De_lgamma_r.c160 int n,ix; in sin_pi() local
162 GET_HIGH_WORD(ix,x); in sin_pi()
163 ix &= 0x7fffffff; in sin_pi()
165 if(ix<0x3fd00000) return __sin(pi*x); in sin_pi()
178 if(ix>=0x43400000) { in sin_pi()
181 if(ix<0x43300000) z = y+two52; /* exact */ in sin_pi()
206 int i,hx,lx,ix; in __ieee754_lgamma_r() local
212 ix = hx&0x7fffffff; in __ieee754_lgamma_r()
213 if(__builtin_expect(ix>=0x7ff00000, 0)) return x*x; in __ieee754_lgamma_r()
214 if(__builtin_expect((ix|lx)==0, 0)) in __ieee754_lgamma_r()
[all …]
Ds_setpayload_main.c33 uint64_t ix; in FUNC() local
34 EXTRACT_WORDS64 (ix, payload); in FUNC()
35 int exponent = ix >> EXPLICIT_MANT_DIG; in FUNC()
39 || (exponent < BIAS && !(SET_HIGH_BIT && ix == 0)) in FUNC()
40 || (ix & ((1ULL << (BIAS + EXPLICIT_MANT_DIG - exponent)) - 1)) != 0) in FUNC()
45 if (ix != 0) in FUNC()
47 ix &= (1ULL << EXPLICIT_MANT_DIG) - 1; in FUNC()
48 ix |= 1ULL << EXPLICIT_MANT_DIG; in FUNC()
49 ix >>= BIAS + EXPLICIT_MANT_DIG - exponent; in FUNC()
51 ix |= 0x7ff0000000000000ULL | (SET_HIGH_BIT ? 0x8000000000000ULL : 0); in FUNC()
[all …]
Ds_fromfp_main.c43 uint64_t ix; in FUNC() local
44 EXTRACT_WORDS64 (ix, x); in FUNC()
45 bool negative = (ix & 0x8000000000000000ULL) != 0; in FUNC()
48 ix &= 0x7fffffffffffffffULL; in FUNC()
49 if (ix == 0) in FUNC()
51 int exponent = ix >> (MANT_DIG - 1); in FUNC()
57 ix &= ((1ULL << (MANT_DIG - 1)) - 1); in FUNC()
58 ix |= 1ULL << (MANT_DIG - 1); in FUNC()
63 uret = ix; in FUNC()
71 half_bit = (ix & h) != 0; in FUNC()
[all …]
Ds_roundeven.c36 uint64_t ix, ux; in __roundeven()
37 EXTRACT_WORDS64 (ix, x); in __roundeven()
38 ux = ix & 0x7fffffffffffffffULL; in __roundeven()
59 if ((ix & (int_bit | (half_bit - 1))) != 0) in __roundeven()
62 ix += half_bit; in __roundeven()
63 ix &= ~(int_bit - 1); in __roundeven()
67 ix = (ix & 0x8000000000000000ULL) | 0x3ff0000000000000ULL; in __roundeven()
70 ix &= 0x8000000000000000ULL; in __roundeven()
71 INSERT_WORDS64 (x, ix); in __roundeven()
/glibc-2.36/sysdeps/ieee754/ldbl-96/
De_j0l.c112 int32_t ix; in __ieee754_j0l() local
116 ix = se & 0x7fff; in __ieee754_j0l()
117 if (__glibc_unlikely (ix >= 0x7fff)) in __ieee754_j0l()
120 if (ix >= 0x4000) /* |x| >= 2.0 */ in __ieee754_j0l()
125 if (ix < 0x7ffe) in __ieee754_j0l()
137 if (__glibc_unlikely (ix > 0x408e)) /* 2^143 */ in __ieee754_j0l()
147 if (__glibc_unlikely (ix < 0x3fef)) /* |x| < 2**-16 */ in __ieee754_j0l()
151 if (ix < 0x3fde) /* |x| < 2^-33 */ in __ieee754_j0l()
159 if (ix < 0x3fff) in __ieee754_j0l()
198 int32_t ix; in __ieee754_y0l() local
[all …]
De_lgammal_r.c211 int n, ix; in sin_pi() local
215 ix = se & 0x7fff; in sin_pi()
216 ix = (ix << 16) | (i0 >> 16); in sin_pi()
217 if (ix < 0x3ffd8000) /* 0.25 */ in sin_pi()
234 if (ix >= 0x403f8000) /* 2^64 */ in sin_pi()
240 if (ix < 0x403e8000) /* 2^63 */ in sin_pi()
278 int i, ix; in __ieee754_lgammal_r() local
283 ix = se & 0x7fff; in __ieee754_lgammal_r()
285 if (__builtin_expect((ix | i0 | i1) == 0, 0)) in __ieee754_lgammal_r()
292 ix = (ix << 16) | (i0 >> 16); in __ieee754_lgammal_r()
[all …]
De_j1l.c116 int32_t ix; in __ieee754_j1l() local
120 ix = se & 0x7fff; in __ieee754_j1l()
121 if (__glibc_unlikely (ix >= 0x7fff)) in __ieee754_j1l()
124 if (ix >= 0x4000) in __ieee754_j1l()
129 if (ix < 0x7ffe) in __ieee754_j1l()
141 if (__glibc_unlikely (ix > 0x408e)) in __ieee754_j1l()
154 if (__glibc_unlikely (ix < 0x3fde)) /* |x| < 2^-33 */ in __ieee754_j1l()
199 int32_t ix; in __ieee754_y1l() local
203 ix = se & 0x7fff; in __ieee754_y1l()
207 if (__glibc_unlikely (ix >= 0x7fff)) in __ieee754_y1l()
[all …]
Ds_erfl.c258 int32_t ix, i; in __erfl() local
262 ix = se & 0x7fff; in __erfl()
264 if (ix >= 0x7fff) in __erfl()
270 ix = (ix << 16) | (i0 >> 16); in __erfl()
271 if (ix < 0x3ffed800) /* |x|<0.84375 */ in __erfl()
273 if (ix < 0x3fde8000) /* |x|<2**-33 */ in __erfl()
275 if (ix < 0x00080000) in __erfl()
292 if (ix < 0x3fffa000) /* 1.25 */ in __erfl()
304 if (ix >= 0x4001d555) /* 6.6666259765625 */ in __erfl()
313 if (ix < 0x4000b6db) /* 2.85711669921875 */ in __erfl()
[all …]
De_sinhl.c49 uint32_t jx,ix,i0,i1; in __ieee754_sinhl() local
53 ix = jx&0x7fff; in __ieee754_sinhl()
56 if(__builtin_expect(ix==0x7fff, 0)) return x+x; in __ieee754_sinhl()
61 if (ix < 0x4003 || (ix == 0x4003 && i0 <= 0xc8000000)) { /* |x|<25 */ in __ieee754_sinhl()
62 if (ix<0x3fdf) { /* |x|<2**-32 */ in __ieee754_sinhl()
67 if(ix<0x3fff) return h*(2.0*t-t*t/(t+one)); in __ieee754_sinhl()
72 if (ix < 0x400c || (ix == 0x400c && i0 < 0xb17217f7)) in __ieee754_sinhl()
76 if (ix<0x400c || (ix == 0x400c && (i0 < 0xb174ddc0 in __ieee754_sinhl()
/glibc-2.36/sysdeps/ieee754/ldbl-128ibm/
Ds_frexpl.c35 uint64_t hx, lx, ix, ixl; in __frexpl() local
43 ix = 0x7fffffffffffffffULL & hx; in __frexpl()
45 if (ix >= 0x7ff0000000000000ULL || ix == 0) in __frexpl()
51 expon = ix >> 52; in __frexpl()
58 if (sizeof (ix) == sizeof (long)) in __frexpl()
59 cnt = __builtin_clzl (ix); in __frexpl()
60 else if ((ix >> 32) != 0) in __frexpl()
61 cnt = __builtin_clzl ((long) (ix >> 32)); in __frexpl()
63 cnt = __builtin_clzl ((long) ix) + 32; in __frexpl()
66 ix <<= cnt + 1; in __frexpl()
[all …]
De_powl.c156 uint32_t ix, iy; in __ieee754_powl() local
163 ix = hx & 0x7fffffff; in __ieee754_powl()
180 if ((ix >= 0x7ff00000 && ((ix - 0x7ff00000) | lx) != 0) in __ieee754_powl()
246 if (ix == 0x7ff00000 || ix == 0 || (ix == 0x3ff00000 && xlo == 0.0)) in __ieee754_powl()
253 if (((ix - 0x3ff00000) | yisint) == 0) in __ieee754_powl()
281 if (ix <= 0x3fefffff) in __ieee754_powl()
283 if (ix >= 0x3ff00000) in __ieee754_powl()
287 if (ix < 0x3fefffff) in __ieee754_powl()
289 if (ix > 0x3ff00000) in __ieee754_powl()
299 if (ix < 0x00100000) in __ieee754_powl()
[all …]

1234567