Subroutine */
| 7307 | } /* zheevd_ */ |
| 7308 | |
| 7309 | /* Subroutine */ int zhetd2_(char *uplo, integer *n, doublecomplex *a, |
| 7310 | integer *lda, doublereal *d__, doublereal *e, doublecomplex *tau, |
| 7311 | integer *info) |
| 7312 | { |
| 7313 | /* System generated locals */ |
| 7314 | integer a_dim1, a_offset, i__1, i__2, i__3; |
| 7315 | doublereal d__1; |
| 7316 | doublecomplex z__1, z__2, z__3, z__4; |
| 7317 | |
| 7318 | /* Local variables */ |
| 7319 | static integer i__; |
| 7320 | static doublecomplex taui; |
| 7321 | extern /* Subroutine */ int zher2_(char *, integer *, doublecomplex *, |
| 7322 | doublecomplex *, integer *, doublecomplex *, integer *, |
| 7323 | doublecomplex *, integer *); |
| 7324 | static doublecomplex alpha; |
| 7325 | extern logical lsame_(char *, char *); |
| 7326 | extern /* Double Complex */ VOID zdotc_(doublecomplex *, integer *, |
| 7327 | doublecomplex *, integer *, doublecomplex *, integer *); |
| 7328 | extern /* Subroutine */ int zhemv_(char *, integer *, doublecomplex *, |
| 7329 | doublecomplex *, integer *, doublecomplex *, integer *, |
| 7330 | doublecomplex *, doublecomplex *, integer *); |
| 7331 | static logical upper; |
| 7332 | extern /* Subroutine */ int zaxpy_(integer *, doublecomplex *, |
| 7333 | doublecomplex *, integer *, doublecomplex *, integer *), xerbla_( |
| 7334 | char *, integer *), zlarfg_(integer *, doublecomplex *, |
| 7335 | doublecomplex *, integer *, doublecomplex *); |
| 7336 | |
| 7337 | |
| 7338 | /* |
| 7339 | -- LAPACK routine (version 3.2) -- |
| 7340 | -- LAPACK is a software package provided by Univ. of Tennessee, -- |
| 7341 | -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
| 7342 | November 2006 |
| 7343 | |
| 7344 | |
| 7345 | Purpose |
| 7346 | ======= |
| 7347 | |
| 7348 | ZHETD2 reduces a complex Hermitian matrix A to real symmetric |
| 7349 | tridiagonal form T by a unitary similarity transformation: |
| 7350 | Q' * A * Q = T. |
| 7351 | |
| 7352 | Arguments |
| 7353 | ========= |
| 7354 | |
| 7355 | UPLO (input) CHARACTER*1 |
| 7356 | Specifies whether the upper or lower triangular part of the |
| 7357 | Hermitian matrix A is stored: |
| 7358 | = 'U': Upper triangular |
| 7359 | = 'L': Lower triangular |
| 7360 | |
| 7361 | N (input) INTEGER |
| 7362 | The order of the matrix A. N >= 0. |
| 7363 | |
| 7364 | A (input/output) COMPLEX*16 array, dimension (LDA,N) |
| 7365 | On entry, the Hermitian matrix A. If UPLO = 'U', the leading |
| 7366 | n-by-n upper triangular part of A contains the upper |