zrqt01()

subroutine zrqt01	(	integer	M,
		integer	N,
		complex*16, dimension( lda, * )	A,
		complex*16, dimension( lda, * )	AF,
		complex*16, dimension( lda, * )	Q,
		complex*16, dimension( lda, * )	R,
		integer	LDA,
		complex*16, dimension( * )	TAU,
		complex*16, dimension( lwork )	WORK,
		integer	LWORK,
		double precision, dimension( * )	RWORK,
		double precision, dimension( * )	RESULT
	)

ZRQT01

Purpose:

 ZRQT01 tests ZGERQF, which computes the RQ factorization of an m-by-n
 matrix A, and partially tests ZUNGRQ which forms the n-by-n
 orthogonal matrix Q.

 ZRQT01 compares R with A*Q', and checks that Q is orthogonal.

Parameters

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

Author

Univ. of Tennessee

Univ. of California Berkeley

Univ. of Colorado Denver

NAG Ltd.