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

NAME

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

SYNOPSIS

10       SUBROUTINE CUPMTR( 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        AP( * ), C( LDC, * ), TAU( * ), WORK( * )
18

PURPOSE

20       CUPMTR 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 CHPTRD 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               CHPTRD; = 'L': Lower triangular packed storage used in previous
40               call to CHPTRD.
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 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 CHPTRD.
55               AP is modified by the routine but restored on exit.
56
57       TAU     (input) COMPLEX 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 CHPTRD.
60
61       C       (input/output) COMPLEX 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 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                       CUPMTR(1)
Impressum