1DLAED6(1)               LAPACK routine (version 3.1.1)               DLAED6(1)
2
3
4

NAME

6       DLAED6 - the positive or negative root (closest to the origin) of  z(1)
7       z(2) z(3) f(x) = rho +  ---------  +  ----------  +  ---------   d(1)-x
8       d(2)-x d(3)-x  It is assumed that   if ORGATI = .true
9

SYNOPSIS

11       SUBROUTINE DLAED6( KNITER, ORGATI, RHO, D, Z, FINIT, TAU, INFO )
12
13           LOGICAL        ORGATI
14
15           INTEGER        INFO, KNITER
16
17           DOUBLE         PRECISION FINIT, RHO, TAU
18
19           DOUBLE         PRECISION D( 3 ), Z( 3 )
20

PURPOSE

22       DLAED6  computes  the positive or negative root (closest to the origin)
23       of
24                        z(1)        z(2)        z(3) f(x) =   rho +  ---------
25       + ---------- + ---------
26                       d(1)-x      d(2)-x      d(3)-x
27             otherwise it is between d(1) and d(2)
28
29       This  routine  will  be called by DLAED4 when necessary. In most cases,
30       the root sought is the smallest in magnitude, though it might not be in
31       some extremely rare situations.
32
33

ARGUMENTS

35       KNITER       (input) INTEGER
36                    Refer to DLAED4 for its significance.
37
38       ORGATI       (input) LOGICAL
39                    If  ORGATI  is  true,  the needed root is between d(2) and
40                    d(3); otherwise it is between d(1) and d(2).   See  DLAED4
41                    for further details.
42
43       RHO          (input) DOUBLE PRECISION
44                    Refer to the equation f(x) above.
45
46       D            (input) DOUBLE PRECISION array, dimension (3)
47                    D satisfies d(1) < d(2) < d(3).
48
49       Z            (input) DOUBLE PRECISION array, dimension (3)
50                    Each of the elements in z must be positive.
51
52       FINIT        (input) DOUBLE PRECISION
53                    The  value  of  f  at  0. It is more accurate than the one
54                    evaluated inside this routine (if someone wants to do so).
55
56       TAU          (output) DOUBLE PRECISION
57                    The root of the equation f(x).
58
59       INFO         (output) INTEGER
60                    = 0: successful exit
61                    > 0: if INFO = 1, failure to converge
62

FURTHER DETAILS

64       30/06/99: Based on contributions by
65          Ren-Cang Li, Computer Science Division, University of California
66          at Berkeley, USA
67
68       10/02/03: This version has a few statements commented  out  for  thread
69       safety (machine parameters are computed on each entry). SJH.
70
71       05/10/06: Modified from a new version of Ren-Cang Li, use
72          Gragg-Thornton-Warner cubic convergent scheme for better stability.
73
74
75
76
77 LAPACK routine (version 3.1.1)  February 2007                       DLAED6(1)
Impressum