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

NAME

6       ZUNGTR  - a complex unitary matrix Q which is defined as the product of
7       n-1 elementary reflectors of order N, as returned by ZHETRD
8

SYNOPSIS

10       SUBROUTINE ZUNGTR( UPLO, N, A, LDA, TAU, WORK, LWORK, INFO )
11
12           CHARACTER      UPLO
13
14           INTEGER        INFO, LDA, LWORK, N
15
16           COMPLEX*16     A( LDA, * ), TAU( * ), WORK( * )
17

PURPOSE

19       ZUNGTR generates a complex unitary matrix Q which  is  defined  as  the
20       product of n-1 elementary reflectors of order N, as returned by ZHETRD:
21
22       if UPLO = 'U', Q = H(n-1) . . . H(2) H(1),
23
24       if UPLO = 'L', Q = H(1) H(2) . . . H(n-1).
25
26

ARGUMENTS

28       UPLO    (input) CHARACTER*1
29               =  'U': Upper triangle of A contains elementary reflectors from
30               ZHETRD; = 'L': Lower triangle of A contains elementary  reflec‐
31               tors from ZHETRD.
32
33       N       (input) INTEGER
34               The order of the matrix Q. N >= 0.
35
36       A       (input/output) COMPLEX*16 array, dimension (LDA,N)
37               On  entry,  the vectors which define the elementary reflectors,
38               as returned by ZHETRD.  On exit, the N-by-N unitary matrix Q.
39
40       LDA     (input) INTEGER
41               The leading dimension of the array A. LDA >= N.
42
43       TAU     (input) COMPLEX*16 array, dimension (N-1)
44               TAU(i) must contain the scalar factor of the elementary reflec‐
45               tor H(i), as returned by ZHETRD.
46
47       WORK    (workspace/output) COMPLEX*16 array, dimension (MAX(1,LWORK))
48               On exit, if INFO = 0, WORK(1) returns the optimal LWORK.
49
50       LWORK   (input) INTEGER
51               The  dimension  of  the  array WORK. LWORK >= N-1.  For optimum
52               performance LWORK >= (N-1)*NB, where NB is the  optimal  block‐
53               size.
54
55               If  LWORK  = -1, then a workspace query is assumed; the routine
56               only calculates the optimal size of  the  WORK  array,  returns
57               this  value  as the first entry of the WORK array, and no error
58               message related to LWORK is issued by XERBLA.
59
60       INFO    (output) INTEGER
61               = 0:  successful exit
62               < 0:  if INFO = -i, the i-th argument had an illegal value
63
64
65
66 LAPACK routine (version 3.1)    November 2006                       ZUNGTR(1)
Impressum