def calc_log_Z(r, s, nu): return ( ((nu + 1.) / 2.) * LOG2 + .5 * LOGPI - .5 * log(r) - (nu/2.) * log(s) + lgamma(nu/2.))