SRC/zgelqf.f(3)
π SRC/zgelqf.f π
Hii bestie π welcome to the SRC/zgelqf.f(3) page from the Library Functions Manual β this is the ZGELQF moment, and yes, weβre doing the whole thing, no skips, no vibes-only summary β¨
πΈ NAME πΈ
SRC/zgelqf.f
β¨ SYNOPSIS β¨
π§ Functions/Subroutines π§
subroutine zgelqf (m, n, a, lda, tau, work, lwork, info)βZGELQF
π«Ά Function/Subroutine Documentation π«Ά
subroutine zgelqf (integer m, integer n, complex*16, dimension( lda, * ) a,
integer lda, complex*16, dimension( * ) tau, complex*16, dimension( * )
work, integer lwork, integer info)
ZGELQF β main character energy, literally π
π Purpose π
Bestie, ZGELQF computes an LQ factorization of a complex M-by-N matrix A:
A = ( L 0 ) * Q
where:
Q is a N-by-N orthogonal matrix;
L is a lower-triangular M-by-M matrix;
0 is a M-by-(N-M) zero matrix, if M < N.
π Parameters π
MβMisINTEGER. Itβs the number of rows of the matrixA. Condition:M >= 0.NβNisINTEGER. Itβs the number of columns of the matrixA. Condition:N >= 0.AβAisCOMPLEX*16array, dimension(LDA,N). On entry, itβs the M-by-N matrixA. On exit, the elements on and below the diagonal of the array contain the m-by-min(m,n) lower trapezoidal matrixL(Lis lower triangular if m <= n); the elements above the diagonal, with the arrayTAU, represent the unitary matrixQas a product of elementary reflectors (see Further Details).LDAβLDAisINTEGER. Itβs the leading dimension of the arrayA. Condition:LDA >= max(1,M).TAUβTAUisCOMPLEX*16array, dimension(min(M,N)). These are the scalar factors of the elementary reflectors (see Further Details).WORKβWORKisCOMPLEX*16array, dimension(MAX(1,LWORK)). On exit, ifINFO = 0,WORK(1)returns the optimalLWORK.LWORKβLWORKisINTEGER. Itβs the dimension of the arrayWORK.LWORK >= 1, ifMIN(M,N) = 0, andLWORK >= M, otherwise. For optimum performanceLWORK >= M*NB, whereNBis the optimal blocksize. IfLWORK = -1, then a workspace query is assumed; the routine only calculates the optimal size of theWORKarray, returns this value as the first entry of theWORKarray, and no error message related toLWORKis issued byXERBLA.INFOβINFOisINTEGER.= 0: successful exit.< 0: ifINFO = -i, the i-th argument had an illegal value.
π Author π
- Univ. of Tennessee
- Univ. of California Berkeley
- Univ. of Colorado Denver
- NAG Ltd.
π₯ Further Details π₯
The matrix Q is represented as a product of elementary reflectors:
Q = H(k)**H . . . H(2)**H H(1)**H, where k = min(m,n).
Each H(i) has the form:
H(i) = I - tau * v * v**H
where tau is a complex scalar, and v is a complex vector with v(1:i-1) = 0 and v(i) = 1; conjg(v(i+1:n)) is stored on exit in A(i,i+1:n), and tau in TAU(i).
Definition at line 143 of file zgelqf.f.
π Author π
Generated automatically by Doxygen for LAPACK from the source code.
Thatβs the whole page, bestie π«Ά source: SRC/zgelqf.f(3) Β· version: LAPACK Version 3.12.0 Β· slay forever πβ¨
π¦ 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