breeze.numerics.lgamma
log Incomplete gamma function = \log \int_0x \exp(-t)pow(t,a-1) dt May store lgamma(a) in lgam(0) if it's non-null and needs to be computed. Based on NR
log Incomplete gamma function = \log \int_0x \exp(-t)pow(t,a-1) dt May store lgamma(a) in lgam(0) if it's non-null and needs to be computed. Based on NR