diff options
author | Tim Hopkins, University of Kent <@> | 2016-06-17 17:10:39 -0400 |
---|---|---|
committer | Julien Langou <julien.langou@ucdenver.edu> | 2016-06-17 17:10:39 -0400 |
commit | d735f09838dcae48c8860266680b4b95d93f2c2b (patch) | |
tree | 78b4bc55a9535ab469031ad638b272b806640d79 /TESTING | |
parent | 3ec6c79d228388e8ccfadfa10e49e6c8cc2596e0 (diff) | |
download | lapack-d735f09838dcae48c8860266680b4b95d93f2c2b.tar.gz lapack-d735f09838dcae48c8860266680b4b95d93f2c2b.tar.bz2 lapack-d735f09838dcae48c8860266680b4b95d93f2c2b.zip |
Patch from Tim Hopkins sent by email on March 14th, 2016
(merge from Julien on June 17th, 2016)
Diffstat (limited to 'TESTING')
-rw-r--r-- | TESTING/EIG/cchkee.f | 4 | ||||
-rw-r--r-- | TESTING/EIG/cerrgg.f | 40 | ||||
-rw-r--r-- | TESTING/EIG/dchkee.f | 3 | ||||
-rw-r--r-- | TESTING/EIG/ddrvbd.f | 2 | ||||
-rw-r--r-- | TESTING/EIG/derrgg.f | 24 | ||||
-rw-r--r-- | TESTING/EIG/schkee.f | 5 | ||||
-rw-r--r-- | TESTING/EIG/serrgg.f | 26 | ||||
-rw-r--r-- | TESTING/EIG/zchkee.f | 3 | ||||
-rw-r--r-- | TESTING/EIG/zerrgg.f | 40 |
9 files changed, 80 insertions, 67 deletions
diff --git a/TESTING/EIG/cchkee.f b/TESTING/EIG/cchkee.f index 91214b0c..74329653 100644 --- a/TESTING/EIG/cchkee.f +++ b/TESTING/EIG/cchkee.f @@ -2098,6 +2098,7 @@ MAXTYP = 26 NTYPES = MIN( MAXTYP, NTYPES ) CALL ALAREQ( C3, NTYPES, DOTYPE, MAXTYP, NIN, NOUT ) + CALL XLAENV(1,1) IF( TSTCHK .AND. TSTERR ) $ CALL CERRGG( C3, NOUT ) DO 350 I = 1, NPARMS @@ -2157,6 +2158,7 @@ * * Blocked version * + CALL XLAENV(16,2) CALL CDRGES3( NN, NVAL, MAXTYP, DOTYPE, ISEED, THRESH, NOUT, $ A( 1, 1 ), NMAX, A( 1, 2 ), A( 1, 3 ), $ A( 1, 4 ), A( 1, 7 ), NMAX, A( 1, 8 ), @@ -2223,6 +2225,7 @@ * * Blocked version * + CALL XLAENV(16,2) CALL CDRGEV3( NN, NVAL, MAXTYP, DOTYPE, ISEED, THRESH, NOUT, $ A( 1, 1 ), NMAX, A( 1, 2 ), A( 1, 3 ), $ A( 1, 4 ), A( 1, 7 ), NMAX, A( 1, 8 ), @@ -2348,6 +2351,7 @@ * GSV: Generalized Singular Value Decomposition * ---------------------------------------------- * + CALL XLAENV(1,1) IF( TSTERR ) $ CALL CERRGG( 'GSV', NOUT ) CALL CCKGSV( NN, MVAL, PVAL, NVAL, NTYPES, ISEED, THRESH, NMAX, diff --git a/TESTING/EIG/cerrgg.f b/TESTING/EIG/cerrgg.f index a60c4eb2..0901af94 100644 --- a/TESTING/EIG/cerrgg.f +++ b/TESTING/EIG/cerrgg.f @@ -83,7 +83,7 @@ * .. * .. Local Arrays .. LOGICAL BW( NMAX ), SEL( NMAX ) - INTEGER IW( LW ) + INTEGER IW( LW ), IDUM(NMAX) REAL LS( NMAX ), R1( NMAX ), R2( NMAX ), $ RCE( NMAX ), RCV( NMAX ), RS( NMAX ), RW( LW ) COMPLEX A( NMAX, NMAX ), ALPHA( NMAX ), @@ -306,57 +306,57 @@ SRNAMT = 'CGGSVD3' INFOT = 1 CALL CGGSVD3( '/', 'N', 'N', 0, 0, 0, DUMMYK, DUMMYL, A, 1, B, - $ 1, R1, R2, U, 1, V, 1, Q, 1, W, RW, IW, LWORK, + $ 1, R1, R2, U, 1, V, 1, Q, 1, W, LWORK, RW, IDUM, $ INFO ) CALL CHKXER( 'CGGSVD3', INFOT, NOUT, LERR, OK ) INFOT = 2 CALL CGGSVD3( 'N', '/', 'N', 0, 0, 0, DUMMYK, DUMMYL, A, 1, B, - $ 1, R1, R2, U, 1, V, 1, Q, 1, W, RW, IW, LWORK, + $ 1, R1, R2, U, 1, V, 1, Q, 1, W, LWORK, RW, IDUM, $ INFO ) CALL CHKXER( 'CGGSVD3', INFOT, NOUT, LERR, OK ) INFOT = 3 CALL CGGSVD3( 'N', 'N', '/', 0, 0, 0, DUMMYK, DUMMYL, A, 1, B, - $ 1, R1, R2, U, 1, V, 1, Q, 1, W, RW, IW, LWORK, + $ 1, R1, R2, U, 1, V, 1, Q, 1, W, LWORK, RW, IDUM, $ INFO ) CALL CHKXER( 'CGGSVD3', INFOT, NOUT, LERR, OK ) INFOT = 4 CALL CGGSVD3( 'N', 'N', 'N', -1, 0, 0, DUMMYK, DUMMYL, A, 1, B, - $ 1, R1, R2, U, 1, V, 1, Q, 1, W, RW, IW, LWORK, + $ 1, R1, R2, U, 1, V, 1, Q, 1, W, LWORK, RW, IDUM, $ INFO ) CALL CHKXER( 'CGGSVD3', INFOT, NOUT, LERR, OK ) INFOT = 5 CALL CGGSVD3( 'N', 'N', 'N', 0, -1, 0, DUMMYK, DUMMYL, A, 1, B, - $ 1, R1, R2, U, 1, V, 1, Q, 1, W, RW, IW, LWORK, + $ 1, R1, R2, U, 1, V, 1, Q, 1, W, LWORK, RW, IDUM, $ INFO ) CALL CHKXER( 'CGGSVD3', INFOT, NOUT, LERR, OK ) INFOT = 6 CALL CGGSVD3( 'N', 'N', 'N', 0, 0, -1, DUMMYK, DUMMYL, A, 1, B, - $ 1, R1, R2, U, 1, V, 1, Q, 1, W, RW, IW, LWORK, + $ 1, R1, R2, U, 1, V, 1, Q, 1, W, LWORK, RW, IDUM, $ INFO ) CALL CHKXER( 'CGGSVD3', INFOT, NOUT, LERR, OK ) INFOT = 10 CALL CGGSVD3( 'N', 'N', 'N', 2, 1, 1, DUMMYK, DUMMYL, A, 1, B, - $ 1, R1, R2, U, 1, V, 1, Q, 1, W, RW, IW, LWORK, + $ 1, R1, R2, U, 1, V, 1, Q, 1, W, LWORK, RW, IDUM, $ INFO ) CALL CHKXER( 'CGGSVD3', INFOT, NOUT, LERR, OK ) INFOT = 12 CALL CGGSVD3( 'N', 'N', 'N', 1, 1, 2, DUMMYK, DUMMYL, A, 1, B, - $ 1, R1, R2, U, 1, V, 1, Q, 1, W, RW, IW, LWORK, + $ 1, R1, R2, U, 1, V, 1, Q, 1, W, LWORK, RW, IDUM, $ INFO ) CALL CHKXER( 'CGGSVD3', INFOT, NOUT, LERR, OK ) INFOT = 16 CALL CGGSVD3( 'U', 'N', 'N', 2, 2, 2, DUMMYK, DUMMYL, A, 2, B, - $ 2, R1, R2, U, 1, V, 1, Q, 1, W, RW, IW, LWORK, + $ 2, R1, R2, U, 1, V, 1, Q, 1, W, LWORK, RW, IDUM, $ INFO ) CALL CHKXER( 'CGGSVD3', INFOT, NOUT, LERR, OK ) INFOT = 18 CALL CGGSVD3( 'N', 'V', 'N', 2, 2, 2, DUMMYK, DUMMYL, A, 2, B, - $ 2, R1, R2, U, 2, V, 1, Q, 1, W, RW, IW, LWORK, + $ 2, R1, R2, U, 2, V, 1, Q, 1, W, LWORK, RW, IDUM, $ INFO ) CALL CHKXER( 'CGGSVD3', INFOT, NOUT, LERR, OK ) INFOT = 20 CALL CGGSVD3( 'N', 'N', 'Q', 2, 2, 2, DUMMYK, DUMMYL, A, 2, B, - $ 2, R1, R2, U, 2, V, 2, Q, 1, W, RW, IW, LWORK, + $ 2, R1, R2, U, 2, V, 2, Q, 1, W, LWORK, RW, IDUM, $ INFO ) CALL CHKXER( 'CGGSVD3', INFOT, NOUT, LERR, OK ) NT = NT + 11 @@ -573,56 +573,56 @@ INFOT = 7 CALL CUNCSD( 'Y', 'Y', 'Y', 'Y', 'N', 'N', $ -1, 0, 0, A, 1, A, - $ 1, A, 1, A, 1, A, + $ 1, A, 1, A, 1, RS, $ A, 1, A, 1, A, 1, A, $ 1, W, LW, RW, LW, IW, INFO ) CALL CHKXER( 'CUNCSD', INFOT, NOUT, LERR, OK ) INFOT = 8 CALL CUNCSD( 'Y', 'Y', 'Y', 'Y', 'N', 'N', $ 1, -1, 0, A, 1, A, - $ 1, A, 1, A, 1, A, + $ 1, A, 1, A, 1, RS, $ A, 1, A, 1, A, 1, A, $ 1, W, LW, RW, LW, IW, INFO ) CALL CHKXER( 'CUNCSD', INFOT, NOUT, LERR, OK ) INFOT = 9 CALL CUNCSD( 'Y', 'Y', 'Y', 'Y', 'N', 'N', $ 1, 1, -1, A, 1, A, - $ 1, A, 1, A, 1, A, + $ 1, A, 1, A, 1, RS, $ A, 1, A, 1, A, 1, A, $ 1, W, LW, RW, LW, IW, INFO ) CALL CHKXER( 'CUNCSD', INFOT, NOUT, LERR, OK ) INFOT = 11 CALL CUNCSD( 'Y', 'Y', 'Y', 'Y', 'N', 'N', $ 1, 1, 1, A, -1, A, - $ 1, A, 1, A, 1, A, + $ 1, A, 1, A, 1, RS, $ A, 1, A, 1, A, 1, A, $ 1, W, LW, RW, LW, IW, INFO ) CALL CHKXER( 'CUNCSD', INFOT, NOUT, LERR, OK ) INFOT = 20 CALL CUNCSD( 'Y', 'Y', 'Y', 'Y', 'N', 'N', $ 1, 1, 1, A, 1, A, - $ 1, A, 1, A, 1, A, + $ 1, A, 1, A, 1, RS, $ A, -1, A, 1, A, 1, A, $ 1, W, LW, RW, LW, IW, INFO ) CALL CHKXER( 'CUNCSD', INFOT, NOUT, LERR, OK ) INFOT = 22 CALL CUNCSD( 'Y', 'Y', 'Y', 'Y', 'N', 'N', $ 1, 1, 1, A, 1, A, - $ 1, A, 1, A, 1, A, + $ 1, A, 1, A, 1, RS, $ A, 1, A, -1, A, 1, A, $ 1, W, LW, RW, LW, IW, INFO ) CALL CHKXER( 'CUNCSD', INFOT, NOUT, LERR, OK ) INFOT = 24 CALL CUNCSD( 'Y', 'Y', 'Y', 'Y', 'N', 'N', $ 1, 1, 1, A, 1, A, - $ 1, A, 1, A, 1, A, + $ 1, A, 1, A, 1, RS, $ A, 1, A, 1, A, -1, A, $ 1, W, LW, RW, LW, IW, INFO ) CALL CHKXER( 'CUNCSD', INFOT, NOUT, LERR, OK ) INFOT = 26 CALL CUNCSD( 'Y', 'Y', 'Y', 'Y', 'N', 'N', $ 1, 1, 1, A, 1, A, - $ 1, A, 1, A, 1, A, + $ 1, A, 1, A, 1, RS, $ A, 1, A, 1, A, 1, A, $ -1, W, LW, RW, LW, IW, INFO ) CALL CHKXER( 'CUNCSD', INFOT, NOUT, LERR, OK ) diff --git a/TESTING/EIG/dchkee.f b/TESTING/EIG/dchkee.f index 14272bc4..c344bfea 100644 --- a/TESTING/EIG/dchkee.f +++ b/TESTING/EIG/dchkee.f @@ -2105,6 +2105,7 @@ MAXTYP = 26 NTYPES = MIN( MAXTYP, NTYPES ) CALL ALAREQ( C3, NTYPES, DOTYPE, MAXTYP, NIN, NOUT ) + CALL XLAENV(1,1) IF( TSTCHK .AND. TSTERR ) $ CALL DERRGG( C3, NOUT ) DO 350 I = 1, NPARMS @@ -2164,6 +2165,7 @@ * * Blocked version * + CALL XLAENV(16, 2) CALL DDRGES3( NN, NVAL, MAXTYP, DOTYPE, ISEED, THRESH, NOUT, $ A( 1, 1 ), NMAX, A( 1, 2 ), A( 1, 3 ), $ A( 1, 4 ), A( 1, 7 ), NMAX, A( 1, 8 ), @@ -2351,6 +2353,7 @@ * GSV: Generalized Singular Value Decomposition * ---------------------------------------------- * + CALL XLAENV(1,1) IF( TSTERR ) $ CALL DERRGG( 'GSV', NOUT ) CALL DCKGSV( NN, MVAL, PVAL, NVAL, NTYPES, ISEED, THRESH, NMAX, diff --git a/TESTING/EIG/ddrvbd.f b/TESTING/EIG/ddrvbd.f index 6d671f3b..dc552a98 100644 --- a/TESTING/EIG/ddrvbd.f +++ b/TESTING/EIG/ddrvbd.f @@ -1141,7 +1141,7 @@ $ / '21 = | S - Spartial | / ( min(M,N) ulp |S| )', $ / '22 = 0 if S contains min(M,N) nonnegative values in', $ ' decreasing order, else 1/ulp', - $ / '23 = | A - U diag(S) VT | / ( |A| max(M,N) ulp ),' + $ / '23 = | A - U diag(S) VT | / ( |A| max(M,N) ulp ),', $ ' DGESVDX(V,V,A) ', $ / '24 = | I - U**T U | / ( M ulp ) ', $ / '25 = | I - VT VT**T | / ( N ulp ) ', diff --git a/TESTING/EIG/derrgg.f b/TESTING/EIG/derrgg.f index 917cf048..10a5fe78 100644 --- a/TESTING/EIG/derrgg.f +++ b/TESTING/EIG/derrgg.f @@ -83,7 +83,7 @@ * .. * .. Local Arrays .. LOGICAL BW( NMAX ), SEL( NMAX ) - INTEGER IW( NMAX ) + INTEGER IW( NMAX ), IDUM(NMAX) DOUBLE PRECISION A( NMAX, NMAX ), B( NMAX, NMAX ), LS( NMAX ), $ Q( NMAX, NMAX ), R1( NMAX ), R2( NMAX ), $ R3( NMAX ), RCE( 2 ), RCV( 2 ), RS( NMAX ), @@ -305,47 +305,47 @@ SRNAMT = 'DGGSVD3' INFOT = 1 CALL DGGSVD3( '/', 'N', 'N', 0, 0, 0, DUMMYK, DUMMYL, A, 1, B, - $ 1, R1, R2, U, 1, V, 1, Q, 1, W, IW, LWORK, INFO ) + $ 1, R1, R2, U, 1, V, 1, Q, 1, W, LWORK, IDUM, INFO ) CALL CHKXER( 'DGGSVD3', INFOT, NOUT, LERR, OK ) INFOT = 2 CALL DGGSVD3( 'N', '/', 'N', 0, 0, 0, DUMMYK, DUMMYL, A, 1, B, - $ 1, R1, R2, U, 1, V, 1, Q, 1, W, IW, LWORK, INFO ) + $ 1, R1, R2, U, 1, V, 1, Q, 1, W, LWORK, IDUM, INFO ) CALL CHKXER( 'DGGSVD3', INFOT, NOUT, LERR, OK ) INFOT = 3 CALL DGGSVD3( 'N', 'N', '/', 0, 0, 0, DUMMYK, DUMMYL, A, 1, B, - $ 1, R1, R2, U, 1, V, 1, Q, 1, W, IW, LWORK, INFO ) + $ 1, R1, R2, U, 1, V, 1, Q, 1, W, LWORK, IDUM, INFO ) CALL CHKXER( 'DGGSVD3', INFOT, NOUT, LERR, OK ) INFOT = 4 CALL DGGSVD3( 'N', 'N', 'N', -1, 0, 0, DUMMYK, DUMMYL, A, 1, B, - $ 1, R1, R2, U, 1, V, 1, Q, 1, W, IW, LWORK, INFO ) + $ 1, R1, R2, U, 1, V, 1, Q, 1, W, LWORK, IDUM, INFO ) CALL CHKXER( 'DGGSVD3', INFOT, NOUT, LERR, OK ) INFOT = 5 CALL DGGSVD3( 'N', 'N', 'N', 0, -1, 0, DUMMYK, DUMMYL, A, 1, B, - $ 1, R1, R2, U, 1, V, 1, Q, 1, W, IW, LWORK, INFO ) + $ 1, R1, R2, U, 1, V, 1, Q, 1, W, LWORK, IDUM, INFO ) CALL CHKXER( 'DGGSVD3', INFOT, NOUT, LERR, OK ) INFOT = 6 CALL DGGSVD3( 'N', 'N', 'N', 0, 0, -1, DUMMYK, DUMMYL, A, 1, B, - $ 1, R1, R2, U, 1, V, 1, Q, 1, W, IW, LWORK, INFO ) + $ 1, R1, R2, U, 1, V, 1, Q, 1, W, LWORK, IDUM, INFO ) CALL CHKXER( 'DGGSVD3', INFOT, NOUT, LERR, OK ) INFOT = 10 CALL DGGSVD3( 'N', 'N', 'N', 2, 1, 1, DUMMYK, DUMMYL, A, 1, B, - $ 1, R1, R2, U, 1, V, 1, Q, 1, W, IW, LWORK, INFO ) + $ 1, R1, R2, U, 1, V, 1, Q, 1, W, LWORK, IDUM, INFO ) CALL CHKXER( 'DGGSVD3', INFOT, NOUT, LERR, OK ) INFOT = 12 CALL DGGSVD3( 'N', 'N', 'N', 1, 1, 2, DUMMYK, DUMMYL, A, 1, B, - $ 1, R1, R2, U, 1, V, 1, Q, 1, W, IW, LWORK, INFO ) + $ 1, R1, R2, U, 1, V, 1, Q, 1, W, LWORK, IDUM, INFO ) CALL CHKXER( 'DGGSVD3', INFOT, NOUT, LERR, OK ) INFOT = 16 CALL DGGSVD3( 'U', 'N', 'N', 2, 2, 2, DUMMYK, DUMMYL, A, 2, B, - $ 2, R1, R2, U, 1, V, 1, Q, 1, W, IW, LWORK, INFO ) + $ 2, R1, R2, U, 1, V, 1, Q, 1, W, LWORK, IDUM, INFO ) CALL CHKXER( 'DGGSVD3', INFOT, NOUT, LERR, OK ) INFOT = 18 CALL DGGSVD3( 'N', 'V', 'N', 1, 1, 2, DUMMYK, DUMMYL, A, 1, B, - $ 2, R1, R2, U, 1, V, 1, Q, 1, W, IW, LWORK, INFO ) + $ 2, R1, R2, U, 1, V, 1, Q, 1, W, LWORK, IDUM, INFO ) CALL CHKXER( 'DGGSVD3', INFOT, NOUT, LERR, OK ) INFOT = 20 CALL DGGSVD3( 'N', 'N', 'Q', 1, 2, 1, DUMMYK, DUMMYL, A, 1, B, - $ 1, R1, R2, U, 1, V, 1, Q, 1, W, IW, LWORK, INFO ) + $ 1, R1, R2, U, 1, V, 1, Q, 1, W, LWORK, IDUM, INFO ) CALL CHKXER( 'DGGSVD3', INFOT, NOUT, LERR, OK ) NT = NT + 11 * diff --git a/TESTING/EIG/schkee.f b/TESTING/EIG/schkee.f index 0bd994cc..b235c320 100644 --- a/TESTING/EIG/schkee.f +++ b/TESTING/EIG/schkee.f @@ -2105,8 +2105,9 @@ MAXTYP = 26 NTYPES = MIN( MAXTYP, NTYPES ) CALL ALAREQ( C3, NTYPES, DOTYPE, MAXTYP, NIN, NOUT ) + CALL XLAENV(1,1) IF( TSTCHK .AND. TSTERR ) - $ CALL SERRGG( C3, NOUT ) + & CALL SERRGG( C3, NOUT ) DO 350 I = 1, NPARMS CALL XLAENV( 1, NBVAL( I ) ) CALL XLAENV( 2, NBMIN( I ) ) @@ -2165,6 +2166,7 @@ * * Blocked version * + CALL XLAENV(16,1) CALL SDRGES3( NN, NVAL, MAXTYP, DOTYPE, ISEED, THRESH, NOUT, $ A( 1, 1 ), NMAX, A( 1, 2 ), A( 1, 3 ), $ A( 1, 4 ), A( 1, 7 ), NMAX, A( 1, 8 ), @@ -2353,6 +2355,7 @@ * GSV: Generalized Singular Value Decomposition * ---------------------------------------------- * + CALL XLAENV( 1, 1 ) IF( TSTERR ) $ CALL SERRGG( 'GSV', NOUT ) CALL SCKGSV( NN, MVAL, PVAL, NVAL, NTYPES, ISEED, THRESH, NMAX, diff --git a/TESTING/EIG/serrgg.f b/TESTING/EIG/serrgg.f index aede817b..5b9a8956 100644 --- a/TESTING/EIG/serrgg.f +++ b/TESTING/EIG/serrgg.f @@ -78,12 +78,12 @@ * .. Local Scalars .. CHARACTER*2 C2 INTEGER DUMMYK, DUMMYL, I, IFST, ILO, IHI, ILST, INFO, - $ J, M, NCYCLE, NT, SDIM, LWORK + $ J, M, NCYCLE, NT, SDIM, LWORK, JDUM REAL ANRM, BNRM, DIF, SCALE, TOLA, TOLB * .. * .. Local Arrays .. LOGICAL BW( NMAX ), SEL( NMAX ) - INTEGER IW( NMAX ) + INTEGER IW( NMAX ), IDUM(NMAX) REAL A( NMAX, NMAX ), B( NMAX, NMAX ), LS( NMAX ), $ Q( NMAX, NMAX ), R1( NMAX ), R2( NMAX ), $ R3( NMAX ), RCE( 2 ), RCV( 2 ), RS( NMAX ), @@ -305,47 +305,47 @@ SRNAMT = 'SGGSVD3' INFOT = 1 CALL SGGSVD3( '/', 'N', 'N', 0, 0, 0, DUMMYK, DUMMYL, A, 1, B, - $ 1, R1, R2, U, 1, V, 1, Q, 1, W, IW, LWORK, INFO ) + $ 1, R1, R2, U, 1, V, 1, Q, 1, W, LWORK, IDUM, INFO ) CALL CHKXER( 'SGGSVD3', INFOT, NOUT, LERR, OK ) INFOT = 2 CALL SGGSVD3( 'N', '/', 'N', 0, 0, 0, DUMMYK, DUMMYL, A, 1, B, - $ 1, R1, R2, U, 1, V, 1, Q, 1, W, IW, LWORK, INFO ) + $ 1, R1, R2, U, 1, V, 1, Q, 1, W, LWORK, IDUM, INFO ) CALL CHKXER( 'SGGSVD3', INFOT, NOUT, LERR, OK ) INFOT = 3 CALL SGGSVD3( 'N', 'N', '/', 0, 0, 0, DUMMYK, DUMMYL, A, 1, B, - $ 1, R1, R2, U, 1, V, 1, Q, 1, W, IW, LWORK, INFO ) + $ 1, R1, R2, U, 1, V, 1, Q, 1, W, LWORK, IDUM, INFO ) CALL CHKXER( 'SGGSVD3', INFOT, NOUT, LERR, OK ) INFOT = 4 CALL SGGSVD3( 'N', 'N', 'N', -1, 0, 0, DUMMYK, DUMMYL, A, 1, B, - $ 1, R1, R2, U, 1, V, 1, Q, 1, W, IW, LWORK, INFO ) + $ 1, R1, R2, U, 1, V, 1, Q, 1, W, LWORK, IDUM, INFO ) CALL CHKXER( 'SGGSVD3', INFOT, NOUT, LERR, OK ) INFOT = 5 CALL SGGSVD3( 'N', 'N', 'N', 0, -1, 0, DUMMYK, DUMMYL, A, 1, B, - $ 1, R1, R2, U, 1, V, 1, Q, 1, W, IW, LWORK, INFO ) + $ 1, R1, R2, U, 1, V, 1, Q, 1, W, LWORK, IDUM, INFO ) CALL CHKXER( 'SGGSVD3', INFOT, NOUT, LERR, OK ) INFOT = 6 CALL SGGSVD3( 'N', 'N', 'N', 0, 0, -1, DUMMYK, DUMMYL, A, 1, B, - $ 1, R1, R2, U, 1, V, 1, Q, 1, W, IW, LWORK, INFO ) + $ 1, R1, R2, U, 1, V, 1, Q, 1, W, LWORK, IDUM, INFO ) CALL CHKXER( 'SGGSVD3', INFOT, NOUT, LERR, OK ) INFOT = 10 CALL SGGSVD3( 'N', 'N', 'N', 2, 1, 1, DUMMYK, DUMMYL, A, 1, B, - $ 1, R1, R2, U, 1, V, 1, Q, 1, W, IW, LWORK, INFO ) + $ 1, R1, R2, U, 1, V, 1, Q, 1, W, LWORK, IDUM, INFO ) CALL CHKXER( 'SGGSVD3', INFOT, NOUT, LERR, OK ) INFOT = 12 CALL SGGSVD3( 'N', 'N', 'N', 1, 1, 2, DUMMYK, DUMMYL, A, 1, B, - $ 1, R1, R2, U, 1, V, 1, Q, 1, W, IW, LWORK, INFO ) + $ 1, R1, R2, U, 1, V, 1, Q, 1, W, LWORK, IDUM, INFO ) CALL CHKXER( 'SGGSVD3', INFOT, NOUT, LERR, OK ) INFOT = 16 CALL SGGSVD3( 'U', 'N', 'N', 2, 2, 2, DUMMYK, DUMMYL, A, 2, B, - $ 2, R1, R2, U, 1, V, 1, Q, 1, W, IW, LWORK, INFO ) + $ 2, R1, R2, U, 1, V, 1, Q, 1, W, LWORK, IDUM, INFO ) CALL CHKXER( 'SGGSVD3', INFOT, NOUT, LERR, OK ) INFOT = 18 CALL SGGSVD3( 'N', 'V', 'N', 1, 1, 2, DUMMYK, DUMMYL, A, 1, B, - $ 2, R1, R2, U, 1, V, 1, Q, 1, W, IW, LWORK, INFO ) + $ 2, R1, R2, U, 1, V, 1, Q, 1, W, LWORK, IDUM, INFO ) CALL CHKXER( 'SGGSVD3', INFOT, NOUT, LERR, OK ) INFOT = 20 CALL SGGSVD3( 'N', 'N', 'Q', 1, 2, 1, DUMMYK, DUMMYL, A, 1, B, - $ 1, R1, R2, U, 1, V, 1, Q, 1, W, IW, LWORK, INFO ) + $ 1, R1, R2, U, 1, V, 1, Q, 1, W, LWORK, IDUM, INFO ) CALL CHKXER( 'SGGSVD3', INFOT, NOUT, LERR, OK ) NT = NT + 11 * diff --git a/TESTING/EIG/zchkee.f b/TESTING/EIG/zchkee.f index 67221276..bcfecf80 100644 --- a/TESTING/EIG/zchkee.f +++ b/TESTING/EIG/zchkee.f @@ -2098,6 +2098,7 @@ MAXTYP = 26 NTYPES = MIN( MAXTYP, NTYPES ) CALL ALAREQ( C3, NTYPES, DOTYPE, MAXTYP, NIN, NOUT ) + CALL XLAENV(1,1) IF( TSTCHK .AND. TSTERR ) $ CALL ZERRGG( C3, NOUT ) DO 350 I = 1, NPARMS @@ -2222,6 +2223,7 @@ * * Blocked version * + CALL XLAENV(16,2) CALL ZDRGEV3( NN, NVAL, MAXTYP, DOTYPE, ISEED, THRESH, NOUT, $ A( 1, 1 ), NMAX, A( 1, 2 ), A( 1, 3 ), $ A( 1, 4 ), A( 1, 7 ), NMAX, A( 1, 8 ), @@ -2347,6 +2349,7 @@ * GSV: Generalized Singular Value Decomposition * ---------------------------------------------- * + CALL XLAENV(1,1) IF( TSTERR ) $ CALL ZERRGG( 'GSV', NOUT ) CALL ZCKGSV( NN, MVAL, PVAL, NVAL, NTYPES, ISEED, THRESH, NMAX, diff --git a/TESTING/EIG/zerrgg.f b/TESTING/EIG/zerrgg.f index 3366c771..e0589348 100644 --- a/TESTING/EIG/zerrgg.f +++ b/TESTING/EIG/zerrgg.f @@ -83,7 +83,7 @@ * .. * .. Local Arrays .. LOGICAL BW( NMAX ), SEL( NMAX ) - INTEGER IW( LW ) + INTEGER IW( LW ), IDUM(NMAX) DOUBLE PRECISION LS( NMAX ), R1( NMAX ), R2( NMAX ), $ RCE( NMAX ), RCV( NMAX ), RS( NMAX ), RW( LW ) COMPLEX*16 A( NMAX, NMAX ), ALPHA( NMAX ), @@ -306,57 +306,57 @@ SRNAMT = 'ZGGSVD3' INFOT = 1 CALL ZGGSVD3( '/', 'N', 'N', 0, 0, 0, DUMMYK, DUMMYL, A, 1, B, - $ 1, R1, R2, U, 1, V, 1, Q, 1, W, RW, IW, LWORK, + $ 1, R1, R2, U, 1, V, 1, Q, 1, W, LWORK, RW, IDUM, $ INFO ) CALL CHKXER( 'ZGGSVD3', INFOT, NOUT, LERR, OK ) INFOT = 2 CALL ZGGSVD3( 'N', '/', 'N', 0, 0, 0, DUMMYK, DUMMYL, A, 1, B, - $ 1, R1, R2, U, 1, V, 1, Q, 1, W, RW, IW, LWORK, + $ 1, R1, R2, U, 1, V, 1, Q, 1, W, LWORK, RW, IDUM, $ INFO ) CALL CHKXER( 'ZGGSVD3', INFOT, NOUT, LERR, OK ) INFOT = 3 CALL ZGGSVD3( 'N', 'N', '/', 0, 0, 0, DUMMYK, DUMMYL, A, 1, B, - $ 1, R1, R2, U, 1, V, 1, Q, 1, W, RW, IW, LWORK, + $ 1, R1, R2, U, 1, V, 1, Q, 1, W, LWORK, RW, IDUM, $ INFO ) CALL CHKXER( 'ZGGSVD3', INFOT, NOUT, LERR, OK ) INFOT = 4 CALL ZGGSVD3( 'N', 'N', 'N', -1, 0, 0, DUMMYK, DUMMYL, A, 1, B, - $ 1, R1, R2, U, 1, V, 1, Q, 1, W, RW, IW, LWORK, + $ 1, R1, R2, U, 1, V, 1, Q, 1, W, LWORK, RW, IDUM, $ INFO ) CALL CHKXER( 'ZGGSVD3', INFOT, NOUT, LERR, OK ) INFOT = 5 CALL ZGGSVD3( 'N', 'N', 'N', 0, -1, 0, DUMMYK, DUMMYL, A, 1, B, - $ 1, R1, R2, U, 1, V, 1, Q, 1, W, RW, IW, LWORK, + $ 1, R1, R2, U, 1, V, 1, Q, 1, W, LWORK, RW, IDUM, $ INFO ) CALL CHKXER( 'ZGGSVD3', INFOT, NOUT, LERR, OK ) INFOT = 6 CALL ZGGSVD3( 'N', 'N', 'N', 0, 0, -1, DUMMYK, DUMMYL, A, 1, B, - $ 1, R1, R2, U, 1, V, 1, Q, 1, W, RW, IW, LWORK, + $ 1, R1, R2, U, 1, V, 1, Q, 1, W, LWORK, RW, IDUM, $ INFO ) CALL CHKXER( 'ZGGSVD3', INFOT, NOUT, LERR, OK ) INFOT = 10 CALL ZGGSVD3( 'N', 'N', 'N', 2, 1, 1, DUMMYK, DUMMYL, A, 1, B, - $ 1, R1, R2, U, 1, V, 1, Q, 1, W, RW, IW, LWORK, + $ 1, R1, R2, U, 1, V, 1, Q, 1, W, LWORK, RW, IDUM, $ INFO ) CALL CHKXER( 'ZGGSVD3', INFOT, NOUT, LERR, OK ) INFOT = 12 CALL ZGGSVD3( 'N', 'N', 'N', 1, 1, 2, DUMMYK, DUMMYL, A, 1, B, - $ 1, R1, R2, U, 1, V, 1, Q, 1, W, RW, IW, LWORK, + $ 1, R1, R2, U, 1, V, 1, Q, 1, W, LWORK, RW, IDUM, $ INFO ) CALL CHKXER( 'ZGGSVD3', INFOT, NOUT, LERR, OK ) INFOT = 16 CALL ZGGSVD3( 'U', 'N', 'N', 2, 2, 2, DUMMYK, DUMMYL, A, 2, B, - $ 2, R1, R2, U, 1, V, 1, Q, 1, W, RW, IW, LWORK, + $ 2, R1, R2, U, 1, V, 1, Q, 1, W, LWORK, RW, IDUM, $ INFO ) CALL CHKXER( 'ZGGSVD3', INFOT, NOUT, LERR, OK ) INFOT = 18 CALL ZGGSVD3( 'N', 'V', 'N', 2, 2, 2, DUMMYK, DUMMYL, A, 2, B, - $ 2, R1, R2, U, 2, V, 1, Q, 1, W, RW, IW, LWORK, + $ 2, R1, R2, U, 2, V, 1, Q, 1, W, LWORK, RW, IDUM, $ INFO ) CALL CHKXER( 'ZGGSVD3', INFOT, NOUT, LERR, OK ) INFOT = 20 CALL ZGGSVD3( 'N', 'N', 'Q', 2, 2, 2, DUMMYK, DUMMYL, A, 2, B, - $ 2, R1, R2, U, 2, V, 2, Q, 1, W, RW, IW, LWORK, + $ 2, R1, R2, U, 2, V, 2, Q, 1, W, LWORK, RW, IDUM, $ INFO ) CALL CHKXER( 'ZGGSVD3', INFOT, NOUT, LERR, OK ) NT = NT + 11 @@ -573,56 +573,56 @@ INFOT = 7 CALL ZUNCSD( 'Y', 'Y', 'Y', 'Y', 'N', 'N', $ -1, 0, 0, A, 1, A, - $ 1, A, 1, A, 1, A, + $ 1, A, 1, A, 1, RS, $ A, 1, A, 1, A, 1, A, $ 1, W, LW, RW, LW, IW, INFO ) CALL CHKXER( 'ZUNCSD', INFOT, NOUT, LERR, OK ) INFOT = 8 CALL ZUNCSD( 'Y', 'Y', 'Y', 'Y', 'N', 'N', $ 1, -1, 0, A, 1, A, - $ 1, A, 1, A, 1, A, + $ 1, A, 1, A, 1, RS, $ A, 1, A, 1, A, 1, A, $ 1, W, LW, RW, LW, IW, INFO ) CALL CHKXER( 'ZUNCSD', INFOT, NOUT, LERR, OK ) INFOT = 9 CALL ZUNCSD( 'Y', 'Y', 'Y', 'Y', 'N', 'N', $ 1, 1, -1, A, 1, A, - $ 1, A, 1, A, 1, A, + $ 1, A, 1, A, 1, RS, $ A, 1, A, 1, A, 1, A, $ 1, W, LW, RW, LW, IW, INFO ) CALL CHKXER( 'ZUNCSD', INFOT, NOUT, LERR, OK ) INFOT = 11 CALL ZUNCSD( 'Y', 'Y', 'Y', 'Y', 'N', 'N', $ 1, 1, 1, A, -1, A, - $ 1, A, 1, A, 1, A, + $ 1, A, 1, A, 1, RS, $ A, 1, A, 1, A, 1, A, $ 1, W, LW, RW, LW, IW, INFO ) CALL CHKXER( 'ZUNCSD', INFOT, NOUT, LERR, OK ) INFOT = 20 CALL ZUNCSD( 'Y', 'Y', 'Y', 'Y', 'N', 'N', $ 1, 1, 1, A, 1, A, - $ 1, A, 1, A, 1, A, + $ 1, A, 1, A, 1, RS, $ A, -1, A, 1, A, 1, A, $ 1, W, LW, RW, LW, IW, INFO ) CALL CHKXER( 'ZUNCSD', INFOT, NOUT, LERR, OK ) INFOT = 22 CALL ZUNCSD( 'Y', 'Y', 'Y', 'Y', 'N', 'N', $ 1, 1, 1, A, 1, A, - $ 1, A, 1, A, 1, A, + $ 1, A, 1, A, 1, RS, $ A, 1, A, -1, A, 1, A, $ 1, W, LW, RW, LW, IW, INFO ) CALL CHKXER( 'ZUNCSD', INFOT, NOUT, LERR, OK ) INFOT = 24 CALL ZUNCSD( 'Y', 'Y', 'Y', 'Y', 'N', 'N', $ 1, 1, 1, A, 1, A, - $ 1, A, 1, A, 1, A, + $ 1, A, 1, A, 1, RS, $ A, 1, A, 1, A, -1, A, $ 1, W, LW, RW, LW, IW, INFO ) CALL CHKXER( 'ZUNCSD', INFOT, NOUT, LERR, OK ) INFOT = 26 CALL ZUNCSD( 'Y', 'Y', 'Y', 'Y', 'N', 'N', $ 1, 1, 1, A, 1, A, - $ 1, A, 1, A, 1, A, + $ 1, A, 1, A, 1, RS, $ A, 1, A, 1, A, 1, A, $ -1, W, LW, RW, LW, IW, INFO ) CALL CHKXER( 'ZUNCSD', INFOT, NOUT, LERR, OK ) |