!>
!> Solve the least squares problem
!>     min || A*X - B ||
!> using the QL factorization
!>     A = Q*L
!> computed by ZGEQLF.
!> 

M

!>          M is INTEGER
!>          The number of rows of the matrix A.  M >= 0.
!> 

N

!>          N is INTEGER
!>          The number of columns of the matrix A.  M >= N >= 0.
!> 

NRHS

!>          NRHS is INTEGER
!>          The number of columns of B.  NRHS >= 0.
!> 

A

!>          A is COMPLEX*16 array, dimension (LDA,N)
!>          Details of the QL factorization of the original matrix A as
!>          returned by ZGEQLF.
!> 

LDA

!>          LDA is INTEGER
!>          The leading dimension of the array A.  LDA >= M.
!> 

TAU

!>          TAU is COMPLEX*16 array, dimension (N)
!>          Details of the orthogonal matrix Q.
!> 

B

!>          B is COMPLEX*16 array, dimension (LDB,NRHS)
!>          On entry, the m-by-nrhs right hand side matrix B.
!>          On exit, the n-by-nrhs solution matrix X, stored in rows
!>          m-n+1:m.
!> 

LDB

!>          LDB is INTEGER
!>          The leading dimension of the array B. LDB >= M.
!> 

WORK

!>          WORK is COMPLEX*16 array, dimension (LWORK)
!> 

LWORK

!>          LWORK is INTEGER
!>          The length of the array WORK.  LWORK must be at least NRHS,
!>          and should be at least NRHS*NB, where NB is the block size
!>          for this environment.
!> 

INFO

!>          INFO is INTEGER
!>          = 0: successful exit
!>          < 0: if INFO = -i, the i-th argument had an illegal value
!> 

Univ. of Tennessee

Univ. of California Berkeley

Univ. of Colorado Denver

NAG Ltd.