Member PostTheta::eval_log_prior_lat (Vect &theta, Vect &mu, double &val): construct spatial matrix (at the moment this is happening twice. FIX)
Member PostTheta::hess_eval (Vect theta, double eps): not yet parallelised ....