/sysdeps/ieee754/flt-32/ |
A D | e_fmodf.c | 35 hx ^=sx; /* |x| */ in __ieee754_fmodf() 43 if(hx==hy) in __ieee754_fmodf() 58 hx = 0x00800000|(0x007fffff&hx); in __ieee754_fmodf() 61 hx = hx<<n; in __ieee754_fmodf() 73 hz=hx-hy; in __ieee754_fmodf() 74 if(hz<0){hx = hx+hx;} in __ieee754_fmodf() 78 hx = hz+hz; in __ieee754_fmodf() 81 hz=hx-hy; in __ieee754_fmodf() 88 hx = hx+hx; in __ieee754_fmodf() 92 hx = ((hx-0x00800000)|((iy+127)<<23)); in __ieee754_fmodf() [all …]
|
A D | e_ilogbf.c | 25 int32_t hx,ix; in __ieee754_ilogbf() local 27 GET_FLOAT_WORD(hx,x); in __ieee754_ilogbf() 28 hx &= 0x7fffffff; in __ieee754_ilogbf() 29 if(hx<0x00800000) { in __ieee754_ilogbf() 30 if(hx==0) in __ieee754_ilogbf() 33 for (ix = -126,hx<<=8; hx>0; hx<<=1) ix -=1; in __ieee754_ilogbf() 36 else if (hx<0x7f800000) return (hx>>23)-127; in __ieee754_ilogbf() 39 if (hx==0x7f800000) in __ieee754_ilogbf()
|
A D | e_log10f.c | 30 int32_t i,k,hx; in __ieee754_log10f() local 32 GET_FLOAT_WORD(hx,x); in __ieee754_log10f() 35 if (hx < 0x00800000) { /* x < 2**-126 */ in __ieee754_log10f() 36 if (__builtin_expect((hx&0x7fffffff)==0, 0)) in __ieee754_log10f() 38 if (__builtin_expect(hx<0, 0)) in __ieee754_log10f() 41 GET_FLOAT_WORD(hx,x); in __ieee754_log10f() 43 if (__builtin_expect(hx >= 0x7f800000, 0)) return x+x; in __ieee754_log10f() 44 k += (hx>>23)-127; in __ieee754_log10f() 46 hx = (hx&0x007fffff)|((0x7f-i)<<23); in __ieee754_log10f() 50 SET_FLOAT_WORD(x,hx); in __ieee754_log10f()
|
A D | s_nextafterf.c | 28 int32_t hx,hy,ix,iy; in __nextafterf() local 30 GET_FLOAT_WORD(hx,x); in __nextafterf() 32 ix = hx&0x7fffffff; /* |x| */ in __nextafterf() 47 if(hx>=0) { /* x > 0 */ in __nextafterf() 48 if(hx>hy) { /* x > y, x -= ulp */ in __nextafterf() 49 hx -= 1; in __nextafterf() 51 hx += 1; in __nextafterf() 55 hx -= 1; in __nextafterf() 57 hx += 1; in __nextafterf() 60 hy = hx&0x7f800000; in __nextafterf() [all …]
|
A D | s_frexpf.c | 28 int32_t hx,ix; in __frexpf() local 29 GET_FLOAT_WORD(hx,x); in __frexpf() 30 ix = 0x7fffffff&hx; in __frexpf() 35 GET_FLOAT_WORD(hx,x); in __frexpf() 36 ix = hx&0x7fffffff; in __frexpf() 40 hx = (hx&0x807fffff)|0x3f000000; in __frexpf() 41 SET_FLOAT_WORD(x,hx); in __frexpf()
|
A D | s_nextupf.c | 28 int32_t hx, ix; in __nextupf() local 30 GET_FLOAT_WORD (hx, x); in __nextupf() 31 ix = hx & 0x7fffffff; in __nextupf() 36 if (hx >= 0) in __nextupf() 40 hx += 1; in __nextupf() 43 hx -= 1; in __nextupf() 44 SET_FLOAT_WORD (x, hx); in __nextupf()
|
/sysdeps/ieee754/ldbl-96/ |
A D | s_nextupl.c | 28 uint32_t hx, ix; in __nextupl() local 37 if ((ix | hx | lx) == 0) in __nextupl() 46 hx += 1; in __nextupl() 48 if (hx == 0 || (esx == 0 && hx == 0x80000000)) in __nextupl() 50 if (hx == 0) in __nextupl() 54 hx |= 0x80000000; in __nextupl() 66 hx = hx - 1; in __nextupl() 71 hx -= 1; in __nextupl() 74 hx = 0; in __nextupl() 75 if (hx == 0) in __nextupl() [all …]
|
A D | s_roundevenl.c | 34 uint32_t hx, lx; in __roundevenl() local 60 hx++; in __roundevenl() 61 if (hx == 0) in __roundevenl() 63 hx = 0x80000000; in __roundevenl() 79 hx++; in __roundevenl() 80 if (hx == 0) in __roundevenl() 82 hx = 0x80000000; in __roundevenl() 100 hx += half_bit; in __roundevenl() 101 if (hx < half_bit) in __roundevenl() 114 hx = 0x80000000; in __roundevenl() [all …]
|
A D | s_setpayloadl_main.c | 33 uint32_t hx, lx; in FUNC() local 35 GET_LDOUBLE_WORDS (exponent, hx, lx, payload); in FUNC() 40 && exponent == 0 && hx == 0 && lx == 0))) in FUNC() 48 : (lx != 0 || (hx & ((1U << (shift - 32)) - 1)) != 0)) in FUNC() 57 lx = hx >> (shift - 32); in FUNC() 58 hx = 0; in FUNC() 62 lx = (lx >> shift) | (hx << (32 - shift)); in FUNC() 63 hx >>= shift; in FUNC() 66 hx |= 0x80000000 | (SET_HIGH_BIT ? 0x40000000 : 0); in FUNC() 67 SET_LDOUBLE_WORDS (*x, 0x7fff, hx, lx); in FUNC()
|
/sysdeps/ieee754/dbl-64/ |
A D | e_fmod.c | 29 int64_t hx,hy,hz,sx,i; in __ieee754_fmod() local 34 hx ^=sx; /* |x| */ in __ieee754_fmod() 64 hx<<=n; in __ieee754_fmod() 76 hz=hx-hy; in __ieee754_fmod() 77 if(hz<0){hx = hx+hx;} in __ieee754_fmod() 81 hx = hz+hz; in __ieee754_fmod() 84 hz=hx-hy; in __ieee754_fmod() 85 if(hz>=0) {hx=hz;} in __ieee754_fmod() 91 hx = hx+hx; in __ieee754_fmod() 95 hx = ((hx-UINT64_C(0x0010000000000000))|((uint64_t)(iy+1023)<<52)); in __ieee754_fmod() [all …]
|
A D | e_ilogb.c | 31 int32_t hx, lx, ix; in __ieee754_ilogb() local 33 GET_HIGH_WORD (hx, x); in __ieee754_ilogb() 34 hx &= 0x7fffffff; in __ieee754_ilogb() 35 if (hx < 0x00100000) in __ieee754_ilogb() 38 if ((hx | lx) == 0) in __ieee754_ilogb() 41 if (hx == 0) in __ieee754_ilogb() 48 for (ix = -1022, hx <<= 11; hx > 0; hx <<= 1) in __ieee754_ilogb() 53 else if (hx < 0x7ff00000) in __ieee754_ilogb() 54 return (hx >> 20) - 1023; in __ieee754_ilogb() 59 if (((hx ^ 0x7ff00000) | lx) == 0) in __ieee754_ilogb()
|
A D | e_log10.c | 61 int64_t i, hx; in __ieee754_log10() local 64 EXTRACT_WORDS64 (hx, x); in __ieee754_log10() 67 if (hx < INT64_C(0x0010000000000000)) in __ieee754_log10() 69 if (__glibc_unlikely ((hx & UINT64_C(0x7fffffffffffffff)) == 0)) in __ieee754_log10() 71 if (__glibc_unlikely (hx < 0)) in __ieee754_log10() 75 EXTRACT_WORDS64 (hx, x); in __ieee754_log10() 78 if (__glibc_unlikely (hx >= UINT64_C(0x7ff0000000000000))) in __ieee754_log10() 80 k += (hx >> 52) - 1023; in __ieee754_log10() 82 hx = (hx & UINT64_C(0x000fffffffffffff)) | ((0x3ff - i) << 52); in __ieee754_log10() 86 INSERT_WORDS64 (x, hx); in __ieee754_log10()
|
/sysdeps/ieee754/ldbl-128/ |
A D | e_fmodl.c | 42 if(hx<=hy) { in __ieee754_fmodl() 50 if(hx==0) { in __ieee754_fmodl() 68 hx = 0x0001000000000000LL|(0x0000ffffffffffffLL&hx); in __ieee754_fmodl() 72 hx = (hx<<n)|(lx>>(64-n)); in __ieee754_fmodl() 75 hx = lx<<(n-64); in __ieee754_fmodl() 96 if(hz<0){hx = hx+hx+(lx>>63); lx = lx+lx;} in __ieee754_fmodl() 110 hx = hx+hx+(lx>>63); lx = lx+lx; in __ieee754_fmodl() 114 hx = ((hx-0x0001000000000000LL)|((iy+16383)<<48)); in __ieee754_fmodl() 120 hx >>= n; in __ieee754_fmodl() 122 lx = (hx<<(64-n))|(lx>>n); hx = sx; in __ieee754_fmodl() [all …]
|
A D | e_ilogbl.c | 32 int64_t hx,lx; in __ieee754_ilogbl() local 35 GET_LDOUBLE_WORDS64(hx,lx,x); in __ieee754_ilogbl() 36 hx &= 0x7fffffffffffffffLL; in __ieee754_ilogbl() 37 if(hx <= 0x0001000000000000LL) { in __ieee754_ilogbl() 38 if((hx|lx)==0) in __ieee754_ilogbl() 41 if(hx==0) { in __ieee754_ilogbl() 44 for (ix = -16382, hx<<=15; hx>0; hx<<=1) ix -=1; in __ieee754_ilogbl() 48 else if (hx<0x7fff000000000000LL) return (hx>>48)-0x3fff; in __ieee754_ilogbl() 51 if (((hx^0x7fff000000000000LL)|lx) == 0) in __ieee754_ilogbl()
|
A D | s_getpayloadl.c | 27 uint64_t hx, lx; in __getpayloadl() local 28 GET_LDOUBLE_WORDS64 (hx, lx, *x); in __getpayloadl() 29 if ((hx & 0x7fff000000000000ULL) != 0x7fff000000000000ULL in __getpayloadl() 30 || ((hx & 0xffffffffffffULL) | lx) == 0) in __getpayloadl() 32 hx &= 0x7fffffffffffULL; in __getpayloadl() 36 if (hx == 0) in __getpayloadl() 44 lz = __builtin_clzll (hx); in __getpayloadl() 48 hx = lx << (shift - 64); in __getpayloadl() 54 hx = (hx << shift) | (lx >> (64 - shift)); in __getpayloadl() 57 hx = (hx & 0xffffffffffffULL) | ((0x3fffULL + 127 - lz) << 48); in __getpayloadl() [all …]
|
A D | s_setpayloadl_main.c | 33 uint64_t hx, lx; in FUNC() local 34 GET_LDOUBLE_WORDS64 (hx, lx, payload); in FUNC() 35 int exponent = hx >> (EXPLICIT_MANT_DIG - 64); in FUNC() 47 : (lx != 0 || (hx & ((1ULL << (shift - 64)) - 1)) != 0)) in FUNC() 54 hx &= (1ULL << (EXPLICIT_MANT_DIG - 64)) - 1; in FUNC() 55 hx |= 1ULL << (EXPLICIT_MANT_DIG - 64); in FUNC() 58 lx = hx >> (shift - 64); in FUNC() 59 hx = 0; in FUNC() 63 lx = (lx >> shift) | (hx << (64 - shift)); in FUNC() 64 hx >>= shift; in FUNC() [all …]
|
A D | s_nextafterl.c | 34 int64_t hx,hy,ix,iy; in __nextafterl() local 37 GET_LDOUBLE_WORDS64(hx,lx,x); in __nextafterl() 54 if(hx>=0) { /* x > 0 */ in __nextafterl() 55 if(hx>hy||((hx==hy)&&(lx>ly))) { /* x > y, x -= ulp */ in __nextafterl() 56 if(lx==0) hx--; in __nextafterl() 60 if(lx==0) hx++; in __nextafterl() 63 if(hy>=0||hx>hy||((hx==hy)&&(lx>ly))){/* x < y, x -= ulp */ in __nextafterl() 64 if(lx==0) hx--; in __nextafterl() 68 if(lx==0) hx++; in __nextafterl() 71 hy = hx&0x7fff000000000000LL; in __nextafterl() [all …]
|
A D | s_isnanl.c | 29 int64_t hx,lx; in __isnanl() local 30 GET_LDOUBLE_WORDS64(hx,lx,x); in __isnanl() 31 hx &= 0x7fffffffffffffffLL; in __isnanl() 32 hx |= (uint64_t)(lx|(-lx))>>63; in __isnanl() 33 hx = 0x7fff000000000000LL - hx; in __isnanl() 34 return (int)((uint64_t)hx>>63); in __isnanl()
|
A D | s_logbl.c | 32 int64_t lx, hx, ex; in __logbl() local 34 GET_LDOUBLE_WORDS64 (hx, lx, x); in __logbl() 35 hx &= 0x7fffffffffffffffLL; /* high |x| */ in __logbl() 36 if ((hx | lx) == 0) in __logbl() 38 if (hx >= 0x7fff000000000000LL) in __logbl() 40 if ((ex = hx >> 48) == 0) /* IEEE 754 logb */ in __logbl() 45 if (hx == 0) in __logbl() 48 ma = __builtin_clzll (hx); in __logbl()
|
/sysdeps/i386/fpu/ |
A D | s_nextafterl.c | 59 if(esx>esy||((esx==esy) && (hx>hy||((hx==hy)&&(lx>ly))))) { in __nextafterl() 64 --hx; in __nextafterl() 67 hx = hx - 1; in __nextafterl() 72 hx -= 1; in __nextafterl() 78 hx += 1; in __nextafterl() 79 if (hx==0 || (esx == 0 && hx == 0x80000000)) { in __nextafterl() 81 hx |= 0x80000000; in __nextafterl() 91 hx = hx - 1; in __nextafterl() 95 hx -= 1; in __nextafterl() 101 hx += 1; in __nextafterl() [all …]
|
/sysdeps/m68k/m680x0/fpu/ |
A D | s_nextafterl.c | 52 if(esx>esy||((esx==esy) && (hx>hy||((hx==hy)&&(lx>ly))))) { in __nextafterl() 55 if (ix != 0 && hx == 0x80000000) hx = 0; in __nextafterl() 57 hx -= 1; in __nextafterl() 63 hx += 1; in __nextafterl() 64 if (hx==0) { in __nextafterl() 65 hx = 0x80000000; in __nextafterl() 74 if (ix != 0 && hx == 0x80000000) hx = 0; in __nextafterl() 76 hx -= 1; in __nextafterl() 82 hx += 1; in __nextafterl() 83 if (hx==0) { in __nextafterl() [all …]
|
/sysdeps/ieee754/ldbl-128ibm/ |
A D | e_ilogbl.c | 33 int64_t hx, hxs; in __ieee754_ilogbl() local 38 EXTRACT_WORDS64 (hx, xhi); in __ieee754_ilogbl() 39 hxs = hx; in __ieee754_ilogbl() 40 hx &= 0x7fffffffffffffffLL; in __ieee754_ilogbl() 41 if(hx <= 0x0010000000000000LL) { in __ieee754_ilogbl() 42 if(hx==0) in __ieee754_ilogbl() 45 for (ix = -1022, hx<<=11; hx>0; hx<<=1) ix -=1; in __ieee754_ilogbl() 48 else if (hx < 0x7ff0000000000000LL) in __ieee754_ilogbl() 50 int hexp = (hx >> 52) - 0x3ff; in __ieee754_ilogbl() 54 if ((hx & 0x000fffffffffffffLL) == 0) in __ieee754_ilogbl() [all …]
|
A D | e_fmodl.c | 30 int64_t hx, hy, hz, sx, sy; in __ieee754_fmodl() local 36 EXTRACT_WORDS64 (hx, xhi); in __ieee754_fmodl() 42 hx ^= sx; /* |x| */ in __ieee754_fmodl() 51 if (__glibc_unlikely (hx <= hy)) in __ieee754_fmodl() 54 if (hx < hy) in __ieee754_fmodl() 98 hx = (hx << 1) | (lx >> 63); in __ieee754_fmodl() 119 if(hz<0){hx = hx+hx+(lx>>63); lx = lx+lx;} in __ieee754_fmodl() 123 hx = hz+hz+(lz>>63); lx = lz+lz; in __ieee754_fmodl() 127 if(hz>=0) {hx=hz;lx=lz;} in __ieee754_fmodl() 133 hx = hx+hx+(lx>>63); lx = lx+lx; in __ieee754_fmodl() [all …]
|
A D | s_setpayloadl_main.c | 33 uint64_t hx, lx; in FUNC() local 36 EXTRACT_WORDS64 (hx, hi); in FUNC() 38 int exponent = hx >> EXPLICIT_MANT_DIG; in FUNC() 44 || (exponent < BIAS && !(SET_HIGH_BIT && hx == 0)) in FUNC() 45 || (hx & ((1ULL << (BIAS + EXPLICIT_MANT_DIG - exponent)) - 1)) != 0) in FUNC() 50 if (hx != 0) in FUNC() 52 hx &= (1ULL << EXPLICIT_MANT_DIG) - 1; in FUNC() 53 hx |= 1ULL << EXPLICIT_MANT_DIG; in FUNC() 54 hx >>= BIAS + EXPLICIT_MANT_DIG - exponent; in FUNC() 56 hx |= 0x7ff0000000000000ULL | (SET_HIGH_BIT ? 0x8000000000000ULL : 0); in FUNC() [all …]
|
A D | s_logbl.c | 29 int64_t hx, hxs, rhx; in __logbl() local 33 EXTRACT_WORDS64 (hx, xhi); in __logbl() 34 hxs = hx; in __logbl() 35 hx &= 0x7fffffffffffffffLL; /* high |x| */ in __logbl() 36 if (hx == 0) in __logbl() 38 if (hx >= 0x7ff0000000000000LL) in __logbl() 40 if (__glibc_unlikely ((rhx = hx >> 52) == 0)) in __logbl() 44 rhx -= __builtin_clzll (hx) - 12; in __logbl() 46 else if ((hx & 0x000fffffffffffffLL) == 0) in __logbl()
|