SRC/DEPRECATED/cgeqpf.f(3)
π SRC/DEPRECATED/cgeqpf.f(3) π
okay bestie, this page is from the Library Functions Manual, and itβs all about SRC/DEPRECATED/cgeqpf.f β yes, the path is the title, we donβt make the rules. π
β¨ NAME β¨
SRC/DEPRECATED/cgeqpf.f
πΈ SYNOPSIS πΈ
Functions/Subroutines
subroutine cgeqpf (m, n, a, lda, jpvt, tau, work, rwork, info)
CGEQPF
π¦ Function/Subroutine Documentation π¦
subroutine cgeqpf (integer m, integer n, complex, dimension( lda, * ) a,
integer lda, integer, dimension( * ) jpvt, complex, dimension( * ) tau,
complex, dimension( * ) work, real, dimension( * ) rwork, integer info)
CGEQPF
Purpose
This routine is deprecated and has been replaced by routine CGEQP3. π«Ά
CGEQPF computes a QR factorization with column pivoting of a complex M-by-N matrix A: A*P = Q*R. β¨
Parameters
MβM is INTEGER. The number of rows of the matrixA. Constraint:M >= 0. πNβN is INTEGER. The number of columns of the matrixA. Constraint:N >= 0. π¦AβA is COMPLEX array, dimension (LDA,N). On entry, theM-by-NmatrixA. On exit, the upper triangle of the array contains themin(M,N)-by-Nupper triangular matrixR; the elements below the diagonal, together with the arrayTAU, represent the unitary matrixQas a product ofmin(m,n)elementary reflectors. πLDAβLDA is INTEGER. The leading dimension of the arrayA. Constraint:LDA >= max(1,M). πJPVTβJPVT is INTEGER array, dimension (N). On entry, ifJPVT(i) .ne. 0, thei-th column ofAis permuted to the front ofA*P(a leading column); ifJPVT(i) = 0, thei-th column ofAis a free column. On exit, ifJPVT(i) = k, then thei-th column ofA*Pwas thek-th column ofA. π₯TAUβTAU is COMPLEX array, dimension (min(M,N)). The scalar factors of the elementary reflectors. βWORKβWORK is COMPLEX array, dimension (N). πRWORKβRWORK is REAL array, dimension (2*N). πΈ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(2) . . . H(n)
Each H(i) has the form:
H = 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). π§
The matrix P is represented in jpvt as follows: If
jpvt(j) = i
then the jth column of P is the ith canonical unit vector. π
Partial column norm updating strategy modified by Z. Drmac and Z. Bujanovic, Dept. of Mathematics, University of Zagreb, Croatia. -- April 2011 -- For more details see LAPACK Working Note 176.
Definition at line 147 of file cgeqpf.f. π«Ά
π Author π
Generated automatically by Doxygen for LAPACK from the source code.
Thatβs the whole page, bestie β mwah π. This one is LAPACK Version 3.12.0, from SRC/DEPRECATED/cgeqpf.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