Subroutine */
| 8458 | } /* zhseqr_ */ |
| 8459 | |
| 8460 | /* Subroutine */ int zlabrd_(integer *m, integer *n, integer *nb, |
| 8461 | doublecomplex *a, integer *lda, doublereal *d__, doublereal *e, |
| 8462 | doublecomplex *tauq, doublecomplex *taup, doublecomplex *x, integer * |
| 8463 | ldx, doublecomplex *y, integer *ldy) |
| 8464 | { |
| 8465 | /* System generated locals */ |
| 8466 | integer a_dim1, a_offset, x_dim1, x_offset, y_dim1, y_offset, i__1, i__2, |
| 8467 | i__3; |
| 8468 | doublecomplex z__1; |
| 8469 | |
| 8470 | /* Local variables */ |
| 8471 | static integer i__; |
| 8472 | static doublecomplex alpha; |
| 8473 | extern /* Subroutine */ int zscal_(integer *, doublecomplex *, |
| 8474 | doublecomplex *, integer *), zgemv_(char *, integer *, integer *, |
| 8475 | doublecomplex *, doublecomplex *, integer *, doublecomplex *, |
| 8476 | integer *, doublecomplex *, doublecomplex *, integer *), |
| 8477 | zlarfg_(integer *, doublecomplex *, doublecomplex *, integer *, |
| 8478 | doublecomplex *), zlacgv_(integer *, doublecomplex *, integer *); |
| 8479 | |
| 8480 | |
| 8481 | /* |
| 8482 | -- LAPACK auxiliary routine (version 3.2) -- |
| 8483 | -- LAPACK is a software package provided by Univ. of Tennessee, -- |
| 8484 | -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
| 8485 | November 2006 |
| 8486 | |
| 8487 | |
| 8488 | Purpose |
| 8489 | ======= |
| 8490 | |
| 8491 | ZLABRD reduces the first NB rows and columns of a complex general |
| 8492 | m by n matrix A to upper or lower real bidiagonal form by a unitary |
| 8493 | transformation Q' * A * P, and returns the matrices X and Y which |
| 8494 | are needed to apply the transformation to the unreduced part of A. |
| 8495 | |
| 8496 | If m >= n, A is reduced to upper bidiagonal form; if m < n, to lower |
| 8497 | bidiagonal form. |
| 8498 | |
| 8499 | This is an auxiliary routine called by ZGEBRD |
| 8500 | |
| 8501 | Arguments |
| 8502 | ========= |
| 8503 | |
| 8504 | M (input) INTEGER |
| 8505 | The number of rows in the matrix A. |
| 8506 | |
| 8507 | N (input) INTEGER |
| 8508 | The number of columns in the matrix A. |
| 8509 | |
| 8510 | NB (input) INTEGER |
| 8511 | The number of leading rows and columns of A to be reduced. |
| 8512 | |
| 8513 | A (input/output) COMPLEX*16 array, dimension (LDA,N) |
| 8514 | On entry, the m by n general matrix to be reduced. |
| 8515 | On exit, the first NB rows and columns of the matrix are |
| 8516 | overwritten; the rest of the array is unchanged. |
| 8517 | If m >= n, elements on and below the diagonal in the first NB |