SRC/dormlq.f(3)
π NAME π
okay bestie, this is SRC/dormlq.f straight out of the Library Functions Manual, and it's literally the LAPACK routine page for DORMLQ. slay.
π¦ SYNOPSIS π¦
β¨ Functions/Subroutines β¨
subroutine dormlq (side, trans, m, n, k, a, lda, tau, c, ldc, work, lwork, info)βDORMLQ
subroutine dormlq (side, trans, m, n, k, a, lda, tau, c, ldc, work,
lwork, info)
DORMLQ
π Function/Subroutine Documentation π
subroutine dormlq (character side, character trans, integer m, integer n,
integer k, double precision, dimension( lda, * ) a, integer lda, double
precision, dimension( * ) tau, double precision, dimension( ldc, * ) c,
integer ldc, double precision, dimension( * ) work, integer lwork,
integer info)
DORMLQ
π Purpose π
okay so DORMLQ overwrites the general real M-by-N matrix C with this exact energy:
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 DGELQF. Q is of order M if SIDE = 'L' and of order N if SIDE = 'R'. periodt.
π Parameters π
SIDE
SIDEβSIDEisCHARACTER*1= 'L': applyQorQ**Tfrom the Left;= 'R': applyQorQ**Tfrom the Right.
TRANS
TRANSβTRANSisCHARACTER*1= 'N': No transpose, applyQ;= 'T': Transpose, applyQ**T.
M
MβMisINTEGER- The number of rows of the matrix
C.M >= 0.
- The number of rows of the matrix
N
NβNisINTEGER- The number of columns of the matrix
C.N >= 0.
- The number of columns of the matrix
K
KβKisINTEGER- The number of elementary reflectors whose product defines the matrix
Q. - If
SIDE = 'L',M >= K >= 0; - if
SIDE = 'R',N >= K >= 0.
- The number of elementary reflectors whose product defines the matrix
A
AβAisDOUBLE PRECISION array, dimension
(LDA,M) if SIDE = 'L',
(LDA,N) if SIDE = 'R'
- The i-th row must contain the vector which defines the elementary reflector
H(i), fori = 1,2,...,k, as returned byDGELQFin the firstkrows of its array argumentA.
LDA
LDAβLDAisINTEGER- The leading dimension of the array
A.LDA >= max(1,K).
- The leading dimension of the array
TAU
TAUβTAUisDOUBLE PRECISION array, dimension(K)TAU(i)must contain the scalar factor of the elementary reflectorH(i), as returned byDGELQF.
C
CβCisDOUBLE PRECISION array, dimension(LDC,N)- On entry, the M-by-N matrix
C. - On exit,
Cis overwritten byQ*CorQ**T*CorC*Q**TorC*Q.
- On entry, the M-by-N matrix
LDC
LDCβLDCisINTEGER- The leading dimension of the array
C.LDC >= max(1,M).
- The leading dimension of the array
WORK
WORKβWORKisDOUBLE PRECISION array, dimension(MAX(1,LWORK))- On exit, if
INFO = 0,WORK(1)returns the optimalLWORK.
- On exit, if
LWORK
LWORKβLWORKisINTEGER- The dimension of the array
WORK. - If
SIDE = 'L',LWORK >= max(1,N); - if
SIDE = 'R',LWORK >= max(1,M). - For good performance,
LWORKshould generally be larger.
- The dimension of the array
bestie, this next bit is a whole mood: 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
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 165 of file dormlq.f.
π Author π
Generated automatically by Doxygen for LAPACK from the source code.
final sign-off, bestie: LAPACK, Version 3.12.0, page SRC/dormlq.f(3) β that's the footer tea, and there's no date line in the original so I'm not inventing one π«Άπβ¨
π¦ 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