CGEQRT(3)
π
CGEQRT(3) β Library Functions Manual π
okay so bestie, welcome to the group chat for CGEQRT(3). The manual name is Library Functions Manual, and we are keeping every single technical crumb while making it cute. Periodt. π
πΈ NAME πΈ
/home/abuild/rpmbuild/BUILD/lapack-static-3.12.0-build/lapack-3.12.0/SRC/cgeqrt.f
no cap, the NAME section is literally just this file path. It's giving "I have one job" energy. π
π SYNOPSIS π
Functions/Subroutines
subroutine CGEQRT (m, n, nb, a, lda, t, ldt, work, info)
CGEQRT
π Function/Subroutine Documentation π
π subroutine CGEQRT π
subroutine CGEQRT (integer m, integer n, integer nb, complex, dimension(
lda, * ) a, integer lda, complex, dimension( ldt, * ) t, integer ldt,
complex, dimension( * ) work, integer info)
CGEQRT
π Purpose π
CGEQRT computes a blocked QR factorization of a complex M-by-N matrix A using the compact WY representation of Q. slay. β¨
π§Έ Parameters π§Έ
MβM is INTEGER; it's the number of rows of matrix A, andM >= 0.NβN is INTEGER; it's the number of columns of matrix A, andN >= 0.NBβNB is INTEGER; this is the block size for the blocked QR, withMIN(M,N) >= NB >= 1.AβA is COMPLEX array, dimension (LDA,N). On entry, it's the M-by-N matrix A. On exit, the elements on and above the diagonal of the array contain themin(M,N)-by-Nupper trapezoidal matrix R (R is upper triangular ifM >= N); the elements below the diagonal are the columns of V.LDAβLDA is INTEGER; the leading dimension of array A, andLDA >= max(1,M).TβT is COMPLEX array, dimension (LDT,MIN(M,N)). It holds the upper triangular block reflectors in compact form, as a sequence of upper triangular blocks. See below for the deets.LDTβLDT is INTEGER; the leading dimension of array T, andLDT >= NB.WORKβWORK is COMPLEX array, dimension (NB*N).INFOβINFO is INTEGER= 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 V stores the elementary reflectors H(i) in the i-th column below the diagonal. For example, if M=5 and N=3, the matrix V is
V = ( 1 )
( v1 1 )
( v1 v2 1 )
( v1 v2 v3 )
( v1 v2 v3 )
where the vi's represent the vectors which define H(i), which are returned in the matrix A. The 1's along the diagonal of V are not stored in A.
Let K=MIN(M,N). The number of blocks is B = ceiling(K/NB), where each block is of order NB except for the last block, which is of order IB = K - (B-1)*NB. For each of the B blocks, an upper triangular block reflector factor is computed: T1, T2, ..., TB. The NB-by-NB (and IB-by-IB for the last block) T's are stored in the NB-by-K matrix T as
T = (T1 T2 ... TB).
Definition at line 140 of file cgeqrt.f.
π€ Author π€
Generated automatically by Doxygen for LAPACK from the source code.
π Sign-off π
Footer tea, since we're not copying the footer line: the source/version is LAPACK and Version 3.12.0, and the date/path field is /home/abui...RC/cgeqrt.f(3). Mwah bestie, go slay that blocked QR. πβ¨π
π¦ where this page comes from
- opensuse-leap 16.0: lapack-man 3.12.0-160000.2.2 Β· BSD-3-Clause
- opensuse-leap 16.0: lapack-man 3.12.0-160000.2.2 Β· BSD-3-Clause
- opensuse-leap 16.0: lapack-man 3.12.0-160000.2.2 Β· BSD-3-Clause
- opensuse-leap 16.0: lapack-man 3.12.0-160000.2.2 Β· BSD-3-Clause
- π girlified by deepseek-ai/DeepSeek-V4.1-Flash-1017 on : a language model's rewrite of the original, under its terms