Variational Inference for Gaussian Process Models with Linear Complexity