Fast Ip solution of large, sparse, linear systems: Application to seismic travel time tomography