zlarz (l) - Linux Manuals

zlarz: applies a complex elementary reflector H to a complex M-by-N matrix C, from either the left or the right

NAME

ZLARZ - applies a complex elementary reflector H to a complex M-by-N matrix C, from either the left or the right

SYNOPSIS

SUBROUTINE ZLARZ(
SIDE, M, N, L, V, INCV, TAU, C, LDC, WORK )

    
CHARACTER SIDE

    
INTEGER INCV, L, LDC, M, N

    
COMPLEX*16 TAU

    
COMPLEX*16 C( LDC, * ), V( * ), WORK( * )

PURPOSE

ZLARZ applies a complex elementary reflector H to a complex M-by-N matrix C, from either the left or the right. H is represented in the form

H = I - tau * v * vaq
where tau is a complex scalar and v is a complex vector.
If tau = 0, then H is taken to be the unit matrix.
To apply Haq (the conjugate transpose of H), supply conjg(tau) instead tau.
H is a product of k elementary reflectors as returned by ZTZRZF.

ARGUMENTS

SIDE (input) CHARACTER*1
= aqLaq: form H * C
= aqRaq: form C * H
M (input) INTEGER
The number of rows of the matrix C.
N (input) INTEGER
The number of columns of the matrix C.
L (input) INTEGER
The number of entries of the vector V containing the meaningful part of the Householder vectors. If SIDE = aqLaq, M >= L >= 0, if SIDE = aqRaq, N >= L >= 0.
V (input) COMPLEX*16 array, dimension (1+(L-1)*abs(INCV))
The vector v in the representation of H as returned by ZTZRZF. V is not used if TAU = 0.
INCV (input) INTEGER
The increment between elements of v. INCV <> 0.
TAU (input) COMPLEX*16
The value tau in the representation of H.
C (input/output) COMPLEX*16 array, dimension (LDC,N)
On entry, the M-by-N matrix C. On exit, C is overwritten by the matrix H * C if SIDE = aqLaq, or C * H if SIDE = aqRaq.
LDC (input) INTEGER
The leading dimension of the array C. LDC >= max(1,M).
WORK (workspace) COMPLEX*16 array, dimension
(N) if SIDE = aqLaq or (M) if SIDE = aqRaq

FURTHER DETAILS

Based on contributions by

  A. Petitet, Computer Science Dept., Univ. of Tenn., Knoxville, USA

Pages related to zlarz

  • zlarz (3)
  • zlarzb (l) - applies a complex block reflector H or its transpose H**H to a complex distributed M-by-N C from the left or the right
  • zlarzt (l) - forms the triangular factor T of a complex block reflector H of order > n, which is defined as a product of k elementary reflectors
  • zlar1v (l) - computes the (scaled) r-th column of the inverse of the sumbmatrix in rows B1 through BN of the tridiagonal matrix L D L^T - sigma I
  • zlar2v (l) - applies a vector of complex plane rotations with real cosines from both sides to a sequence of 2-by-2 complex Hermitian matrices,
  • zlarcm (l) - performs a very simple matrix-matrix multiplication
  • zlarf (l) - applies a complex elementary reflector H to a complex M-by-N matrix C, from either the left or the right
  • zlarfb (l) - applies a complex block reflector H or its transpose Haq to a complex M-by-N matrix C, from either the left or the right
  • zlarfg (l) - generates a complex elementary reflector H of order n, such that Haq * ( alpha ) = ( beta ), Haq * H = I
  • zlarfp (l) - generates a complex elementary reflector H of order n, such that Haq * ( alpha ) = ( beta ), Haq * H = I
  • zlarft (l) - forms the triangular factor T of a complex block reflector H of order n, which is defined as a product of k elementary reflectors