/picolibc-3.7.0-3.6.0/newlib/libm/math/ |
D | sf_fmod.c | 32 __int32_t n, hx, hy, hz, ix, iy, sx, i; in fmodf() local 35 GET_FLOAT_WORD(hy, y); in fmodf() 38 hy &= 0x7fffffff; /* |y| */ in fmodf() 47 if (hy == 0) /* y=0, domain error */ in fmodf() 50 if (hx < hy) in fmodf() 52 if (hx == hy) in fmodf() 65 if (FLT_UWORD_IS_SUBNORMAL(hy)) { /* subnormal y */ in fmodf() 66 for (iy = -126, i = (hy << 8); i >= 0; i <<= 1) in fmodf() 69 iy = (hy >> 23) - 127; in fmodf() 79 hy = 0x00800000 | (0x007fffff & hy); in fmodf() [all …]
|
D | s_fmod.c | 31 __int32_t n, hx, hy, hz, ix, iy, sx, i; in fmod64() local 35 EXTRACT_WORDS(hy, ly, y); in fmod64() 38 hy &= 0x7fffffff; /* |y| */ in fmod64() 47 if ((hy | ly) == 0) /* y=0, domain error */ in fmod64() 50 if (hx <= hy) { in fmod64() 51 if ((hx < hy) || (lx < ly)) in fmod64() 70 if (hy < 0x00100000) { /* subnormal y */ in fmod64() 71 if (hy == 0) { in fmod64() 75 for (iy = -1022, i = (hy << 11); i > 0; i <<= 1) in fmod64() 79 iy = (hy >> 20) - 1023; in fmod64() [all …]
|
D | sf_atan2.c | 29 __int32_t k, m, hx, hy, ix, iy; in atan2f() local 33 GET_FLOAT_WORD(hy, y); in atan2f() 34 iy = hy & 0x7fffffff; in atan2f() 39 m = ((hy >> 31) & 1) | ((hx >> 30) & 2); /* 2*sign(x)+sign(y) */ in atan2f() 55 return (hy < 0) ? -pi_o_2 - tiny : pi_o_2 + tiny; in atan2f() 85 return (hy < 0) ? -pi_o_2 - tiny : pi_o_2 + tiny; in atan2f()
|
D | sf_exp.c | 87 __uint32_t hy; in expf() local 88 GET_FLOAT_WORD(hy, y); in expf() 89 SET_FLOAT_WORD(y, hy + (k << 23)); /* add k to y's exponent */ in expf() 92 __uint32_t hy; in expf() local 93 GET_FLOAT_WORD(hy, y); in expf() 94 SET_FLOAT_WORD(y, hy + ((k + 100) << 23)); /* add k to y's exponent */ in expf()
|
D | s_atan2.c | 57 __int32_t k, m, hx, hy, ix, iy; in atan264() local 62 EXTRACT_WORDS(hy, ly, y); in atan264() 63 iy = hy & 0x7fffffff; in atan264() 69 m = ((hy >> 31) & 1) | ((hx >> 30) & 2); /* 2*sign(x)+sign(y) */ in atan264() 85 return (hy < 0) ? -pi_o_2 - tiny : pi_o_2 + tiny; in atan264() 115 return (hy < 0) ? -pi_o_2 - tiny : pi_o_2 + tiny; in atan264()
|
D | s_exp.c | 154 __uint32_t hy; in exp64() local 155 GET_HIGH_WORD(hy, y); in exp64() 156 SET_HIGH_WORD(y, hy + (k << 20)); /* add k to y's exponent */ in exp64() 159 __uint32_t hy; in exp64() local 160 GET_HIGH_WORD(hy, y); in exp64() 161 SET_HIGH_WORD(y, hy + ((k + 1000) << 20)); /* add k to y's exponent */ in exp64()
|
D | s_pow.c | 103 __int32_t hx, hy, ix, iy; in pow64() local 107 EXTRACT_WORDS(hy, ly, y); in pow64() 109 iy = hy & 0x7fffffff; in pow64() 156 return (hy >= 0) ? y : zero; in pow64() 158 return (hy < 0) ? -y : zero; in pow64() 161 if (hy < 0) { in pow64() 168 if (hy == 0x40000000 && ix < 0x5ff00000 && ix > 0x1e500000) in pow64() 170 if (hy == 0x3fe00000) { /* y is 0.5 */ in pow64() 181 if (hy < 0) in pow64() 193 if (hy < 0) in pow64() [all …]
|
D | sf_pow.c | 62 __int32_t hx, hy, ix, iy, is; in powf() local 65 GET_FLOAT_WORD(hy, y); in powf() 67 iy = hy & 0x7fffffff; in powf() 106 return (hy >= 0) ? y : zero; in powf() 108 return (hy < 0) ? -y : zero; in powf() 111 if (hy < 0) { in powf() 118 if (hy == 0x40000000 && ix < 0x5f800000 && ix > 0x1a000000) in powf() 120 if (hy == 0x3f000000) { /* y is 0.5 */ in powf() 129 if (hy < 0) in powf() 141 if (hy < 0) in powf() [all …]
|
/picolibc-3.7.0-3.6.0/newlib/libm/common/ |
D | s_nextafter.c | 58 __int32_t hx,hy,ix,iy; in nextafter64() local 62 EXTRACT_WORDS(hy,ly,y); in nextafter64() 64 iy = hy&0x7fffffff; /* |y| */ in nextafter64() 71 INSERT_WORDS(x,hy&0x80000000,1); /* return +-minsubnormal */ in nextafter64() 76 if(hx>hy||((hx==hy)&&(lx>ly))) { /* x > y, x -= ulp */ in nextafter64() 84 if(hy>=0||hx>hy||((hx==hy)&&(lx>ly))){/* x < y, x -= ulp */ in nextafter64() 92 hy = hx&0x7ff00000; in nextafter64() 93 if(hy>=0x7ff00000) in nextafter64() 96 if(hy<0x00100000) /* underflow */ in nextafter64()
|
D | sf_remquo.c | 39 __int32_t n,hx,hy,hz,ix,iy,sx,i; in remquof() local 43 GET_FLOAT_WORD(hy,y); in remquof() 44 sxy = (hx ^ hy) & 0x80000000; in remquof() 47 hy &= 0x7fffffff; /* |y| */ in remquof() 50 if(hy==0||hx>=0x7f800000||hy>0x7f800000) { /* y=0,NaN;or x not finite */ in remquof() 54 if(hx<hy) { in remquof() 57 } else if(hx==hy) { in remquof() 68 if(hy<0x00800000) { /* subnormal y */ in remquof() 69 for (iy = -126,i=(hy<<8); i>0; i<<=1) iy -=1; in remquof() 70 } else iy = (hy>>23)-127; in remquof() [all …]
|
D | sf_nextafter.c | 20 __int32_t hx,hy,ix,iy; in nextafterf() local 23 GET_FLOAT_WORD(hy,y); in nextafterf() 25 iy = hy&0x7fffffff; /* |y| */ in nextafterf() 32 SET_FLOAT_WORD(x,(hy&0x80000000)|FLT_UWORD_MIN); in nextafterf() 37 if(hx>hy) { /* x > y, x -= ulp */ in nextafterf() 43 if(hy>=0||hx>hy){ /* x < y, x -= ulp */ in nextafterf() 49 hy = hx&0x7f800000; in nextafterf() 50 if(hy>FLT_UWORD_MAX) return check_oflowf(x+x); /* overflow */ in nextafterf() 52 if(hy<0x00800000) /* underflow */ in nextafterf()
|
D | s_remquo.c | 88 __int32_t n,hx,hy,hz,ix,iy,sx,i; in remquo64() local 92 EXTRACT_WORDS(hy,ly,y); in remquo64() 93 sxy = (hx ^ hy) & 0x80000000; in remquo64() 96 hy &= 0x7fffffff; /* |y| */ in remquo64() 99 if((hy|ly)==0||(hx>=0x7ff00000)|| /* y=0,or x not finite */ in remquo64() 100 ((hy|((ly|-ly)>>31))>0x7ff00000)) { /* or y is NaN */ in remquo64() 104 if(hx<=hy) { in remquo64() 105 if((hx<hy)||(lx<ly)) { in remquo64() 125 if(hy<0x00100000) { /* subnormal y */ in remquo64() 126 if(hy==0) { in remquo64() [all …]
|
D | s_copysign.c | 61 __uint32_t hx,hy; in copysign64() local 63 GET_HIGH_WORD(hy,y); in copysign64() 64 SET_HIGH_WORD(x,(hx&0x7fffffff)|(hy&0x80000000)); in copysign64()
|
/picolibc-3.7.0-3.6.0/newlib/libm/ld/ld128/ |
D | s_remquol.c | 53 int64_t hx,hz,hy,_hx; in remquol() local 59 GET_LDOUBLE_WORDS64(hy,ly,y); in remquol() 61 sxy = sx ^ ((hy>>48)&0x8000); in remquol() 63 hy &= 0x7fffffffffffffffLL; /* |y| */ in remquol() 65 SET_LDOUBLE_WORDS64(y,hy,ly); in remquol() 68 if((hy|ly)==0 || /* y=0 */ in remquol() 70 ((hy>>48) == BIAS + LDBL_MAX_EXP && in remquol() 71 (((hy&0x0000ffffffffffffLL)&~LDBL_NBIT)|ly)!=0)) /* or y is NaN */ in remquol() 73 if((hx>>48)<=(hy>>48)) { in remquol() 74 if(((hx>>48)<(hy>>48)) || in remquol() [all …]
|
D | s_nextafterl.c | 25 int64_t hx,hy,ix,iy; in nextafterl() local 29 GET_LDOUBLE_WORDS64(hy,ly,y); in nextafterl() 31 iy = hy&0x7fffffffffffffffLL; /* |y| */ in nextafterl() 38 SET_LDOUBLE_WORDS64(x,hy&0x8000000000000000ULL,1);/* return +-minsubnormal */ in nextafterl() 43 if(hx>hy||((hx==hy)&&(lx>ly))) { /* x > y, x -= ulp */ in nextafterl() 51 if(hy>=0||hx>hy||((hx==hy)&&(lx>ly))){/* x < y, x -= ulp */ in nextafterl() 61 return __math_oflowl(hy<0); in nextafterl()
|
D | e_fmodl.c | 26 int64_t n,hx,hy,hz,ix,iy,sx,i; in fmodl() local 30 GET_LDOUBLE_WORDS64(hy,ly,y); in fmodl() 33 hy &= 0x7fffffffffffffffLL; /* |y| */ in fmodl() 44 if(hx<=hy) { in fmodl() 45 if((hx<hy)||(lx<ly)) return x; /* |x|<|y| return x */ in fmodl() 60 if(hy<0x0001000000000000LL) { /* subnormal y */ in fmodl() 61 if(hy==0) { in fmodl() 64 for (iy = -16382, i=hy<<15; i>0; i<<=1) iy -=1; in fmodl() 66 } else iy = (hy>>48)-0x3fff; in fmodl() 82 hy = 0x0001000000000000LL|(0x0000ffffffffffffLL&hy); in fmodl() [all …]
|
D | s_nexttoward.c | 26 int64_t hy,iy; in nexttoward64() local 31 GET_LDOUBLE_WORDS64(hy,ly,y); in nexttoward64() 33 iy = hy&0x7fffffffffffffffLL; /* |y| */ in nexttoward64() 44 INSERT_WORDS(x,(u_int32_t)((hy>>32)&0x80000000),1);/* return +-minsub */ in nexttoward64() 67 return __math_oflow(hy < 0); in nexttoward64()
|
D | s_copysignl.c | 57 int64_t hx, hy; in copysignl() local 59 GET_LDOUBLE_MSW64(hy, y); in copysignl() 60 SET_LDOUBLE_MSW64(x, (hx & 0x7fffffffffffffffLL)|(hy & 0x8000000000000000LL)); in copysignl()
|
D | s_nexttowardf.c | 19 int64_t hy,iy; in nexttowardf() local 23 GET_LDOUBLE_WORDS64(hy,ly,y); in nexttowardf() 25 iy = hy&0x7fffffffffffffffLL; /* |y| */ in nexttowardf() 36 SET_FLOAT_WORD(x,(u_int32_t)((hy>>32)&0x80000000)|1);/* return +-minsub*/ in nexttowardf()
|
D | e_powl.c | 150 int32_t hx, hy; in powl() local 158 hy = q.parts32.mswhi; in powl() 159 iy = hy & 0x7fffffff; in powl() 218 if (hy < 0) in powl() 247 return (hy >= 0) ? y : zero; in powl() 249 return (hy < 0) ? -y : zero; in powl() 253 if (hy < 0) in powl() 258 if (hy == 0x40000000) in powl() 260 if (hy == 0x3ffe0000) in powl() 281 return (hy < 0) ? __math_oflowl(neg) : __math_uflowl(neg); in powl() [all …]
|
/picolibc-3.7.0-3.6.0/newlib/libm/ld/ld80/ |
D | s_nexttowardf.c | 19 u_int32_t hy,ly,esy; in nexttowardf() local 22 GET_LDOUBLE_WORDS(esy,hy,ly,y); in nexttowardf() 25 hy &= 0x7fffffff; /* mask off leading 1 */ in nexttowardf() 28 (iy>=0x7fff&&((hy|ly)!=0))) /* y is nan */ in nexttowardf() 49 hy = hx&0x7f800000; in nexttowardf() 50 if(hy>=0x7f800000) in nexttowardf() 53 if(hy<0x00800000) in nexttowardf()
|
D | s_nexttoward.c | 26 u_int32_t lx,hy,ly,esy; in nexttoward() local 29 GET_LDOUBLE_WORDS(esy,hy,ly,y); in nexttoward() 32 hy &= 0x7fffffff; /* mask off leading 1 */ in nexttoward() 35 ((iy>=0x7fff)&&(hy|ly)!=0)) /* y is nan */ in nexttoward() 60 hy = hx&0x7ff00000; in nexttoward() 61 if(hy>=0x7ff00000) in nexttoward() 64 if(hy<0x00100000) in nexttoward()
|
D | s_nextafterl.c | 25 u_int32_t hx,hy,ix,iy; in nextafterl() local 30 GET_LDOUBLE_WORDS(esy,hy,ly,y); in nextafterl() 35 ((iy==0x7fff)&&(((hy&0x7fffffff)|ly)!=0))) /* y is nan */ in nextafterl() 44 if(esy<0||(ix>iy||((ix==iy) && (hx>hy||((hx==hy)&&(lx>ly)))))) { in nextafterl() 59 if(esy>=0||(ix>iy||((ix==iy)&&(hx>hy||((hx==hy)&&(lx>ly)))))){ in nextafterl()
|
D | s_remquol.c | 52 uint32_t hy; in remquol() local 58 GET_LDOUBLE_WORDS(esy,hy,ly,y); in remquol() 67 if((esy|hy|ly)==0 || /* y=0 */ in remquol() 70 ((hy&~LDBL_NBIT)|ly)!=0)) { /* or y is NaN */ in remquol() 76 (hx<=hy && in remquol() 77 (hx<hy || in remquol() 82 if(hx==hy && lx==ly) { in remquol() 100 GET_LDOUBLE_WORDS(esy,hy,ly,y); in remquol() 115 hz=hx-hy;lz=lx-ly; if(lx<ly) hz -= 1; in remquol() 120 hz=hx-hy;lz=lx-ly; if(lx<ly) hz -= 1; in remquol()
|
D | e_fmodl.c | 55 uint32_t hy; in fmodl() local 103 hy = SET_NBIT(uy.bits.ext_frach); in fmodl() 111 hz=hx-hy;lz=lx-ly; if(lx<ly) hz -= 1; in fmodl() 119 hz=hx-hy;lz=lx-ly; if(lx<ly) hz -= 1; in fmodl()
|