Massively parallel approximate Gaussian process regression