__ldexp_cexp
double complex __ldexp_cexp(double complex,int);
z = __ldexp_cexp(CMPLX(fabs(x), y), -1);
return (__ldexp_cexp(z, 0));