CUNGHR(l)		LAPACK routine (version	1.1)		    CUNGHR(l)

NAME
  CUNGHR - generate a complex unitary matrix Q which is	defined	as the pro-
  duct of IHI-ILO elementary reflectors	of order N, as returned	by CGEHRD

SYNOPSIS

  SUBROUTINE CUNGHR( N,	ILO, IHI, A, LDA, TAU, WORK, LWORK, INFO )

      INTEGER	     IHI, ILO, INFO, LDA, LWORK, N

      COMPLEX	     A(	LDA, * ), TAU( * ), WORK( LWORK	)

PURPOSE
  CUNGHR generates a complex unitary matrix Q which is defined as the product
  of IHI-ILO elementary	reflectors of order N, as returned by CGEHRD:

  Q = H(ilo) H(ilo+1) .	. . H(ihi-1).

ARGUMENTS

  N	  (input) INTEGER
	  The order of the matrix Q. N >= 0.

  ILO	  (input) INTEGER
	  IHI	  (input) INTEGER ILO and IHI must have	the same values	as in
	  the previous call of CGEHRD. Q is equal to the unit matrix except
	  in the submatrix Q(ilo+1:ihi,ilo+1:ihi).  If N > 0, 1	<= ILO <= IHI
	  <= N;	otherwise ILO =	1 and IHI = N.

  A	  (input/output) COMPLEX array,	dimension (LDA,N)
	  On entry, the	vectors	which define the elementary reflectors,	as
	  returned by CGEHRD.  On exit,	the N-by-N unitary matrix Q.

  LDA	  (input) INTEGER
	  The leading dimension	of the array A.	LDA >= max(1,N).

  TAU	  (input) COMPLEX array, dimension (N-1)
	  TAU(i) must contain the scalar factor	of the elementary reflector
	  H(i),	as returned by CGEHRD.

  WORK	  (workspace) COMPLEX array, dimension (LWORK)
	  On exit, if INFO = 0,	WORK(1)	returns	the optimal LWORK.

  LWORK	  (input) INTEGER
	  The dimension	of the array WORK. LWORK >= IHI-ILO.  For optimum
	  performance LWORK >= (IHI-ILO)*NB, where NB is the optimal block-
	  size.

  INFO	  (output) INTEGER
	  = 0:	successful exit
	  < 0:	if INFO	= -i, the i-th argument	had an illegal value