SRC/cunmr3.f(3)
π NAME
SRC/cunmr3.f β that's the name, bestie. This page is from the Library Functions Manual, and it's giving very specific LAPACK energy.
β¨ SYNOPSIS
Functions/Subroutines
subroutine cunmr3 (side, trans, m, n, k, l, a, lda, tau, c, ldc, work,
info)
CUNMR3 multiplies a general matrix by the unitary matrix from a RZ
factorization determined by ctzrzf (unblocked algorithm).
π¦ Function/Subroutine Documentation
subroutine cunmr3 (character side, character trans, integer m, integer n,
integer k, integer l, complex, dimension( lda, * ) a, integer lda,
complex, dimension( * ) tau, complex, dimension( ldc, * ) c, integer
ldc, complex, dimension( * ) work, integer info)
CUNMR3 multiplies a general matrix by the unitary matrix from a RZ factorization determined by ctzrzf (unblocked algorithm). Iconic behavior, honestly.
π Purpose
Okay so CUNMR3 overwrites the general complex m by n matrix C with:
Q * C if SIDE = 'L' and TRANS = 'N', or
Q**H* C if SIDE = 'L' and TRANS = 'C', or
C * Q if SIDE = 'R' and TRANS = 'N', or
C * Q**H if SIDE = 'R' and TRANS = 'C',
where Q is a complex unitary matrix defined as the product of k elementary reflectors:
Q = H(1) H(2) . . . H(k)
as returned by CTZRZF. Q is of order m if SIDE = 'L' and of order n if SIDE = 'R'. Slay.
πΈ Parameters
SIDEβSIDE is CHARACTER*1;= 'L': applyQorQ**Hfrom the Left;= 'R': applyQorQ**Hfrom the Right. Pick a side, bestie.TRANSβTRANS is CHARACTER*1;= 'N': applyQ(No transpose);= 'C': applyQ**H(Conjugate transpose). TheCis for conjugate, not for "casual."MβM is INTEGER; the number of rows of the matrixC.M >= 0.NβN is INTEGER; the number of columns of the matrixC.N >= 0.KβK is INTEGER; the number of elementary reflectors whose product defines the matrixQ. IfSIDE = 'L',M >= K >= 0; ifSIDE = 'R',N >= K >= 0.LβL is INTEGER; the number of columns of the matrixAcontaining the meaningful part of the Householder reflectors. IfSIDE = 'L',M >= L >= 0, ifSIDE = 'R',N >= L >= 0.AβA is COMPLEX array, dimension(LDA,M)ifSIDE = 'L',(LDA,N)ifSIDE = 'R'. The i-th row must contain the vector which defines the elementary reflectorH(i), fori = 1,2,...,k, as returned byCTZRZFin the lastkrows of its array argumentA.Ais modified by the routine but restored on exit β so it gets touched, but it comes back, no cap.LDAβLDA is INTEGER; the leading dimension of the arrayA.LDA >= max(1,K).TAUβTAU is COMPLEX array, dimension(K);TAU(i)must contain the scalar factor of the elementary reflectorH(i), as returned byCTZRZF.CβC is COMPLEX array, dimension(LDC,N). On entry, them-by-nmatrixC. On exit,Cis overwritten byQ*CorQ**H*CorC*Q**HorC*Q.LDCβLDC is INTEGER; the leading dimension of the arrayC.LDC >= max(1,M).WORKβWORK is COMPLEX array, dimension(N)ifSIDE = 'L',(M)ifSIDE = 'R'.INFOβINFO is INTEGER;= 0: successful exit;< 0: ifINFO = -i, the i-th argument had an illegal value. Rude little error code, but at least it's clear.
π Author
- Univ. of Tennessee
- Univ. of California Berkeley
- Univ. of Colorado Denver
- NAG Ltd.
π«Ά Contributors:
A. Petitet, Computer Science Dept., Univ. of Tenn., Knoxville, USA
π§ Further Details
Literally just !> and then nothing β the section is empty, bestie, no extra lore.
Definition at line 176 of file cunmr3.f.
π Author
Generated automatically by Doxygen for LAPACK from the source code.
okay bye bestie β this whole page is from LAPACK, Version 3.12.0, for SRC/cunmr3.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