Bowei "Bobbie" Wu

I agree with @putianyi889 The reciprocal function `1/gamma` frequently appears in applications. Since `1/gamma` is an entire function with simple zeros, it makes sense that `gamma` returns `Inf` at its...

An update: I notice that `exp(loggamma(n))` gives `Inf` for any negative integer `n`, is there a reason why this behavior is not used to compute `gamma(n) == exp(loggamma(n))`? Also, I...

Hi @stevengj, thanks for the explanation. I thought that all the negative inputs to `loggamma` would be promoted to complex, but apparently they are not, i.e. I have to manually...

Hi @stevengj, I see that type stability is important. But back to the original question of giving a value to `gamma(-n)`, I think defining `gamma(-n)=Inf` is good because: 1. It...