Subroutine */
| 26591 | } /* cungbr_ */ |
| 26592 | |
| 26593 | /* Subroutine */ int cunghr_(integer *n, integer *ilo, integer *ihi, complex * |
| 26594 | a, integer *lda, complex *tau, complex *work, integer *lwork, integer |
| 26595 | *info) |
| 26596 | { |
| 26597 | /* System generated locals */ |
| 26598 | integer a_dim1, a_offset, i__1, i__2, i__3, i__4; |
| 26599 | |
| 26600 | /* Local variables */ |
| 26601 | static integer i__, j, nb, nh, iinfo; |
| 26602 | extern /* Subroutine */ int xerbla_(char *, integer *); |
| 26603 | extern integer ilaenv_(integer *, char *, char *, integer *, integer *, |
| 26604 | integer *, integer *, ftnlen, ftnlen); |
| 26605 | extern /* Subroutine */ int cungqr_(integer *, integer *, integer *, |
| 26606 | complex *, integer *, complex *, complex *, integer *, integer *); |
| 26607 | static integer lwkopt; |
| 26608 | static logical lquery; |
| 26609 | |
| 26610 | |
| 26611 | /* |
| 26612 | -- LAPACK routine (version 3.2) -- |
| 26613 | -- LAPACK is a software package provided by Univ. of Tennessee, -- |
| 26614 | -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
| 26615 | November 2006 |
| 26616 | |
| 26617 | |
| 26618 | Purpose |
| 26619 | ======= |
| 26620 | |
| 26621 | CUNGHR generates a complex unitary matrix Q which is defined as the |
| 26622 | product of IHI-ILO elementary reflectors of order N, as returned by |
| 26623 | CGEHRD: |
| 26624 | |
| 26625 | Q = H(ilo) H(ilo+1) . . . H(ihi-1). |
| 26626 | |
| 26627 | Arguments |
| 26628 | ========= |
| 26629 | |
| 26630 | N (input) INTEGER |
| 26631 | The order of the matrix Q. N >= 0. |
| 26632 | |
| 26633 | ILO (input) INTEGER |
| 26634 | IHI (input) INTEGER |
| 26635 | ILO and IHI must have the same values as in the previous call |
| 26636 | of CGEHRD. Q is equal to the unit matrix except in the |
| 26637 | submatrix Q(ilo+1:ihi,ilo+1:ihi). |
| 26638 | 1 <= ILO <= IHI <= N, if N > 0; ILO=1 and IHI=0, if N=0. |
| 26639 | |
| 26640 | A (input/output) COMPLEX array, dimension (LDA,N) |
| 26641 | On entry, the vectors which define the elementary reflectors, |
| 26642 | as returned by CGEHRD. |
| 26643 | On exit, the N-by-N unitary matrix Q. |
| 26644 | |
| 26645 | LDA (input) INTEGER |
| 26646 | The leading dimension of the array A. LDA >= max(1,N). |
| 26647 | |
| 26648 | TAU (input) COMPLEX array, dimension (N-1) |
| 26649 | TAU(i) must contain the scalar factor of the elementary |
| 26650 | reflector H(i), as returned by CGEHRD. |