The previous code did not take into account the 2nd term of the Taylor
expansion (t^3/3) which is only t^2/2 relatively to the main term, thus
for y/x near 2^-27 it gave a contribution of the order of 1 ulp.