version 1.2, 2017/06/17 11:06:47
|
version 1.6, 2023/08/07 08:39:23
|
Line 81
|
Line 81
|
*> The desired accuracy of the output can be specified by the input |
*> The desired accuracy of the output can be specified by the input |
*> parameter ABSTOL. |
*> parameter ABSTOL. |
*> |
*> |
*> For more details, see DSTEMR's documentation and: |
*> For more details, see ZSTEMR's documentation and: |
*> - Inderjit S. Dhillon and Beresford N. Parlett: "Multiple representations |
*> - Inderjit S. Dhillon and Beresford N. Parlett: "Multiple representations |
*> to compute orthogonal eigenvectors of symmetric tridiagonal matrices," |
*> to compute orthogonal eigenvectors of symmetric tridiagonal matrices," |
*> Linear Algebra and its Applications, 387(1), pp. 1-28, August 2004. |
*> Linear Algebra and its Applications, 387(1), pp. 1-28, August 2004. |
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 354
|
Line 354
|
*> \author Univ. of Colorado Denver |
*> \author Univ. of Colorado Denver |
*> \author NAG Ltd. |
*> \author NAG Ltd. |
* |
* |
*> \date June 2016 |
|
* |
|
*> \ingroup complex16HEeigen |
*> \ingroup complex16HEeigen |
* |
* |
*> \par Contributors: |
*> \par Contributors: |
Line 408
|
Line 406
|
* |
* |
IMPLICIT NONE |
IMPLICIT NONE |
* |
* |
* -- LAPACK driver routine (version 3.7.0) -- |
* -- LAPACK driver routine -- |
* -- 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 |
|
* |
* |
* .. Scalar Arguments .. |
* .. Scalar Arguments .. |
CHARACTER JOBZ, RANGE, UPLO |
CHARACTER JOBZ, RANGE, UPLO |
Line 445
|
Line 442
|
* .. |
* .. |
* .. 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 468
|
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 ) |