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

NAME

6       CUNGQL - an M-by-N complex matrix Q with orthonormal columns,
7

SYNOPSIS

9       SUBROUTINE CUNGQL( M, N, K, A, LDA, TAU, WORK, LWORK, INFO )
10
11           INTEGER        INFO, K, LDA, LWORK, M, N
12
13           COMPLEX        A( LDA, * ), TAU( * ), WORK( * )
14

PURPOSE

16       CUNGQL  generates  an M-by-N complex matrix Q with orthonormal columns,
17       which is defined as the last N columns of a  product  of  K  elementary
18       reflectors of order M
19
20             Q  =  H(k) . . . H(2) H(1)
21
22       as returned by CGEQLF.
23
24

ARGUMENTS

26       M       (input) INTEGER
27               The number of rows of the matrix Q. M >= 0.
28
29       N       (input) INTEGER
30               The number of columns of the matrix Q. M >= N >= 0.
31
32       K       (input) INTEGER
33               The  number  of elementary reflectors whose product defines the
34               matrix Q. N >= K >= 0.
35
36       A       (input/output) COMPLEX array, dimension (LDA,N)
37               On entry, the (n-k+i)-th column must contain the  vector  which
38               defines  the  elementary  reflector H(i), for i = 1,2,...,k, as
39               returned by CGEQLF in the last k columns of its array  argument
40               A.  On exit, the M-by-N matrix Q.
41
42       LDA     (input) INTEGER
43               The first dimension of the array A. LDA >= max(1,M).
44
45       TAU     (input) COMPLEX array, dimension (K)
46               TAU(i) must contain the scalar factor of the elementary reflec‐
47               tor H(i), as returned by CGEQLF.
48
49       WORK    (workspace/output) COMPLEX array, dimension (MAX(1,LWORK))
50               On exit, if INFO = 0, WORK(1) returns the optimal LWORK.
51
52       LWORK   (input) INTEGER
53               The dimension of the array WORK. LWORK >= max(1,N).  For  opti‐
54               mum  performance  LWORK >= N*NB, where NB is the optimal block‐
55               size.
56
57               If LWORK = -1, then a workspace query is assumed;  the  routine
58               only  calculates  the  optimal  size of the WORK array, returns
59               this value as the first entry of the WORK array, and  no  error
60               message related to LWORK is issued by XERBLA.
61
62       INFO    (output) INTEGER
63               = 0:  successful exit
64               < 0:  if INFO = -i, the i-th argument has an illegal value
65
66
67
68 LAPACK routine (version 3.1)    November 2006                       CUNGQL(1)
Impressum