double remainder(double x, double y) /* wrapper remainder */ { #ifdef CYGSEM_LIBM_COMPAT_IEEE_ONLY return __ieee754_remainder(x,y); #else double z; z = __ieee754_remainder(x,y); if(cyg_libm_get_compat_mode() == CYGNUM_LIBM_COMPAT_IEEE || isnan(y)) return z; if(y==0.0) return __kernel_standard(x,y,28); /* remainder(x,0) */ else return z; #endif }
double remainder(double x, double y) /* wrapper remainder */ { #ifdef _IEEE_LIBM return __ieee754_remainder(x,y); #else double z; z = __ieee754_remainder(x,y); if(_LIB_VERSION == _IEEE_ || isnan(y)) return z; if(y==0.0) return __kernel_standard(x,y,28); /* remainder(x,0) */ else return z; #endif }
void Math_remainder(void *fp) { F_Math_remainder *f; f = fp; *f->ret = __ieee754_remainder(f->x, f->p); }
/* wrapper remainder */ double __remainder (double x, double y) { if (((__builtin_expect (y == 0.0, 0) && ! __isnan (x)) || (__builtin_expect (__isinf_ns (x), 0) && ! __isnan (y))) && _LIB_VERSION != _IEEE_) return __kernel_standard (x, y, 28); /* remainder domain */ return __ieee754_remainder (x, y); }
/* wrapper remainderf */ float remainderf (float x, float y) { #if defined(__UCLIBC_HAS_FENV__) if (((__builtin_expect (y == 0.0f, 0) && ! isnan (x)) || (__builtin_expect (isinf (x), 0) && ! isnan (y))) && _LIB_VERSION != _IEEE_) return __kernel_standard_f (x, y, 128); /* remainder domain */ #endif return (float) __ieee754_remainder ((double) x,(double) y); }
double __ieee754_remainder(double x, double y) { double z,d,xx; #if 0 double yy; #endif int4 kx,ky,n,nn,n1,m1,l; #if 0 int4 m; #endif mynumber u,t,w={{0,0}},v={{0,0}},ww={{0,0}},r; u.x=x; t.x=y; kx=u.i[HIGH_HALF]&0x7fffffff; /* no sign for x*/ t.i[HIGH_HALF]&=0x7fffffff; /*no sign for y */ ky=t.i[HIGH_HALF]; /*------ |x| < 2^1023 and 2^-970 < |y| < 2^1024 ------------------*/ if (kx<0x7fe00000 && ky<0x7ff00000 && ky>=0x03500000) { if (kx+0x00100000<ky) return x; if ((kx-0x01500000)<ky) { z=x/t.x; v.i[HIGH_HALF]=t.i[HIGH_HALF]; d=(z+big.x)-big.x; xx=(x-d*v.x)-d*(t.x-v.x); if (d-z!=0.5&&d-z!=-0.5) return (xx!=0)?xx:((x>0)?ZERO.x:nZERO.x); else { if (ABS(xx)>0.5*t.x) return (z>d)?xx-t.x:xx+t.x; else return xx; } } /* (kx<(ky+0x01500000)) */ else { r.x=1.0/t.x; n=t.i[HIGH_HALF]; nn=(n&0x7ff00000)+0x01400000; w.i[HIGH_HALF]=n; ww.x=t.x-w.x; l=(kx-nn)&0xfff00000; n1=ww.i[HIGH_HALF]; m1=r.i[HIGH_HALF]; while (l>0) { r.i[HIGH_HALF]=m1-l; z=u.x*r.x; w.i[HIGH_HALF]=n+l; ww.i[HIGH_HALF]=(n1)?n1+l:n1; d=(z+big.x)-big.x; u.x=(u.x-d*w.x)-d*ww.x; l=(u.i[HIGH_HALF]&0x7ff00000)-nn; } r.i[HIGH_HALF]=m1; w.i[HIGH_HALF]=n; ww.i[HIGH_HALF]=n1; z=u.x*r.x; d=(z+big.x)-big.x; u.x=(u.x-d*w.x)-d*ww.x; if (ABS(u.x)<0.5*t.x) return (u.x!=0)?u.x:((x>0)?ZERO.x:nZERO.x); else if (ABS(u.x)>0.5*t.x) return (d>z)?u.x+t.x:u.x-t.x; else {z=u.x/t.x; d=(z+big.x)-big.x; return ((u.x-d*w.x)-d*ww.x);} } } /* (kx<0x7fe00000&&ky<0x7ff00000&&ky>=0x03500000) */ else { if (kx<0x7fe00000&&ky<0x7ff00000&&(ky>0||t.i[LOW_HALF]!=0)) { y=ABS(y)*t128.x; z=__ieee754_remainder(x,y)*t128.x; z=__ieee754_remainder(z,y)*tm128.x; return z; } else { if ((kx&0x7ff00000)==0x7fe00000&&ky<0x7ff00000&&(ky>0||t.i[LOW_HALF]!=0)) { y=ABS(y); z=2.0*__ieee754_remainder(0.5*x,y); d = ABS(z); if (d <= ABS(d-y)) return z; else return (z>0)?z-y:z+y; } else { /* if x is too big */ if (kx == 0x7ff00000 && u.i[LOW_HALF] == 0 && y == 1.0) return x / x; if (kx>=0x7ff00000||(ky==0&&t.i[LOW_HALF]==0)||ky>0x7ff00000|| (ky==0x7ff00000&&t.i[LOW_HALF]!=0)) return (u.i[HIGH_HALF]&0x80000000)?nNAN.x:NAN.x; else return x; } } } }
Err mathlib_remainder(UInt16 refnum, double x, double y, double *result) { #pragma unused(refnum) *result = __ieee754_remainder(x, y); return mlErrNone; }
double __ieee754_remainder (double x, double y) { double z, d, xx; int4 kx, ky, n, nn, n1, m1, l; mynumber u, t, w = { { 0, 0 } }, v = { { 0, 0 } }, ww = { { 0, 0 } }, r; u.x = x; t.x = y; kx = u.i[HIGH_HALF] & 0x7fffffff; /* no sign for x*/ t.i[HIGH_HALF] &= 0x7fffffff; /*no sign for y */ ky = t.i[HIGH_HALF]; /*------ |x| < 2^1023 and 2^-970 < |y| < 2^1024 ------------------*/ if (kx < 0x7fe00000 && ky < 0x7ff00000 && ky >= 0x03500000) { SET_RESTORE_ROUND_NOEX (FE_TONEAREST); if (kx + 0x00100000 < ky) return x; if ((kx - 0x01500000) < ky) { z = x / t.x; v.i[HIGH_HALF] = t.i[HIGH_HALF]; d = (z + big.x) - big.x; xx = (x - d * v.x) - d * (t.x - v.x); if (d - z != 0.5 && d - z != -0.5) return (xx != 0) ? xx : ((x > 0) ? ZERO.x : nZERO.x); else { if (fabs (xx) > 0.5 * t.x) return (z > d) ? xx - t.x : xx + t.x; else return xx; } } /* (kx<(ky+0x01500000)) */ else { r.x = 1.0 / t.x; n = t.i[HIGH_HALF]; nn = (n & 0x7ff00000) + 0x01400000; w.i[HIGH_HALF] = n; ww.x = t.x - w.x; l = (kx - nn) & 0xfff00000; n1 = ww.i[HIGH_HALF]; m1 = r.i[HIGH_HALF]; while (l > 0) { r.i[HIGH_HALF] = m1 - l; z = u.x * r.x; w.i[HIGH_HALF] = n + l; ww.i[HIGH_HALF] = (n1) ? n1 + l : n1; d = (z + big.x) - big.x; u.x = (u.x - d * w.x) - d * ww.x; l = (u.i[HIGH_HALF] & 0x7ff00000) - nn; } r.i[HIGH_HALF] = m1; w.i[HIGH_HALF] = n; ww.i[HIGH_HALF] = n1; z = u.x * r.x; d = (z + big.x) - big.x; u.x = (u.x - d * w.x) - d * ww.x; if (fabs (u.x) < 0.5 * t.x) return (u.x != 0) ? u.x : ((x > 0) ? ZERO.x : nZERO.x); else if (fabs (u.x) > 0.5 * t.x) return (d > z) ? u.x + t.x : u.x - t.x; else { z = u.x / t.x; d = (z + big.x) - big.x; return ((u.x - d * w.x) - d * ww.x); } } } /* (kx<0x7fe00000&&ky<0x7ff00000&&ky>=0x03500000) */ else { if (kx < 0x7fe00000 && ky < 0x7ff00000 && (ky > 0 || t.i[LOW_HALF] != 0)) { y = fabs (y) * t128.x; z = __ieee754_remainder (x, y) * t128.x; z = __ieee754_remainder (z, y) * tm128.x; return z; } else { if ((kx & 0x7ff00000) == 0x7fe00000 && ky < 0x7ff00000 && (ky > 0 || t.i[LOW_HALF] != 0)) { y = fabs (y); z = 2.0 * __ieee754_remainder (0.5 * x, y); d = fabs (z); if (d <= fabs (d - y)) return z; else if (d == y) return 0.0 * x; else return (z > 0) ? z - y : z + y; } else /* if x is too big */ { if (ky == 0 && t.i[LOW_HALF] == 0) /* y = 0 */ return (x * y) / (x * y); else if (kx >= 0x7ff00000 /* x not finite */ || (ky > 0x7ff00000 /* y is NaN */ || (ky == 0x7ff00000 && t.i[LOW_HALF] != 0))) return (x * y) / (x * y); else return x; } } } }