ln10
return (scalbn(exp(ln10 * ((x - k * logt2hi) - k * logt2lo)), k));
ln10 = 2.302585092994045684017991454684364207601e+0000L,
return (scalbnl(expl(ln10 * ((x - t * logt2hi) -