/sysdeps/ieee754/flt-32/ |
A D | e_sqrtf.c | 40 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() 58 ix += ix; in __ieee754_sqrtf() 64 if(t<=ix) { in __ieee754_sqrtf() 66 ix -= t; in __ieee754_sqrtf() 69 ix += ix; in __ieee754_sqrtf() 74 if(ix!=0) { in __ieee754_sqrtf() [all …]
|
A D | s_setpayloadf_main.c | 33 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 …]
|
A D | s_getpayloadf.c | 28 uint32_t ix; in __getpayloadf() local 29 GET_FLOAT_WORD (ix, *x); in __getpayloadf() 30 if ((ix & 0x7f800000) != 0x7f800000 in __getpayloadf() 31 || (ix & 0x7fffff) == 0) in __getpayloadf() 33 ix &= 0x3fffff; in __getpayloadf() 34 if (FIX_INT_FP_CONVERT_ZERO && ix == 0) in __getpayloadf() 36 return (float) ix; in __getpayloadf()
|
A D | s_fromfpf_main.c | 43 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 …]
|
A D | s_logbf.c | 23 int32_t ix, rix; in __logbf() local 25 GET_FLOAT_WORD (ix, x); in __logbf() 26 ix &= 0x7fffffff; /* high |x| */ in __logbf() 27 if (ix == 0) in __logbf() 29 if (ix >= 0x7f800000) in __logbf() 31 if (__glibc_unlikely ((rix = ix >> 23) == 0)) in __logbf() 35 rix -= __builtin_clz (ix) - 9; in __logbf()
|
A D | s_roundevenf.c | 37 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()
|
A D | e_coshf.c | 27 int32_t ix; in __ieee754_coshf() local 29 GET_FLOAT_WORD(ix,x); in __ieee754_coshf() 30 ix &= 0x7fffffff; in __ieee754_coshf() 33 if (ix < 0x41b00000) { in __ieee754_coshf() 35 if(ix<0x3eb17218) { in __ieee754_coshf() 36 if (ix<0x24000000) return one; /* cosh(tiny) = 1 */ in __ieee754_coshf() 48 if (ix < 0x42b17180) return half*__ieee754_expf(fabsf(x)); in __ieee754_coshf() 51 if (ix<=0x42b2d4fc) { in __ieee754_coshf() 58 if(ix>=0x7f800000) return x*x; in __ieee754_coshf()
|
A D | s_isnanf.c | 30 int32_t ix; in __isnanf() local 31 GET_FLOAT_WORD(ix,x); in __isnanf() 32 ix &= 0x7fffffff; in __isnanf() 33 ix = 0x7f800000 - ix; in __isnanf() 34 return (int)(((uint32_t)(ix))>>31); in __isnanf()
|
A D | s_frexpf.c | 28 int32_t hx,ix; in __frexpf() local 30 ix = 0x7fffffff&hx; in __frexpf() 32 if(ix>=0x7f800000||(ix==0)) return x + x; /* 0,inf,nan */ in __frexpf() 33 if (ix<0x00800000) { /* subnormal */ in __frexpf() 36 ix = hx&0x7fffffff; in __frexpf() 39 *eptr += (ix>>23)-126; in __frexpf()
|
A D | e_log2f.c | 43 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 …]
|
A D | e_logf.c | 44 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 …]
|
A D | s_scalbnf.c | 27 int32_t k,ix; in __scalbnf() local 28 GET_FLOAT_WORD(ix,x); in __scalbnf() 29 k = (ix&0x7f800000)>>23; /* extract exponent */ in __scalbnf() 31 if ((ix&0x7fffffff)==0) return x; /* +-0 */ in __scalbnf() 33 GET_FLOAT_WORD(ix,x); in __scalbnf() 34 k = ((ix&0x7f800000)>>23) - 25; in __scalbnf() 45 {SET_FLOAT_WORD(x,(ix&0x807fffff)|(k<<23)); return x;} in __scalbnf() 49 SET_FLOAT_WORD(x,(ix&0x807fffff)|(k<<23)); in __scalbnf()
|
A D | s_scalblnf.c | 27 int32_t k,ix; in __scalblnf() local 28 GET_FLOAT_WORD(ix,x); in __scalblnf() 29 k = (ix&0x7f800000)>>23; /* extract exponent */ in __scalblnf() 31 if ((ix&0x7fffffff)==0) return x; /* +-0 */ in __scalblnf() 33 GET_FLOAT_WORD(ix,x); in __scalblnf() 34 k = ((ix&0x7f800000)>>23) - 25; in __scalblnf() 45 {SET_FLOAT_WORD(x,(ix&0x807fffff)|(k<<23)); return x;} in __scalblnf() 49 SET_FLOAT_WORD(x,(ix&0x807fffff)|(k<<23)); in __scalblnf()
|
/sysdeps/ieee754/dbl-64/ |
A D | s_setpayload_main.c | 33 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 …]
|
A D | s_getpayload.c | 28 uint64_t ix; in __getpayload() local 29 EXTRACT_WORDS64 (ix, *x); in __getpayload() 30 if ((ix & 0x7ff0000000000000ULL) != 0x7ff0000000000000ULL in __getpayload() 31 || (ix & 0xfffffffffffffULL) == 0) in __getpayload() 33 ix &= 0x7ffffffffffffULL; in __getpayload() 34 if (FIX_INT_FP_CONVERT_ZERO && ix == 0) in __getpayload() 36 return (double) ix; in __getpayload()
|
A D | s_fromfp_main.c | 43 uint64_t ix; in FUNC() local 44 EXTRACT_WORDS64 (ix, x); 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() 72 more_bits = (ix & (h - 1)) != 0; in FUNC() [all …]
|
A D | s_roundeven.c | 36 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()
|
A D | s_frexp.c | 37 int64_t ix; in __frexp() local 38 EXTRACT_WORDS64 (ix, x); in __frexp() 39 int32_t ex = 0x7ff & (ix >> 52); in __frexp() 50 EXTRACT_WORDS64 (ix, x); in __frexp() 51 ex = 0x7ff & (ix >> 52); in __frexp() 55 ix = (ix & INT64_C (0x800fffffffffffff)) | INT64_C (0x3fe0000000000000); in __frexp() 56 INSERT_WORDS64 (x, ix); in __frexp()
|
A D | s_logb.c | 28 int64_t ix, ex; in __logb() local 30 EXTRACT_WORDS64 (ix, x); in __logb() 31 ix &= UINT64_C(0x7fffffffffffffff); in __logb() 32 if (ix == 0) in __logb() 34 ex = ix >> 52; in __logb() 39 int m = __builtin_clzll (ix); in __logb()
|
/sysdeps/ieee754/ldbl-128ibm/ |
A D | s_getpayloadl.c | 28 uint64_t ix; in __getpayloadl() local 29 EXTRACT_WORDS64 (ix, xhi); in __getpayloadl() 30 if ((ix & 0x7ff0000000000000ULL) != 0x7ff0000000000000ULL in __getpayloadl() 31 || (ix & 0xfffffffffffffULL) == 0) in __getpayloadl() 33 ix &= 0x7ffffffffffffULL; in __getpayloadl() 34 if (FIX_INT_FP_CONVERT_ZERO && ix == 0) in __getpayloadl() 36 return (long double) ix; in __getpayloadl()
|
A D | s_frexpl.c | 35 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() 66 ix <<= cnt + 1; in __frexpl() 69 ix &= 0x000fffffffffffffULL; in __frexpl() 71 hx |= (1022LL << 52) | ix; in __frexpl() [all …]
|
A D | e_coshl.c | 44 int64_t ix; in __ieee754_coshl() local 49 EXTRACT_WORDS64 (ix, xhi); in __ieee754_coshl() 50 ix &= 0x7fffffffffffffffLL; in __ieee754_coshl() 53 if(ix>=0x7ff0000000000000LL) return x*x; in __ieee754_coshl() 56 if(ix<0x3fd62e42fefa39efLL) { in __ieee754_coshl() 57 if (ix<0x3c80000000000000LL) return one; /* cosh(tiny) = 1 */ in __ieee754_coshl() 64 if (ix < 0x4044000000000000LL) { in __ieee754_coshl() 70 if (ix < 0x40862e42fefa39efLL) return half*__ieee754_expl(fabsl(x)); in __ieee754_coshl() 73 if (ix < 0x408633ce8fb9f87fLL) { in __ieee754_coshl()
|
/sysdeps/ieee754/ldbl-96/ |
A D | e_sinhl.c | 49 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()
|
A D | s_frexpl.c | 44 uint32_t se, hx, ix, lx; in __frexpl() local 46 ix = 0x7fff&se; in __frexpl() 48 if(ix==0x7fff||((ix|hx|lx)==0)) return x + x; /* 0,inf,nan */ in __frexpl() 49 if (ix==0x0000) { /* subnormal */ in __frexpl() 52 ix = se&0x7fff; in __frexpl() 55 *eptr += ix-16382; in __frexpl()
|
/sysdeps/ieee754/ldbl-128/ |
A D | s_frexpl.c | 38 uint64_t hx, lx, ix; in __frexpl() local 40 ix = 0x7fffffffffffffffULL&hx; in __frexpl() 42 if(ix>=0x7fff000000000000ULL||((ix|lx)==0)) return x + x;/* 0,inf,nan */ in __frexpl() 43 if (ix<0x0001000000000000ULL) { /* subnormal */ in __frexpl() 46 ix = hx&0x7fffffffffffffffULL; in __frexpl() 49 *eptr += (ix>>48)-16382; in __frexpl()
|