It has no remaining callers (kept only for BC) so nothing else exercises it. Values checked against mpmath's gammainc(a, 0, x); the deep-tail case (a=4, x=80) reproduces the convergence bug this PR fixes.