FORTRAN subroutines for updating the QR decomposition