version 1.1, 2017/06/17 11:02:54
|
version 1.5, 2020/05/21 21:46:05
|
Line 217
|
Line 217
|
*> eigenvalues are computed to high relative accuracy when |
*> eigenvalues are computed to high relative accuracy when |
*> possible in future releases. The current code does not |
*> possible in future releases. The current code does not |
*> make any guarantees about high relative accuracy, but |
*> make any guarantees about high relative accuracy, but |
*> furutre releases will. See J. Barlow and J. Demmel, |
*> future releases will. See J. Barlow and J. Demmel, |
*> "Computing Accurate Eigensystems of Scaled Diagonally |
*> "Computing Accurate Eigensystems of Scaled Diagonally |
*> Dominant Matrices", LAPACK Working Note #7, for a discussion |
*> Dominant Matrices", LAPACK Working Note #7, for a discussion |
*> of which matrices define their eigenvalues to high relative |
*> of which matrices define their eigenvalues to high relative |
Line 408
|
Line 408
|
* |
* |
IMPLICIT NONE |
IMPLICIT NONE |
* |
* |
* -- LAPACK driver routine (version 3.7.0) -- |
* -- LAPACK driver routine (version 3.8.0) -- |
* -- 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..-- |
* June 2016 |
* June 2016 |
Line 445
|
Line 445
|
* .. |
* .. |
* .. External Functions .. |
* .. External Functions .. |
LOGICAL LSAME |
LOGICAL LSAME |
INTEGER ILAENV |
INTEGER ILAENV, ILAENV2STAGE |
DOUBLE PRECISION DLAMCH, ZLANSY |
DOUBLE PRECISION DLAMCH, ZLANSY |
EXTERNAL LSAME, ILAENV, DLAMCH, ZLANSY |
EXTERNAL LSAME, DLAMCH, ZLANSY, ILAENV, ILAENV2STAGE |
* .. |
* .. |
* .. External Subroutines .. |
* .. External Subroutines .. |
EXTERNAL DCOPY, DSCAL, DSTEBZ, DSTERF, XERBLA, ZDSCAL, |
EXTERNAL DCOPY, DSCAL, DSTEBZ, DSTERF, XERBLA, ZDSCAL, |
Line 471
|
Line 471
|
LQUERY = ( ( LWORK.EQ.-1 ) .OR. ( LRWORK.EQ.-1 ) .OR. |
LQUERY = ( ( LWORK.EQ.-1 ) .OR. ( LRWORK.EQ.-1 ) .OR. |
$ ( LIWORK.EQ.-1 ) ) |
$ ( LIWORK.EQ.-1 ) ) |
* |
* |
KD = ILAENV( 17, 'DSYTRD_2STAGE', JOBZ, N, -1, -1, -1 ) |
KD = ILAENV2STAGE( 1, 'ZHETRD_2STAGE', JOBZ, N, -1, -1, -1 ) |
IB = ILAENV( 18, 'DSYTRD_2STAGE', JOBZ, N, KD, -1, -1 ) |
IB = ILAENV2STAGE( 2, 'ZHETRD_2STAGE', JOBZ, N, KD, -1, -1 ) |
LHTRD = ILAENV( 19, 'DSYTRD_2STAGE', JOBZ, N, KD, IB, -1 ) |
LHTRD = ILAENV2STAGE( 3, 'ZHETRD_2STAGE', JOBZ, N, KD, IB, -1 ) |
LWTRD = ILAENV( 20, 'DSYTRD_2STAGE', JOBZ, N, KD, IB, -1 ) |
LWTRD = ILAENV2STAGE( 4, 'ZHETRD_2STAGE', JOBZ, N, KD, IB, -1 ) |
LWMIN = N + LHTRD + LWTRD |
LWMIN = N + LHTRD + LWTRD |
LRWMIN = MAX( 1, 24*N ) |
LRWMIN = MAX( 1, 24*N ) |
LIWMIN = MAX( 1, 10*N ) |
LIWMIN = MAX( 1, 10*N ) |