Fast Gaussian process codes in Fortran