!>
!> DQRT11 computes the test ratio
!>
!>       || Q'*Q - I || / (eps * m)
!>
!> where the orthogonal matrix Q is represented as a product of
!> elementary transformations.  Each transformation has the form
!>
!>    H(k) = I - tau(k) v(k) v(k)'
!>
!> where tau(k) is stored in TAU(k) and v(k) is an m-vector of the form
!> [ 0 ... 0 1 x(k) ]', where x(k) is a vector of length m-k stored
!> in A(k+1:m,k).
!> 

M

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

K

!>          K is INTEGER
!>          The number of columns of A whose subdiagonal entries
!>          contain information about orthogonal transformations.
!> 

A

!>          A is DOUBLE PRECISION array, dimension (LDA,K)
!>          The (possibly partial) output of a QR reduction routine.
!> 

LDA

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

TAU

!>          TAU is DOUBLE PRECISION array, dimension (K)
!>          The scaling factors tau for the elementary transformations as
!>          computed by the QR factorization routine.
!> 

WORK

!>          WORK is DOUBLE PRECISION array, dimension (LWORK)
!> 

LWORK

!>          LWORK is INTEGER
!>          The length of the array WORK.  LWORK >= M*M + M.
!> 

Univ. of Tennessee

Univ. of California Berkeley

Univ. of Colorado Denver

NAG Ltd.