Using a design matrix of dimension 512 by 256000 and 256000 observations, we compare 9 different methods of solving the least squares problem.