def METHOD_NAME(x0):
    """Multivariate normal log-likelihood."""
    x = x0[:ndim_gau]
    return -0.5 * np.dot((x - mean_gau), np.dot(cov_inv_gau,
                                                (x - mean_gau))) + lnorm_gau