Is there such a method in Julia? Right now I'm multiplying both sides by the transpose of the coefficient matrix to form and solve the normal equations using \. I suspect that this is not very efficient compared to a sparse least squares solver.
Thanks, --Peter
