Paul Zimmermann <Paul.Zimmermann@inria.fr> writes:

> thank you Jeff for your answer.
>
> In return I found I believe an issue with newlib's tgammaf function, for
> input -0, where it gives +inf instead of -inf:

Yeah, it's easy to fix the underlying __ieee754_lgamma functions to
adjust the sign of the return: