Scalable Inference for Gaussian Process Models with Black-Box Likelihoods