Subroutine */
| 12369 | } /* dlaed8_ */ |
| 12370 | |
| 12371 | /* Subroutine */ int dlaed9_(integer *k, integer *kstart, integer *kstop, |
| 12372 | integer *n, doublereal *d__, doublereal *q, integer *ldq, doublereal * |
| 12373 | rho, doublereal *dlamda, doublereal *w, doublereal *s, integer *lds, |
| 12374 | integer *info) |
| 12375 | { |
| 12376 | /* System generated locals */ |
| 12377 | integer q_dim1, q_offset, s_dim1, s_offset, i__1, i__2; |
| 12378 | doublereal d__1; |
| 12379 | |
| 12380 | /* Local variables */ |
| 12381 | static integer i__, j; |
| 12382 | static doublereal temp; |
| 12383 | extern doublereal dnrm2_(integer *, doublereal *, integer *); |
| 12384 | extern /* Subroutine */ int dcopy_(integer *, doublereal *, integer *, |
| 12385 | doublereal *, integer *), dlaed4_(integer *, integer *, |
| 12386 | doublereal *, doublereal *, doublereal *, doublereal *, |
| 12387 | doublereal *, integer *); |
| 12388 | extern doublereal dlamc3_(doublereal *, doublereal *); |
| 12389 | extern /* Subroutine */ int xerbla_(char *, integer *); |
| 12390 | |
| 12391 | |
| 12392 | /* |
| 12393 | -- LAPACK routine (version 3.2) -- |
| 12394 | -- LAPACK is a software package provided by Univ. of Tennessee, -- |
| 12395 | -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
| 12396 | November 2006 |
| 12397 | |
| 12398 | |
| 12399 | Purpose |
| 12400 | ======= |
| 12401 | |
| 12402 | DLAED9 finds the roots of the secular equation, as defined by the |
| 12403 | values in D, Z, and RHO, between KSTART and KSTOP. It makes the |
| 12404 | appropriate calls to DLAED4 and then stores the new matrix of |
| 12405 | eigenvectors for use in calculating the next level of Z vectors. |
| 12406 | |
| 12407 | Arguments |
| 12408 | ========= |
| 12409 | |
| 12410 | K (input) INTEGER |
| 12411 | The number of terms in the rational function to be solved by |
| 12412 | DLAED4. K >= 0. |
| 12413 | |
| 12414 | KSTART (input) INTEGER |
| 12415 | KSTOP (input) INTEGER |
| 12416 | The updated eigenvalues Lambda(I), KSTART <= I <= KSTOP |
| 12417 | are to be computed. 1 <= KSTART <= KSTOP <= K. |
| 12418 | |
| 12419 | N (input) INTEGER |
| 12420 | The number of rows and columns in the Q matrix. |
| 12421 | N >= K (delation may result in N > K). |
| 12422 | |
| 12423 | D (output) DOUBLE PRECISION array, dimension (N) |
| 12424 | D(I) contains the updated eigenvalues |
| 12425 | for KSTART <= I <= KSTOP. |
| 12426 | |
| 12427 | Q (workspace) DOUBLE PRECISION array, dimension (LDQ,N) |
| 12428 |