SRC/sormrq.f(3)
πΈ NAME πΈ
Bestie, welcome to SRC/sormrq.f(3) from the Library Functions Manual β the source file is SRC/sormrq.f, and the subroutine living here is SORMRQ. Itβs giving iconic linear-algebra behavior. π
β¨ SYNOPSIS β¨
π¦ Functions/Subroutines π¦
subroutine sormrq (side, trans, m, n, k, a, lda, tau, c, ldc, work,
lwork, info)
SORMRQ
π Function/Subroutine Documentation π
subroutine sormrq (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)
SORMRQ
π Purpose π
Okay so SORMRQ overwrites the general real M-by-N matrix C with one of these four looks:
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)
as returned by SGERQF. Q is of order M if SIDE = 'L' and of order N if SIDE = 'R'. No cap. π
π§ 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βAisREAL array, dimension(LDA,M)ifSIDE = 'L',(LDA,N)ifSIDE = 'R'.
The i-th row must contain the vector which defines the elementary reflector H(i), for i = 1,2,...,k, as returned by SGERQF in the last k rows of its array argument A.
LDAβLDAisINTEGER. The leading dimension of the arrayA.LDA >= max(1,K).TAUβTAUisREAL array, dimension(K).TAU(i)must contain the scalar factor of the elementary reflectorH(i), as returned bySGERQF.CβCisREAL array, 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βWORKisREAL array, 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.
Heads up, bestie: 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 sormrq.f.
π Author π
Generated automatically by Doxygen for LAPACK from the source code.
Thatβs the whole vibe, bestie β source SRC/sormrq.f(3), LAPACK Version 3.12.0. 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