computing the laplace approximation of the model evidence

[ToennisStef]

I wanted to ask how one would compute the model evidence as given here:

in particular the $\Sigma^{-1}$ part. Computing the hessian is not so straightforward to me. Is there already an implementation for this?

To do this one would need to register priors for every hyperparameter of the model first, right?