The current code calculates the Hessian contribution of the spatial smoother, and this seems to work fine. Need to do a better job with the Hessian calculations of the ImageGP part, but that's work for when i grow up...