SRC/sgemlqt.f(3)
π
SRC/sgemlqt.f(3) β Library Functions Manual π
π NAME π
SRC/sgemlqt.f β bestie, that's the file this page is about. The routine inside is sgemlqt, and it's giving compact-WY reflector realness.
β¨ SYNOPSIS β¨
Functions/Subroutines
subroutine sgemlqt (side, trans, m, n, k, mb, v, ldv, t, ldt, c, ldc,
work, info)
SGEMLQT
π Function/Subroutine Documentation π
subroutine sgemlqt (character side, character trans, integer m, integer n,
integer k, integer mb, real, dimension( ldv, * ) v, integer ldv, real,
dimension( ldt, * ) t, integer ldt, real, dimension( ldc, * ) c, integer
ldc, real, dimension( * ) work, integer info)
SGEMLQT
πΈ Purpose
okay so SGEMLQT overwrites the general real M-by-N matrix C with:
SIDE = 'L' SIDE = 'R'
TRANS = 'N': Q C C Q
TRANS = 'T': Q**T C C Q**T
where Q is a real orthogonal matrix defined as the product of K elementary reflectors:
Q = H(1) H(2) . . . H(K) = I - V T V**T
generated using the compact WY representation as returned by SGELQT. π
Q is of order M if SIDE = 'L' and of order N if SIDE = 'R'.
π¦ Parameters
SIDEβSIDEisCHARACTER*1;= 'L': apply Q or QT from the Left;= 'R': apply Q or QT from the Right.TRANSβTRANSisCHARACTER*1;= 'N': No transpose, apply Q;= 'C': Transpose, apply Q**T.MβMisINTEGER; the number of rows of the matrix C.M >= 0.NβNisINTEGER; the number of columns of the matrix C.N >= 0.KβKisINTEGER; the number of elementary reflectors whose product defines the matrix Q. IfSIDE = 'L',M >= K >= 0; ifSIDE = 'R',N >= K >= 0.MBβMBisINTEGER; the block size used for the storage of T.K >= MB >= 1. This must be the same value of MB used to generate T inSGELQT.VβVisREALarray, dimension(LDV,M)ifSIDE = 'L',(LDV,N)ifSIDE = 'R'. The i-th row must contain the vector which defines the elementary reflectorH(i), fori = 1,2,...,k, as returned bySGELQTin the first K rows of its array argumentA.LDVβLDVisINTEGER; the leading dimension of the array V.LDV >= max(1,K).TβTisREALarray, dimension(LDT,K). The upper triangular factors of the block reflectors as returned bySGELQT, stored as aMB-by-Kmatrix.LDTβLDTisINTEGER; the leading dimension of the array T.LDT >= MB.CβCisREALarray, dimension(LDC,N). On entry, the M-by-N matrix C. On exit, C is overwritten byQ C,Q**T C,C Q**TorC Q.LDCβLDCisINTEGER; the leading dimension of the array C.LDC >= max(1,M).WORKβWORKisREALarray. The dimension ofWORKisN*MBifSIDE = 'L', orM*MBifSIDE = 'R'.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.
Definition at line 151 of file sgemlqt.f.
π Author π
Generated automatically by Doxygen for LAPACK from the source code.
okay bestie, that's the whole page π β it's from LAPACK Version 3.12.0, source SRC/sgemlqt.f(3), and the footer doesn't print a separate date, so no secret date lore here. slay and bye! π¦β¨
π¦ 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