Subroutine */
| 3577 | } /* zgeqr2_ */ |
| 3578 | |
| 3579 | /* Subroutine */ int zgeqrf_(integer *m, integer *n, doublecomplex *a, |
| 3580 | integer *lda, doublecomplex *tau, doublecomplex *work, integer *lwork, |
| 3581 | integer *info) |
| 3582 | { |
| 3583 | /* System generated locals */ |
| 3584 | integer a_dim1, a_offset, i__1, i__2, i__3, i__4; |
| 3585 | |
| 3586 | /* Local variables */ |
| 3587 | static integer i__, k, ib, nb, nx, iws, nbmin, iinfo; |
| 3588 | extern /* Subroutine */ int zgeqr2_(integer *, integer *, doublecomplex *, |
| 3589 | integer *, doublecomplex *, doublecomplex *, integer *), xerbla_( |
| 3590 | char *, integer *); |
| 3591 | extern integer ilaenv_(integer *, char *, char *, integer *, integer *, |
| 3592 | integer *, integer *, ftnlen, ftnlen); |
| 3593 | extern /* Subroutine */ int zlarfb_(char *, char *, char *, char *, |
| 3594 | integer *, integer *, integer *, doublecomplex *, integer *, |
| 3595 | doublecomplex *, integer *, doublecomplex *, integer *, |
| 3596 | doublecomplex *, integer *); |
| 3597 | static integer ldwork; |
| 3598 | extern /* Subroutine */ int zlarft_(char *, char *, integer *, integer *, |
| 3599 | doublecomplex *, integer *, doublecomplex *, doublecomplex *, |
| 3600 | integer *); |
| 3601 | static integer lwkopt; |
| 3602 | static logical lquery; |
| 3603 | |
| 3604 | |
| 3605 | /* |
| 3606 | -- LAPACK routine (version 3.2) -- |
| 3607 | -- LAPACK is a software package provided by Univ. of Tennessee, -- |
| 3608 | -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
| 3609 | November 2006 |
| 3610 | |
| 3611 | |
| 3612 | Purpose |
| 3613 | ======= |
| 3614 | |
| 3615 | ZGEQRF computes a QR factorization of a complex M-by-N matrix A: |
| 3616 | A = Q * R. |
| 3617 | |
| 3618 | Arguments |
| 3619 | ========= |
| 3620 | |
| 3621 | M (input) INTEGER |
| 3622 | The number of rows of the matrix A. M >= 0. |
| 3623 | |
| 3624 | N (input) INTEGER |
| 3625 | The number of columns of the matrix A. N >= 0. |
| 3626 | |
| 3627 | A (input/output) COMPLEX*16 array, dimension (LDA,N) |
| 3628 | On entry, the M-by-N matrix A. |
| 3629 | On exit, the elements on and above the diagonal of the array |
| 3630 | contain the min(M,N)-by-N upper trapezoidal matrix R (R is |
| 3631 | upper triangular if m >= n); the elements below the diagonal, |
| 3632 | with the array TAU, represent the unitary matrix Q as a |
| 3633 | product of min(m,n) elementary reflectors (see Further |
| 3634 | Details). |
| 3635 | |
| 3636 | LDA (input) INTEGER |