Subroutine */
| 3555 | } /* cgeqr2_ */ |
| 3556 | |
| 3557 | /* Subroutine */ int cgeqrf_(integer *m, integer *n, complex *a, integer *lda, |
| 3558 | complex *tau, complex *work, integer *lwork, integer *info) |
| 3559 | { |
| 3560 | /* System generated locals */ |
| 3561 | integer a_dim1, a_offset, i__1, i__2, i__3, i__4; |
| 3562 | |
| 3563 | /* Local variables */ |
| 3564 | static integer i__, k, ib, nb, nx, iws, nbmin, iinfo; |
| 3565 | extern /* Subroutine */ int cgeqr2_(integer *, integer *, complex *, |
| 3566 | integer *, complex *, complex *, integer *), clarfb_(char *, char |
| 3567 | *, char *, char *, integer *, integer *, integer *, complex *, |
| 3568 | integer *, complex *, integer *, complex *, integer *, complex *, |
| 3569 | integer *), clarft_(char *, char * |
| 3570 | , integer *, integer *, complex *, integer *, complex *, complex * |
| 3571 | , integer *), xerbla_(char *, integer *); |
| 3572 | extern integer ilaenv_(integer *, char *, char *, integer *, integer *, |
| 3573 | integer *, integer *, ftnlen, ftnlen); |
| 3574 | static integer ldwork, lwkopt; |
| 3575 | static logical lquery; |
| 3576 | |
| 3577 | |
| 3578 | /* |
| 3579 | -- LAPACK routine (version 3.2) -- |
| 3580 | -- LAPACK is a software package provided by Univ. of Tennessee, -- |
| 3581 | -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
| 3582 | November 2006 |
| 3583 | |
| 3584 | |
| 3585 | Purpose |
| 3586 | ======= |
| 3587 | |
| 3588 | CGEQRF computes a QR factorization of a complex M-by-N matrix A: |
| 3589 | A = Q * R. |
| 3590 | |
| 3591 | Arguments |
| 3592 | ========= |
| 3593 | |
| 3594 | M (input) INTEGER |
| 3595 | The number of rows of the matrix A. M >= 0. |
| 3596 | |
| 3597 | N (input) INTEGER |
| 3598 | The number of columns of the matrix A. N >= 0. |
| 3599 | |
| 3600 | A (input/output) COMPLEX array, dimension (LDA,N) |
| 3601 | On entry, the M-by-N matrix A. |
| 3602 | On exit, the elements on and above the diagonal of the array |
| 3603 | contain the min(M,N)-by-N upper trapezoidal matrix R (R is |
| 3604 | upper triangular if m >= n); the elements below the diagonal, |
| 3605 | with the array TAU, represent the unitary matrix Q as a |
| 3606 | product of min(m,n) elementary reflectors (see Further |
| 3607 | Details). |
| 3608 | |
| 3609 | LDA (input) INTEGER |
| 3610 | The leading dimension of the array A. LDA >= max(1,M). |
| 3611 | |
| 3612 | TAU (output) COMPLEX array, dimension (min(M,N)) |
| 3613 | The scalar factors of the elementary reflectors (see Further |
| 3614 | Details). |