version 1.4, 2010/08/06 15:32:29
|
version 1.5, 2010/08/07 13:18:07
|
Line 4
|
Line 4
|
$ IBLOCK, INDEXW, GERS, Z, LDZ, ISUPPZ, |
$ IBLOCK, INDEXW, GERS, Z, LDZ, ISUPPZ, |
$ WORK, IWORK, INFO ) |
$ WORK, IWORK, INFO ) |
* |
* |
* -- LAPACK auxiliary routine (version 3.2) -- |
* -- LAPACK auxiliary routine (version 3.2.2) -- |
* -- 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 |
* June 2010 |
* |
* |
* .. Scalar Arguments .. |
* .. Scalar Arguments .. |
INTEGER DOL, DOU, INFO, LDZ, M, N |
INTEGER DOL, DOU, INFO, LDZ, M, N |
Line 51
|
Line 51
|
* is stored the corresponding shift as given by DLARRE. |
* is stored the corresponding shift as given by DLARRE. |
* On exit, L is overwritten. |
* On exit, L is overwritten. |
* |
* |
* PIVMIN (in) DOUBLE PRECISION |
* PIVMIN (input) DOUBLE PRECISION |
* The minimum pivot allowed in the Sturm sequence. |
* The minimum pivot allowed in the Sturm sequence. |
* |
* |
* ISPLIT (input) INTEGER array, dimension (N) |
* ISPLIT (input) INTEGER array, dimension (N) |
Line 429
|
Line 429
|
* within the current block |
* within the current block |
P = INDEXW( WBEGIN-1+OLDFST ) |
P = INDEXW( WBEGIN-1+OLDFST ) |
Q = INDEXW( WBEGIN-1+OLDLST ) |
Q = INDEXW( WBEGIN-1+OLDLST ) |
* Offset for the arrays WORK, WGAP and WERR, i.e., th P-OFFSET |
* Offset for the arrays WORK, WGAP and WERR, i.e., the P-OFFSET |
* thru' Q-OFFSET elements of these arrays are to be used. |
* through the Q-OFFSET elements of these arrays are to be used. |
C OFFSET = P-OLDFST |
C OFFSET = P-OLDFST |
OFFSET = INDEXW( WBEGIN ) - 1 |
OFFSET = INDEXW( WBEGIN ) - 1 |
* perform limited bisection (if necessary) to get approximate |
* perform limited bisection (if necessary) to get approximate |