1DLAED9(1)                LAPACK routine (version 3.2)                DLAED9(1)
2
3
4

NAME

6       DLAED9  -  finds  the  roots of the secular equation, as defined by the
7       values in D, Z, and RHO, between KSTART and KSTOP
8

SYNOPSIS

10       SUBROUTINE DLAED9( K, KSTART, KSTOP, N, D, Q, LDQ, RHO, DLAMDA,  W,  S,
11                          LDS, INFO )
12
13           INTEGER        INFO, K, KSTART, KSTOP, LDQ, LDS, N
14
15           DOUBLE         PRECISION RHO
16
17           DOUBLE         PRECISION  D( * ), DLAMDA( * ), Q( LDQ, * ), S( LDS,
18                          * ), W( * )
19

PURPOSE

21       DLAED9 finds the roots of the secular equation, as defined by the  val‐
22       ues in D, Z, and RHO, between KSTART and KSTOP.  It makes the appropri‐
23       ate calls to DLAED4 and then stores the new matrix of eigenvectors  for
24       use in calculating the next level of Z vectors.
25

ARGUMENTS

27       K       (input) INTEGER
28               The  number  of  terms in the rational function to be solved by
29               DLAED4.  K >= 0.
30
31       KSTART  (input) INTEGER
32               KSTOP   (input)  INTEGER  The  updated  eigenvalues  Lambda(I),
33               KSTART  <= I <= KSTOP are to be computed.  1 <= KSTART <= KSTOP
34               <= K.
35
36       N       (input) INTEGER
37               The number of rows and columns in the Q matrix.  N >= K  (dela‐
38               tion may result in N > K).
39
40       D       (output) DOUBLE PRECISION array, dimension (N)
41               D(I) contains the updated eigenvalues for KSTART <= I <= KSTOP.
42
43       Q       (workspace) DOUBLE PRECISION array, dimension (LDQ,N)
44
45       LDQ     (input) INTEGER
46               The leading dimension of the array Q.  LDQ >= max( 1, N ).
47
48       RHO     (input) DOUBLE PRECISION
49               The  value  of  the  parameter in the rank one update equation.
50               RHO >= 0 required.
51
52       DLAMDA  (input) DOUBLE PRECISION array, dimension (K)
53               The first K elements of this array contain the old roots of the
54               deflated  updating problem.  These are the poles of the secular
55               equation.
56
57       W       (input) DOUBLE PRECISION array, dimension (K)
58               The first K elements of this array contain  the  components  of
59               the deflation-adjusted updating vector.
60
61       S       (output) DOUBLE PRECISION array, dimension (LDS, K)
62               Will contain the eigenvectors of the repaired matrix which will
63               be stored for subsequent Z vector calculation and multiplied by
64               the previously accumulated eigenvectors to update the system.
65
66       LDS     (input) INTEGER
67               The leading dimension of S.  LDS >= max( 1, K ).
68
69       INFO    (output) INTEGER
70               = 0:  successful exit.
71               < 0:  if INFO = -i, the i-th argument had an illegal value.
72               > 0:  if INFO = 1, an eigenvalue did not converge
73

FURTHER DETAILS

75       Based on contributions by
76          Jeff Rutter, Computer Science Division, University of California
77          at Berkeley, USA
78
79
80
81 LAPACK routine (version 3.2)    November 2008                       DLAED9(1)