There are actually specialized implementations for these cases (in gammainc()). Why not open these private functions? For instance, in SymPy they could be added to the MPMATH_TRANSLATIONS map (instead of using custom _eval_evalf() helpers).