Subroutine */
| 7646 | } /* zhetd2_ */ |
| 7647 | |
| 7648 | /* Subroutine */ int zhetrd_(char *uplo, integer *n, doublecomplex *a, |
| 7649 | integer *lda, doublereal *d__, doublereal *e, doublecomplex *tau, |
| 7650 | doublecomplex *work, integer *lwork, integer *info) |
| 7651 | { |
| 7652 | /* System generated locals */ |
| 7653 | integer a_dim1, a_offset, i__1, i__2, i__3, i__4, i__5; |
| 7654 | doublecomplex z__1; |
| 7655 | |
| 7656 | /* Local variables */ |
| 7657 | static integer i__, j, nb, kk, nx, iws; |
| 7658 | extern logical lsame_(char *, char *); |
| 7659 | static integer nbmin, iinfo; |
| 7660 | static logical upper; |
| 7661 | extern /* Subroutine */ int zhetd2_(char *, integer *, doublecomplex *, |
| 7662 | integer *, doublereal *, doublereal *, doublecomplex *, integer *), zher2k_(char *, char *, integer *, integer *, |
| 7663 | doublecomplex *, doublecomplex *, integer *, doublecomplex *, |
| 7664 | integer *, doublereal *, doublecomplex *, integer *), xerbla_(char *, integer *); |
| 7665 | extern integer ilaenv_(integer *, char *, char *, integer *, integer *, |
| 7666 | integer *, integer *, ftnlen, ftnlen); |
| 7667 | extern /* Subroutine */ int zlatrd_(char *, integer *, integer *, |
| 7668 | doublecomplex *, integer *, doublereal *, doublecomplex *, |
| 7669 | doublecomplex *, integer *); |
| 7670 | static integer ldwork, lwkopt; |
| 7671 | static logical lquery; |
| 7672 | |
| 7673 | |
| 7674 | /* |
| 7675 | -- LAPACK routine (version 3.2) -- |
| 7676 | -- LAPACK is a software package provided by Univ. of Tennessee, -- |
| 7677 | -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
| 7678 | November 2006 |
| 7679 | |
| 7680 | |
| 7681 | Purpose |
| 7682 | ======= |
| 7683 | |
| 7684 | ZHETRD reduces a complex Hermitian matrix A to real symmetric |
| 7685 | tridiagonal form T by a unitary similarity transformation: |
| 7686 | Q**H * A * Q = T. |
| 7687 | |
| 7688 | Arguments |
| 7689 | ========= |
| 7690 | |
| 7691 | UPLO (input) CHARACTER*1 |
| 7692 | = 'U': Upper triangle of A is stored; |
| 7693 | = 'L': Lower triangle of A is stored. |
| 7694 | |
| 7695 | N (input) INTEGER |
| 7696 | The order of the matrix A. N >= 0. |
| 7697 | |
| 7698 | A (input/output) COMPLEX*16 array, dimension (LDA,N) |
| 7699 | On entry, the Hermitian matrix A. If UPLO = 'U', the leading |
| 7700 | N-by-N upper triangular part of A contains the upper |
| 7701 | triangular part of the matrix A, and the strictly lower |
| 7702 | triangular part of A is not referenced. If UPLO = 'L', the |
| 7703 | leading N-by-N lower triangular part of A contains the lower |
| 7704 | triangular part of the matrix A, and the strictly upper |
| 7705 | triangular part of A is not referenced. |