Subroutine */
| 7597 | } /* chetd2_ */ |
| 7598 | |
| 7599 | /* Subroutine */ int chetrd_(char *uplo, integer *n, complex *a, integer *lda, |
| 7600 | real *d__, real *e, complex *tau, complex *work, integer *lwork, |
| 7601 | integer *info) |
| 7602 | { |
| 7603 | /* System generated locals */ |
| 7604 | integer a_dim1, a_offset, i__1, i__2, i__3, i__4, i__5; |
| 7605 | complex q__1; |
| 7606 | |
| 7607 | /* Local variables */ |
| 7608 | static integer i__, j, nb, kk, nx, iws; |
| 7609 | extern logical lsame_(char *, char *); |
| 7610 | static integer nbmin, iinfo; |
| 7611 | static logical upper; |
| 7612 | extern /* Subroutine */ int chetd2_(char *, integer *, complex *, integer |
| 7613 | *, real *, real *, complex *, integer *), cher2k_(char *, |
| 7614 | char *, integer *, integer *, complex *, complex *, integer *, |
| 7615 | complex *, integer *, real *, complex *, integer *), clatrd_(char *, integer *, integer *, complex *, integer |
| 7616 | *, real *, complex *, complex *, integer *), xerbla_(char |
| 7617 | *, integer *); |
| 7618 | extern integer ilaenv_(integer *, char *, char *, integer *, integer *, |
| 7619 | integer *, integer *, ftnlen, ftnlen); |
| 7620 | static integer ldwork, lwkopt; |
| 7621 | static logical lquery; |
| 7622 | |
| 7623 | |
| 7624 | /* |
| 7625 | -- LAPACK routine (version 3.2) -- |
| 7626 | -- LAPACK is a software package provided by Univ. of Tennessee, -- |
| 7627 | -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
| 7628 | November 2006 |
| 7629 | |
| 7630 | |
| 7631 | Purpose |
| 7632 | ======= |
| 7633 | |
| 7634 | CHETRD reduces a complex Hermitian matrix A to real symmetric |
| 7635 | tridiagonal form T by a unitary similarity transformation: |
| 7636 | Q**H * A * Q = T. |
| 7637 | |
| 7638 | Arguments |
| 7639 | ========= |
| 7640 | |
| 7641 | UPLO (input) CHARACTER*1 |
| 7642 | = 'U': Upper triangle of A is stored; |
| 7643 | = 'L': Lower triangle of A is stored. |
| 7644 | |
| 7645 | N (input) INTEGER |
| 7646 | The order of the matrix A. N >= 0. |
| 7647 | |
| 7648 | A (input/output) COMPLEX array, dimension (LDA,N) |
| 7649 | On entry, the Hermitian matrix A. If UPLO = 'U', the leading |
| 7650 | N-by-N upper triangular part of A contains the upper |
| 7651 | triangular part of the matrix A, and the strictly lower |
| 7652 | triangular part of A is not referenced. If UPLO = 'L', the |
| 7653 | leading N-by-N lower triangular part of A contains the lower |
| 7654 | triangular part of the matrix A, and the strictly upper |
| 7655 | triangular part of A is not referenced. |
| 7656 | On exit, if UPLO = 'U', the diagonal and first superdiagonal |