!>
!> Compute a minimum-norm solution
!>     min || A*X - B ||
!> using the LQ factorization
!>     A = L*Q
!> computed by ZGELQF.
!> 

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.  N >= M >= 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 LQ factorization of the original matrix A as
!>          returned by ZGELQF.
!> 

LDA

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

TAU

!>          TAU is COMPLEX*16 array, dimension (M)
!>          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.
!> 

LDB

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

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.