SRC/sgeql2.f(3)
๐ SRC/sgeql2.f(3) โ Library Functions Manual ๐
okay so welcome to the SRC/sgeql2.f(3) page from the Library Functions Manual, bestie. This is the manual entry for the sgeql2 subroutine, and yes, we're gonna get into the QL factorization tea. โจ
๐ธ NAME ๐ธ
SRC/sgeql2.f โ that's the page name, no cap. ๐ซถ
๐ SYNOPSIS ๐
โจ Functions/Subroutines โจ
subroutine sgeql2 (m, n, a, lda, tau, work, info)โ SGEQL2 computes the QL factorization of a general rectangular matrix using an unblocked algorithm. ๐ฅ
๐ฆ Function/Subroutine Documentation ๐ฆ
๐ subroutine sgeql2 (integer m, integer n, real, dimension( lda, * ) a, integer lda, real, dimension( * ) tau, real, dimension( * ) work, integer info) ๐
SGEQL2 computes the QL factorization of a general rectangular matrix using an unblocked algorithm. That's the whole vibe of this routine. ๐
๐ง Purpose ๐ง
The source comment says it plain: SGEQL2 computes a QL factorization of a real m by n matrix A:
!>
!> SGEQL2 computes a QL factorization of a real m by n matrix A:
!> A = Q * L.
!>
๐ Parameters ๐
MโMisINTEGER. The number of rows of the matrixA.M >= 0. ๐NโNisINTEGER. The number of columns of the matrixA.N >= 0. ๐AโAisREAL 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โLDAisINTEGER. The leading dimension of the arrayA.LDA >= max(1,M). ๐งฎTAUโTAUisREAL array, dimension (min(M,N)). The scalar factors of the elementary reflectors (see Further Details). ๐WORKโWORKisREAL array, dimension (N). ๐ ๏ธINFOโINFOisINTEGER.= 0: successful exit.< 0: ifINFO = -i, thei-th argument had an illegal value. ๐จ
๐ Author ๐
- Univ. of Tennessee
- Univ. of California Berkeley
- Univ. of Colorado Denver
- NAG Ltd.
๐ Further Details ๐
Here's the exact source-comment block for the reflector math, because we do not paraphrase equations in this group chat:
!>
!> 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 sgeql2.f. ๐
๐ Author ๐
Generated automatically by Doxygen for LAPACK from the source code. ๐ค๐
Cute sign-off, bestie: source/version is LAPACK Version 3.12.0, the page is SRC/sgeql2.f(3), and the footer shows no date, so I'm not inventing one โ periodt. ๐
โจ
๐ฆ 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