?unmrz
?unmrz
Multiplies a complex matrix by the unitary matrix defined from the factorization formed by
?tzrzf
.Syntax
lapack_int
LAPACKE_cunmrz
(
int
matrix_layout
,
char
side
,
char
trans
,
lapack_int
m
,
lapack_int
n
,
lapack_int
k
,
lapack_int
l
,
const
lapack_complex_float
*
a
,
lapack_int
lda
,
const
lapack_complex_float
*
tau
,
lapack_complex_float
*
c
,
lapack_int
ldc
);
lapack_int
LAPACKE_zunmrz
(
int
matrix_layout
,
char
side
,
char
trans
,
lapack_int
m
,
lapack_int
n
,
lapack_int
k
,
lapack_int
l
,
const
lapack_complex_double
*
a
,
lapack_int
lda
,
const
lapack_complex_double
*
tau
,
lapack_complex_double
*
c
,
lapack_int
ldc
);
Include Files
- mkl.h
Description
The routine multiplies a complex , where :
m
-by-n
matrix C
by Q
or Q
H
Q
is the unitary matrix defined as a product of k
elementary reflectors H
(i
)Q
= H
(1)H
H
(2)
H
H
(k
)H
Depending on the parameters (overwriting the result over
side
and trans
, the routine can form one of the matrix products Q*C
, Q
H
*C
, C*Q
, or C*Q
H
C
).The matrix and of order
.
Q
is of order m
if side
= 'L'
n
if side
= 'R'
Input Parameters
- matrix_layout
- Specifies whether matrix storage layout is row major (LAPACK_ROW_MAJOR) or column major (LAPACK_COL_MAJOR).
- side
- Must be either'L'or'R'.If,side='L'QorQis applied toHCfrom the left.If,side='R'QorQis applied toHCfrom the right.
- trans
- Must be either'N'or'C'.If, the routine multipliestrans='N'CbyQ.If, the routine multipliestrans='C'CbyQ.H
- m
- The number of rows in the matrixC().m≥0
- n
- The number of columns inC().n≥0
- k
- The number of elementary reflectors whose product defines the matrixQ. Constraints:0, if≤k≤m;side='L'0, if≤k≤n.side='R'
- l
- The number of columns of the matrixAcontaining the meaningful part of the Householder reflectors. Constraints:0, if≤l≤m;side='L'0, if≤l≤n.side='R'
- a,tau,c
- Arrays:a(size forside= 'L': max(1,lda*m) for column major layout and max(1,lda*k) for row major layout; forside= 'R': max(1,lda*b) for column major layout and max(1,lda*k) for row major layout),tau,c(size max(1,ldc*n) for column major layout and max(1,ldc*m) for row major layout).On entry, theith row ofamust contain the vector which defines the elementary reflector, for i = 1,2,...,H(i)k, as returned byctzrzf/ztzrzfin the lastkrows of its array argumenta.tau[must contain the scalar factor of the elementary reflectori- 1], as returned byH(i)ctzrzf/ztzrzf.The size oftaumust be at least max(1,k).ccontains them-by-nmatrixC.
- lda
- The leading dimension ofa;lda≥max(1,k)for column major layout. For row major layout,.iflda≥max(1,m)side= 'L', andiflda≥max(1,n)side= 'R'
- ldc
- The leading dimension ofc;ldc≥max(1,m)for column major layout and max(1,.n) for row major layout
Output Parameters
- c
- Overwritten by the productQ*C,QH*C,C*Q, orC*Q(as specified byHsideandtrans).
Return Values
This function returns a value
info
.If , the execution is successful.
info
=0If , the
info
= -i
i
-th parameter had an illegal value.Application Notes
The real counterpart of this routine is ormrz.