version 1.2, 2017/06/17 11:06:34
|
version 1.5, 2020/05/21 21:46:01
|
Line 143
|
Line 143
|
*> \author Univ. of Colorado Denver |
*> \author Univ. of Colorado Denver |
*> \author NAG Ltd. |
*> \author NAG Ltd. |
* |
* |
*> \date December 2016 |
*> \date November 2017 |
* |
* |
*> \ingroup doubleSYeigen |
*> \ingroup doubleSYeigen |
* |
* |
Line 185
|
Line 185
|
* |
* |
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..-- |
* December 2016 |
* November 2017 |
* |
* |
* .. Scalar Arguments .. |
* .. Scalar Arguments .. |
CHARACTER JOBZ, UPLO |
CHARACTER JOBZ, UPLO |
Line 213
|
Line 213
|
* .. |
* .. |
* .. External Functions .. |
* .. External Functions .. |
LOGICAL LSAME |
LOGICAL LSAME |
INTEGER ILAENV |
INTEGER ILAENV2STAGE |
DOUBLE PRECISION DLAMCH, DLANSY |
DOUBLE PRECISION DLAMCH, DLANSY |
EXTERNAL LSAME, ILAENV, DLAMCH, DLANSY |
EXTERNAL LSAME, DLAMCH, DLANSY, ILAENV2STAGE |
* .. |
* .. |
* .. External Subroutines .. |
* .. External Subroutines .. |
EXTERNAL DLASCL, DORGTR, DSCAL, DSTEQR, DSTERF, |
EXTERNAL DLASCL, DORGTR, DSCAL, DSTEQR, DSTERF, |
Line 244
|
Line 244
|
END IF |
END IF |
* |
* |
IF( INFO.EQ.0 ) THEN |
IF( INFO.EQ.0 ) THEN |
KD = ILAENV( 17, 'DSYTRD_2STAGE', JOBZ, N, -1, -1, -1 ) |
KD = ILAENV2STAGE( 1, 'DSYTRD_2STAGE', JOBZ, N, -1, -1, -1 ) |
IB = ILAENV( 18, 'DSYTRD_2STAGE', JOBZ, N, KD, -1, -1 ) |
IB = ILAENV2STAGE( 2, 'DSYTRD_2STAGE', JOBZ, N, KD, -1, -1 ) |
LHTRD = ILAENV( 19, 'DSYTRD_2STAGE', JOBZ, N, KD, IB, -1 ) |
LHTRD = ILAENV2STAGE( 3, 'DSYTRD_2STAGE', JOBZ, N, KD, IB, -1 ) |
LWTRD = ILAENV( 20, 'DSYTRD_2STAGE', JOBZ, N, KD, IB, -1 ) |
LWTRD = ILAENV2STAGE( 4, 'DSYTRD_2STAGE', JOBZ, N, KD, IB, -1 ) |
LWMIN = 2*N + LHTRD + LWTRD |
LWMIN = 2*N + LHTRD + LWTRD |
WORK( 1 ) = LWMIN |
WORK( 1 ) = LWMIN |
* |
* |
Line 317
|
Line 317
|
IF( .NOT.WANTZ ) THEN |
IF( .NOT.WANTZ ) THEN |
CALL DSTERF( N, W, WORK( INDE ), INFO ) |
CALL DSTERF( N, W, WORK( INDE ), INFO ) |
ELSE |
ELSE |
* Not available in this release, and agrument checking should not |
* Not available in this release, and argument checking should not |
* let it getting here |
* let it getting here |
RETURN |
RETURN |
CALL DORGTR( UPLO, N, A, LDA, WORK( INDTAU ), WORK( INDWRK ), |
CALL DORGTR( UPLO, N, A, LDA, WORK( INDTAU ), WORK( INDWRK ), |