Subroutine */
| 17339 | } /* dlapy3_ */ |
| 17340 | |
| 17341 | /* Subroutine */ int dlaqr0_(logical *wantt, logical *wantz, integer *n, |
| 17342 | integer *ilo, integer *ihi, doublereal *h__, integer *ldh, doublereal |
| 17343 | *wr, doublereal *wi, integer *iloz, integer *ihiz, doublereal *z__, |
| 17344 | integer *ldz, doublereal *work, integer *lwork, integer *info) |
| 17345 | { |
| 17346 | /* System generated locals */ |
| 17347 | integer h_dim1, h_offset, z_dim1, z_offset, i__1, i__2, i__3, i__4, i__5; |
| 17348 | doublereal d__1, d__2, d__3, d__4; |
| 17349 | |
| 17350 | /* Local variables */ |
| 17351 | static integer i__, k; |
| 17352 | static doublereal aa, bb, cc, dd; |
| 17353 | static integer ld; |
| 17354 | static doublereal cs; |
| 17355 | static integer nh, it, ks, kt; |
| 17356 | static doublereal sn; |
| 17357 | static integer ku, kv, ls, ns; |
| 17358 | static doublereal ss; |
| 17359 | static integer nw, inf, kdu, nho, nve, kwh, nsr, nwr, kwv, ndec, ndfl, |
| 17360 | kbot, nmin; |
| 17361 | static doublereal swap; |
| 17362 | static integer ktop; |
| 17363 | static doublereal zdum[1] /* was [1][1] */; |
| 17364 | static integer kacc22, itmax, nsmax, nwmax, kwtop; |
| 17365 | extern /* Subroutine */ int dlanv2_(doublereal *, doublereal *, |
| 17366 | doublereal *, doublereal *, doublereal *, doublereal *, |
| 17367 | doublereal *, doublereal *, doublereal *, doublereal *), dlaqr3_( |
| 17368 | logical *, logical *, integer *, integer *, integer *, integer *, |
| 17369 | doublereal *, integer *, integer *, integer *, doublereal *, |
| 17370 | integer *, integer *, integer *, doublereal *, doublereal *, |
| 17371 | doublereal *, integer *, integer *, doublereal *, integer *, |
| 17372 | integer *, doublereal *, integer *, doublereal *, integer *), |
| 17373 | dlaqr4_(logical *, logical *, integer *, integer *, integer *, |
| 17374 | doublereal *, integer *, doublereal *, doublereal *, integer *, |
| 17375 | integer *, doublereal *, integer *, doublereal *, integer *, |
| 17376 | integer *), dlaqr5_(logical *, logical *, integer *, integer *, |
| 17377 | integer *, integer *, integer *, doublereal *, doublereal *, |
| 17378 | doublereal *, integer *, integer *, integer *, doublereal *, |
| 17379 | integer *, doublereal *, integer *, doublereal *, integer *, |
| 17380 | integer *, doublereal *, integer *, integer *, doublereal *, |
| 17381 | integer *); |
| 17382 | static integer nibble; |
| 17383 | extern /* Subroutine */ int dlahqr_(logical *, logical *, integer *, |
| 17384 | integer *, integer *, doublereal *, integer *, doublereal *, |
| 17385 | doublereal *, integer *, integer *, doublereal *, integer *, |
| 17386 | integer *), dlacpy_(char *, integer *, integer *, doublereal *, |
| 17387 | integer *, doublereal *, integer *); |
| 17388 | extern integer ilaenv_(integer *, char *, char *, integer *, integer *, |
| 17389 | integer *, integer *, ftnlen, ftnlen); |
| 17390 | static char jbcmpz[2]; |
| 17391 | static integer nwupbd; |
| 17392 | static logical sorted; |
| 17393 | static integer lwkopt; |
| 17394 | |
| 17395 | |
| 17396 | /* |
| 17397 | -- LAPACK auxiliary routine (version 3.2) -- |
| 17398 | Univ. of Tennessee, Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd.. |