/l4re-core-master/uclibc/lib/contrib/uclibc/libm/ |
A D | e_fmod.c | 31 hx ^=sx; /* |x| */ in __ieee754_fmod() 38 if(hx<=hy) { in __ieee754_fmod() 46 if(hx==0) { in __ieee754_fmod() 64 hx = 0x00100000|(0x000fffff&hx); in __ieee754_fmod() 68 hx = (hx<<n)|(lx>>(32-n)); in __ieee754_fmod() 92 if(hz<0){hx = hx+hx+(lx>>31); lx = lx+lx;} in __ieee754_fmod() 106 hx = hx+hx+(lx>>31); lx = lx+lx; in __ieee754_fmod() 110 hx = ((hx-0x00100000)|((iy+1023)<<20)); in __ieee754_fmod() 116 hx >>= n; in __ieee754_fmod() 118 lx = (hx<<(32-n))|(lx>>n); hx = sx; in __ieee754_fmod() [all …]
|
A D | s_nextafter.c | 24 int32_t hx,hy,ix,iy; in nextafter() local 27 EXTRACT_WORDS(hx,lx,x); in nextafter() 41 if(hx>=0) { /* x > 0 */ in nextafter() 42 if(hx>hy||((hx==hy)&&(lx>ly))) { /* x > y, x -= ulp */ in nextafter() 43 if(lx==0) hx -= 1; in nextafter() 47 if(lx==0) hx += 1; in nextafter() 50 if(hy>=0||hx>hy||((hx==hy)&&(lx>ly))){/* x < y, x -= ulp */ in nextafter() 51 if(lx==0) hx -= 1; in nextafter() 55 if(lx==0) hx += 1; in nextafter() 58 hy = hx&0x7ff00000; in nextafter() [all …]
|
A D | e_log2.c | 76 int32_t k,hx,i,j; in __ieee754_log2() local 79 EXTRACT_WORDS(hx,lx,x); in __ieee754_log2() 82 if (hx < 0x00100000) { /* x < 2**-1022 */ in __ieee754_log2() 83 if (((hx&0x7fffffff)|lx)==0) in __ieee754_log2() 87 GET_HIGH_WORD(hx,x); in __ieee754_log2() 89 if (hx >= 0x7ff00000) return x+x; in __ieee754_log2() 90 k += (hx>>20)-1023; in __ieee754_log2() 91 hx &= 0x000fffff; in __ieee754_log2() 92 i = (hx+0x95f64)&0x100000; in __ieee754_log2() 104 i = hx-0x6147a; in __ieee754_log2() [all …]
|
A D | s_nextafterf.c | 21 int32_t hx, hy, ix, iy; in nextafterf() local 23 GET_FLOAT_WORD(hx, x); in nextafterf() 25 ix = hx & 0x7fffffff; /* |x| */ in nextafterf() 45 if (hx >= 0) { /* x > 0 */ in nextafterf() 46 if (hx > hy) { /* x > y: x -= ulp */ in nextafterf() 47 hx -= 1; in nextafterf() 49 hx += 1; in nextafterf() 53 hx -= 1; in nextafterf() 55 hx += 1; in nextafterf() 58 hy = hx & 0x7f800000; in nextafterf() [all …]
|
A D | s_ilogb.c | 24 int32_t hx,lx,ix; in ilogb() local 26 GET_HIGH_WORD(hx, x); in ilogb() 27 hx &= 0x7fffffff; in ilogb() 29 if (hx < 0x00100000) { in ilogb() 31 if ((hx|lx)==0) /* +-0, ilogb(0) = FP_ILOGB0 */ in ilogb() 35 if (hx != 0) { in ilogb() 37 lx = (hx << 11); in ilogb() 45 if (hx < 0x7ff00000) /* normal x */ in ilogb() 46 return (hx>>20) - 1023; in ilogb() 50 if (hx == 0x7ff00000 && lx == 0) /* +-inf */ in ilogb()
|
A D | e_log10.c | 59 int32_t i,k,hx; in __ieee754_log10() local 62 EXTRACT_WORDS(hx,lx,x); in __ieee754_log10() 65 if (hx < 0x00100000) { /* x < 2**-1022 */ in __ieee754_log10() 66 if (((hx&0x7fffffff)|lx)==0) in __ieee754_log10() 68 if (hx<0) return (x-x)/zero; /* log(-#) = NaN */ in __ieee754_log10() 70 GET_HIGH_WORD(hx,x); in __ieee754_log10() 72 if (hx >= 0x7ff00000) return x+x; in __ieee754_log10() 73 k += (hx>>20)-1023; in __ieee754_log10() 75 hx = (hx&0x000fffff)|((0x3ff-i)<<20); in __ieee754_log10() 77 SET_HIGH_WORD(x,hx); in __ieee754_log10()
|
A D | e_log.c | 83 int32_t k,hx,i,j; in __ieee754_log() local 86 EXTRACT_WORDS(hx,lx,x); in __ieee754_log() 89 if (hx < 0x00100000) { /* x < 2**-1022 */ in __ieee754_log() 90 if (((hx&0x7fffffff)|lx)==0) in __ieee754_log() 94 GET_HIGH_WORD(hx,x); in __ieee754_log() 96 if (hx >= 0x7ff00000) return x+x; in __ieee754_log() 97 k += (hx>>20)-1023; in __ieee754_log() 98 hx &= 0x000fffff; in __ieee754_log() 99 i = (hx+0x95f64)&0x100000; in __ieee754_log() 114 i = hx-0x6147a; in __ieee754_log() [all …]
|
A D | s_frexp.c | 30 int32_t hx, ix, lx; in frexp() local 31 EXTRACT_WORDS(hx,lx,x); in frexp() 32 ix = 0x7fffffff&hx; in frexp() 37 GET_HIGH_WORD(hx,x); in frexp() 38 ix = hx&0x7fffffff; in frexp() 42 hx = (hx&0x800fffff)|0x3fe00000; in frexp() 43 SET_HIGH_WORD(x,hx); in frexp()
|
A D | s_isnan.c | 22 int32_t hx,lx; in __isnan() local 23 EXTRACT_WORDS(hx,lx,x); in __isnan() 24 hx &= 0x7fffffff; in __isnan() 25 hx |= (u_int32_t)(lx|(-lx))>>31; in __isnan() 26 hx = 0x7ff00000 - hx; in __isnan() 27 return (int)(((u_int32_t)hx)>>31); in __isnan()
|
A D | s_fpclassify.c | 26 u_int32_t hx, lx; in __fpclassify() local 29 EXTRACT_WORDS (hx, lx, x); in __fpclassify() 30 lx |= hx & 0xfffff; in __fpclassify() 31 hx &= 0x7ff00000; in __fpclassify() 32 if ((hx | lx) == 0) in __fpclassify() 34 else if (hx == 0) in __fpclassify() 36 else if (hx == 0x7ff00000) in __fpclassify()
|
A D | s_cbrt.c | 31 int32_t hx; in cbrt() local 36 GET_HIGH_WORD(hx,x); in cbrt() 37 sign=hx&0x80000000; /* sign= sign(x) */ in cbrt() 38 hx ^=sign; in cbrt() 39 if(hx>=0x7ff00000) return(x+x); /* cbrt(NaN,INF) is itself */ in cbrt() 41 if((hx|low)==0) in cbrt() 44 SET_HIGH_WORD(x,hx); /* x <- |x| */ in cbrt() 46 if(hx<0x00100000) /* subnormal number */ in cbrt() 51 SET_HIGH_WORD(t,hx/3+B1); in cbrt()
|
A D | s_scalbn.c | 31 int32_t k, hx, lx; in scalbln() local 33 EXTRACT_WORDS(hx, lx, x); in scalbln() 34 k = (hx & 0x7ff00000) >> 20; /* extract exponent */ in scalbln() 36 if ((lx | (hx & 0x7fffffff)) == 0) in scalbln() 39 GET_HIGH_WORD(hx, x); in scalbln() 40 k = ((hx & 0x7ff00000) >> 20) - 54; in scalbln() 50 SET_HIGH_WORD(x, (hx & 0x800fffff) | (k << 20)); in scalbln() 59 SET_HIGH_WORD(x, (hx & 0x800fffff) | (k << 20)); in scalbln()
|
A D | s_expm1.c | 129 u_int32_t hx; in expm1() local 131 GET_HIGH_WORD(hx,x); in expm1() 132 xsb = hx&0x80000000; /* sign bit of x */ in expm1() 134 hx &= 0x7fffffff; /* high word of |x| */ in expm1() 137 if(hx >= 0x4043687A) { /* if |x|>=56*ln2 */ in expm1() 138 if(hx >= 0x40862E42) { /* if |x|>=709.78... */ in expm1() 139 if(hx>=0x7ff00000) { in expm1() 142 if(((hx&0xfffff)|low)!=0) in expm1() 155 if(hx > 0x3fd62e42) { /* if |x| > 0.5 ln2 */ in expm1() 156 if(hx < 0x3FF0A2B2) { /* and |x| < 1.5 ln2 */ in expm1() [all …]
|
A D | e_exp.c | 105 u_int32_t hx; in __ieee754_exp() local 107 GET_HIGH_WORD(hx,x); in __ieee754_exp() 108 xsb = (hx>>31)&1; /* sign bit of x */ in __ieee754_exp() 109 hx &= 0x7fffffff; /* high word of |x| */ in __ieee754_exp() 112 if(hx >= 0x40862E42) { /* if |x|>=709.78... */ in __ieee754_exp() 113 if(hx>=0x7ff00000) { in __ieee754_exp() 116 if(((hx&0xfffff)|lx)!=0) in __ieee754_exp() 125 if(hx > 0x3fd62e42) { /* if |x| > 0.5 ln2 */ in __ieee754_exp() 126 if(hx < 0x3FF0A2B2) { /* and |x| < 1.5 ln2 */ in __ieee754_exp() 136 else if(hx < 0x3e300000) { /* when |x|<2**-28 */ in __ieee754_exp()
|
A D | e_acosh.c | 36 int32_t hx; in __ieee754_acosh() local 38 EXTRACT_WORDS(hx,lx,x); in __ieee754_acosh() 39 if(hx<0x3ff00000) { /* x < 1 */ in __ieee754_acosh() 41 } else if(hx >=0x41b00000) { /* x > 2**28 */ in __ieee754_acosh() 42 if(hx >=0x7ff00000) { /* x is inf of NaN */ in __ieee754_acosh() 46 } else if(((hx-0x3ff00000)|lx)==0) { in __ieee754_acosh() 48 } else if (hx > 0x40000000) { /* 2**28 > x > 2 */ in __ieee754_acosh()
|
A D | e_remainder.c | 28 int32_t hx,hp; in __ieee754_remainder() local 32 EXTRACT_WORDS(hx,lx,x); in __ieee754_remainder() 34 sx = hx&0x80000000; in __ieee754_remainder() 36 hx &= 0x7fffffff; in __ieee754_remainder() 40 if((hx>=0x7ff00000)|| /* x not finite */ in __ieee754_remainder() 47 if (((hx-hp)|(lx-lp))==0) return zero*x; in __ieee754_remainder() 62 GET_HIGH_WORD(hx,x); in __ieee754_remainder() 63 SET_HIGH_WORD(x,hx^sx); in __ieee754_remainder()
|
A D | s_erf.c | 186 int32_t hx,ix,i; in erf() local 188 GET_HIGH_WORD(hx,x); in erf() 189 ix = hx&0x7fffffff; in erf() 191 i = ((u_int32_t)hx>>31)<<1; in erf() 232 if(hx>=0) return one-r/x; else return r/x-one; in erf() 238 int32_t hx,ix; in libm_hidden_def() local 240 GET_HIGH_WORD(hx,x); in libm_hidden_def() 241 ix = hx&0x7fffffff; in libm_hidden_def() 254 if(hx < 0x3fd00000) { /* x<1/4 */ in libm_hidden_def() 266 if(hx>=0) { in libm_hidden_def() [all …]
|
A D | s_log1p.c | 97 int32_t k,hx,hu=0,ax; in log1p() local 99 GET_HIGH_WORD(hx,x); in log1p() 100 ax = hx&0x7fffffff; in log1p() 103 if (hx < 0x3FDA827A) { /* x < 0.41422 */ in log1p() 115 if(hx>0||hx<=((int32_t)0xbfd2bec3)) { in log1p() 118 if (hx >= 0x7ff00000) return x+x; in log1p() 120 if(hx<0x43400000) { in log1p()
|
A D | s_isinf.c | 17 int32_t hx,lx; in __isinf() local 18 EXTRACT_WORDS(hx,lx,x); in __isinf() 19 lx |= (hx & 0x7fffffff) ^ 0x7ff00000; in __isinf() 21 return ~(lx >> 31) & (hx >> 30); in __isinf()
|
A D | e_atan2.c | 53 int32_t k,m,hx,hy,ix,iy; in __ieee754_atan2() local 56 EXTRACT_WORDS(hx,lx,x); in __ieee754_atan2() 57 ix = hx&0x7fffffff; in __ieee754_atan2() 63 if(((hx-0x3ff00000)|lx)==0) return atan(y); /* x=1.0 */ in __ieee754_atan2() 64 m = ((hy>>31)&1)|((hx>>30)&2); /* 2*sign(x)+sign(y) */ in __ieee754_atan2() 102 else if(hx<0&&k<-60) z=0.0; /* |y|/x < -2**60 */ in __ieee754_atan2()
|
A D | e_lgamma_r.c | 203 int i,hx,lx,ix; in __ieee754_lgamma_r() local 205 EXTRACT_WORDS(hx,lx,x); in __ieee754_lgamma_r() 209 ix = hx&0x7fffffff; in __ieee754_lgamma_r() 217 if(hx<0) { in __ieee754_lgamma_r() 222 if(hx<0) { in __ieee754_lgamma_r() 294 if(hx<0) r = nadj - r; in __ieee754_lgamma_r() 371 int32_t hx; in strong_alias() local 378 EXTRACT_WORDS(hx, lx, x); in strong_alias() 380 if (((hx & 0x7fffffff) | lx) == 0) { in strong_alias() 384 if (hx < 0 && (u_int32_t)hx < 0xfff00000 && rint(x) == x) { in strong_alias() [all …]
|
A D | s_atan.c | 70 int32_t ix,hx,id; in atan() local 72 GET_HIGH_WORD(hx,x); in atan() 73 ix = hx&0x7fffffff; in atan() 80 if(hx>0) return atanhi[3]+atanlo[3]; in atan() 111 return (hx<0)? -z:z; in atan()
|
A D | e_jn.c | 50 int32_t i,hx,ix,lx, sgn; in __ieee754_jn() local 57 EXTRACT_WORDS(hx,lx,x); in __ieee754_jn() 58 ix = 0x7fffffff&hx; in __ieee754_jn() 64 hx ^= 0x80000000; in __ieee754_jn() 68 sgn = (n&1)&(hx>>31); /* even n -- 0, odd n -- sign(x) */ in __ieee754_jn() 222 int32_t i,hx,ix,lx; in strong_alias() local 226 EXTRACT_WORDS(hx,lx,x); in strong_alias() 227 ix = 0x7fffffff&hx; in strong_alias() 231 if(hx<0) return zero/zero; in strong_alias()
|
A D | s_asinh.c | 34 int32_t hx,ix; in asinh() local 35 GET_HIGH_WORD(hx,x); in asinh() 36 ix = hx&0x7fffffff; in asinh() 50 if(hx>0) return w; else return -w; in asinh()
|
A D | k_tan.c | 72 int32_t ix,hx; in __kernel_tan() local 73 GET_HIGH_WORD(hx,x); in __kernel_tan() 74 ix = hx&0x7fffffff; /* high word of |x| */ in __kernel_tan() 84 if(hx<0) {x = -x; y = -y;} in __kernel_tan() 103 return (double)(1-((hx>>30)&2))*(v-2.0*(x-(w*w/(w+v)-r))); in __kernel_tan()
|