W. Morven Gentleman, Algorithm AS 75: Basic Procedures for Large, Sparse or Weighted Linear Least Problems, Journal of the Royal Statistical Society. Series C (Applied Statistics), Vol. 23, No. 3 ...