Iterative Methods for Full-Scale Gaussian Process Approximations for Large Spatial Data