Subroutine */
| 7261 | } /* cheevd_ */ |
| 7262 | |
| 7263 | /* Subroutine */ int chetd2_(char *uplo, integer *n, complex *a, integer *lda, |
| 7264 | real *d__, real *e, complex *tau, integer *info) |
| 7265 | { |
| 7266 | /* System generated locals */ |
| 7267 | integer a_dim1, a_offset, i__1, i__2, i__3; |
| 7268 | real r__1; |
| 7269 | complex q__1, q__2, q__3, q__4; |
| 7270 | |
| 7271 | /* Local variables */ |
| 7272 | static integer i__; |
| 7273 | static complex taui; |
| 7274 | extern /* Subroutine */ int cher2_(char *, integer *, complex *, complex * |
| 7275 | , integer *, complex *, integer *, complex *, integer *); |
| 7276 | static complex alpha; |
| 7277 | extern /* Complex */ VOID cdotc_(complex *, integer *, complex *, integer |
| 7278 | *, complex *, integer *); |
| 7279 | extern logical lsame_(char *, char *); |
| 7280 | extern /* Subroutine */ int chemv_(char *, integer *, complex *, complex * |
| 7281 | , integer *, complex *, integer *, complex *, complex *, integer * |
| 7282 | ), caxpy_(integer *, complex *, complex *, integer *, |
| 7283 | complex *, integer *); |
| 7284 | static logical upper; |
| 7285 | extern /* Subroutine */ int clarfg_(integer *, complex *, complex *, |
| 7286 | integer *, complex *), xerbla_(char *, integer *); |
| 7287 | |
| 7288 | |
| 7289 | /* |
| 7290 | -- LAPACK routine (version 3.2) -- |
| 7291 | -- LAPACK is a software package provided by Univ. of Tennessee, -- |
| 7292 | -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
| 7293 | November 2006 |
| 7294 | |
| 7295 | |
| 7296 | Purpose |
| 7297 | ======= |
| 7298 | |
| 7299 | CHETD2 reduces a complex Hermitian matrix A to real symmetric |
| 7300 | tridiagonal form T by a unitary similarity transformation: |
| 7301 | Q' * A * Q = T. |
| 7302 | |
| 7303 | Arguments |
| 7304 | ========= |
| 7305 | |
| 7306 | UPLO (input) CHARACTER*1 |
| 7307 | Specifies whether the upper or lower triangular part of the |
| 7308 | Hermitian matrix A is stored: |
| 7309 | = 'U': Upper triangular |
| 7310 | = 'L': Lower triangular |
| 7311 | |
| 7312 | N (input) INTEGER |
| 7313 | The order of the matrix A. N >= 0. |
| 7314 | |
| 7315 | A (input/output) COMPLEX array, dimension (LDA,N) |
| 7316 | On entry, the Hermitian matrix A. If UPLO = 'U', the leading |
| 7317 | n-by-n upper triangular part of A contains the upper |
| 7318 | triangular part of the matrix A, and the strictly lower |
| 7319 | triangular part of A is not referenced. If UPLO = 'L', the |
| 7320 | leading n-by-n lower triangular part of A contains the lower |