MEP1
return (scalbnl(expl((w - MEP1 * ln2hi) - MEP1 * ln2lo), ME));
r = copysignl(expl((r - MEP1 * ln2hi) - MEP1 * ln2lo), x);