SRC/zgeqrfp.f(3)
π
SRC/zgeqrfp.f π
okay so bestie, this is the SRC/zgeqrfp.f(3) page from the Library Functions Manual β¨. Itβs giving very βLAPACK routine, but make it group chat.β We are not dropping a single fact. Periodt.
πΈ NAME πΈ
SRC/zgeqrfp.f β thatβs the name of the page/file, and yes, itβs short. Iconic.
β¨ SYNOPSIS β¨
π§ Functions/Subroutines π§
subroutine zgeqrfp (m, n, a, lda, tau, work, lwork, info)βZGEQRFP.
subroutine zgeqrfp (integer m, integer n, complex*16, dimension( lda, * ) a,
integer lda, complex*16, dimension( * ) tau, complex*16, dimension( * )
work, integer lwork, integer info)
π Function/Subroutine Documentation π
subroutine zgeqrfp (integer m, integer n, complex*16, dimension( lda, * ) a, integer lda, complex*16, dimension( * ) tau, complex*16, dimension( * ) work, integer lwork, integer info) β label ZGEQRFP.
π― Purpose π―
The purpose block says: ZGEQR2P computes a QR factorization of a complex M-by-N matrix A:
A = Q * ( R ),
( 0 )
where:
Qis a M-by-M orthogonal matrix;Ris an upper-triangular N-by-N matrix with nonnegative diagonal entries;0is a (M-N)-by-N zero matrix, if M > N.
π Parameters π
MβMisINTEGER; the number of rows of the matrixA.M >= 0.NβNisINTEGER; the number of columns of the matrixA.N >= 0.AβAisCOMPLEX*16array, dimension(LDA,N). On entry, the M-by-N matrixA. On exit, the elements on and above the diagonal of the array contain themin(M,N)-by-N upper trapezoidal matrixR(Ris upper triangular ifm >= n). The diagonal entries ofRare real and nonnegative; the elements below the diagonal, with the arrayTAU, represent the unitary matrixQas a product ofmin(m,n)elementary reflectors (see Further Details).LDAβLDAisINTEGER; the leading dimension of the arrayA.LDA >= max(1,M).TAUβTAUisCOMPLEX*16array, dimension(min(M,N)); 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; the dimension of the arrayWORK.LWORK >= 1, if MIN(M,N) = 0, and LWORK >= N, otherwise.For optimum performanceLWORK >= N*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, thei-th argument had an illegal value.
π Author π
Univ. of TennesseeUniv. of California BerkeleyUniv. of Colorado DenverNAG Ltd.
π Further Details π
The matrix Q is represented as a product of elementary reflectors:
Q = H(1) H(2) . . . H(k), 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; v(i+1:m) is stored on exit in A(i+1:m,i), and tau in TAU(i).
See Lapack Working Note 203 for details.
Definition at line 149 of file zgeqrfp.f.
π€ Author π€
Generated automatically by Doxygen for LAPACK from the source code.
byeee bestie, thatβs the whole page β LAPACK Version 3.12.0, SRC/zgeqrfp.f(3) π«Άβ¨
π¦ 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