Sequential Preprocessing of Linear Least-Squares Data
These routines are intended for linear least squares problems where there
are many more rows than columns and one would like to solve the problem
without storing all the rows.
This software is part of the MATH77/mathc90 collection here.