SRC/dgeql2.f(3)
πΈ NAME πΈ
okay bestie, youβre looking at SRC/dgeql2.f from the Library Functions Manual β manual name folded in near the top, as requested, and yes weβre being cute and correct. This page is all about dgeql2, a LAPACK subroutine situation. π
β¨
β¨ SYNOPSIS β¨
Functions/Subroutines
subroutine dgeql2 (m, n, a, lda, tau, work, info)β DGEQL2 computes the QL factorization of a general rectangular matrix using an unblocked algorithm.
π Function/Subroutine Documentation π
Hereβs the full signature, copied exactly because we do not play with code:
subroutine dgeql2 (integer m, integer n, double precision, dimension( lda, *
) a, integer lda, double precision, dimension( * ) tau, double
precision, dimension( * ) work, integer info)
and the vibe is: DGEQL2 computes the QL factorization of a general rectangular matrix using an unblocked algorithm. No cap. π¦
π Purpose π
fr fr, hereβs the purpose tea: DGEQL2 computes a QL factorization of a real m by n matrix A:
A = Q * L.
Thatβs it, thatβs the whole purpose. β¨π
π«Ά Parameters π«Ά
Okay, parameter parade, bestie β every single one is here:
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). On entry, thembynmatrixA. On exit, ifm >= n, the lower triangle of the subarrayA(m-n+1:m,1:n)contains thenbynlower triangular matrixL; ifm <= n, the elements on and below the(n-m)-th superdiagonal contain thembynlower trapezoidal matrixL; the remaining elements, with the arrayTAU, represent the orthogonal matrixQas a product of elementary reflectors (see Further Details).LDAβLDA is INTEGER. The leading dimension of the arrayA.LDA >= max(1,M).TAUβTAU is DOUBLE PRECISION array, dimension (min(M,N)). The scalar factors of the elementary reflectors (see Further Details).WORKβWORK is DOUBLE PRECISION array, dimension (N).INFOβINFO is INTEGER.= 0: successful exit.< 0: ifINFO = -i, thei-th argument had an illegal value.
π§ Author π§
Shoutout to the author squad:
- Univ. of Tennessee
- Univ. of California Berkeley
- Univ. of Colorado Denver
- NAG Ltd.
π Further Details π
Okay, this is the details section and itβs lowkey iconic:
The matrix Q is represented as a product of elementary reflectors:
Q = H(k) . . . H(2) H(1), where k = min(m,n).
Each H(i) has the form:
H(i) = I - tau * v * v**T
where tau is a real scalar, and v is a real vector with v(m-k+i+1:m) = 0 and v(m-k+i) = 1; v(1:m-k+i-1) is stored on exit in A(1:m-k+i-1,n-k+i), and tau in TAU(i).
Definition at line 122 of file dgeql2.f.
π§ Author (Doxygenβs Version) π§
Generated automatically by Doxygen for LAPACK from the source code.
Bye bestie πβ¨ β this masterpiece is from LAPACK Version 3.12.0, page SRC/dgeql2.f(3). Thatβs the footer tea, and weβre out. π«Άπ
π¦ 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