Basic Procedures for Large, Sparse or Weighted Linear Least Squares Problems