pints
pints copied to clipboard
Add evaluateS1 method for HighDimensionalGaussianLogPDF
Why? Because it'll be very useful for functionally testing gradient-based methods. I think this can just be copied straight over from that for GaussianLogPDF
.
Whilst we're at it, it might be good to update the docstrings as they're not very informative.
Also, the sigma_inv
calculation in the method in GaussianLogPDF
can be cached, which can be addressed in same PR, I reckon.
Also, add note about result being from chapter 2 of matrix cookbook (plus maybe https://stats.stackexchange.com/questions/27436/how-to-take-derivative-of-multivariate-normal-density).