Subroutine */
| 29125 | } /* zunmlq_ */ |
| 29126 | |
| 29127 | /* Subroutine */ int zunmql_(char *side, char *trans, integer *m, integer *n, |
| 29128 | integer *k, doublecomplex *a, integer *lda, doublecomplex *tau, |
| 29129 | doublecomplex *c__, integer *ldc, doublecomplex *work, integer *lwork, |
| 29130 | integer *info) |
| 29131 | { |
| 29132 | /* System generated locals */ |
| 29133 | address a__1[2]; |
| 29134 | integer a_dim1, a_offset, c_dim1, c_offset, i__1, i__2, i__3[2], i__4, |
| 29135 | i__5; |
| 29136 | char ch__1[2]; |
| 29137 | |
| 29138 | /* Local variables */ |
| 29139 | static integer i__; |
| 29140 | static doublecomplex t[4160] /* was [65][64] */; |
| 29141 | static integer i1, i2, i3, ib, nb, mi, ni, nq, nw, iws; |
| 29142 | static logical left; |
| 29143 | extern logical lsame_(char *, char *); |
| 29144 | static integer nbmin, iinfo; |
| 29145 | extern /* Subroutine */ int zunm2l_(char *, char *, integer *, integer *, |
| 29146 | integer *, doublecomplex *, integer *, doublecomplex *, |
| 29147 | doublecomplex *, integer *, doublecomplex *, integer *), xerbla_(char *, integer *); |
| 29148 | extern integer ilaenv_(integer *, char *, char *, integer *, integer *, |
| 29149 | integer *, integer *, ftnlen, ftnlen); |
| 29150 | extern /* Subroutine */ int zlarfb_(char *, char *, char *, char *, |
| 29151 | integer *, integer *, integer *, doublecomplex *, integer *, |
| 29152 | doublecomplex *, integer *, doublecomplex *, integer *, |
| 29153 | doublecomplex *, integer *); |
| 29154 | static logical notran; |
| 29155 | static integer ldwork; |
| 29156 | extern /* Subroutine */ int zlarft_(char *, char *, integer *, integer *, |
| 29157 | doublecomplex *, integer *, doublecomplex *, doublecomplex *, |
| 29158 | integer *); |
| 29159 | static integer lwkopt; |
| 29160 | static logical lquery; |
| 29161 | |
| 29162 | |
| 29163 | /* |
| 29164 | -- LAPACK routine (version 3.2) -- |
| 29165 | -- LAPACK is a software package provided by Univ. of Tennessee, -- |
| 29166 | -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
| 29167 | November 2006 |
| 29168 | |
| 29169 | |
| 29170 | Purpose |
| 29171 | ======= |
| 29172 | |
| 29173 | ZUNMQL overwrites the general complex M-by-N matrix C with |
| 29174 | |
| 29175 | SIDE = 'L' SIDE = 'R' |
| 29176 | TRANS = 'N': Q * C C * Q |
| 29177 | TRANS = 'C': Q**H * C C * Q**H |
| 29178 | |
| 29179 | where Q is a complex unitary matrix defined as the product of k |
| 29180 | elementary reflectors |
| 29181 | |
| 29182 | Q = H(k) . . . H(2) H(1) |
| 29183 | |
| 29184 | as returned by ZGEQLF. Q is of order M if SIDE = 'L' and of order N |