Subroutine */
| 20346 | } /* dlaqr4_ */ |
| 20347 | |
| 20348 | /* Subroutine */ int dlaqr5_(logical *wantt, logical *wantz, integer *kacc22, |
| 20349 | integer *n, integer *ktop, integer *kbot, integer *nshfts, doublereal |
| 20350 | *sr, doublereal *si, doublereal *h__, integer *ldh, integer *iloz, |
| 20351 | integer *ihiz, doublereal *z__, integer *ldz, doublereal *v, integer * |
| 20352 | ldv, doublereal *u, integer *ldu, integer *nv, doublereal *wv, |
| 20353 | integer *ldwv, integer *nh, doublereal *wh, integer *ldwh) |
| 20354 | { |
| 20355 | /* System generated locals */ |
| 20356 | integer h_dim1, h_offset, u_dim1, u_offset, v_dim1, v_offset, wh_dim1, |
| 20357 | wh_offset, wv_dim1, wv_offset, z_dim1, z_offset, i__1, i__2, i__3, |
| 20358 | i__4, i__5, i__6, i__7; |
| 20359 | doublereal d__1, d__2, d__3, d__4, d__5; |
| 20360 | |
| 20361 | /* Local variables */ |
| 20362 | static integer i__, j, k, m, i2, j2, i4, j4, k1; |
| 20363 | static doublereal h11, h12, h21, h22; |
| 20364 | static integer m22, ns, nu; |
| 20365 | static doublereal vt[3], scl; |
| 20366 | static integer kdu, kms; |
| 20367 | static doublereal ulp; |
| 20368 | static integer knz, kzs; |
| 20369 | static doublereal tst1, tst2, beta; |
| 20370 | static logical blk22, bmp22; |
| 20371 | static integer mend, jcol, jlen, jbot, mbot; |
| 20372 | static doublereal swap; |
| 20373 | static integer jtop, jrow, mtop; |
| 20374 | static doublereal alpha; |
| 20375 | static logical accum; |
| 20376 | extern /* Subroutine */ int dgemm_(char *, char *, integer *, integer *, |
| 20377 | integer *, doublereal *, doublereal *, integer *, doublereal *, |
| 20378 | integer *, doublereal *, doublereal *, integer *); |
| 20379 | static integer ndcol, incol, krcol, nbmps; |
| 20380 | extern /* Subroutine */ int dtrmm_(char *, char *, char *, char *, |
| 20381 | integer *, integer *, doublereal *, doublereal *, integer *, |
| 20382 | doublereal *, integer *), dlaqr1_( |
| 20383 | integer *, doublereal *, integer *, doublereal *, doublereal *, |
| 20384 | doublereal *, doublereal *, doublereal *), dlabad_(doublereal *, |
| 20385 | doublereal *); |
| 20386 | |
| 20387 | extern /* Subroutine */ int dlarfg_(integer *, doublereal *, doublereal *, |
| 20388 | integer *, doublereal *), dlacpy_(char *, integer *, integer *, |
| 20389 | doublereal *, integer *, doublereal *, integer *); |
| 20390 | static doublereal safmin; |
| 20391 | extern /* Subroutine */ int dlaset_(char *, integer *, integer *, |
| 20392 | doublereal *, doublereal *, doublereal *, integer *); |
| 20393 | static doublereal safmax, refsum; |
| 20394 | static integer mstart; |
| 20395 | static doublereal smlnum; |
| 20396 | |
| 20397 | |
| 20398 | /* |
| 20399 | -- LAPACK auxiliary routine (version 3.2) -- |
| 20400 | -- LAPACK is a software package provided by Univ. of Tennessee, -- |
| 20401 | -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
| 20402 | November 2006 |
| 20403 | |
| 20404 | |
| 20405 | This auxiliary subroutine called by DLAQR0 performs a |