Sequential Solution of a Banded Least-Squares Problem
These routines are intended for banded 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.