1CUNGTR(1) LAPACK routine (version 3.1) CUNGTR(1)
2
3
4
6 CUNGTR - a complex unitary matrix Q which is defined as the product of
7 n-1 elementary reflectors of order N, as returned by CHETRD
8
10 SUBROUTINE CUNGTR( UPLO, N, A, LDA, TAU, WORK, LWORK, INFO )
11
12 CHARACTER UPLO
13
14 INTEGER INFO, LDA, LWORK, N
15
16 COMPLEX A( LDA, * ), TAU( * ), WORK( * )
17
19 CUNGTR generates a complex unitary matrix Q which is defined as the
20 product of n-1 elementary reflectors of order N, as returned by CHETRD:
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
28 UPLO (input) CHARACTER*1
29 = 'U': Upper triangle of A contains elementary reflectors from
30 CHETRD; = 'L': Lower triangle of A contains elementary reflec‐
31 tors from CHETRD.
32
33 N (input) INTEGER
34 The order of the matrix Q. N >= 0.
35
36 A (input/output) COMPLEX array, dimension (LDA,N)
37 On entry, the vectors which define the elementary reflectors,
38 as returned by CHETRD. 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 array, dimension (N-1)
44 TAU(i) must contain the scalar factor of the elementary reflec‐
45 tor H(i), as returned by CHETRD.
46
47 WORK (workspace/output) COMPLEX 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 CUNGTR(1)