Scaling Gaussian Processes with Derivative Information Using Variational Inference