SRC/sormlq.f(3)
πΈ NAME πΈ
okay so this is the page for SRC/sormlq.f from the Library Functions Manual, bestie. It's literally all about SORMLQ, the real orthogonal-matrix-applying subroutine, and we're gonna get into every single detail. π
β¨
β¨ SYNOPSIS β¨
π Functions/Subroutines π
subroutine sormlq (side, trans, m, n, k, a, lda, tau, c, ldc, work,
lwork, info)
SORMLQ
π Function/Subroutine Documentation π
subroutine sormlq (character side, character trans, integer m, integer n,
integer k, real, dimension( lda, * ) a, integer lda, real, dimension( *
) tau, real, dimension( ldc, * ) c, integer ldc, real, dimension( * )
work, integer lwork, integer info)
SORMLQ
π¦ Purpose π¦
okay so SORMLQ overwrites the general real M-by-N matrix C with the following, and yes this is the moment:
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(k) . . . H(2) H(1)
as returned by SGELQF. Q is of order M if SIDE = 'L' and of order N if SIDE = 'R'. periodt. π
π Parameters π
SIDEβSIDEisCHARACTER*1;= 'L': applyQorQ**Tfrom the Left;= 'R': applyQorQ**Tfrom the Right. π«ΆTRANSβTRANSisCHARACTER*1;= 'N': No transpose, applyQ;= 'T': Transpose, applyQ**T. β¨MβMisINTEGER; the number of rows of the matrixC.M >= 0. π«NβNisINTEGER; the number of columns of the matrixC.N >= 0. πΈKβKisINTEGER; the number of elementary reflectors whose product defines the matrixQ. IfSIDE = 'L',M >= K >= 0; ifSIDE = 'R',N >= K >= 0. πAβAisREALarray, dimension(LDA,M)ifSIDE = 'L',(LDA,N)ifSIDE = 'R'. Thei-th row must contain the vector which defines the elementary reflectorH(i), fori = 1,2,...,k, as returned bySGELQFin the firstkrows of its array argumentA. π¦LDAβLDAisINTEGER; the leading dimension of the arrayA.LDA >= max(1,K). πTAUβTAUisREALarray, dimension(K);TAU(i)must contain the scalar factor of the elementary reflectorH(i), as returned bySGELQF. πCβCisREALarray, dimension(LDC,N). On entry, theM-by-NmatrixC. On exit,Cis overwritten byQ*CorQ**T*CorC*Q**TorC*Q. π₯LDCβLDCisINTEGER; the leading dimension of the arrayC.LDC >= max(1,M). πWORKβWORKisREALarray, dimension(MAX(1,LWORK)). On exit, ifINFO = 0,WORK(1)returns the optimalLWORK. β¨LWORKβLWORKisINTEGER; the dimension of the arrayWORK. IfSIDE = 'L',LWORK >= max(1,N); ifSIDE = 'R',LWORK >= max(1,M). For good performance,LWORKshould generally be larger. 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 Tennessee
- Univ. of California Berkeley
- Univ. of Colorado Denver
- NAG Ltd.
Definition at line 166 of file sormlq.f. π―
π Author π
Generated automatically by Doxygen for LAPACK from the source code. π«Ά
okay bestie, that's the whole thing β this cutie is from LAPACK, Version 3.12.0, and the page source is SRC/sormlq.f(3). mwah! πβ¨π
π¦ 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