version 1.18, 2018/05/29 07:17:51
|
version 1.20, 2023/08/07 08:38:48
|
Line 272
|
Line 272
|
*> \author Univ. of Colorado Denver |
*> \author Univ. of Colorado Denver |
*> \author NAG Ltd. |
*> \author NAG Ltd. |
* |
* |
*> \date June 2016 |
|
* |
|
*> \ingroup doubleGEeigen |
*> \ingroup doubleGEeigen |
* |
* |
* ===================================================================== |
* ===================================================================== |
Line 281
|
Line 279
|
$ WR, WI, VS, LDVS, RCONDE, RCONDV, WORK, LWORK, |
$ WR, WI, VS, LDVS, RCONDE, RCONDV, WORK, LWORK, |
$ IWORK, LIWORK, BWORK, INFO ) |
$ IWORK, LIWORK, BWORK, INFO ) |
* |
* |
* -- 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 JOBVS, SENSE, SORT |
CHARACTER JOBVS, SENSE, SORT |
Line 385
|
Line 382
|
* |
* |
CALL DHSEQR( 'S', JOBVS, N, 1, N, A, LDA, WR, WI, VS, LDVS, |
CALL DHSEQR( 'S', JOBVS, N, 1, N, A, LDA, WR, WI, VS, LDVS, |
$ WORK, -1, IEVAL ) |
$ WORK, -1, IEVAL ) |
HSWORK = WORK( 1 ) |
HSWORK = INT( WORK( 1 ) ) |
* |
* |
IF( .NOT.WANTVS ) THEN |
IF( .NOT.WANTVS ) THEN |
MAXWRK = MAX( MAXWRK, N + HSWORK ) |
MAXWRK = MAX( MAXWRK, N + HSWORK ) |
Line 583
|
Line 580
|
IF( N.GT.I+1 ) |
IF( N.GT.I+1 ) |
$ CALL DSWAP( N-I-1, A( I, I+2 ), LDA, |
$ CALL DSWAP( N-I-1, A( I, I+2 ), LDA, |
$ A( I+1, I+2 ), LDA ) |
$ A( I+1, I+2 ), LDA ) |
CALL DSWAP( N, VS( 1, I ), 1, VS( 1, I+1 ), 1 ) |
IF( WANTVS ) THEN |
|
CALL DSWAP( N, VS( 1, I ), 1, VS( 1, I+1 ), 1 ) |
|
END IF |
A( I, I+1 ) = A( I+1, I ) |
A( I, I+1 ) = A( I+1, I ) |
A( I+1, I ) = ZERO |
A( I+1, I ) = ZERO |
END IF |
END IF |