SRC/sgemlqt.f(3)

πŸ’… SRC/sgemlqt.f(3) β€” Library Functions Manual πŸ’…

πŸŽ€ NAME πŸŽ€

SRC/sgemlqt.f β€” bestie, that's the file this page is about. The routine inside is sgemlqt, and it's giving compact-WY reflector realness.

✨ SYNOPSIS ✨

Functions/Subroutines

subroutine sgemlqt (side, trans, m, n, k, mb, v, ldv, t, ldt, c, ldc,
     work, info)

SGEMLQT

πŸ’– Function/Subroutine Documentation πŸ’–

subroutine sgemlqt (character side, character trans, integer m, integer n,
     integer k, integer mb, real, dimension( ldv, * ) v, integer ldv, real,
     dimension( ldt, * ) t, integer ldt, real, dimension( ldc, * ) c, integer
     ldc, real, dimension( * ) work, integer info)

SGEMLQT

🌸 Purpose

okay so SGEMLQT overwrites the general real M-by-N matrix C with:

                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) = I - V T V**T

generated using the compact WY representation as returned by SGELQT. πŸ‘‘

Q is of order M if SIDE = 'L' and of order N if SIDE = 'R'.

πŸ¦‹ Parameters

  • SIDE β€” SIDE is CHARACTER*1; = 'L': apply Q or QT from the Left; = 'R': apply Q or QT from the Right.
  • TRANS β€” TRANS is CHARACTER*1; = 'N': No transpose, apply Q; = 'C': Transpose, apply Q**T.
  • M β€” M is INTEGER; the number of rows of the matrix C. M >= 0.
  • N β€” N is INTEGER; the number of columns of the matrix C. N >= 0.
  • K β€” K is INTEGER; the number of elementary reflectors whose product defines the matrix Q. If SIDE = 'L', M >= K >= 0; if SIDE = 'R', N >= K >= 0.
  • MB β€” MB is INTEGER; the block size used for the storage of T. K >= MB >= 1. This must be the same value of MB used to generate T in SGELQT.
  • V β€” V is REAL array, dimension (LDV,M) if SIDE = 'L', (LDV,N) if SIDE = 'R'. The i-th row must contain the vector which defines the elementary reflector H(i), for i = 1,2,...,k, as returned by SGELQT in the first K rows of its array argument A.
  • LDV β€” LDV is INTEGER; the leading dimension of the array V. LDV >= max(1,K).
  • T β€” T is REAL array, dimension (LDT,K). The upper triangular factors of the block reflectors as returned by SGELQT, stored as a MB-by-K matrix.
  • LDT β€” LDT is INTEGER; the leading dimension of the array T. LDT >= MB.
  • C β€” C is REAL array, dimension (LDC,N). On entry, the M-by-N matrix C. On exit, C is overwritten by Q C, Q**T C, C Q**T or C Q.
  • LDC β€” LDC is INTEGER; the leading dimension of the array C. LDC >= max(1,M).
  • WORK β€” WORK is REAL array. The dimension of WORK is N*MB if SIDE = 'L', or M*MB if SIDE = 'R'.
  • INFO β€” INFO is INTEGER; = 0: successful exit; < 0: if INFO = -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 151 of file sgemlqt.f.

🌈 Author 🌈

Generated automatically by Doxygen for LAPACK from the source code.


okay bestie, that's the whole page πŸ’– β€” it's from LAPACK Version 3.12.0, source SRC/sgemlqt.f(3), and the footer doesn't print a separate date, so no secret date lore here. slay and bye! πŸ¦‹βœ¨

πŸ“¦ where this page comes from