1ZUPMTR(1)                LAPACK routine (version 3.1)                ZUPMTR(1)
2
3
4

NAME

6       ZUPMTR  -  the general complex M-by-N matrix C with   SIDE = 'L' SIDE =
7       'R' TRANS = 'N'
8

SYNOPSIS

10       SUBROUTINE ZUPMTR( SIDE, UPLO, TRANS, M, N, AP, TAU, C, LDC, WORK, INFO
11                          )
12
13           CHARACTER      SIDE, TRANS, UPLO
14
15           INTEGER        INFO, LDC, M, N
16
17           COMPLEX*16     AP( * ), C( LDC, * ), TAU( * ), WORK( * )
18

PURPOSE

20       ZUPMTR overwrites the general complex M-by-N matrix C with TRANS = 'C':
21       Q**H * C       C * Q**H
22
23       where Q is a complex unitary matrix of order nq, with nq = m if SIDE  =
24       'L'  and nq = n if SIDE = 'R'. Q is defined as the product of nq-1 ele‐
25       mentary reflectors, as returned by ZHPTRD using packed storage:
26
27       if UPLO = 'U', Q = H(nq-1) . . . H(2) H(1);
28
29       if UPLO = 'L', Q = H(1) H(2) . . . H(nq-1).
30
31

ARGUMENTS

33       SIDE    (input) CHARACTER*1
34               = 'L': apply Q or Q**H from the Left;
35               = 'R': apply Q or Q**H from the Right.
36
37       UPLO    (input) CHARACTER*1
38               = 'U': Upper triangular packed storage used in previous call to
39               ZHPTRD; = 'L': Lower triangular packed storage used in previous
40               call to ZHPTRD.
41
42       TRANS   (input) CHARACTER*1
43               = 'N':  No transpose, apply Q;
44               = 'C':  Conjugate transpose, apply Q**H.
45
46       M       (input) INTEGER
47               The number of rows of the matrix C. M >= 0.
48
49       N       (input) INTEGER
50               The number of columns of the matrix C. N >= 0.
51
52       AP      (input) COMPLEX*16 array, dimension
53               (M*(M+1)/2) if SIDE = 'L' (N*(N+1)/2) if SIDE = 'R' The vectors
54               which  define the elementary reflectors, as returned by ZHPTRD.
55               AP is modified by the routine but restored on exit.
56
57       TAU     (input) COMPLEX*16 array, dimension (M-1) if SIDE = 'L'
58               or (N-1) if SIDE = 'R' TAU(i) must contain the scalar factor of
59               the elementary reflector H(i), as returned by ZHPTRD.
60
61       C       (input/output) COMPLEX*16 array, dimension (LDC,N)
62               On  entry,  the  M-by-N matrix C.  On exit, C is overwritten by
63               Q*C or Q**H*C or C*Q**H or C*Q.
64
65       LDC     (input) INTEGER
66               The leading dimension of the array C. LDC >= max(1,M).
67
68       WORK    (workspace) COMPLEX*16 array, dimension
69               (N) if SIDE = 'L' (M) if SIDE = 'R'
70
71       INFO    (output) INTEGER
72               = 0:  successful exit
73               < 0:  if INFO = -i, the i-th argument had an illegal value
74
75
76
77 LAPACK routine (version 3.1)    November 2006                       ZUPMTR(1)
Impressum