Looks like we need an actual gamma function? I think gammaln is relatively slow.