scalbn
extern double scalbn(double x, int y);
z = scalbn(z,q0); /* actual value of z */
z -= scalbn(1.0,q0);
z = scalbn(z,-q0);
fw = scalbn(1.0,q0);
return scalbn(r, e);
return scalbn(x, n);
if ( fn > 65000.0) return scalbn(x, 65000);
if (-fn > 65000.0) return scalbn(x,-65000);
return scalbn(x,(int)fn);
return scalbn(x, n);
return scalbn(x, n);
return scalbn(x, -ilogb(x));