regression in tgamma?

[Paul Zimmermann]

       Hi again,

I have checked both tgammaf and tgamma with commit b2f3d59 (on x86_64).

tgammaf (exhaustive test of all the 2^32 inputs) no longer has the issue
of tgammaf(-0) = +inf instead of -inf, and similarly for -2^-128 <= x < 0.
The larger ulp error is 239 ulps as in 3.3.0 (for x=-0x1.01743cp+5).

For tgamma, by random testing I get the same maximal error as before:
2244.45 ulps for x=-0x1.53f198fe3b278p+7.

Thus the regression we had in 4.0.0 has gone.

Thanks Jeff,
Paul