SRC/stplqt2.f(3)
π SRC/stplqt2.f(3) π
okay so this page is from the Library Functions Manual, and it is literally all about SRC/stplqt2.f. Bestie, it is giving Fortran LAPACK main-character energy, and we are going to get through every single detail together. no cap. π«Ά
πΈ NAME πΈ
SRC/stplqt2.f β that is the whole name, simple and iconic. β¨
β¨ SYNOPSIS β¨
π§ Functions/Subroutines π§
subroutine stplqt2 (m, n, l, a, lda, b, ldb, t, ldt, info)
STPLQT2 computes a LQ factorization of a real or complex triangular-
pentagonal"matrix,whichiscomposedofatriangularblockandapentagonalblock,usingthecompactWYrepresentationforQ.
π Function/Subroutine Documentation π
subroutine stplqt2 (integer m, integer n, integer l, real, dimension( lda, *
) a, integer lda, real, dimension( ldb, * ) b, integer ldb, real,
dimension( ldt, * ) t, integer ldt, integer info)
STPLQT2 computes a LQ factorization of a real or complex "triangular-pentagonal" matrix, which is composed of a triangular block and a pentagonal block, using the compact WY representation for Q. π
π Purpose π
The purpose, fr fr: STPLQT2 computes a LQ a factorization of a real matrix C, which is composed of a triangular block A and pentagonal block B, using the compact WY representation for Q. π¦
π Parameters π
MβM is INTEGER; the total number of rows of the matrixB;M >= 0. πNβN is INTEGER; the number of columns of the matrixB, and the order of the triangular matrixA;N >= 0. β¨LβL is INTEGER; the number of rows of the lower trapezoidal part ofB;MIN(M,N) >= L >= 0; see Further Details. π§AβA is REAL array, dimension (LDA,M); on entry, the lower triangularM-by-MmatrixA; on exit, the elements on and below the diagonal of the array contain the lower triangular matrixL. πLDAβLDA is INTEGER; the leading dimension of the arrayA;LDA >= max(1,M). πΈBβB is REAL array, dimension (LDB,N); on entry, the pentagonalM-by-NmatrixB; the firstN-Lcolumns are rectangular, and the lastLcolumns are lower trapezoidal; on exit,Bcontains the pentagonal matrixV; see Further Details. π¦LDBβLDB is INTEGER; the leading dimension of the arrayB;LDB >= max(1,M). βTβT is REAL array, dimension (LDT,M); theN-by-Nupper triangular factorTof the block reflector; see Further Details. π₯LDTβLDT is INTEGER; the leading dimension of the arrayT;LDT >= max(1,M). π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 π¦
The input matrix C is a M-by-(M+N) matrix:
C = [ A ][ B ]
where A is an lower triangular M-by-M matrix, and B is M-by-N pentagonal matrix consisting of a M-by-(N-L) rectangular matrix B1 left of a M-by-L upper trapezoidal matrix B2:
B = [ B1 ][ B2 ]
[ B1 ] <- M-by-(N-L) rectangular
[ B2 ] <- M-by-L lower trapezoidal.
The lower trapezoidal matrix B2 consists of the first L columns of a N-by-N lower triangular matrix, where 0 <= L <= MIN(M,N). If L=0, B is rectangular M-by-N; if M=L=N, B is lower triangular. π
The matrix W stores the elementary reflectors H(i) in the i-th row above the diagonal (of A) in the M-by-(M+N) input matrix C:
C = [ A ][ B ]
[ A ] <- lower triangular M-by-M
[ B ] <- M-by-N pentagonal
so that W can be represented as:
W = [ I ][ V ]
[ I ] <- identity, M-by-M
[ V ] <- M-by-N, same form as B.
Thus, all of information needed for W is contained on exit in B, which we call V above. Note that V has the same form as B; that is:
W = [ V1 ][ V2 ]
[ V1 ] <- M-by-(N-L) rectangular
[ V2 ] <- M-by-L lower trapezoidal.
The rows of V represent the vectors which define the H(i)'s. The (M+N)-by-(M+N) block reflector H is then given by:
H = I - W**T * T * W
where W^H is the conjugate transpose of W and T is the upper triangular factor of the block reflector. πΈ
Definition at line 176 of file stplqt2.f. β
πΈ Author πΈ
Generated automatically by Doxygen for LAPACK from the source code. π«Ά
signed off with LAPACK Version 3.12.0 β source SRC/stplqt2.f(3). periodt. π
β¨
π¦ 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