Subroutine */
| 3417 | } /* cgelsd_ */ |
| 3418 | |
| 3419 | /* Subroutine */ int cgeqr2_(integer *m, integer *n, complex *a, integer *lda, |
| 3420 | complex *tau, complex *work, integer *info) |
| 3421 | { |
| 3422 | /* System generated locals */ |
| 3423 | integer a_dim1, a_offset, i__1, i__2, i__3; |
| 3424 | complex q__1; |
| 3425 | |
| 3426 | /* Local variables */ |
| 3427 | static integer i__, k; |
| 3428 | static complex alpha; |
| 3429 | extern /* Subroutine */ int clarf_(char *, integer *, integer *, complex * |
| 3430 | , integer *, complex *, complex *, integer *, complex *), |
| 3431 | clarfg_(integer *, complex *, complex *, integer *, complex *), |
| 3432 | xerbla_(char *, integer *); |
| 3433 | |
| 3434 | |
| 3435 | /* |
| 3436 | -- LAPACK routine (version 3.2.2) -- |
| 3437 | -- LAPACK is a software package provided by Univ. of Tennessee, -- |
| 3438 | -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
| 3439 | June 2010 |
| 3440 | |
| 3441 | |
| 3442 | Purpose |
| 3443 | ======= |
| 3444 | |
| 3445 | CGEQR2 computes a QR factorization of a complex m by n matrix A: |
| 3446 | A = Q * R. |
| 3447 | |
| 3448 | Arguments |
| 3449 | ========= |
| 3450 | |
| 3451 | M (input) INTEGER |
| 3452 | The number of rows of the matrix A. M >= 0. |
| 3453 | |
| 3454 | N (input) INTEGER |
| 3455 | The number of columns of the matrix A. N >= 0. |
| 3456 | |
| 3457 | A (input/output) COMPLEX array, dimension (LDA,N) |
| 3458 | On entry, the m by n matrix A. |
| 3459 | On exit, the elements on and above the diagonal of the array |
| 3460 | contain the min(m,n) by n upper trapezoidal matrix R (R is |
| 3461 | upper triangular if m >= n); the elements below the diagonal, |
| 3462 | with the array TAU, represent the unitary matrix Q as a |
| 3463 | product of elementary reflectors (see Further Details). |
| 3464 | |
| 3465 | LDA (input) INTEGER |
| 3466 | The leading dimension of the array A. LDA >= max(1,M). |
| 3467 | |
| 3468 | TAU (output) COMPLEX array, dimension (min(M,N)) |
| 3469 | The scalar factors of the elementary reflectors (see Further |
| 3470 | Details). |
| 3471 | |
| 3472 | WORK (workspace) COMPLEX array, dimension (N) |
| 3473 | |
| 3474 | INFO (output) INTEGER |
| 3475 | = 0: successful exit |
| 3476 | < 0: if INFO = -i, the i-th argument had an illegal value |