Scaling Gaussian Process Regression with Derivatives