DQLT01

   Univ. of Tennessee, Univ. of California Berkeley and NAG Ltd..
   November 2006

Purpose

DQLT01 tests DGEQLF, which computes the QL factorization of an m-by-n
matrix A, and partially tests DORGQL which forms the m-by-m
orthogonal matrix Q.

DQLT01 compares L with Q'*A, and checks that Q is orthogonal.

Arguments

M
(input) INTEGER
The number of rows of the matrix A.  M >= 0.
N
(input) INTEGER
The number of columns of the matrix A.  N >= 0.
A
(input) DOUBLE PRECISION array, dimension (LDA,N)
The m-by-n matrix A.
AF
(output) DOUBLE PRECISION array, dimension (LDA,N)
Details of the QL factorization of A, as returned by DGEQLF.
See DGEQLF for further details.
Q
(output) DOUBLE PRECISION array, dimension (LDA,M)
The m-by-m orthogonal matrix Q.
L
(workspace) DOUBLE PRECISION array, dimension (LDA,max(M,N))
LDA
(input) INTEGER
The leading dimension of the arrays A, AF, Q and R.
LDA >= max(M,N).
TAU
(output) DOUBLE PRECISION array, dimension (min(M,N))
The scalar factors of the elementary reflectors, as returned
by DGEQLF.
WORK
(workspace) DOUBLE PRECISION array, dimension (LWORK)
LWORK
(input) INTEGER
The dimension of the array WORK.
RWORK
(workspace) DOUBLE PRECISION array, dimension (M)
RESULT
(output) DOUBLE PRECISION array, dimension (2)
The test ratios:
RESULT(1) = norm( L - Q'*A ) / ( M * norm(A) * EPS )
RESULT(2) = norm( I - Q'*Q ) / ( M * EPS )

Call Graph

Caller Graph