frexpl
extern long double frexpl(long double num, int *_exponent);
xs = frexpl(x, &ex);
ys = frexpl(y, &ey);
zs = frexpl(z, &ez);
x = frexpl(x*0x1p120, e);
x = frexpl(x, &e);
x = frexpl(x, &e);
x = frexpl(x, &e);
x = frexpl(x, &e);
x = frexpl(x, &i);
s = frexpl( s, &lx);