ZGRQTS(3)
π ZGRQTS(3) π
Okay bestie, welcome to the ZGRQTS(3) page from the Library Functions Manual. The header/footer shorten the file to /home/abui...IG/zgrqts.f(3), and the full path is waiting in NAME below. This is LAPACK test-subroutine territory, and yes, it's giving "I test the factorization so you don't have to" energy. π¦
π NAME π
/home/abuild/rpmbuild/BUILD/lapack-static-3.12.0-build/lapack-3.12.0/TESTING/EIG/zgrqts.f
That's the whole NAME, no cap. The file is zgrqts.f, and the short header/footer path is /home/abui...IG/zgrqts.f(3). π
β¨ SYNOPSIS β¨
π§ Functions/Subroutines π§
subroutine ZGRQTS (m, p, n, a, af, q, r, lda, taua, b, bf, z, t, bwk,
ldb, taub, work, lwork, rwork, result)
ZGRQTS
π¦ Function/Subroutine Documentation π¦
subroutine ZGRQTS (integer m, integer p, 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( * ) taua, complex*16, dimension( ldb, * ) b, complex*16,
dimension( ldb, * ) bf, complex*16, dimension( ldb, * ) z, complex*16,
dimension( ldb, * ) t, complex*16, dimension( ldb, * ) bwk, integer ldb,
complex*16, dimension( * ) taub, complex*16, dimension( lwork ) work,
integer lwork, double precision, dimension( * ) rwork, double precision,
dimension( 4 ) result)
ZGRQTS
π Purpose π
ZGRQTS is the test bestie for ZGGRQF β the routine that computes the GRQ factorization of an M-by-N matrix A and a P-by-N matrix B, specifically:
A = R*Q and B = Z*T*Q.
It's literally giving "I check the factorization" main-character energy. π
π Parameters π
MβMisINTEGER; the number of rows of the matrixA. Constraint:M >= 0.PβPisINTEGER; the number of rows of the matrixB. Constraint:P >= 0.NβNisINTEGER; the number of columns of the matricesAandB. Constraint:N >= 0.AβAisCOMPLEX*16array, dimension(LDA,N); theM-by-NmatrixA.AFβAFisCOMPLEX*16array, dimension(LDA,N); details of the GRQ factorization ofAandB, as returned byZGGRQF, seeCGGRQFfor further details.QβQisCOMPLEX*16array, dimension(LDA,N); theN-by-Nunitary matrixQ.RβRisCOMPLEX*16array, dimension(LDA,MAX(M,N)).LDAβLDAisINTEGER; the leading dimension of the arraysA,AF,RandQ. Constraint:LDA >= max(M,N).TAUAβTAUAisCOMPLEX*16array, dimension(min(M,N)); the scalar factors of the elementary reflectors, as returned byDGGQRC.BβBisCOMPLEX*16array, dimension(LDB,N); on entry, theP-by-NmatrixA(yes, the source saysAhere, bestie β keeping it exactly as written).BFβBFisCOMPLEX*16array, dimension(LDB,N); details of the GQR factorization ofAandB, as returned byZGGRQF, seeCGGRQFfor further details.ZβZisDOUBLE PRECISIONarray, dimension(LDB,P); theP-by-Punitary matrixZ(the parameter text saysDOUBLE PRECISION, so I'm keeping that exactly).TβTisCOMPLEX*16array, dimension(LDB,max(P,N)).BWKβBWKisCOMPLEX*16array, dimension(LDB,N).LDBβLDBisINTEGER; the leading dimension of the arraysB,BF,ZandT. Constraint:LDB >= max(P,N).TAUBβTAUBisCOMPLEX*16array, dimension(min(P,N)); the scalar factors of the elementary reflectors, as returned byDGGRQF.WORKβWORKisCOMPLEX*16array, dimension(LWORK).LWORKβLWORKisINTEGER; the dimension of the arrayWORK, withLWORK >= max(M,P,N)**2.RWORKβRWORKisDOUBLE PRECISIONarray, dimension(M).RESULTβRESULTisDOUBLE PRECISIONarray, dimension(4). The test ratios are:
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 zgrqts.f.
π Author π
Generated automatically by Doxygen for LAPACK from the source code.
π That's the whole page, bestie! Footer check: LAPACK, Version 3.12.0, and the page path /home/abui...IG/zgrqts.f(3). There's no separate date line in this one, so no invented date β just the source/version vibe. 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