Numerical stability of inverse Lévy tail mass for gamma process

The inverse transform returns negative values for large inputs when the rate parameter is small. E.g.

```julia
p = GammaProcess(2.1, 0.5)
inverse_levy_tail_mass(p, 6.0)
#> -0.06249282411901989
```