TESTING/LIN/dqlt03.f(3)
π NAME π
TESTING/LIN/dqlt03.f
okay so bestie, this is the TESTING/LIN/dqlt03.f(3) page from the Library Functions Manual β the one and only DQLT03. Itβs giving βFortran test routine thatβs lowkey the momentβ π
.
β¨ SYNOPSIS β¨
Functions/Subroutines
Functions/Subroutines
subroutine dqlt03 (m, n, k, af, c, cc, q, lda, tau, work, lwork, rwork,
result)
DQLT03
π¦ FUNCTION/SUBROUTINE DOCUMENTATION π¦
subroutine dqlt03 (integer m, integer n, integer k, double precision,
dimension( lda, * ) af, double precision, dimension( lda, * ) c, double
precision, dimension( lda, * ) cc, double precision, dimension( lda, * )
q, integer lda, double precision, dimension( * ) tau, double precision,
dimension( lwork ) work, integer lwork, double precision, dimension( * )
rwork, double precision, dimension( * ) result)
DQLT03
π PURPOSE π
DQLT03 tests DORMQL, which computes Q*C, Q'*C, C*Q or C*Q'.
DQLT03 compares the results of a call to DORMQL with the results of forming Q explicitly by a call to DORGQL and then performing matrix multiplication by a call to DGEMM.
πΈ PARAMETERS πΈ
MβM is INTEGER; itβs the order of the orthogonal matrixQ.M >= 0.NβN is INTEGER; the number of rows or columns of the matrixC;Cis m-by-n ifQis applied from the left, or n-by-m ifQis applied from the right.N >= 0.KβK is INTEGER; the number of elementary reflectors whose product defines the orthogonal matrixQ.M >= K >= 0.AFβAF is DOUBLE PRECISION array, dimension (LDA,N); details of the QL factorization of an m-by-n matrix, as returned byDGEQLF. SeeSGEQLFfor further details.CβC is DOUBLE PRECISION array, dimension (LDA,N)β thatβs the whole description, no cap.CCβCC is DOUBLE PRECISION array, dimension (LDA,N)β same energy asC, bestie.QβQ is DOUBLE PRECISION array, dimension (LDA,M).LDAβLDA is INTEGER; the leading dimension of the arraysAF,C,CC, andQ.TAUβTAU is DOUBLE PRECISION array, dimension (min(M,N)); the scalar factors of the elementary reflectors corresponding to the QL factorization inAF.WORKβWORK is DOUBLE PRECISION array, dimension (LWORK).LWORKβLWORK is INTEGER; the length ofWORK.LWORKmust be at leastM, and should beM*NB, whereNBis the blocksize for this environment.RWORKβRWORK is DOUBLE PRECISION array, dimension (M).RESULTβRESULT is DOUBLE PRECISION array, dimension (4); the test ratios compare two techniques for multiplying a random matrixCby an m-by-m orthogonal matrixQ.
RESULT(1) = norm( Q*C - Q*C ) / ( M * norm(C) * EPS )
RESULT(2) = norm( C*Q - C*Q ) / ( M * norm(C) * EPS )
RESULT(3) = norm( Q'*C - Q'*C )/ ( M * norm(C) * EPS )
RESULT(4) = norm( C*Q' - C*Q' )/ ( M * norm(C) * EPS )
π AUTHOR π
- Univ. of Tennessee
- Univ. of California Berkeley
- Univ. of Colorado Denver
- NAG Ltd.
Definition at line 134 of file dqlt03.f.
ποΈ AUTHOR ποΈ
Generated automatically by Doxygen for LAPACK from the source code.
π thatβs the whole page, bestie β LAPACK Version 3.12.0, source/version vibes from TESTING/LIN/dqlt03.f(3), and a little periodt for the road. π«Ά
π¦ where this page comes from
- arch rolling: lapack-doc 3.12.1-2 Β· custom
- π girlified by deepseek-ai/DeepSeek-V4.1-Flash-1017 on : a language model's rewrite of the original, under its terms