Subroutine */
| 8405 | } /* chseqr_ */ |
| 8406 | |
| 8407 | /* Subroutine */ int clabrd_(integer *m, integer *n, integer *nb, complex *a, |
| 8408 | integer *lda, real *d__, real *e, complex *tauq, complex *taup, |
| 8409 | complex *x, integer *ldx, complex *y, integer *ldy) |
| 8410 | { |
| 8411 | /* System generated locals */ |
| 8412 | integer a_dim1, a_offset, x_dim1, x_offset, y_dim1, y_offset, i__1, i__2, |
| 8413 | i__3; |
| 8414 | complex q__1; |
| 8415 | |
| 8416 | /* Local variables */ |
| 8417 | static integer i__; |
| 8418 | static complex alpha; |
| 8419 | extern /* Subroutine */ int cscal_(integer *, complex *, complex *, |
| 8420 | integer *), cgemv_(char *, integer *, integer *, complex *, |
| 8421 | complex *, integer *, complex *, integer *, complex *, complex *, |
| 8422 | integer *), clarfg_(integer *, complex *, complex *, |
| 8423 | integer *, complex *), clacgv_(integer *, complex *, integer *); |
| 8424 | |
| 8425 | |
| 8426 | /* |
| 8427 | -- LAPACK auxiliary routine (version 3.2) -- |
| 8428 | -- LAPACK is a software package provided by Univ. of Tennessee, -- |
| 8429 | -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
| 8430 | November 2006 |
| 8431 | |
| 8432 | |
| 8433 | Purpose |
| 8434 | ======= |
| 8435 | |
| 8436 | CLABRD reduces the first NB rows and columns of a complex general |
| 8437 | m by n matrix A to upper or lower real bidiagonal form by a unitary |
| 8438 | transformation Q' * A * P, and returns the matrices X and Y which |
| 8439 | are needed to apply the transformation to the unreduced part of A. |
| 8440 | |
| 8441 | If m >= n, A is reduced to upper bidiagonal form; if m < n, to lower |
| 8442 | bidiagonal form. |
| 8443 | |
| 8444 | This is an auxiliary routine called by CGEBRD |
| 8445 | |
| 8446 | Arguments |
| 8447 | ========= |
| 8448 | |
| 8449 | M (input) INTEGER |
| 8450 | The number of rows in the matrix A. |
| 8451 | |
| 8452 | N (input) INTEGER |
| 8453 | The number of columns in the matrix A. |
| 8454 | |
| 8455 | NB (input) INTEGER |
| 8456 | The number of leading rows and columns of A to be reduced. |
| 8457 | |
| 8458 | A (input/output) COMPLEX array, dimension (LDA,N) |
| 8459 | On entry, the m by n general matrix to be reduced. |
| 8460 | On exit, the first NB rows and columns of the matrix are |
| 8461 | overwritten; the rest of the array is unchanged. |
| 8462 | If m >= n, elements on and below the diagonal in the first NB |
| 8463 | columns, with the array TAUQ, represent the unitary |
| 8464 | matrix Q as a product of elementary reflectors; and |