zunmr3 - torization determined by ctzrzf (unblocked algorithm)
SUBROUTINE ZUNMR3(SIDE, TRANS, M, N, K, L, A, LDA, TAU, C, LDC, WORK, INFO) CHARACTER*1 SIDE, TRANS INTEGER INFO, K, L, LDA, LDC, M, N DOUBLE COMPLEX A(LDA,*), C(LDC,*), TAU(*),WORK(*) SUBROUTINE ZUNMR3_64(SIDE, TRANS, M, N, K, L, A, LDA, TAU, C, LDC, WORK, INFO) CHARACTER*1 SIDE, TRANS INTEGER*8 INFO, K, L, LDA, LDC, M, N DOUBLE COMPLEX A(LDA,*), C(LDC,*), TAU(*),WORK(*) F95 INTERFACE SUBROUTINE UNMR3(SIDE, TRANS, M, N, K, L, A, LDA, TAU, C, LDC, WORK, INFO) INTEGER :: M, N, K, L, LDA, LDC, INFO CHARACTER(LEN=1) :: SIDE, TRANS COMPLEX(8), DIMENSION(:,:) :: A, C COMPLEX(8), DIMENSION(:) :: TAU, WORK SUBROUTINE UNMR3_64(SIDE, TRANS, M, N, K, L, A, LDA, TAU, C, LDC, WORK, INFO) INTEGER(8) :: M, N, K, L, LDA, LDC, INFO CHARACTER(LEN=1) :: SIDE, TRANS COMPLEX(8), DIMENSION(:,:) :: A, C COMPLEX(8), DIMENSION(:) :: TAU, WORK C INTERFACE #include <sunperf.h> void zunmr3 (char side, char trans, int m, int n, int k, int l, double- complex *a, int lda, doublecomplex *tau, doublecomplex *c, int ldc, int *info); void zunmr3_64 (char side, char trans, long m, long n, long k, long l, doublecomplex *a, long lda, doublecomplex *tau, doublecomplex *c, long ldc, long *info);
Oracle Solaris Studio Performance Library zunmr3(3P)
NAME
zunmr3 - multiply a general matrix by the unitary matrix from a RZ fac-
torization determined by ctzrzf (unblocked algorithm)
SYNOPSIS
SUBROUTINE ZUNMR3(SIDE, TRANS, M, N, K, L, A, LDA, TAU, C, LDC, WORK,
INFO)
CHARACTER*1 SIDE, TRANS
INTEGER INFO, K, L, LDA, LDC, M, N
DOUBLE COMPLEX A(LDA,*), C(LDC,*), TAU(*),WORK(*)
SUBROUTINE ZUNMR3_64(SIDE, TRANS, M, N, K, L, A, LDA, TAU, C, LDC,
WORK, INFO)
CHARACTER*1 SIDE, TRANS
INTEGER*8 INFO, K, L, LDA, LDC, M, N
DOUBLE COMPLEX A(LDA,*), C(LDC,*), TAU(*),WORK(*)
F95 INTERFACE
SUBROUTINE UNMR3(SIDE, TRANS, M, N, K, L, A, LDA, TAU, C, LDC, WORK,
INFO)
INTEGER :: M, N, K, L, LDA, LDC, INFO
CHARACTER(LEN=1) :: SIDE, TRANS
COMPLEX(8), DIMENSION(:,:) :: A, C
COMPLEX(8), DIMENSION(:) :: TAU, WORK
SUBROUTINE UNMR3_64(SIDE, TRANS, M, N, K, L, A, LDA, TAU, C, LDC, WORK,
INFO)
INTEGER(8) :: M, N, K, L, LDA, LDC, INFO
CHARACTER(LEN=1) :: SIDE, TRANS
COMPLEX(8), DIMENSION(:,:) :: A, C
COMPLEX(8), DIMENSION(:) :: TAU, WORK
C INTERFACE
#include <sunperf.h>
void zunmr3 (char side, char trans, int m, int n, int k, int l, double-
complex *a, int lda, doublecomplex *tau, doublecomplex *c,
int ldc, int *info);
void zunmr3_64 (char side, char trans, long m, long n, long k, long l,
doublecomplex *a, long lda, doublecomplex *tau, doublecomplex
*c, long ldc, long *info);
PURPOSE
zunmr3 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 elemen-
tary reflectors
Q = H(1) H(2) . . . H(k)
as returned by ZTZRZF. Q is of order m if SIDE = 'L' and of order n if
SIDE = 'R'.
ARGUMENTS
SIDE (input)
SIDE is CHARACTER*1
= 'L': apply Q or Q**H from the Left
= 'R': apply Q or Q**H from the Right
TRANS (input)
TRANS is CHARACTER*1
= 'N': apply Q (No transpose)
= 'C': apply Q**H (Conjugate transpose)
M (input)
M is INTEGER
The number of rows of the matrix C. M >= 0.
N (input)
N is INTEGER
The number of columns of the matrix C. N >= 0.
K (input)
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.
L (input)
L is INTEGER
The number of columns of the matrix A containing
the meaningful part of the Householder reflectors.
If SIDE = 'L', M >= L >= 0, if SIDE = 'R', N >= L >= 0.
A (input)
A is COMPLEX*16 array, dimension
(LDA,M) if SIDE = 'L',
(LDA,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
ZTZRZF in the last k rows of its array argument A.
A is modified by the routine but restored on exit.
LDA (input)
LDA is INTEGER
The leading dimension of the array A. LDA >= max(1,K).
TAU (input)
TAU is COMPLEX*16 array, dimension (K)
TAU(i) must contain the scalar factor of the elementary
reflector H(i), as returned by ZTZRZF.
C (input/output)
C is COMPLEX*16 array, dimension (LDC,N)
On entry, the m-by-n matrix C.
On exit, C is overwritten by Q*C or Q**H*C or C*Q**H or C*Q.
LDC (input)
LDC is INTEGER
The leading dimension of the array C. LDC >= max(1,M).
WORK (output)
WORK is COMPLEX*16 array, dimension
(N) if SIDE = 'L',
(M) if SIDE = 'R'
INFO (output)
INFO is INTEGER
= 0: successful exit;
< 0: if INFO = -i, the i-th argument had an illegal value.
7 Nov 2015 zunmr3(3P)