Iterative Methods for Vecchia-Laplace Approximations for Latent Gaussian Process Models