Variational Gaussian Processes For Linear Inverse Problems