![]() |
LAPACK
3.7.0
LAPACK: Linear Algebra PACKage
|
subroutine dgqrts | ( | integer | N, |
integer | M, | ||
integer | P, | ||
double precision, dimension( lda, * ) | A, | ||
double precision, dimension( lda, * ) | AF, | ||
double precision, dimension( lda, * ) | Q, | ||
double precision, dimension( lda, * ) | R, | ||
integer | LDA, | ||
double precision, dimension( * ) | TAUA, | ||
double precision, dimension( ldb, * ) | B, | ||
double precision, dimension( ldb, * ) | BF, | ||
double precision, dimension( ldb, * ) | Z, | ||
double precision, dimension( ldb, * ) | T, | ||
double precision, dimension( ldb, * ) | BWK, | ||
integer | LDB, | ||
double precision, dimension( * ) | TAUB, | ||
double precision, dimension( lwork ) | WORK, | ||
integer | LWORK, | ||
double precision, dimension( * ) | RWORK, | ||
double precision, dimension( 4 ) | RESULT | ||
) |
DGQRTS
DGQRTS tests DGGQRF, which computes the GQR factorization of an N-by-M matrix A and a N-by-P matrix B: A = Q*R and B = Q*T*Z.
[in] | N | N is INTEGER The number of rows of the matrices A and B. N >= 0. |
[in] | M | M is INTEGER The number of columns of the matrix A. M >= 0. |
[in] | P | P is INTEGER The number of columns of the matrix B. P >= 0. |
[in] | A | A is DOUBLE PRECISION array, dimension (LDA,M) The N-by-M matrix A. |
[out] | AF | AF is DOUBLE PRECISION array, dimension (LDA,N) Details of the GQR factorization of A and B, as returned by DGGQRF, see SGGQRF for further details. |
[out] | Q | Q is DOUBLE PRECISION array, dimension (LDA,N) The M-by-M orthogonal matrix Q. |
[out] | R | R is DOUBLE PRECISION array, dimension (LDA,MAX(M,N)) |
[in] | LDA | LDA is INTEGER The leading dimension of the arrays A, AF, R and Q. LDA >= max(M,N). |
[out] | TAUA | TAUA is DOUBLE PRECISION array, dimension (min(M,N)) The scalar factors of the elementary reflectors, as returned by DGGQRF. |
[in] | B | B is DOUBLE PRECISION array, dimension (LDB,P) On entry, the N-by-P matrix A. |
[out] | BF | BF is DOUBLE PRECISION array, dimension (LDB,N) Details of the GQR factorization of A and B, as returned by DGGQRF, see SGGQRF for further details. |
[out] | Z | Z is DOUBLE PRECISION array, dimension (LDB,P) The P-by-P orthogonal matrix Z. |
[out] | T | T is DOUBLE PRECISION array, dimension (LDB,max(P,N)) |
[out] | BWK | BWK is DOUBLE PRECISION array, dimension (LDB,N) |
[in] | LDB | LDB is INTEGER The leading dimension of the arrays B, BF, Z and T. LDB >= max(P,N). |
[out] | TAUB | TAUB is DOUBLE PRECISION array, dimension (min(P,N)) The scalar factors of the elementary reflectors, as returned by DGGRQF. |
[out] | WORK | WORK is DOUBLE PRECISION array, dimension (LWORK) |
[in] | LWORK | LWORK is INTEGER The dimension of the array WORK, LWORK >= max(N,M,P)**2. |
[out] | RWORK | RWORK is DOUBLE PRECISION array, dimension (max(N,M,P)) |
[out] | RESULT | RESULT is DOUBLE PRECISION array, dimension (4) The test ratios: RESULT(1) = norm( R - Q'*A ) / ( MAX(M,N)*norm(A)*ULP) RESULT(2) = norm( T*Z - Q'*B ) / (MAX(P,N)*norm(B)*ULP) RESULT(3) = norm( I - Q'*Q ) / ( M*ULP ) RESULT(4) = norm( I - Z'*Z ) / ( P*ULP ) |
Definition at line 178 of file dgqrts.f.