s += 0.5 * (fast_logdet(sigma_) + n_samples * log(alpha_) + np.sum(np.log(lambda_)))
s += alpha_1 * log(alpha_) - alpha_2 * alpha_