CGRQTS(3)
π NAME π
okay so bestie, welcome to cgrqts.f(3), straight out of the Library Functions Manual. This page is literally about the file at:
/home/abuild/rpmbuild/BUILD/lapack-static-3.12.0-build/lapack-3.12.0/TESTING/EIG/cgrqts.f
β¨ SYNOPSIS β¨
Functions/Subroutines:
Functions/Subroutines
subroutine CGRQTS (m, p, n, a, af, q, r, lda, taua, b, bf, z, t, bwk,
ldb, taub, work, lwork, rwork, result)
CGRQTS
π¦ Function/Subroutine Documentation π¦
The signature, exactly as the page gives it:
subroutine CGRQTS (integer m, integer p, integer n, complex, dimension( lda,
* ) a, complex, dimension( lda, * ) af, complex, dimension( lda, * ) q,
complex, dimension( lda, * ) r, integer lda, complex, dimension( * )
taua, complex, dimension( ldb, * ) b, complex, dimension( ldb, * ) bf,
complex, dimension( ldb, * ) z, complex, dimension( ldb, * ) t, complex,
dimension( ldb, * ) bwk, integer ldb, complex, dimension( * ) taub,
complex, dimension( lwork ) work, integer lwork, real, dimension( * )
rwork, real, dimension( 4 ) result)
CGRQTS
πΈ Purpose πΈ
No cap, CGRQTS tests CGGRQF, which computes the GRQ factorization of an M-by-N matrix A and a P-by-N matrix B: A = R*Q and B = Z*T*Q. That's the whole vibe of this subroutine. π
π Parameters π
MβMis INTEGER. It's the number of rows of the matrixA. Constraint:M >= 0.PβPis INTEGER. It's the number of rows of the matrixB. Constraint:P >= 0.NβNis INTEGER. It's the number of columns of the matricesAandB. Constraint:N >= 0.AβAis COMPLEX array, dimension(LDA,N). It's theM-by-NmatrixA.AFβAFis COMPLEX array, dimension(LDA,N). It holds details of the GRQ factorization ofAandB, as returned byCGGRQF; seeCGGRQFfor further details.QβQis COMPLEX array, dimension(LDA,N). It's theN-by-Nunitary matrixQ.RβRis COMPLEX array, dimension(LDA,MAX(M,N)).LDAβLDAis INTEGER. It's the leading dimension of the arraysA,AF,RandQ. Constraint:LDA >= max(M,N).TAUAβTAUAis COMPLEX array, dimension(min(M,N)). These are the scalar factors of the elementary reflectors, as returned bySGGQRC.BβBis COMPLEX array, dimension(LDB,N). On entry, theP-by-NmatrixA. (yes, it saysAright there in the source β I'm just relaying the exact words, bestie π)BFβBFis COMPLEX array, dimension(LDB,N). It holds details of the GQR factorization ofAandB, as returned byCGGRQF; seeCGGRQFfor further details.ZβZis REAL array, dimension(LDB,P). It's theP-by-Punitary matrixZ.TβTis COMPLEX array, dimension(LDB,max(P,N)).BWKβBWKis COMPLEX array, dimension(LDB,N).LDBβLDBis INTEGER. It's the leading dimension of the arraysB,BF,ZandT. Constraint:LDB >= max(P,N).TAUBβTAUBis COMPLEX array, dimension(min(P,N)). These are the scalar factors of the elementary reflectors, as returned bySGGRQF.WORKβWORKis COMPLEX array, dimension(LWORK).LWORKβLWORKis INTEGER. It's the dimension of the arrayWORK. Constraint:LWORK >= max(M,P,N)**2.RWORKβRWORKis REAL array, dimension(M).RESULTβRESULTis REAL array, dimension(4). These are the test ratios:
RESULT(1) = norm( R - A*Q' ) / ( MAX(M,N)*norm(A)*ULP)
RESULT(2) = norm( T*Q - Z'*B ) / (MAX(P,N)*norm(B)*ULP)
RESULT(3) = norm( I - Q'*Q ) / ( N*ULP )
RESULT(4) = norm( I - Z'*Z ) / ( P*ULP )
π Author π
- Univ. of Tennessee
- Univ. of California Berkeley
- Univ. of Colorado Denver
- NAG Ltd.
Definition at line 174 of file cgrqts.f.
π Author π
Generated automatically by Doxygen for LAPACK from the source code.
Cute sign-off, bestie: source LAPACK, version Version 3.12.0, and the footer date line β which in the source page is the path-ish bit /home/abui...IG/cgrqts.f(3) with no separate calendar date shown β mwah π
ππ«Ά
π¦ where this page comes from
- opensuse-leap 16.0: lapack-man 3.12.0-160000.2.2 Β· BSD-3-Clause
- opensuse-leap 16.0: lapack-man 3.12.0-160000.2.2 Β· BSD-3-Clause
- opensuse-leap 16.0: lapack-man 3.12.0-160000.2.2 Β· BSD-3-Clause
- opensuse-leap 16.0: lapack-man 3.12.0-160000.2.2 Β· BSD-3-Clause
- π girlified by deepseek-ai/DeepSeek-V4.1-Flash-1017 on : a language model's rewrite of the original, under its terms