double_t
fn = rnint((double_t)x*invpio2);
double_t hi, lo;
double_t hi, lo;
double_t ax, ax2h, ax2l, axh, axl, ay, ay2h, ay2l, ayh, ayl, sh, sl, t;
double_t y = x;
double_t hi, lo, t;