Fast Gaussian process inference by exact Matérn kernel decomposition