| /netbsd-src/external/gpl3/gcc.old/dist/libquadmath/math/ |
| H A D | fmodq.c | 28 int64_t n,hx,hy,hz,ix,iy,sx,i; in fmodq() local 31 GET_FLT128_WORDS64(hx,lx,x); in fmodq() 33 sx = hx&0x8000000000000000ULL; /* sign of x */ in fmodq() 34 hx ^=sx; /* |x| */ in fmodq() 38 if((hy|ly)==0||(hx>=0x7fff000000000000LL)|| /* y=0,or x not finite */ in fmodq() 41 if(hx<=hy) { in fmodq() 42 if((hx<hy)||(lx<ly)) return x; /* |x|<|y| return x */ in fmodq() 48 if(hx<0x0001000000000000LL) { /* subnormal x */ in fmodq() 49 if(hx==0) { in fmodq() 52 for (ix = -16382, i=hx<<15; i>0; i<<=1) ix -=1; in fmodq() [all …]
|
| H A D | nextafterq.c | 31 int64_t hx,hy,ix,iy; in nextafterq() local 34 GET_FLT128_WORDS64(hx,lx,x); in nextafterq() 36 ix = hx&0x7fffffffffffffffLL; /* |x| */ in nextafterq() 51 if(hx>=0) { /* x > 0 */ in nextafterq() 52 if(hx>hy||((hx==hy)&&(lx>ly))) { /* x > y, x -= ulp */ in nextafterq() 53 if(lx==0) hx--; in nextafterq() 57 if(lx==0) hx++; in nextafterq() 60 if(hy>=0||hx>hy||((hx==hy)&&(lx>ly))){/* x < y, x -= ulp */ in nextafterq() 61 if(lx==0) hx--; in nextafterq() 65 if(lx==0) hx++; in nextafterq() [all …]
|
| H A D | ilogbq.c | 38 int64_t hx,lx; in ilogbq() local 41 GET_FLT128_WORDS64(hx,lx,x); in ilogbq() 42 hx &= 0x7fffffffffffffffLL; in ilogbq() 43 if(hx <= 0x0001000000000000LL) { in ilogbq() 44 if((hx|lx)==0) in ilogbq() 47 if(hx==0) { in ilogbq() 50 for (ix = -16382, hx<<=15; hx>0; hx<<=1) ix -=1; in ilogbq() 54 else if (hx<0x7fff000000000000LL) return (hx>>48)-0x3fff; in ilogbq() 57 if (((hx^0x7fff000000000000LL)|lx) == 0) in ilogbq()
|
| H A D | logbq.c | 31 int64_t lx, hx, ex; in logbq() local 33 GET_FLT128_WORDS64 (hx, lx, x); in logbq() 34 hx &= 0x7fffffffffffffffLL; /* high |x| */ in logbq() 35 if ((hx | lx) == 0) in logbq() 37 if (hx >= 0x7fff000000000000LL) in logbq() 39 if ((ex = hx >> 48) == 0) /* IEEE 754 logb */ in logbq() 44 if (hx == 0) in logbq() 47 ma = __builtin_clzll (hx); in logbq()
|
| /netbsd-src/external/gpl3/gcc/dist/libquadmath/math/ |
| H A D | fmodq.c | 28 int64_t n,hx,hy,hz,ix,iy,sx,i; in fmodq() local 31 GET_FLT128_WORDS64(hx,lx,x); in fmodq() 33 sx = hx&0x8000000000000000ULL; /* sign of x */ in fmodq() 34 hx ^=sx; /* |x| */ in fmodq() 38 if((hy|ly)==0||(hx>=0x7fff000000000000LL)|| /* y=0,or x not finite */ in fmodq() 41 if(hx<=hy) { in fmodq() 42 if((hx<hy)||(lx<ly)) return x; /* |x|<|y| return x */ in fmodq() 48 if(hx<0x0001000000000000LL) { /* subnormal x */ in fmodq() 49 if(hx==0) { in fmodq() 52 for (ix = -16382, i=hx<<15; i>0; i<<=1) ix -=1; in fmodq() [all …]
|
| H A D | nextafterq.c | 31 int64_t hx,hy,ix,iy; in nextafterq() local 34 GET_FLT128_WORDS64(hx,lx,x); in nextafterq() 36 ix = hx&0x7fffffffffffffffLL; /* |x| */ in nextafterq() 51 if(hx>=0) { /* x > 0 */ in nextafterq() 52 if(hx>hy||((hx==hy)&&(lx>ly))) { /* x > y, x -= ulp */ in nextafterq() 53 if(lx==0) hx--; in nextafterq() 57 if(lx==0) hx++; in nextafterq() 60 if(hy>=0||hx>hy||((hx==hy)&&(lx>ly))){/* x < y, x -= ulp */ in nextafterq() 61 if(lx==0) hx--; in nextafterq() 65 if(lx==0) hx++; in nextafterq() [all …]
|
| H A D | ilogbq.c | 38 int64_t hx,lx; in ilogbq() local 41 GET_FLT128_WORDS64(hx,lx,x); in ilogbq() 42 hx &= 0x7fffffffffffffffLL; in ilogbq() 43 if(hx <= 0x0001000000000000LL) { in ilogbq() 44 if((hx|lx)==0) in ilogbq() 47 if(hx==0) { in ilogbq() 50 for (ix = -16382, hx<<=15; hx>0; hx<<=1) ix -=1; in ilogbq() 54 else if (hx<0x7fff000000000000LL) return (hx>>48)-0x3fff; in ilogbq() 57 if (((hx^0x7fff000000000000LL)|lx) == 0) in ilogbq()
|
| H A D | logbq.c | 31 int64_t lx, hx, ex; in logbq() local 33 GET_FLT128_WORDS64 (hx, lx, x); in logbq() 34 hx &= 0x7fffffffffffffffLL; /* high |x| */ in logbq() 35 if ((hx | lx) == 0) in logbq() 37 if (hx >= 0x7fff000000000000LL) in logbq() 39 if ((ex = hx >> 48) == 0) /* IEEE 754 logb */ in logbq() 44 if (hx == 0) in logbq() 47 ma = __builtin_clzll (hx); in logbq()
|
| /netbsd-src/lib/libm/src/ |
| H A D | e_fmod.c | 36 int32_t n,hx,hy,hz,ix,iy,sx,i; in __strong_alias() local 39 EXTRACT_WORDS(hx,lx,x); in __strong_alias() 41 sx = hx&0x80000000; /* sign of x */ in __strong_alias() 42 hx ^=sx; /* |x| */ in __strong_alias() 46 if((hy|ly)==0||(hx>=0x7ff00000)|| /* y=0,or x not finite */ in __strong_alias() 49 if(hx<=hy) { in __strong_alias() 50 if((hx<hy)||(lx<ly)) return x; /* |x|<|y| return x */ in __strong_alias() 56 if(hx<0x00100000) { /* subnormal x */ in __strong_alias() 57 if(hx==0) { in __strong_alias() 60 for (ix = -1022,i=(hx<<11); i>0; i<<=1) ix -=1; in __strong_alias() [all …]
|
| H A D | s_remquo.c | 39 int32_t n,hx,hy,hz,ix,iy,sx,i; in remquo() local 42 EXTRACT_WORDS(hx,lx,x); in remquo() 44 sxy = (hx ^ hy) & 0x80000000; in remquo() 45 sx = hx&0x80000000; /* sign of x */ in remquo() 46 hx ^=sx; /* |x| */ in remquo() 50 if((hy|ly)==0||(hx>=0x7ff00000)|| /* y=0,or x not finite */ in remquo() 53 if(hx<=hy) { in remquo() 54 if((hx<hy)||(lx<ly)) { in remquo() 65 if(hx<0x00100000) { /* subnormal x */ in remquo() 66 if(hx in remquo() [all...] |
| H A D | e_fmodf.c | 35 int32_t n,hx,hy,hz,ix,iy,sx,i; in __ieee754_fmodf() local 37 GET_FLOAT_WORD(hx,x); in __ieee754_fmodf() 39 sx = hx&0x80000000; /* sign of x */ in __ieee754_fmodf() 40 hx ^=sx; /* |x| */ in __ieee754_fmodf() 44 if(hy==0||(hx>=0x7f800000)|| /* y=0,or x not finite */ in __ieee754_fmodf() 47 if(hx<hy) return x; /* |x|<|y| return x */ in __ieee754_fmodf() 48 if(hx==hy) in __ieee754_fmodf() 52 if(hx<0x00800000) { /* subnormal x */ in __ieee754_fmodf() 53 for (ix = -126,i=(hx<<8); i>0; i<<=1) ix -=1; in __ieee754_fmodf() 54 } else ix = (hx>>23)-127; in __ieee754_fmodf() [all …]
|
| H A D | s_remquof.c | 37 int32_t n,hx,hy,hz,ix,iy,sx,i; in remquof() local 40 GET_FLOAT_WORD(hx,x); in remquof() 42 sxy = (hx ^ hy) & 0x80000000; in remquof() 43 sx = hx&0x80000000; /* sign of x */ in remquof() 44 hx ^=sx; /* |x| */ in remquof() 48 if(hy==0||hx>=0x7f800000||hy>0x7f800000) /* y=0,NaN;or x not finite */ in remquof() 50 if(hx<hy) { in remquof() 53 } else if(hx==hy) { in remquof() 59 if(hx<0x00800000) { /* subnormal x */ in remquof() 60 for (ix = -126,i=(hx<<8); i>0; i<<=1) ix -=1; in remquof() [all …]
|
| H A D | s_nextafter.c | 37 int32_t hx,hy,ix,iy; in __strong_alias() local 40 EXTRACT_WORDS(hx,lx,x); in __strong_alias() 42 ix = hx&0x7fffffff; /* |x| */ in __strong_alias() 54 if(hx>=0) { /* x > 0 */ in __strong_alias() 55 if(hx>hy||((hx==hy)&&(lx>ly))) { /* x > y, x -= ulp */ in __strong_alias() 56 if(lx==0) hx -= 1; in __strong_alias() 60 if(lx==0) hx += 1; in __strong_alias() 63 if(hy>=0||hx>hy||((hx==hy)&&(lx>ly))){/* x < y, x -= ulp */ in __strong_alias() 64 if(lx==0) hx -= 1; in __strong_alias() 68 if(lx==0) hx += 1; in __strong_alias() [all …]
|
| H A D | e_fmodl.c | 50 #define SET_NBIT(hx) ((hx) | (1ULL << EXT_FRACHBITS)) argument 54 #define SET_NBIT(hx) (hx) argument 77 int64_t hx,hz; /* We need a carry bit even if EXT_FRACHBITS is 32. */ in __ieee754_fmodl() local 118 /* set up {hx,lx}, {hy,ly} and align y to x */ in __ieee754_fmodl() 119 hx = SET_NBIT(ux.extu_frach); in __ieee754_fmodl() 128 hz=hx-hy;lz=lx-ly; if(lx<ly) hz -= 1; in __ieee754_fmodl() 129 if(hz<0){hx = hx in __ieee754_fmodl() [all...] |
| H A D | s_nextafterf.c | 27 int32_t hx,hy,ix,iy; in nextafterf() local 29 GET_FLOAT_WORD(hx,x); in nextafterf() 31 ix = hx&0x7fffffff; /* |x| */ in nextafterf() 43 if(hx>=0) { /* x > 0 */ in nextafterf() 44 if(hx>hy) { /* x > y, x -= ulp */ in nextafterf() 45 hx -= 1; in nextafterf() 47 hx += 1; in nextafterf() 50 if(hy>=0||hx>hy){ /* x < y, x -= ulp */ in nextafterf() 51 hx -= 1; in nextafterf() 53 hx += 1; in nextafterf() [all …]
|
| H A D | s_remquol.c | 53 #define SET_NBIT(hx) ((hx) | (1ULL << EXT_FRACHBITS)) argument 56 #define SET_NBIT(hx) (hx) argument 81 int64_t hx,hz; /* We need a carry bit even if EXT_FRACHBITS is 32. */ in remquol() local 129 /* set up {hx,lx}, {hy,ly} and align y to x */ in remquol() 130 hx = SET_NBIT(ux.extu_frach); in remquol() 139 hz=hx-hy;lz=lx-ly; if(lx<ly) hz -= 1; in remquol() 140 if(hz<0){hx = hx in remquol() [all...] |
| H A D | e_log2.c | 39 int32_t k,hx,i,j; in __ieee754_log2() local 42 EXTRACT_WORDS(hx,lx,x); in __ieee754_log2() 45 if (hx < 0x00100000) { /* x < 2**-1022 */ in __ieee754_log2() 46 if (((hx&0x7fffffff)|lx)==0) in __ieee754_log2() 48 if (hx<0) return (x-x)/zero; /* log(-#) = NaN */ in __ieee754_log2() 50 GET_HIGH_WORD(hx,x); in __ieee754_log2() 52 if (hx >= 0x7ff00000) return x+x; in __ieee754_log2() 53 k += (hx>>20)-1023; in __ieee754_log2() 54 hx &= 0x000fffff; in __ieee754_log2() 55 i = (hx+0x95f64)&0x100000; in __ieee754_log2() [all …]
|
| H A D | s_ilogb.c | 39 int32_t hx, lx, ix; in __strong_alias() local 41 GET_HIGH_WORD(hx, x); in __strong_alias() 42 hx &= 0x7fffffff; in __strong_alias() 43 if (hx < 0x00100000) { in __strong_alias() 45 if ((hx | lx) == 0) { in __strong_alias() 49 if (hx == 0) { in __strong_alias() 52 for (ix = -1022, hx <<= 11; hx > 0; hx <<= 1) ix -= 1; in __strong_alias() 57 if (hx < 0x7ff00000) { in __strong_alias() 58 return (hx >> 20) - 1023; in __strong_alias()
|
| H A D | e_log.c | 90 int32_t k,hx,i,j; in __ieee754_log() local 93 EXTRACT_WORDS(hx,lx,x); in __ieee754_log() 96 if (hx < 0x00100000) { /* x < 2**-1022 */ in __ieee754_log() 97 if (((hx&0x7fffffff)|lx)==0) in __ieee754_log() 99 if (hx<0) return (x-x)/zero; /* log(-#) = NaN */ in __ieee754_log() 101 GET_HIGH_WORD(hx,x); in __ieee754_log() 103 if (hx >= 0x7ff00000) return x+x; in __ieee754_log() 104 k += (hx>>20)-1023; in __ieee754_log() 105 hx &= 0x000fffff; in __ieee754_log() 106 i = (hx+0x95f64)&0x100000; in __ieee754_log() [all …]
|
| H A D | e_log10f.c | 36 int32_t i,k,hx; in __ieee754_log10f() local 38 GET_FLOAT_WORD(hx,x); in __ieee754_log10f() 41 if (hx < 0x00800000) { /* x < 2**-126 */ in __ieee754_log10f() 42 if ((hx&0x7fffffff)==0) in __ieee754_log10f() 44 if (hx<0) return (x-x)/zero; /* log(-#) = NaN */ in __ieee754_log10f() 46 GET_FLOAT_WORD(hx,x); in __ieee754_log10f() 48 if (hx >= 0x7f800000) return x+x; in __ieee754_log10f() 49 k += (hx>>23)-127; in __ieee754_log10f() 51 hx = (hx&0x007fffff)|((0x7f-i)<<23); in __ieee754_log10f() 53 SET_FLOAT_WORD(x,hx); in __ieee754_log10f()
|
| H A D | e_log10.c | 66 int32_t i,k,hx; in __ieee754_log10() local 69 EXTRACT_WORDS(hx,lx,x); in __ieee754_log10() 72 if (hx < 0x00100000) { /* x < 2**-1022 */ in __ieee754_log10() 73 if (((hx&0x7fffffff)|lx)==0) in __ieee754_log10() 75 if (hx<0) return (x-x)/zero; /* log(-#) = NaN */ in __ieee754_log10() 77 GET_HIGH_WORD(hx,x); in __ieee754_log10() 79 if (hx >= 0x7ff00000) return x+x; in __ieee754_log10() 80 k += (hx>>20)-1023; in __ieee754_log10() 82 hx = (hx&0x000fffff)|((0x3ff-i)<<20); in __ieee754_log10() 84 SET_HIGH_WORD(x,hx); in __ieee754_log10()
|
| H A D | s_ilogbf.c | 32 int32_t hx, ix; in ilogbf() local 34 GET_FLOAT_WORD(hx, x); in ilogbf() 35 hx &= 0x7fffffff; in ilogbf() 36 if (hx < 0x00800000) { in ilogbf() 37 if (hx == 0) { in ilogbf() 41 for (ix = -126, hx <<= 8; hx > 0; hx <<= 1) ix -= 1; in ilogbf() 45 if (hx < 0x7f800000) { in ilogbf() 46 return (hx >> 23) - 127; in ilogbf()
|
| H A D | s_cbrtf.c | 42 int32_t hx; in cbrtf() local 46 GET_FLOAT_WORD(hx,x); in cbrtf() 47 sign=hx&0x80000000; /* sign= sign(x) */ in cbrtf() 48 hx ^=sign; in cbrtf() 49 if(hx>=0x7f800000) return(x+x); /* cbrt(NaN,INF) is itself */ in cbrtf() 50 if(hx==0) in cbrtf() 53 SET_FLOAT_WORD(x,hx); /* x <- |x| */ in cbrtf() 55 if(hx<0x00800000) /* subnormal number */ in cbrtf() 60 SET_FLOAT_WORD(t,hx/3+B1); in cbrtf()
|
| H A D | s_frexpf.c | 30 int32_t hx,ix; in frexpf() local 31 GET_FLOAT_WORD(hx,x); in frexpf() 32 ix = 0x7fffffff&hx; in frexpf() 37 GET_FLOAT_WORD(hx,x); in frexpf() 38 ix = hx&0x7fffffff; in frexpf() 42 hx = (hx&0x807fffff)|0x3f000000; in frexpf() 43 SET_FLOAT_WORD(x,hx); in frexpf()
|
| H A D | s_nexttoward.c | 50 int32_t hx, ix; in nexttoward() local 53 EXTRACT_WORDS(hx, lx, x); in nexttoward() 54 ix = hx & 0x7fffffff; /* |x| */ in nexttoward() 74 if ((hx >= 0) ^ (x < y)) { /* x -= ulp */ in nexttoward() 75 if (lx == 0) hx -= 1; in nexttoward() 79 if (lx == 0) hx += 1; in nexttoward() 81 ix = hx & 0x7ff00000; in nexttoward() 86 INSERT_WORDS(y, hx, lx); in nexttoward() 90 INSERT_WORDS(x, hx, lx); in nexttoward()
|