SRC/zgerqf.f(3)
πΈ NAME πΈ
SRC/zgerqf.f β okay bestie, this is the manual page for that file, and the manual-name energy is SRC/zgerqf.f(3) from the Library Functions Manual π
β¨ SYNOPSIS β¨
Functions/Subroutines
subroutine zgerqf (m, n, a, lda, tau, work, lwork, info)
ZGERQF
π Function/Subroutine Documentation π
subroutine zgerqf (integer m, integer n, complex*16, dimension( lda, * ) a,
integer lda, complex*16, dimension( * ) tau, complex*16, dimension( * )
work, integer lwork, integer info)
π ZGERQF π
β¨ Purpose β¨
The way that ZGERQF computes an RQ factorization of a complex M-by-N matrix A: A = R * Q. That's the whole mission, no cap.
π¦ 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 COMPLEX*16 array, dimension (LDA,N). On entry, theM-by-NmatrixA. On exit:- if
m <= n, the upper triangle of the subarrayA(1:m,n-m+1:n)contains theM-by-Mupper triangular matrixR; - if
m >= n, the elements on and above the(m-n)-th subdiagonal contain theM-by-Nupper trapezoidal matrixR; - the remaining elements, with the array
TAU, represent the unitary matrixQas a product ofmin(m,n)elementary reflectors (see Further Details).
- if
LDAβLDA is INTEGER. The leading dimension of the arrayA.LDA >= max(1,M).TAUβTAU is COMPLEX*16 array, dimension (min(M,N)). The scalar factors of the elementary reflectors (see Further Details).WORKβWORK is COMPLEX*16 array, dimension (MAX(1,LWORK)). On exit, ifINFO = 0,WORK(1)returns the optimalLWORK.LWORKβLWORK is INTEGER. The dimension of the arrayWORK.LWORK >= 1, if MIN(M,N) = 0, and LWORK >= 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βINFO is INTEGER:= 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 π
The matrix Q is represented as a product of elementary reflectors:
Q = H(1)**H H(2)**H . . . H(k)**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(n-k+i+1:n) = 0 and v(n-k+i) = 1; conjg(v(1:n-k+i-1)) is stored on exit in A(m-k+i,1:n-k+i-1), and tau in TAU(i).
Definition at line 138 of file zgerqf.f.
π Author π
Generated automatically by Doxygen for LAPACK from the source code.
k bye bestie! π Footer sign-off: LAPACK Version 3.12.0, page SRC/zgerqf.f(3) β that's the source/version line, and there's no extra date shown in the footer, so we're not inventing one π
π¦ 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