drem
drem + 1) > drem) {
size_t drem;
drem = PAXPATHLEN - dlen;
double drem(double, double);
else if ( (t=drem(y,two)) == zero)
x=drem(x,PI2); /* reduce x into [-PI,PI] */
x=drem(x,PI2); /* reduce x into [-PI,PI] */
p *= b; x = drem(x,p); x *= b; return(drem(x,p)/b);}
{ p /= 2 ; x /= 2; return(drem(x,p)*2);}
x = drem(x,PI); /* reduce x into [-PI/2, PI/2] */