gammaIncomplete(M,M) fails to return when M >= 1e20

[tedgin]

The function std.mathspecial.gammaIncomplete appears to get into a infinite loop for certain arguments. In particular, if M = 1e20, gammaIncomplete(M, M) never returns. When M = 1e10 it does return though.