def calc_log_Z(a, b, V_inv): # Equation 19. return gammaln(a) + log(sqrt(1./det(V_inv))) - a * np.log(b)