Adaptive Gaussian process surrogates for Bayesian inference