https://github.com/ekg/freebayes
Tip revision: 4cde5d32ff745018e0b8ec31e81a71e8bc9ecc6f authored by Erik Garrison on 03 June 2019, 13:29:43 UTC
Merge pull request #545 from ekg/skip-coverage
Merge pull request #545 from ekg/skip-coverage
Tip revision: 4cde5d3
factorialln.py
from scipy.special import gamma, gammaln
def factorialln(n):
if n == 1:
return 0
elif n == 0:
return 0
elif n < 0:
raise Exception("factorial is not defined for n < 0")
else:
return gammaln(n + 1)