TESTING/LIN/dqrt01.f(3)
πΈ TESTING/LIN/dqrt01.f πΈ
Okay bestie, welcome to the Library Functions Manual page for TESTING/LIN/dqrt01.f β this one is lowkey a whole Fortran test-routine moment. π
β¨
β¨ NAME β¨
TESTING/LIN/dqrt01.f
π SYNOPSIS π
π§ Functions/Subroutines π§
subroutine dqrt01 (m, n, a, af, q, r, lda, tau, work, lwork, rwork,
result)
DQRT01
π Function/Subroutine Documentation π
Here's the full signature, fr fr β every type and dimension stays exactly where it lives:
subroutine dqrt01 (integer m, integer n, 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( * ) tau, double precision, dimension(
lwork ) work, integer lwork, double precision, dimension( * ) rwork,
double precision, dimension( * ) result)
DQRT01
π¦ Purpose π¦
Okay so DQRT01 tests DGEQRF, which computes the QR factorization of an m-by-n matrix A, and partially tests DORGQR which forms the m-by-m orthogonal matrix Q. Iconic behavior. β¨
DQRT01 compares R with Q'*A, and checks that Q is orthogonal. Periodt. π
π Parameters π
MβM is INTEGER; the number of rows of the matrixA.M >= 0.NβN is INTEGER; the number of columns of the matrixA.N >= 0.AβA is DOUBLE PRECISION array, dimension (LDA,N); them-by-nmatrixA.AFβAF is DOUBLE PRECISION array, dimension (LDA,N); details of the QR factorization ofA, as returned byDGEQRF. SeeDGEQRFfor further details.QβQ is DOUBLE PRECISION array, dimension (LDA,M); them-by-morthogonal matrixQ.RβR is DOUBLE PRECISION array, dimension (LDA,max(M,N)).LDAβLDA is INTEGER; the leading dimension of the arraysA,AF,QandR.LDA >= max(M,N).TAUβTAU is DOUBLE PRECISION array, dimension (min(M,N)); the scalar factors of the elementary reflectors, as returned byDGEQRF.WORKβWORK is DOUBLE PRECISION array, dimension (LWORK).LWORKβLWORK is INTEGER; the dimension of the arrayWORK.RWORKβRWORK is DOUBLE PRECISION array, dimension (M).RESULTβRESULT is DOUBLE PRECISION array, dimension (2); the test ratios:RESULT(1) = norm( R - Q'*A ) / ( M * norm(A) * EPS )RESULT(2) = norm( I - Q'*Q ) / ( M * EPS )
Highkey important: those RESULT lines are the vibe check for whether the QR stuff is actually behaving. π
π Author π
- Univ. of Tennessee
- Univ. of California Berkeley
- Univ. of Colorado Denver
- NAG Ltd.
Definition at line 124 of file dqrt01.f.
π Author π
Generated automatically by Doxygen for LAPACK from the source code.
That's the tea, bestie! Source/version sign-off: LAPACK Version 3.12.0; page label: TESTING/LIN/dqrt01.f(3). Slay and goodbye. πβ¨π
π¦ 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