Subroutine */
| 38358 | } /* ssteqr_ */ |
| 38359 | |
| 38360 | /* Subroutine */ int ssterf_(integer *n, real *d__, real *e, integer *info) |
| 38361 | { |
| 38362 | /* System generated locals */ |
| 38363 | integer i__1; |
| 38364 | real r__1, r__2, r__3; |
| 38365 | |
| 38366 | /* Local variables */ |
| 38367 | static real c__; |
| 38368 | static integer i__, l, m; |
| 38369 | static real p, r__, s; |
| 38370 | static integer l1; |
| 38371 | static real bb, rt1, rt2, eps, rte; |
| 38372 | static integer lsv; |
| 38373 | static real eps2, oldc; |
| 38374 | static integer lend, jtot; |
| 38375 | extern /* Subroutine */ int slae2_(real *, real *, real *, real *, real *) |
| 38376 | ; |
| 38377 | static real gamma, alpha, sigma, anorm; |
| 38378 | extern doublereal slapy2_(real *, real *); |
| 38379 | static integer iscale; |
| 38380 | static real oldgam; |
| 38381 | extern doublereal slamch_(char *); |
| 38382 | static real safmin; |
| 38383 | extern /* Subroutine */ int xerbla_(char *, integer *); |
| 38384 | static real safmax; |
| 38385 | extern /* Subroutine */ int slascl_(char *, integer *, integer *, real *, |
| 38386 | real *, integer *, integer *, real *, integer *, integer *); |
| 38387 | static integer lendsv; |
| 38388 | static real ssfmin; |
| 38389 | static integer nmaxit; |
| 38390 | static real ssfmax; |
| 38391 | extern doublereal slanst_(char *, integer *, real *, real *); |
| 38392 | extern /* Subroutine */ int slasrt_(char *, integer *, real *, integer *); |
| 38393 | |
| 38394 | |
| 38395 | /* |
| 38396 | -- LAPACK routine (version 3.2) -- |
| 38397 | -- LAPACK is a software package provided by Univ. of Tennessee, -- |
| 38398 | -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
| 38399 | November 2006 |
| 38400 | |
| 38401 | |
| 38402 | Purpose |
| 38403 | ======= |
| 38404 | |
| 38405 | SSTERF computes all eigenvalues of a symmetric tridiagonal matrix |
| 38406 | using the Pal-Walker-Kahan variant of the QL or QR algorithm. |
| 38407 | |
| 38408 | Arguments |
| 38409 | ========= |
| 38410 | |
| 38411 | N (input) INTEGER |
| 38412 | The order of the matrix. N >= 0. |
| 38413 | |
| 38414 | D (input/output) REAL array, dimension (N) |
| 38415 | On entry, the n diagonal elements of the tridiagonal matrix. |
| 38416 | On exit, if INFO = 0, the eigenvalues in ascending order. |
| 38417 |