logt2hi
return (scalbn(exp(ln10 * ((x - k * logt2hi) - k * logt2lo)), k));
logt2hi = LOGT2HI,
return (scalbnl(expl(ln10 * ((x - t * logt2hi) -