float_t
if (sizeof(float_t) > sizeof(float))
return (r - hfsq + f) * ((float_t)ivln10lo + ivln10hi) +
y * ((float_t)log10_2lo + log10_2hi);
if (sizeof(float_t) > sizeof(float))
return (r - hfsq + f) * ((float_t)ivln2lo + ivln2hi) + y;
fn = rnint((float_t)x*invpio2);
float_t ax, ax2h, ax2l, axh, axl, ay, ay2h, ay2l, ayh, ayl, sh, sl, t;
float_t y = x;