typo
creates draws from a GP given a kernel function using 2d interpoalion
hyperameter search, fixed kernel likelihood evaluation to use cho decomposition and slogdet() instead of inversion and log(det())