Scaling up Gauss-Newton methods for expensive least squares problems