I think this comment is a bit misleading with the code since our implementation of `incomplete_gamma` already includes the division by `gamma(k/2)`.