handle zero quickly in get_d
refactor stirling param selection and make more robust for large arguments
add harmonic sums for digamma function argument reduction
fix unnecessary recomputation of bernoulli numbers
work to support the digamma function