Gaussian processes in C++