Fast radial basis function interpolation with Gaussians by localization and iteration