Module 544 in TOMSL2A and L2B: Fortran subroutines for solving weighted least squares problems by modified Gram-Schmidt with iterative refinement. The types of problems solved include overdetermined and underdetermined systems of linear equations, and problems where the solution is subject to linear equality constraints. The covariance matrix of the solution vector is computed. (See R.H. Wampler, ACM TOMS 5 (1979) pp. 494-499.).