Subroutine */
| 39550 | } /* dsytd2_ */ |
| 39551 | |
| 39552 | /* Subroutine */ int dsytrd_(char *uplo, integer *n, doublereal *a, integer * |
| 39553 | lda, doublereal *d__, doublereal *e, doublereal *tau, doublereal * |
| 39554 | work, integer *lwork, integer *info) |
| 39555 | { |
| 39556 | /* System generated locals */ |
| 39557 | integer a_dim1, a_offset, i__1, i__2, i__3; |
| 39558 | |
| 39559 | /* Local variables */ |
| 39560 | static integer i__, j, nb, kk, nx, iws; |
| 39561 | extern logical lsame_(char *, char *); |
| 39562 | static integer nbmin, iinfo; |
| 39563 | static logical upper; |
| 39564 | extern /* Subroutine */ int dsytd2_(char *, integer *, doublereal *, |
| 39565 | integer *, doublereal *, doublereal *, doublereal *, integer *), dsyr2k_(char *, char *, integer *, integer *, doublereal |
| 39566 | *, doublereal *, integer *, doublereal *, integer *, doublereal *, |
| 39567 | doublereal *, integer *), dlatrd_(char *, |
| 39568 | integer *, integer *, doublereal *, integer *, doublereal *, |
| 39569 | doublereal *, doublereal *, integer *), xerbla_(char *, |
| 39570 | integer *); |
| 39571 | extern integer ilaenv_(integer *, char *, char *, integer *, integer *, |
| 39572 | integer *, integer *, ftnlen, ftnlen); |
| 39573 | static integer ldwork, lwkopt; |
| 39574 | static logical lquery; |
| 39575 | |
| 39576 | |
| 39577 | /* |
| 39578 | -- LAPACK routine (version 3.2) -- |
| 39579 | -- LAPACK is a software package provided by Univ. of Tennessee, -- |
| 39580 | -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
| 39581 | November 2006 |
| 39582 | |
| 39583 | |
| 39584 | Purpose |
| 39585 | ======= |
| 39586 | |
| 39587 | DSYTRD reduces a real symmetric matrix A to real symmetric |
| 39588 | tridiagonal form T by an orthogonal similarity transformation: |
| 39589 | Q**T * A * Q = T. |
| 39590 | |
| 39591 | Arguments |
| 39592 | ========= |
| 39593 | |
| 39594 | UPLO (input) CHARACTER*1 |
| 39595 | = 'U': Upper triangle of A is stored; |
| 39596 | = 'L': Lower triangle of A is stored. |
| 39597 | |
| 39598 | N (input) INTEGER |
| 39599 | The order of the matrix A. N >= 0. |
| 39600 | |
| 39601 | A (input/output) DOUBLE PRECISION array, dimension (LDA,N) |
| 39602 | On entry, the symmetric matrix A. If UPLO = 'U', the leading |
| 39603 | N-by-N upper triangular part of A contains the upper |
| 39604 | triangular part of the matrix A, and the strictly lower |
| 39605 | triangular part of A is not referenced. If UPLO = 'L', the |
| 39606 | leading N-by-N lower triangular part of A contains the lower |
| 39607 | triangular part of the matrix A, and the strictly upper |
| 39608 | triangular part of A is not referenced. |
| 39609 | On exit, if UPLO = 'U', the diagonal and first superdiagonal |