version 1.8, 2010/12/21 13:53:46
|
version 1.9, 2011/07/22 07:38:15
|
Line 2
|
Line 2
|
$ ABSTOL, M, W, Z, LDZ, ISUPPZ, WORK, LWORK, |
$ ABSTOL, M, W, Z, LDZ, ISUPPZ, WORK, LWORK, |
$ RWORK, LRWORK, IWORK, LIWORK, INFO ) |
$ RWORK, LRWORK, IWORK, LIWORK, INFO ) |
* |
* |
* -- LAPACK driver routine (version 3.2.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..-- |
* June 2010 |
* -- April 2011 -- |
* |
* |
* .. Scalar Arguments .. |
* .. Scalar Arguments .. |
CHARACTER JOBZ, RANGE, UPLO |
CHARACTER JOBZ, RANGE, UPLO |
Line 92
|
Line 92
|
* = 'V': all eigenvalues in the half-open interval (VL,VU] |
* = 'V': all eigenvalues in the half-open interval (VL,VU] |
* will be found. |
* will be found. |
* = 'I': the IL-th through IU-th eigenvalues will be found. |
* = 'I': the IL-th through IU-th eigenvalues will be found. |
********** For RANGE = 'V' or 'I' and IU - IL < N - 1, DSTEBZ and |
* For RANGE = 'V' or 'I' and IU - IL < N - 1, DSTEBZ and |
********** ZSTEIN are called |
* ZSTEIN are called |
* |
* |
* UPLO (input) CHARACTER*1 |
* UPLO (input) CHARACTER*1 |
* = 'U': Upper triangle of A is stored; |
* = 'U': Upper triangle of A is stored; |
Line 183
|
Line 183
|
* indicating the nonzero elements in Z. The i-th eigenvector |
* indicating the nonzero elements in Z. The i-th eigenvector |
* is nonzero only in elements ISUPPZ( 2*i-1 ) through |
* is nonzero only in elements ISUPPZ( 2*i-1 ) through |
* ISUPPZ( 2*i ). |
* ISUPPZ( 2*i ). |
********** Implemented only for RANGE = 'A' or 'I' and IU - IL = N - 1 |
* Implemented only for RANGE = 'A' or 'I' and IU - IL = N - 1 |
* |
* |
* WORK (workspace/output) COMPLEX*16 array, dimension (MAX(1,LWORK)) |
* WORK (workspace/output) COMPLEX*16 array, dimension (MAX(1,LWORK)) |
* On exit, if INFO = 0, WORK(1) returns the optimal LWORK. |
* On exit, if INFO = 0, WORK(1) returns the optimal LWORK. |
Line 204
|
Line 204
|
* On exit, if INFO = 0, RWORK(1) returns the optimal |
* On exit, if INFO = 0, RWORK(1) returns the optimal |
* (and minimal) LRWORK. |
* (and minimal) LRWORK. |
* |
* |
* LRWORK (input) INTEGER |
* LRWORK (input) INTEGER |
* The length of the array RWORK. LRWORK >= max(1,24*N). |
* The length of the array RWORK. LRWORK >= max(1,24*N). |
* |
* |
* If LRWORK = -1, then a workspace query is assumed; the |
* If LRWORK = -1, then a workspace query is assumed; the |
Line 217
|
Line 217
|
* On exit, if INFO = 0, IWORK(1) returns the optimal |
* On exit, if INFO = 0, IWORK(1) returns the optimal |
* (and minimal) LIWORK. |
* (and minimal) LIWORK. |
* |
* |
* LIWORK (input) INTEGER |
* LIWORK (input) INTEGER |
* The dimension of the array IWORK. LIWORK >= max(1,10*N). |
* The dimension of the array IWORK. LIWORK >= max(1,10*N). |
* |
* |
* If LIWORK = -1, then a workspace query is assumed; the |
* If LIWORK = -1, then a workspace query is assumed; the |
Line 242
|
Line 242
|
* Jason Riedy, Computer Science Division, University of |
* Jason Riedy, Computer Science Division, University of |
* California at Berkeley, USA |
* California at Berkeley, USA |
* |
* |
* ===================================================================== |
* ===================================================================== |
* |
* |
* .. Parameters .. |
* .. Parameters .. |
DOUBLE PRECISION ZERO, ONE, TWO |
DOUBLE PRECISION ZERO, ONE, TWO |