Variational Gaussian processes for linear inverse problems