SRC/zgemlqt.f(3)

πŸŽ€ SRC/zgemlqt.f(3)

okay bestie, this is the Library Functions Manual page for SRC/zgemlqt.f(3) β€” we're folding the manual name in right here instead of copy-pasting that raw header line, because we're classy like that. πŸ’…

πŸ’– NAME

SRC/zgemlqt.f

✨ SYNOPSIS

Functions/Subroutines

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

That's the ZGEMLQT subroutine, no cap.

🌸 Function/Subroutine Documentation

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

ZGEMLQT

πŸ’‘ Purpose

ZGEMLQT overwrites the general complex M-by-N matrix C with:

                SIDE = 'L'     SIDE = 'R'
TRANS = 'N':      Q C            C Q
TRANS = 'C':   Q**H C            C Q**H

where Q is a complex unitary matrix defined as the product of K elementary reflectors:

      Q = H(1) H(2) . . . H(K) = I - V T V**H

generated using the compact WY representation as returned by ZGELQT.

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

πŸ“‹ Parameters

  • SIDE β€” CHARACTER*1
    • = 'L': apply Q or Q**H from the Left;
    • = 'R': apply Q or Q**H from the Right.
  • TRANS β€” CHARACTER*1
    • = 'N': No transpose, apply Q;
    • = 'C': Conjugate transpose, apply Q**H.
  • M β€” INTEGER The number of rows of the matrix C. M >= 0.
  • N β€” INTEGER The number of columns of the matrix C. N >= 0.
  • K β€” 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 β€” 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 ZGELQT.
  • V β€” COMPLEX*16 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 ZGELQT in the first K rows of its array argument A.
  • LDV β€” INTEGER The leading dimension of the array V. LDV >= max(1,K).
  • T β€” COMPLEX*16 array, dimension (LDT,K) The upper triangular factors of the block reflectors as returned by ZGELQT, stored as a MB-by-K matrix.
  • LDT β€” INTEGER The leading dimension of the array T. LDT >= MB.
  • C β€” COMPLEX*16 array, dimension (LDC,N) On entry, the M-by-N matrix C. On exit, C is overwritten by Q C, Q**H C, C Q**H or C Q.
  • LDC β€” INTEGER The leading dimension of the array C. LDC >= max(1,M).
  • WORK β€” COMPLEX*16 array. The dimension of WORK is N*MB if SIDE = 'L', or M*MB if SIDE = 'R'.
  • INFO β€” 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 166 of file zgemlqt.f.

🫢 Author

Generated automatically by Doxygen for LAPACK from the source code.


Cute footer sign-off: this page is SRC/zgemlqt.f(3) from LAPACK, Version 3.12.0; the footer gives us that source/version + page stamp, and no extra date line, so we're honoring exactly that. Mwah πŸ’‹βœ¨

πŸ“¦ where this page comes from