version 1.3, 2010/08/06 15:28:48
|
version 1.9, 2011/07/22 07:38:11
|
Line 2
|
Line 2
|
$ VL, VU, IL, IU, ABSTOL, M, W, Z, LDZ, WORK, |
$ VL, VU, IL, IU, ABSTOL, M, W, Z, LDZ, WORK, |
$ LWORK, IWORK, IFAIL, INFO ) |
$ LWORK, IWORK, IFAIL, INFO ) |
* |
* |
* -- LAPACK driver routine (version 3.2) -- |
* -- LAPACK driver 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 .. |
CHARACTER JOBZ, RANGE, UPLO |
CHARACTER JOBZ, RANGE, UPLO |
Line 68
|
Line 68
|
* LDA (input) INTEGER |
* LDA (input) INTEGER |
* The leading dimension of the array A. LDA >= max(1,N). |
* The leading dimension of the array A. LDA >= max(1,N). |
* |
* |
* B (input/output) DOUBLE PRECISION array, dimension (LDA, N) |
* B (input/output) DOUBLE PRECISION array, dimension (LDB, N) |
* On entry, the symmetric matrix B. If UPLO = 'U', the |
* On entry, the symmetric matrix B. If UPLO = 'U', the |
* leading N-by-N upper triangular part of B contains the |
* leading N-by-N upper triangular part of B contains the |
* upper triangular part of the matrix B. If UPLO = 'L', |
* upper triangular part of the matrix B. If UPLO = 'L', |
Line 296
|
Line 296
|
IF( ITYPE.EQ.1 .OR. ITYPE.EQ.2 ) THEN |
IF( ITYPE.EQ.1 .OR. ITYPE.EQ.2 ) THEN |
* |
* |
* For A*x=(lambda)*B*x and A*B*x=(lambda)*x; |
* For A*x=(lambda)*B*x and A*B*x=(lambda)*x; |
* backtransform eigenvectors: x = inv(L)'*y or inv(U)*y |
* backtransform eigenvectors: x = inv(L)**T*y or inv(U)*y |
* |
* |
IF( UPPER ) THEN |
IF( UPPER ) THEN |
TRANS = 'N' |
TRANS = 'N' |
Line 310
|
Line 310
|
ELSE IF( ITYPE.EQ.3 ) THEN |
ELSE IF( ITYPE.EQ.3 ) THEN |
* |
* |
* For B*A*x=(lambda)*x; |
* For B*A*x=(lambda)*x; |
* backtransform eigenvectors: x = L*y or U'*y |
* backtransform eigenvectors: x = L*y or U**T*y |
* |
* |
IF( UPPER ) THEN |
IF( UPPER ) THEN |
TRANS = 'T' |
TRANS = 'T' |