powered by
likelihood function
log_marginal_likelihood(k, x.tilde_m, y.tilde, D_inv, ydinvy, dinvy, g)
likelihood