Subroutine */
| 11235 | } /* dlaed5_ */ |
| 11236 | |
| 11237 | /* Subroutine */ int dlaed6_(integer *kniter, logical *orgati, doublereal * |
| 11238 | rho, doublereal *d__, doublereal *z__, doublereal *finit, doublereal * |
| 11239 | tau, integer *info) |
| 11240 | { |
| 11241 | /* System generated locals */ |
| 11242 | integer i__1; |
| 11243 | doublereal d__1, d__2, d__3, d__4; |
| 11244 | |
| 11245 | /* Local variables */ |
| 11246 | static doublereal a, b, c__, f; |
| 11247 | static integer i__; |
| 11248 | static doublereal fc, df, ddf, lbd, eta, ubd, eps, base; |
| 11249 | static integer iter; |
| 11250 | static doublereal temp, temp1, temp2, temp3, temp4; |
| 11251 | static logical scale; |
| 11252 | static integer niter; |
| 11253 | static doublereal small1, small2, sminv1, sminv2; |
| 11254 | |
| 11255 | static doublereal dscale[3], sclfac, zscale[3], erretm, sclinv; |
| 11256 | |
| 11257 | |
| 11258 | /* |
| 11259 | -- LAPACK routine (version 3.2) -- |
| 11260 | -- LAPACK is a software package provided by Univ. of Tennessee, -- |
| 11261 | -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
| 11262 | February 2007 |
| 11263 | |
| 11264 | |
| 11265 | Purpose |
| 11266 | ======= |
| 11267 | |
| 11268 | DLAED6 computes the positive or negative root (closest to the origin) |
| 11269 | of |
| 11270 | z(1) z(2) z(3) |
| 11271 | f(x) = rho + --------- + ---------- + --------- |
| 11272 | d(1)-x d(2)-x d(3)-x |
| 11273 | |
| 11274 | It is assumed that |
| 11275 | |
| 11276 | if ORGATI = .true. the root is between d(2) and d(3); |
| 11277 | otherwise it is between d(1) and d(2) |
| 11278 | |
| 11279 | This routine will be called by DLAED4 when necessary. In most cases, |
| 11280 | the root sought is the smallest in magnitude, though it might not be |
| 11281 | in some extremely rare situations. |
| 11282 | |
| 11283 | Arguments |
| 11284 | ========= |
| 11285 | |
| 11286 | KNITER (input) INTEGER |
| 11287 | Refer to DLAED4 for its significance. |
| 11288 | |
| 11289 | ORGATI (input) LOGICAL |
| 11290 | If ORGATI is true, the needed root is between d(2) and |
| 11291 | d(3); otherwise it is between d(1) and d(2). See |
| 11292 | DLAED4 for further details. |
| 11293 | |
| 11294 | RHO (input) DOUBLE PRECISION |