lnovlo
return (scalbnl(expl((w - lnovft) - lnovlo), 16383));
r = copysignl(expl((r - lnovft) - lnovlo), x);