Subroutine */
| 26786 | } /* cunghr_ */ |
| 26787 | |
| 26788 | /* Subroutine */ int cungl2_(integer *m, integer *n, integer *k, complex *a, |
| 26789 | integer *lda, complex *tau, complex *work, integer *info) |
| 26790 | { |
| 26791 | /* System generated locals */ |
| 26792 | integer a_dim1, a_offset, i__1, i__2, i__3; |
| 26793 | complex q__1, q__2; |
| 26794 | |
| 26795 | /* Local variables */ |
| 26796 | static integer i__, j, l; |
| 26797 | extern /* Subroutine */ int cscal_(integer *, complex *, complex *, |
| 26798 | integer *), clarf_(char *, integer *, integer *, complex *, |
| 26799 | integer *, complex *, complex *, integer *, complex *), |
| 26800 | clacgv_(integer *, complex *, integer *), xerbla_(char *, integer |
| 26801 | *); |
| 26802 | |
| 26803 | |
| 26804 | /* |
| 26805 | -- LAPACK routine (version 3.2) -- |
| 26806 | -- LAPACK is a software package provided by Univ. of Tennessee, -- |
| 26807 | -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
| 26808 | November 2006 |
| 26809 | |
| 26810 | |
| 26811 | Purpose |
| 26812 | ======= |
| 26813 | |
| 26814 | CUNGL2 generates an m-by-n complex matrix Q with orthonormal rows, |
| 26815 | which is defined as the first m rows of a product of k elementary |
| 26816 | reflectors of order n |
| 26817 | |
| 26818 | Q = H(k)' . . . H(2)' H(1)' |
| 26819 | |
| 26820 | as returned by CGELQF. |
| 26821 | |
| 26822 | Arguments |
| 26823 | ========= |
| 26824 | |
| 26825 | M (input) INTEGER |
| 26826 | The number of rows of the matrix Q. M >= 0. |
| 26827 | |
| 26828 | N (input) INTEGER |
| 26829 | The number of columns of the matrix Q. N >= M. |
| 26830 | |
| 26831 | K (input) INTEGER |
| 26832 | The number of elementary reflectors whose product defines the |
| 26833 | matrix Q. M >= K >= 0. |
| 26834 | |
| 26835 | A (input/output) COMPLEX array, dimension (LDA,N) |
| 26836 | On entry, the i-th row must contain the vector which defines |
| 26837 | the elementary reflector H(i), for i = 1,2,...,k, as returned |
| 26838 | by CGELQF in the first k rows of its array argument A. |
| 26839 | On exit, the m by n matrix Q. |
| 26840 | |
| 26841 | LDA (input) INTEGER |
| 26842 | The first dimension of the array A. LDA >= max(1,M). |
| 26843 | |
| 26844 | TAU (input) COMPLEX array, dimension (K) |
| 26845 | TAU(i) must contain the scalar factor of the elementary |