Subroutine */
| 7470 | } /* dgetrs_ */ |
| 7471 | |
| 7472 | /* Subroutine */ int dhseqr_(char *job, char *compz, integer *n, integer *ilo, |
| 7473 | integer *ihi, doublereal *h__, integer *ldh, doublereal *wr, |
| 7474 | doublereal *wi, doublereal *z__, integer *ldz, doublereal *work, |
| 7475 | integer *lwork, integer *info) |
| 7476 | { |
| 7477 | /* System generated locals */ |
| 7478 | address a__1[2]; |
| 7479 | integer h_dim1, h_offset, z_dim1, z_offset, i__1, i__2[2], i__3; |
| 7480 | doublereal d__1; |
| 7481 | char ch__1[2]; |
| 7482 | |
| 7483 | /* Local variables */ |
| 7484 | static integer i__; |
| 7485 | static doublereal hl[2401] /* was [49][49] */; |
| 7486 | static integer kbot, nmin; |
| 7487 | extern logical lsame_(char *, char *); |
| 7488 | static logical initz; |
| 7489 | static doublereal workl[49]; |
| 7490 | static logical wantt, wantz; |
| 7491 | extern /* Subroutine */ int dlaqr0_(logical *, logical *, integer *, |
| 7492 | integer *, integer *, doublereal *, integer *, doublereal *, |
| 7493 | doublereal *, integer *, integer *, doublereal *, integer *, |
| 7494 | doublereal *, integer *, integer *), dlahqr_(logical *, logical *, |
| 7495 | integer *, integer *, integer *, doublereal *, integer *, |
| 7496 | doublereal *, doublereal *, integer *, integer *, doublereal *, |
| 7497 | integer *, integer *), dlacpy_(char *, integer *, integer *, |
| 7498 | doublereal *, integer *, doublereal *, integer *), |
| 7499 | dlaset_(char *, integer *, integer *, doublereal *, doublereal *, |
| 7500 | doublereal *, integer *); |
| 7501 | extern integer ilaenv_(integer *, char *, char *, integer *, integer *, |
| 7502 | integer *, integer *, ftnlen, ftnlen); |
| 7503 | extern /* Subroutine */ int xerbla_(char *, integer *); |
| 7504 | static logical lquery; |
| 7505 | |
| 7506 | |
| 7507 | /* |
| 7508 | -- LAPACK computational routine (version 3.2.2) -- |
| 7509 | Univ. of Tennessee, Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd.. |
| 7510 | June 2010 |
| 7511 | |
| 7512 | Purpose |
| 7513 | ======= |
| 7514 | |
| 7515 | DHSEQR computes the eigenvalues of a Hessenberg matrix H |
| 7516 | and, optionally, the matrices T and Z from the Schur decomposition |
| 7517 | H = Z T Z**T, where T is an upper quasi-triangular matrix (the |
| 7518 | Schur form), and Z is the orthogonal matrix of Schur vectors. |
| 7519 | |
| 7520 | Optionally Z may be postmultiplied into an input orthogonal |
| 7521 | matrix Q so that this routine can give the Schur factorization |
| 7522 | of a matrix A which has been reduced to the Hessenberg form H |
| 7523 | by the orthogonal matrix Q: A = Q*H*Q**T = (QZ)*T*(QZ)**T. |
| 7524 | |
| 7525 | Arguments |
| 7526 | ========= |
| 7527 | |
| 7528 | JOB (input) CHARACTER*1 |
| 7529 | = 'E': compute eigenvalues only; |