--- rpl/lapack/lapack/zstemr.f 2012/12/14 12:30:34 1.11 +++ rpl/lapack/lapack/zstemr.f 2015/11/26 11:44:26 1.15 @@ -311,7 +311,7 @@ *> \author Univ. of Colorado Denver *> \author NAG Ltd. * -*> \date September 2012 +*> \date November 2015 * *> \ingroup complex16OTHERcomputational * @@ -329,10 +329,10 @@ $ M, W, Z, LDZ, NZC, ISUPPZ, TRYRAC, WORK, LWORK, $ IWORK, LIWORK, INFO ) * -* -- LAPACK computational routine (version 3.4.2) -- +* -- LAPACK computational routine (version 3.6.0) -- * -- LAPACK is a software package provided by Univ. of Tennessee, -- * -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- -* September 2012 +* November 2015 * * .. Scalar Arguments .. CHARACTER JOBZ, RANGE @@ -408,7 +408,8 @@ WU = ZERO IIL = 0 IIU = 0 - + NSPLIT = 0 + IF( VALEIG ) THEN * We do not reference VL, VU in the cases RANGE = 'I','A' * The interval (WL, WU] contains all the wanted eigenvalues. @@ -525,10 +526,10 @@ IF (SN.NE.ZERO) THEN IF (CS.NE.ZERO) THEN ISUPPZ(2*M-1) = 1 - ISUPPZ(2*M-1) = 2 + ISUPPZ(2*M) = 2 ELSE ISUPPZ(2*M-1) = 1 - ISUPPZ(2*M-1) = 1 + ISUPPZ(2*M) = 1 END IF ELSE ISUPPZ(2*M-1) = 2 @@ -549,10 +550,10 @@ IF (SN.NE.ZERO) THEN IF (CS.NE.ZERO) THEN ISUPPZ(2*M-1) = 1 - ISUPPZ(2*M-1) = 2 + ISUPPZ(2*M) = 2 ELSE ISUPPZ(2*M-1) = 1 - ISUPPZ(2*M-1) = 1 + ISUPPZ(2*M) = 1 END IF ELSE ISUPPZ(2*M-1) = 2