version 1.6, 2010/08/13 21:03:58
|
version 1.8, 2011/07/22 07:38:11
|
Line 1
|
Line 1
|
SUBROUTINE DSTERF( N, D, E, INFO ) |
SUBROUTINE DSTERF( N, D, E, INFO ) |
* |
* |
* -- LAPACK routine (version 3.2) -- |
* -- LAPACK routine (version 3.3.1) -- |
* -- LAPACK is a software package provided by Univ. of Tennessee, -- |
* -- LAPACK is a software package provided by Univ. of Tennessee, -- |
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
* November 2006 |
* -- April 2011 -- |
* |
* |
* .. Scalar Arguments .. |
* .. Scalar Arguments .. |
INTEGER INFO, N |
INTEGER INFO, N |
Line 54
|
Line 54
|
$ NMAXIT |
$ NMAXIT |
DOUBLE PRECISION ALPHA, ANORM, BB, C, EPS, EPS2, GAMMA, OLDC, |
DOUBLE PRECISION ALPHA, ANORM, BB, C, EPS, EPS2, GAMMA, OLDC, |
$ OLDGAM, P, R, RT1, RT2, RTE, S, SAFMAX, SAFMIN, |
$ OLDGAM, P, R, RT1, RT2, RTE, S, SAFMAX, SAFMIN, |
$ SIGMA, SSFMAX, SSFMIN |
$ SIGMA, SSFMAX, SSFMIN, RMAX |
* .. |
* .. |
* .. External Functions .. |
* .. External Functions .. |
DOUBLE PRECISION DLAMCH, DLANST, DLAPY2 |
DOUBLE PRECISION DLAMCH, DLANST, DLAPY2 |
Line 90
|
Line 90
|
SAFMAX = ONE / SAFMIN |
SAFMAX = ONE / SAFMIN |
SSFMAX = SQRT( SAFMAX ) / THREE |
SSFMAX = SQRT( SAFMAX ) / THREE |
SSFMIN = SQRT( SAFMIN ) / EPS2 |
SSFMIN = SQRT( SAFMIN ) / EPS2 |
|
RMAX = DLAMCH( 'O' ) |
* |
* |
* Compute the eigenvalues of the tridiagonal matrix. |
* Compute the eigenvalues of the tridiagonal matrix. |
* |
* |
Line 128
|
Line 129
|
* |
* |
* Scale submatrix in rows and columns L to LEND |
* Scale submatrix in rows and columns L to LEND |
* |
* |
ANORM = DLANST( 'I', LEND-L+1, D( L ), E( L ) ) |
ANORM = DLANST( 'M', LEND-L+1, D( L ), E( L ) ) |
ISCALE = 0 |
ISCALE = 0 |
IF( ANORM.GT.SSFMAX ) THEN |
IF( ANORM.EQ.ZERO ) |
|
$ GO TO 10 |
|
IF( (ANORM.GT.SSFMAX) ) THEN |
ISCALE = 1 |
ISCALE = 1 |
CALL DLASCL( 'G', 0, 0, ANORM, SSFMAX, LEND-L+1, 1, D( L ), N, |
CALL DLASCL( 'G', 0, 0, ANORM, SSFMAX, LEND-L+1, 1, D( L ), N, |
$ INFO ) |
$ INFO ) |