summaryrefslogtreecommitdiff
path: root/TESTING
diff options
context:
space:
mode:
authorTim Hopkins, University of Kent <@>2016-06-17 17:10:39 -0400
committerJulien Langou <julien.langou@ucdenver.edu>2016-06-17 17:10:39 -0400
commitd735f09838dcae48c8860266680b4b95d93f2c2b (patch)
tree78b4bc55a9535ab469031ad638b272b806640d79 /TESTING
parent3ec6c79d228388e8ccfadfa10e49e6c8cc2596e0 (diff)
downloadlapack-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.f4
-rw-r--r--TESTING/EIG/cerrgg.f40
-rw-r--r--TESTING/EIG/dchkee.f3
-rw-r--r--TESTING/EIG/ddrvbd.f2
-rw-r--r--TESTING/EIG/derrgg.f24
-rw-r--r--TESTING/EIG/schkee.f5
-rw-r--r--TESTING/EIG/serrgg.f26
-rw-r--r--TESTING/EIG/zchkee.f3
-rw-r--r--TESTING/EIG/zerrgg.f40
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 )