Subroutine */
| 26906 | } /* zunghr_ */ |
| 26907 | |
| 26908 | /* Subroutine */ int zungl2_(integer *m, integer *n, integer *k, |
| 26909 | doublecomplex *a, integer *lda, doublecomplex *tau, doublecomplex * |
| 26910 | work, integer *info) |
| 26911 | { |
| 26912 | /* System generated locals */ |
| 26913 | integer a_dim1, a_offset, i__1, i__2, i__3; |
| 26914 | doublecomplex z__1, z__2; |
| 26915 | |
| 26916 | /* Local variables */ |
| 26917 | static integer i__, j, l; |
| 26918 | extern /* Subroutine */ int zscal_(integer *, doublecomplex *, |
| 26919 | doublecomplex *, integer *), zlarf_(char *, integer *, integer *, |
| 26920 | doublecomplex *, integer *, doublecomplex *, doublecomplex *, |
| 26921 | integer *, doublecomplex *), xerbla_(char *, integer *), zlacgv_(integer *, doublecomplex *, integer *); |
| 26922 | |
| 26923 | |
| 26924 | /* |
| 26925 | -- LAPACK routine (version 3.2) -- |
| 26926 | -- LAPACK is a software package provided by Univ. of Tennessee, -- |
| 26927 | -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
| 26928 | November 2006 |
| 26929 | |
| 26930 | |
| 26931 | Purpose |
| 26932 | ======= |
| 26933 | |
| 26934 | ZUNGL2 generates an m-by-n complex matrix Q with orthonormal rows, |
| 26935 | which is defined as the first m rows of a product of k elementary |
| 26936 | reflectors of order n |
| 26937 | |
| 26938 | Q = H(k)' . . . H(2)' H(1)' |
| 26939 | |
| 26940 | as returned by ZGELQF. |
| 26941 | |
| 26942 | Arguments |
| 26943 | ========= |
| 26944 | |
| 26945 | M (input) INTEGER |
| 26946 | The number of rows of the matrix Q. M >= 0. |
| 26947 | |
| 26948 | N (input) INTEGER |
| 26949 | The number of columns of the matrix Q. N >= M. |
| 26950 | |
| 26951 | K (input) INTEGER |
| 26952 | The number of elementary reflectors whose product defines the |
| 26953 | matrix Q. M >= K >= 0. |
| 26954 | |
| 26955 | A (input/output) COMPLEX*16 array, dimension (LDA,N) |
| 26956 | On entry, the i-th row must contain the vector which defines |
| 26957 | the elementary reflector H(i), for i = 1,2,...,k, as returned |
| 26958 | by ZGELQF in the first k rows of its array argument A. |
| 26959 | On exit, the m by n matrix Q. |
| 26960 | |
| 26961 | LDA (input) INTEGER |
| 26962 | The first dimension of the array A. LDA >= max(1,M). |
| 26963 | |
| 26964 | TAU (input) COMPLEX*16 array, dimension (K) |
| 26965 | TAU(i) must contain the scalar factor of the elementary |