def __init__(self):
self._mode_finding_tolerance = <NUM_LIT><EOL>self._mode_finding_max_iter = <NUM_LIT:30><EOL>self.bad_fhat = False<EOL>self.first_run = True<EOL>self._previous_Ki_fhat = None<EOL>
Laplace Approximation

Find the moments \hat{f} and the hessian at this point
(using Newton-Raphson) of the unnormalised posterior