Subroutine */
| 29427 | } /* zunmql_ */ |
| 29428 | |
| 29429 | /* Subroutine */ int zunmqr_(char *side, char *trans, integer *m, integer *n, |
| 29430 | integer *k, doublecomplex *a, integer *lda, doublecomplex *tau, |
| 29431 | doublecomplex *c__, integer *ldc, doublecomplex *work, integer *lwork, |
| 29432 | integer *info) |
| 29433 | { |
| 29434 | /* System generated locals */ |
| 29435 | address a__1[2]; |
| 29436 | integer a_dim1, a_offset, c_dim1, c_offset, i__1, i__2, i__3[2], i__4, |
| 29437 | i__5; |
| 29438 | char ch__1[2]; |
| 29439 | |
| 29440 | /* Local variables */ |
| 29441 | static integer i__; |
| 29442 | static doublecomplex t[4160] /* was [65][64] */; |
| 29443 | static integer i1, i2, i3, ib, ic, jc, nb, mi, ni, nq, nw, iws; |
| 29444 | static logical left; |
| 29445 | extern logical lsame_(char *, char *); |
| 29446 | static integer nbmin, iinfo; |
| 29447 | extern /* Subroutine */ int zunm2r_(char *, char *, integer *, integer *, |
| 29448 | integer *, doublecomplex *, integer *, doublecomplex *, |
| 29449 | doublecomplex *, integer *, doublecomplex *, integer *), xerbla_(char *, integer *); |
| 29450 | extern integer ilaenv_(integer *, char *, char *, integer *, integer *, |
| 29451 | integer *, integer *, ftnlen, ftnlen); |
| 29452 | extern /* Subroutine */ int zlarfb_(char *, char *, char *, char *, |
| 29453 | integer *, integer *, integer *, doublecomplex *, integer *, |
| 29454 | doublecomplex *, integer *, doublecomplex *, integer *, |
| 29455 | doublecomplex *, integer *); |
| 29456 | static logical notran; |
| 29457 | static integer ldwork; |
| 29458 | extern /* Subroutine */ int zlarft_(char *, char *, integer *, integer *, |
| 29459 | doublecomplex *, integer *, doublecomplex *, doublecomplex *, |
| 29460 | integer *); |
| 29461 | static integer lwkopt; |
| 29462 | static logical lquery; |
| 29463 | |
| 29464 | |
| 29465 | /* |
| 29466 | -- LAPACK routine (version 3.2) -- |
| 29467 | -- LAPACK is a software package provided by Univ. of Tennessee, -- |
| 29468 | -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
| 29469 | November 2006 |
| 29470 | |
| 29471 | |
| 29472 | Purpose |
| 29473 | ======= |
| 29474 | |
| 29475 | ZUNMQR overwrites the general complex M-by-N matrix C with |
| 29476 | |
| 29477 | SIDE = 'L' SIDE = 'R' |
| 29478 | TRANS = 'N': Q * C C * Q |
| 29479 | TRANS = 'C': Q**H * C C * Q**H |
| 29480 | |
| 29481 | where Q is a complex unitary matrix defined as the product of k |
| 29482 | elementary reflectors |
| 29483 | |
| 29484 | Q = H(1) H(2) . . . H(k) |
| 29485 | |
| 29486 | as returned by ZGEQRF. Q is of order M if SIDE = 'L' and of order N |