/glibc-2.36/sysdeps/ieee754/ldbl-128/ |
D | e_fmodl.c | 29 int64_t n,hx,hy,hz,ix,iy,sx,i; in __ieee754_fmodl() local 32 GET_LDOUBLE_WORDS64(hx,lx,x); in __ieee754_fmodl() 34 sx = hx&0x8000000000000000ULL; /* sign of x */ in __ieee754_fmodl() 35 hx ^=sx; /* |x| */ in __ieee754_fmodl() 39 if((hy|ly)==0||(hx>=0x7fff000000000000LL)|| /* y=0,or x not finite */ in __ieee754_fmodl() 42 if(hx<=hy) { in __ieee754_fmodl() 43 if((hx<hy)||(lx<ly)) return x; /* |x|<|y| return x */ in __ieee754_fmodl() 49 if(hx<0x0001000000000000LL) { /* subnormal x */ in __ieee754_fmodl() 50 if(hx==0) { in __ieee754_fmodl() 53 for (ix = -16382, i=hx<<15; i>0; i<<=1) ix -=1; in __ieee754_fmodl() [all …]
|
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 …]
|
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() 39 || (exponent < BIAS && !(SET_HIGH_BIT && hx == 0 && lx == 0))) 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() [all …]
|
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()
|
D | s_nextafterl.c | 34 int64_t hx,hy,ix,iy; in __nextafterl() local 37 GET_LDOUBLE_WORDS64(hx,lx,x); in __nextafterl() 39 ix = hx&0x7fffffffffffffffLL; /* |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() [all …]
|
D | s_roundevenl.c | 37 uint64_t hx, lx, uhx; in __roundevenl() 38 GET_LDOUBLE_WORDS64 (hx, lx, x); in __roundevenl() 39 uhx = hx & 0x7fffffffffffffffULL; in __roundevenl() 63 hx += lx < half_bit; in __roundevenl() 71 if (((hx & 1) | (lx & 0x7fffffffffffffffULL)) != 0) in __roundevenl() 74 hx += lx < 0x8000000000000000ULL; in __roundevenl() 89 if (((hx & (int_bit | (half_bit - 1))) | lx) != 0) in __roundevenl() 90 hx += half_bit; in __roundevenl() 91 hx &= ~(int_bit - 1); in __roundevenl() 97 hx = (hx & 0x8000000000000000ULL) | 0x3fff000000000000ULL; in __roundevenl() [all …]
|
D | s_fromfpl_main.c | 43 uint64_t hx, lx; in FUNC() local 44 GET_LDOUBLE_WORDS64 (hx, lx, x); in FUNC() 45 bool negative = (hx & 0x8000000000000000ULL) != 0; in FUNC() 48 hx &= 0x7fffffffffffffffULL; in FUNC() 49 if ((hx | lx) == 0) in FUNC() 51 int exponent = hx >> (MANT_DIG - 1 - 64); in FUNC() 57 hx &= ((1ULL << (MANT_DIG - 1 - 64)) - 1); in FUNC() 58 hx |= 1ULL << (MANT_DIG - 1 - 64); in FUNC() 71 uret = hx << (64 - shift); in FUNC() 78 half_bit = (hx & h) != 0; in FUNC() [all …]
|
/glibc-2.36/sysdeps/ieee754/flt-32/ |
D | e_fmodf.c | 30 int32_t n,hx,hy,hz,ix,iy,sx,i; in __ieee754_fmodf() local 32 GET_FLOAT_WORD(hx,x); in __ieee754_fmodf() 34 sx = hx&0x80000000; /* sign of x */ in __ieee754_fmodf() 35 hx ^=sx; /* |x| */ in __ieee754_fmodf() 39 if(hy==0||(hx>=0x7f800000)|| /* y=0,or x not finite */ in __ieee754_fmodf() 42 if(hx<hy) return x; /* |x|<|y| return x */ in __ieee754_fmodf() 43 if(hx==hy) in __ieee754_fmodf() 47 if(hx<0x00800000) { /* subnormal x */ in __ieee754_fmodf() 48 for (ix = -126,i=(hx<<8); i>0; i<<=1) ix -=1; in __ieee754_fmodf() 49 } else ix = (hx>>23)-127; in __ieee754_fmodf() [all …]
|
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()
|
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() 54 if(hy>=0||hx>hy){ /* x < y, x -= ulp */ in __nextafterf() 55 hx -= 1; in __nextafterf() 57 hx += 1; in __nextafterf() [all …]
|
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()
|
/glibc-2.36/sysdeps/ieee754/dbl-64/ |
D | e_fmod.c | 29 int64_t hx,hy,hz,sx,i; in __ieee754_fmod() local 31 EXTRACT_WORDS64(hx,x); in __ieee754_fmod() 33 sx = hx&UINT64_C(0x8000000000000000); /* sign of x */ in __ieee754_fmod() 34 hx ^=sx; /* |x| */ in __ieee754_fmod() 39 || hx >= UINT64_C(0x7ff0000000000000) in __ieee754_fmod() 43 if(__builtin_expect(hx<=hy, 0)) { in __ieee754_fmod() 44 if(hx<hy) return x; /* |x|<|y| return x */ in __ieee754_fmod() 49 if(__builtin_expect(hx<UINT64_C(0x0010000000000000), 0)) { in __ieee754_fmod() 51 for (ix = -1022,i=(hx<<11); i>0; i<<=1) ix -=1; in __ieee754_fmod() 52 } else ix = (hx>>52)-1023; in __ieee754_fmod() [all …]
|
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()
|
/glibc-2.36/sysdeps/i386/fpu/ |
D | s_nextafterl.c | 35 uint32_t hx,hy,ix,iy; in __nextafterl() local 39 GET_LDOUBLE_WORDS(esx,hx,lx,x); in __nextafterl() 46 if(((ix==0x7fff)&&(((hx&0x7fffffff)|lx)!=0)) || /* x is nan */ in __nextafterl() 50 if((ix|hx|lx)==0) { /* x == 0 */ in __nextafterl() 59 if(esx>esy||((esx==esy) && (hx>hy||((hx==hy)&&(lx>ly))))) { in __nextafterl() 62 if (hx <= 0x80000000) { in __nextafterl() 64 --hx; in __nextafterl() 67 hx = hx - 1; in __nextafterl() 69 hx |= 0x80000000; in __nextafterl() 72 hx -= 1; in __nextafterl() [all …]
|
/glibc-2.36/sysdeps/m68k/m680x0/fpu/ |
D | s_nextafterl.c | 34 uint32_t hx,hy,lx,ly; in __nextafterl() local 36 GET_LDOUBLE_WORDS(esx,hx,lx,x); in __nextafterl() 41 if(((ix==0x7fff)&&((hx&0x7fffffff)|lx)!=0) || /* x is nan */ in __nextafterl() 45 if((ix|hx|lx)==0) { /* x == 0 */ in __nextafterl() 52 if(esx>esy||((esx==esy) && (hx>hy||((hx==hy)&&(lx>ly))))) { in __nextafterl() 55 if (ix != 0 && hx == 0x80000000) hx = 0; in __nextafterl() 56 if (hx==0) esx -= 1; in __nextafterl() 57 hx -= 1; in __nextafterl() 63 hx += 1; in __nextafterl() 64 if (hx==0) { in __nextafterl() [all …]
|
/glibc-2.36/sysdeps/ieee754/ldbl-128ibm/ |
D | e_fmodl.c | 30 int64_t hx, hy, hz, sx, sy; in __ieee754_fmodl() local 36 EXTRACT_WORDS64 (hx, xhi); in __ieee754_fmodl() 41 sx = hx&0x8000000000000000ULL; /* sign of x */ in __ieee754_fmodl() 42 hx ^= sx; /* |x| */ in __ieee754_fmodl() 48 (hx>=0x7ff0000000000000LL)|| /* y=0,or x not finite */ in __ieee754_fmodl() 51 if (__glibc_unlikely (hx <= hy)) in __ieee754_fmodl() 54 if (hx < hy) in __ieee754_fmodl() 90 ldbl_extract_mantissa(&hx, &lx, &ix, x); in __ieee754_fmodl() 96 while ((hx & (1LL << 48)) == 0) in __ieee754_fmodl() 98 hx = (hx << 1) | (lx >> 63); in __ieee754_fmodl() [all …]
|
D | s_nextafterl.c | 35 int64_t hx, hy, ihx, ihy, lx; in __nextafterl() local 39 EXTRACT_WORDS64 (hx, xhi); in __nextafterl() 43 ihx = hx&0x7fffffffffffffffLL; /* |hx| */ in __nextafterl() 70 if((hx==0xffefffffffffffffLL)&&(lx==0xfc8ffffffffffffeLL)) { in __nextafterl() 76 if (hx >= 0x7ff0000000000000LL) { in __nextafterl() 84 || (hx > 0 && lx <= 0) in __nextafterl() 85 || (hx < 0 && lx > 1)) { in __nextafterl() 101 if ((hx & 0x000fffffffffffffLL) == 0 in __nextafterl() 102 && ((lx != 0 && (hx ^ lx) < 0) in __nextafterl() 103 || (lx == 0 && hx >= 0))) in __nextafterl() [all …]
|
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 …]
|
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 …]
|
D | s_scalblnl.c | 41 int64_t k,l,hx,lx; in __scalblnl() local 46 EXTRACT_WORDS64 (hx, xhi); in __scalblnl() 48 k = (hx>>52)&0x7ff; /* extract exponent */ in __scalblnl() 51 if ((hx&0x7fffffffffffffffULL)==0) return x; /* +-0 */ in __scalblnl() 52 u.i = hx; in __scalblnl() 54 hx = u.i; in __scalblnl() 55 k = ((hx>>52)&0x7ff) - 54; in __scalblnl() 65 hx = (hx&0x800fffffffffffffULL)|(k<<52); in __scalblnl() 67 INSERT_WORDS64 (xhi, hx); in __scalblnl() 89 INSERT_WORDS64 (xhi, hx); in __scalblnl() [all …]
|
D | s_scalbnl.c | 41 int64_t k,l,hx,lx; in __scalbnl() local 46 EXTRACT_WORDS64 (hx, xhi); in __scalbnl() 48 k = (hx>>52)&0x7ff; /* extract exponent */ in __scalbnl() 51 if ((hx&0x7fffffffffffffffULL)==0) return x; /* +-0 */ in __scalbnl() 52 u.i = hx; in __scalbnl() 54 hx = u.i; in __scalbnl() 55 k = ((hx>>52)&0x7ff) - 54; in __scalbnl() 65 hx = (hx&0x800fffffffffffffULL)|(k<<52); in __scalbnl() 67 INSERT_WORDS64 (xhi, hx); in __scalbnl() 89 INSERT_WORDS64 (xhi, hx); in __scalbnl() [all …]
|
D | s_fromfpl_main.c | 43 uint64_t hx, lx; in FUNC() local 45 EXTRACT_WORDS64 (hx, hi); in FUNC() 47 bool negative = (hx & 0x8000000000000000ULL) != 0; in FUNC() 51 hx &= 0x7fffffffffffffffULL; in FUNC() 53 if ((hx | lx) == 0) in FUNC() 55 int hi_exponent = hx >> (MANT_DIG - 1); in FUNC() 58 hx &= ((1ULL << (MANT_DIG - 1)) - 1); in FUNC() 59 if (hx == 0 && lx != 0 && lo_negative != negative) in FUNC() 68 hx |= 1ULL << (MANT_DIG - 1); in FUNC() 73 uret = hx; in FUNC() [all …]
|
/glibc-2.36/sysdeps/ieee754/ldbl-96/ |
D | s_nextupl.c | 28 uint32_t hx, ix; in __nextupl() local 32 GET_LDOUBLE_WORDS (esx, hx, lx, x); in __nextupl() 35 if (((ix == 0x7fff) && (((hx & 0x7fffffff) | lx) != 0))) /* x is nan. */ in __nextupl() 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() 63 if (hx <= 0x80000000 && esx != 0xffff8000) in __nextupl() 66 hx = hx - 1; in __nextupl() [all …]
|
D | s_roundevenl.c | 34 uint32_t hx, lx; in __roundevenl() local 35 GET_LDOUBLE_WORDS (se, hx, lx, x); in __roundevenl() 60 hx++; in __roundevenl() 61 if (hx == 0) in __roundevenl() 63 hx = 0x80000000; in __roundevenl() 74 if (((hx & 1) | (lx & 0x7fffffff)) != 0) in __roundevenl() 79 hx++; in __roundevenl() 80 if (hx == 0) in __roundevenl() 82 hx = 0x80000000; in __roundevenl() 98 if (((hx & (int_bit | (half_bit - 1))) | lx) != 0) in __roundevenl() [all …]
|
/glibc-2.36/math/ |
D | s_nextafter.c | 37 int32_t hx,hy,ix,iy; in __nextafter() local 40 EXTRACT_WORDS(hx,lx,x); in __nextafter() 42 ix = hx&0x7fffffff; /* |x| */ in __nextafter() 57 if(hx>=0) { /* x > 0 */ in __nextafter() 58 if(hx>hy||((hx==hy)&&(lx>ly))) { /* x > y, x -= ulp */ in __nextafter() 59 if(lx==0) hx -= 1; in __nextafter() 63 if(lx==0) hx += 1; in __nextafter() 66 if(hy>=0||hx>hy||((hx==hy)&&(lx>ly))){/* x < y, x -= ulp */ in __nextafter() 67 if(lx==0) hx -= 1; in __nextafter() 71 if(lx==0) hx += 1; in __nextafter() [all …]
|