SRC/zlaswlq.f(3)
π NAME π
SRC/zlaswlq.f
okay so this is the SRC/zlaswlq.f(3) page from the Library Functions Manual β literally the file, the whole file, and nothing but the file, bestie. π«Ά
β¨ SYNOPSIS β¨
Functions/Subroutines
subroutine zlaswlq (m, n, mb, nb, a, lda, t, ldt, work, lwork, info)
ZLASWLQ
π¦ Function/Subroutine Documentation π¦
subroutine zlaswlq (integer m, integer n, integer mb, integer nb,
complex*16, dimension( lda, * ) a, integer lda, complex*16, dimension(
ldt, * ) t, integer ldt, complex*16, dimension( * ) work, integer lwork,
integer info)
ZLASWLQ
π Purpose π
okay so ZLASWLQ computes a blocked Tall-Skinny LQ factorization of a complexx M-by-N matrix A for M <= N:
A = ( L 0 ) * Q,
where:
Qis an-by-Northogonal matrix, stored on exit in an implicit form in the elements above the diagonal of the arrayAand in the elements of the arrayT;Lis a lower-triangularM-by-Mmatrix stored on exit in the elements on and below the diagonal of the arrayA.0is aM-by-(N-M)zero matrix, ifM < N, and is not stored.
π Parameters π
MβM is INTEGER; the number of rows of the matrixA.M >= 0.NβN is INTEGER; the number of columns of the matrixA.N >= M >= 0.MBβMB is INTEGER; the row block size to be used in the blocked QR.M >= MB >= 1.NBβNB is INTEGER; the column block size to be used in the blocked QR.NB > 0.AβA is COMPLEX*16 array, dimension (LDA,N). On entry, theM-by-NmatrixA. On exit, the elements on and below the diagonal of the array contain theN-by-Nlower triangular matrixL; the elements above the diagonal representQby the rows of blockedV(see Further Details).LDAβLDA is INTEGER; the leading dimension of the arrayA.LDA >= max(1,M).TβT is COMPLEX*16 array, dimension (LDT, N * Number_of_row_blocks)whereNumber_of_row_blocks = CEIL((N-M)/(NB-M)). The blocked upper triangular block reflectors stored in compact form as a sequence of upper triangular blocks. See Further Details below.LDTβLDT is INTEGER; the leading dimension of the arrayT.LDT >= MB.WORKβ(workspace) COMPLEX*16 array, dimension (MAX(1,LWORK)). On exit, ifINFO = 0,WORK(1)returns the minimalLWORK.LWORKβLWORK is INTEGER; the dimension of the arrayWORK.LWORK >= 1, if MIN(M,N) = 0, and LWORK >= MB*M, otherwise.IfLWORK = -1, then a workspace query is assumed; the routine only calculates the minimal size of theWORKarray, returns this value as the first entry of theWORKarray, and no error message related toLWORKis issued byXERBLA.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: π₯
Short-Wide LQ (SWLQ) performs LQ by a sequence of orthogonal transformations, representing Q as a product of other orthogonal matrices
Q = Q(1) * Q(2) * . . . * Q(k)
where each Q(i) zeros out upper diagonal entries of a block of NB rows of A:
Q(1)zeros out the upper diagonal entries of rows1:NBofAQ(2)zeros out the bottomMB-Nrows of rows[1:M,NB+1:2*NB-M]ofAQ(3)zeros out the bottomMB-Nrows of rows[1:M,2*NB-M+1:3*NB-2*M]ofA. . .
Q(1) is computed by GELQT, which represents Q(1) by Householder vectors stored under the diagonal of rows 1:MB of A, and by upper triangular block reflectors, stored in array T(1:LDT,1:N). For more information see Further Details in GELQT.
Q(i) for i>1 is computed by TPLQT, which represents Q(i) by Householder vectors stored in columns [(i-1)*(NB-M)+M+1:i*(NB-M)+M] of A, and by upper triangular block reflectors, stored in array T(1:LDT,(i-1)*M+1:i*M). The last Q(k) may use fewer rows. For more information see Further Details in TPQRT.
For more details of the overall algorithm, see the description of Sequential TSQR in Section 2.2 of [1].
[1] βCommunication-Optimal Parallel and Sequential QR and LU Factorizations,β J. Demmel, L. Grigori, M. Hoemmen, J. Langou, SIAM J. Sci. Comput, vol. 34, no. 1, 2012
Definition at line 166 of file zlaswlq.f.
π Author π
Generated automatically by Doxygen for LAPACK from the source code.
kisses, bestie π β this was SRC/zlaswlq.f(3) from LAPACK, Version 3.12.0. 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