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