Started work on fast least squares with LSQR