*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2011
+*> \date June 2016
*
*> \ingroup complex_eig
*
$ Y, Z, Q, LDQ, PT, LDPT, U, VT, WORK, LWORK,
$ RWORK, NOUT, INFO )
*
-* -- LAPACK test routine (version 3.4.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2011
+* June 2016
*
* .. Scalar Arguments ..
INTEGER INFO, LDA, LDPT, LDQ, LDX, LWORK, NOUT, NRHS,
EXTERNAL SLAMCH, SLARND
* ..
* .. External Subroutines ..
- EXTERNAL ALASUM, CBDSQR, CBDT01, CBDT02, CBDT03, CGEBRD,
- $ CGEMM, CLACPY, CLASET, CLATMR, CLATMS, CUNGBR,
- $ CUNT01, SCOPY, SLABAD, SLAHD2, SSVDCH, XERBLA
+ EXTERNAL ALASUM, CBDSQR, CBDT01, CBDT02, CBDT03,
+ $ CGEBRD, CGEMM, CLACPY, CLASET, CLATMR,
+ $ CLATMS, CUNGBR, CUNT01, SCOPY, SLABAD,
+ $ SLAHD2, SSVDCH, XERBLA
* ..
* .. Intrinsic Functions ..
INTRINSIC ABS, EXP, INT, LOG, MAX, MIN, SQRT
COMMON / SRNAMC / SRNAMT
* ..
* .. Data statements ..
- DATA KTYPE / 1, 2, 5*4, 5*6, 3*9, 10 /
- DATA KMAGN / 2*1, 3*1, 2, 3, 3*1, 2, 3, 1, 2, 3, 0 /
- DATA KMODE / 2*0, 4, 3, 1, 4, 4, 4, 3, 1, 4, 4, 0,
+ DATA KTYPE / 1, 2, 5*4, 5*6, 3*9, 10 /
+ DATA KMAGN / 2*1, 3*1, 2, 3, 3*1, 2, 3, 1, 2, 3, 0 /
+ DATA KMODE / 2*0, 4, 3, 1, 4, 4, 4, 3, 1, 4, 4, 0,
$ 0, 0, 0 /
* ..
* .. Executable Statements ..
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2015
+*> \date June 2016
*
*> \ingroup complex_eig
*
* =====================================================================
PROGRAM CCHKEE
*
-* -- LAPACK test routine (version 3.6.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2015
+* June 2016
*
* =====================================================================
*
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
*
* 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 ),
*
* 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 ),
* 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,
*> and each type of matrix, one matrix will be generated and used
*> to test the nonsymmetric eigenroutines. For each matrix, 13
*> tests will be performed. The first twelve "test ratios" should be
-*> small -- O(1). They will be compared with the threshhold THRESH:
+*> small -- O(1). They will be compared with the threshold THRESH:
*>
*> H
*> (1) | A - U H V | / ( |A| n ulp )
*> |beta(Q,Z computed) - beta(Q,Z not computed)|/|P| ) / ulp
*>
*> In addition, the normalization of L and R are checked, and compared
-*> with the threshhold THRSHN.
+*> with the threshold THRSHN.
*>
*> Test Matrices
*> ---- --------
*> \param[in] THRSHN
*> \verbatim
*> THRSHN is REAL
-*> Threshhold for reporting eigenvector normalization error.
+*> Threshold for reporting eigenvector normalization error.
*> If the normalization of any eigenvector differs from 1 by
*> more than THRSHN*ulp, then a special error message will be
*> printed. (This is handled separately from the other tests,
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2011
+*> \date June 2016
*
*> \ingroup complex_eig
*
$ ALPHA3, BETA3, EVECTL, EVECTR, WORK, LWORK,
$ RWORK, LLWORK, RESULT, INFO )
*
-* -- LAPACK test routine (version 3.4.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2011
+* June 2016
*
* .. Scalar Arguments ..
LOGICAL TSTDIF
*> number of matrix "TYPES" are specified. For each size ("N")
*> and each TYPE of matrix, a pair of matrices (A, B) will be generated
*> and used for testing. For each matrix pair, the following 13 tests
-*> will be performed and compared with the threshhold THRESH except
+*> will be performed and compared with the threshold THRESH except
*> the tests (5), (11) and (13).
*>
*>
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2011
+*> \date June 2016
*
*> \ingroup complex_eig
*
$ NOUNIT, A, LDA, B, S, T, Q, LDQ, Z, ALPHA,
$ BETA, WORK, LWORK, RWORK, RESULT, BWORK, INFO )
*
-* -- LAPACK test routine (version 3.4.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2011
+* June 2016
*
* .. Scalar Arguments ..
INTEGER INFO, LDA, LDQ, LWORK, NOUNIT, NSIZES, NTYPES
*> number of matrix "TYPES" are specified. For each size ("N")
*> and each TYPE of matrix, a pair of matrices (A, B) will be generated
*> and used for testing. For each matrix pair, the following 13 tests
-*> will be performed and compared with the threshhold THRESH except
+*> will be performed and compared with the threshold THRESH except
*> the tests (5), (11) and (13).
*>
*>
$ BETA, WORK, LWORK, RWORK, RESULT, BWORK,
$ INFO )
*
-* -- LAPACK test routine (version 3.6.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
* February 2015
*> number of matrix "types" are specified. For each size ("n")
*> and each type of matrix, a pair of matrices (A, B) will be generated
*> and used for testing. For each matrix pair, the following tests
-*> will be performed and compared with the threshhold THRESH.
+*> will be performed and compared with the threshold THRESH.
*>
*> Results from CGGEV:
*>
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2015
+*> \date June 2016
*
*> \ingroup complex_eig
*
$ ALPHA, BETA, ALPHA1, BETA1, WORK, LWORK, RWORK,
$ RESULT, INFO )
*
-* -- LAPACK test routine (version 3.6.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2015
+* June 2016
*
* .. Scalar Arguments ..
INTEGER INFO, LDA, LDQ, LDQE, LWORK, NOUNIT, NSIZES,
*> number of matrix "types" are specified. For each size ("n")
*> and each type of matrix, a pair of matrices (A, B) will be generated
*> and used for testing. For each matrix pair, the following tests
-*> will be performed and compared with the threshhold THRESH.
+*> will be performed and compared with the threshold THRESH.
*>
*> Results from CGGEV3:
*>
$ ALPHA, BETA, ALPHA1, BETA1, WORK, LWORK,
$ RWORK, RESULT, INFO )
*
-* -- LAPACK test routine (version 3.6.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
* January 2015
*> (need more details on what kind of read-in data are needed).
*>
*> For each matrix pair, the following tests will be performed and
-*> compared with the threshhold THRESH except for the tests (7) and (9):
+*> compared with the threshold THRESH except for the tests (7) and (9):
*>
*> (1) | A - Q S Z' | / ( |A| n ulp )
*>
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2011
+*> \date June 2016
*
*> \ingroup complex_eig
*
$ AI, BI, Z, Q, ALPHA, BETA, C, LDC, S, WORK,
$ LWORK, RWORK, IWORK, LIWORK, BWORK, INFO )
*
-* -- LAPACK test routine (version 3.4.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2011
+* June 2016
*
* .. Scalar Arguments ..
INTEGER INFO, LDA, LDC, LIWORK, LWORK, NCMAX, NIN,
*> corresponding the first and last eigenvalues are also know
*> ``exactly'' (see CLATM6).
*> For each matrix pair, the following tests will be performed and
-*> compared with the threshhold THRESH.
+*> compared with the threshold THRESH.
*>
*> (1) max over all left eigenvalue/-vector pairs (beta/alpha,l) of
*>
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2011
+*> \date June 2016
*
*> \ingroup complex_eig
*
$ S, STRU, DIF, DIFTRU, WORK, LWORK, RWORK,
$ IWORK, LIWORK, RESULT, BWORK, INFO )
*
-* -- LAPACK test routine (version 3.4.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2011
+* June 2016
*
* .. Scalar Arguments ..
INTEGER IHI, ILO, INFO, LDA, LIWORK, LWORK, NIN, NOUT,
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2015
+*> \date June 2016
*
*> \ingroup complex_eig
*
$ SSAV, E, WORK, LWORK, RWORK, IWORK, NOUNIT,
$ INFO )
*
-* -- LAPACK test routine (version 3.6.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2015
+* June 2016
*
* .. Scalar Arguments ..
INTEGER INFO, LDA, LDU, LDVT, LWORK, NOUNIT, NSIZES,
* .. Local Scalars ..
LOGICAL BADMM, BADNN
CHARACTER JOBQ, JOBU, JOBVT, RANGE
- INTEGER I, IINFO, IJQ, IJU, IJVT, IL, IU, ITEMP, IWSPC,
- $ IWTMP, J, JSIZE, JTYPE, LSWORK, M, MINWRK,
- $ MMAX, MNMAX, MNMIN, MTYPES, N, NERRS, NFAIL,
- $ NMAX, NS, NSI, NSV, NTEST, NTESTF, NTESTT,
- $ LRWORK
+ INTEGER I, IINFO, IJQ, IJU, IJVT, IL, IU, ITEMP,
+ $ IWSPC, IWTMP, J, JSIZE, JTYPE, LSWORK, M,
+ $ MINWRK, MMAX, MNMAX, MNMIN, MTYPES, N,
+ $ NERRS, NFAIL, NMAX, NS, NSI, NSV, NTEST,
+ $ NTESTF, NTESTT, LRWORK
REAL ANORM, DIF, DIV, OVFL, RTUNFL, ULP, ULPINV,
$ UNFL, VL, VU
* ..
EXTERNAL SLAMCH, SLARND
* ..
* .. External Subroutines ..
- EXTERNAL ALASVM, XERBLA, CBDT01, CBDT05, CGESDD, CGESVD,
- $ CGESVJ, CGEJSV, CGESVDX, CLACPY, CLASET, CLATMS,
- $ CUNT01, CUNT03
+ EXTERNAL ALASVM, XERBLA, CBDT01, CBDT05, CGESDD,
+ $ CGESVD, CGESVJ, CGEJSV, CGESVDX, CLACPY,
+ $ CLASET, CLATMS, CUNT01, CUNT03
* ..
* .. Intrinsic Functions ..
INTRINSIC ABS, REAL, MAX, MIN
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2011
+*> \date June 2016
*
*> \ingroup complex_eig
*
$ NOUNIT, A, LDA, H, HT, W, WT, VS, LDVS, RESULT,
$ WORK, NWORK, RWORK, IWORK, BWORK, INFO )
*
-* -- LAPACK test routine (version 3.4.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2011
+* June 2016
*
* .. Scalar Arguments ..
INTEGER INFO, LDA, LDVS, NOUNIT, NSIZES, NTYPES, NWORK
$ ' 1/ulp otherwise', /
$ ' 12 = 0 if W same no matter if VS computed (sort),',
$ ' 1/ulp otherwise', /
- $ ' 13 = 0 if sorting succesful, 1/ulp otherwise', / )
+ $ ' 13 = 0 if sorting successful, 1/ulp otherwise', / )
9993 FORMAT( ' N=', I5, ', IWK=', I2, ', seed=', 4( I4, ',' ),
$ ' type ', I2, ', test(', I2, ')=', G10.3 )
9992 FORMAT( ' CDRVES: ', A, ' returned INFO=', I6, '.', / 9X, 'N=',
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2011
+*> \date June 2016
*
*> \ingroup complex_eig
*
$ LDVS, VS1, RESULT, WORK, LWORK, RWORK, BWORK,
$ INFO )
*
-* -- LAPACK test routine (version 3.4.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2011
+* June 2016
*
* .. Scalar Arguments ..
INTEGER INFO, LDA, LDVS, LWORK, NIUNIT, NOUNIT, NSIZES,
$ ' 1/ulp otherwise', /
$ ' 12 = 0 if W same no matter what else computed ',
$ '(sort), 1/ulp otherwise', /
- $ ' 13 = 0 if sorting succesful, 1/ulp otherwise',
+ $ ' 13 = 0 if sorting successful, 1/ulp otherwise',
$ / ' 14 = 0 if RCONDE same no matter what else computed,',
$ ' 1/ulp otherwise', /
$ ' 15 = 0 if RCONDv same no matter what else computed,',
*> \verbatim
*> INFO is INTEGER
*> If 0, then successful exit.
-*> If <0, then input paramter -INFO is incorrect.
+*> If <0, then input parameter -INFO is incorrect.
*> If >0, CLATMR, CLATMS, CLATME or CGET23 returned an error
*> code, and INFO is its absolute value.
*>
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2011
+*> \date June 2016
*
*> \ingroup complex_eig
*
$ RCONDE, RCNDE1, RCDEIN, SCALE, SCALE1, RESULT,
$ WORK, NWORK, RWORK, INFO )
*
-* -- LAPACK test routine (version 3.4.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2011
+* June 2016
*
* .. Scalar Arguments ..
INTEGER INFO, LDA, LDLRE, LDVL, LDVR, NIUNIT, NOUNIT,
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2015
+*> \date June 2016
*
*> \ingroup complex_eig
*
* =====================================================================
SUBROUTINE CERRED( PATH, NUNIT )
*
-* -- LAPACK test routine (version 3.6.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2015
+* June 2016
*
* .. Scalar Arguments ..
CHARACTER*3 PATH
$ VT( NMAX, NMAX ), W( 10*NMAX ), X( NMAX )
* ..
* .. External Subroutines ..
- EXTERNAL CHKXER, CGEES, CGEESX, CGEEV, CGEEVX, CGEJSV
+ EXTERNAL CHKXER, CGEES, CGEESX, CGEEV, CGEEVX, CGEJSV,
$ CGESDD, CGESVD
* ..
* .. External Functions ..
CALL CGESVDX( 'V', 'N', 'A', 2, 2, A, 2, ZERO, ZERO,
$ 0, 0, NS, S, U, 1, VT, 1, W, 1, RW, IW, INFO )
CALL CHKXER( 'CGESVDX', INFOT, NOUT, LERR, OK )
- INFOT = 16
+ INFOT = 17
CALL CGESVDX( 'N', 'V', 'A', 2, 2, A, 2, ZERO, ZERO,
$ 0, 0, NS, S, U, 1, VT, 1, W, 1, RW, IW, INFO )
CALL CHKXER( 'CGESVDX', INFOT, NOUT, LERR, OK )
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2015
+*> \date June 2016
*
*> \ingroup complex_eig
*
* =====================================================================
SUBROUTINE CERRGG( PATH, NUNIT )
*
-* -- LAPACK test routine (version 3.6.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2015
+* June 2016
*
* .. Scalar Arguments ..
CHARACTER*3 PATH
* ..
* .. 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 ),
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
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 )
*> eigenvalue is negative), and otherwise it returns .FALSE..
*>
*> It is used by the test routine CDRGES to test whether the driver
-*> routine CGGES succesfully sorts eigenvalues.
+*> routine CGGES successfully sorts eigenvalues.
*> \endverbatim
*
* Arguments:
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2011
+*> \date June 2016
*
*> \ingroup complex_eig
*
* =====================================================================
LOGICAL FUNCTION CLCTES( Z, D )
*
-* -- LAPACK test routine (version 3.4.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2011
+* June 2016
*
* .. Scalar Arguments ..
COMPLEX D, Z
*>
*> CSLECT returns .TRUE. if the eigenvalue Z is to be selected,
*> otherwise it returns .FALSE.
-*> It is used by CCHK41 to test if CGEES succesfully sorts eigenvalues,
-*> and by CCHK43 to test if CGEESX succesfully sorts eigenvalues.
+*> It is used by CCHK41 to test if CGEES successfully sorts eigenvalues,
+*> and by CCHK43 to test if CGEESX successfully sorts eigenvalues.
*>
*> The common block /SSLCT/ controls how eigenvalues are selected.
*> If SELOPT = 0, then CSLECT return .TRUE. when real(Z) is less than
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2011
+*> \date June 2016
*
*> \ingroup complex_eig
*
* =====================================================================
LOGICAL FUNCTION CSLECT( Z )
*
-* -- LAPACK test routine (version 3.4.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2011
+* June 2016
*
* .. Scalar Arguments ..
COMPLEX Z
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2015
+*> \date June 2016
*
*> \ingroup double_eig
*
$ Y, Z, Q, LDQ, PT, LDPT, U, VT, WORK, LWORK,
$ IWORK, NOUT, INFO )
*
-* -- LAPACK test routine (version 3.6.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2015
+* June 2016
*
* .. Scalar Arguments ..
INTEGER INFO, LDA, LDPT, LDQ, LDX, LWORK, NOUT, NRHS,
EXTERNAL DLAMCH, DLARND, DSXT1
* ..
* .. External Subroutines ..
- EXTERNAL ALASUM, DBDSDC, DBDSQR, DBDSVDX, DBDT01, DBDT02,
- $ DBDT03, DBDT04, DCOPY, DGEBRD, DGEMM, DLABAD,
- $ DLACPY, DLAHD2, DLASET, DLATMR, DLATMS,
- $ DORGBR, DORT01, XERBLA
+ EXTERNAL ALASUM, DBDSDC, DBDSQR, DBDSVDX, DBDT01,
+ $ DBDT02, DBDT03, DBDT04, DCOPY, DGEBRD,
+ $ DGEMM, DLABAD, DLACPY, DLAHD2, DLASET,
+ $ DLATMR, DLATMS, DORGBR, DORT01, XERBLA
* ..
* .. Intrinsic Functions ..
INTRINSIC ABS, EXP, INT, LOG, MAX, MIN, SQRT
COMMON / SRNAMC / SRNAMT
* ..
* .. Data statements ..
- DATA KTYPE / 1, 2, 5*4, 5*6, 3*9, 10 /
- DATA KMAGN / 2*1, 3*1, 2, 3, 3*1, 2, 3, 1, 2, 3, 0 /
- DATA KMODE / 2*0, 4, 3, 1, 4, 4, 4, 3, 1, 4, 4, 0,
- $ 0, 0, 0 /
+ DATA KTYPE / 1, 2, 5*4, 5*6, 3*9, 10 /
+ DATA KMAGN / 2*1, 3*1, 2, 3, 3*1, 2, 3, 1, 2, 3, 0 /
+ DATA KMODE / 2*0, 4, 3, 1, 4, 4, 4, 3, 1, 4, 4, 0,
+ $ 0, 0, 0 /
* ..
* .. Executable Statements ..
*
IWBD = IWBS + MNMIN
IWBE = IWBD + MNMIN
IWBZ = IWBE + MNMIN
- IWWORK = IWBZ + MNMIN*(MNMIN*2+1)
+ IWWORK = IWBZ + 2*MNMIN*(MNMIN+1)
MNMIN2 = MAX( 1,MNMIN*2 )
*
CALL DCOPY( MNMIN, BD, 1, WORK( IWBD ), 1 )
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2015
+*> \date June 2016
*
*> \ingroup double_eig
*
* =====================================================================
PROGRAM DCHKEE
*
-* -- LAPACK test routine (version 3.6.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2015
+* June 2016
*
* =====================================================================
*
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
*
* 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 ),
* 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,
*> and each type of matrix, one matrix will be generated and used
*> to test the nonsymmetric eigenroutines. For each matrix, 15
*> tests will be performed. The first twelve "test ratios" should be
-*> small -- O(1). They will be compared with the threshhold THRESH:
+*> small -- O(1). They will be compared with the threshold THRESH:
*>
*> T
*> (1) | A - U H V | / ( |A| n ulp )
*> |beta(Q,Z computed) - beta(Q,Z not computed)|/|P| ) / ulp
*>
*> In addition, the normalization of L and R are checked, and compared
-*> with the threshhold THRSHN.
+*> with the threshold THRSHN.
*>
*> Test Matrices
*> ---- --------
*> \param[in] THRSHN
*> \verbatim
*> THRSHN is DOUBLE PRECISION
-*> Threshhold for reporting eigenvector normalization error.
+*> Threshold for reporting eigenvector normalization error.
*> If the normalization of any eigenvector differs from 1 by
*> more than THRSHN*ulp, then a special error message will be
*> printed. (This is handled separately from the other tests,
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2011
+*> \date June 2016
*
*> \ingroup double_eig
*
$ BETA1, ALPHR3, ALPHI3, BETA3, EVECTL, EVECTR,
$ WORK, LWORK, LLWORK, RESULT, INFO )
*
-* -- LAPACK test routine (version 3.4.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2011
+* June 2016
*
* .. Scalar Arguments ..
LOGICAL TSTDIF
*> number of matrix "TYPES" are specified. For each size ("N")
*> and each TYPE of matrix, a pair of matrices (A, B) will be generated
*> and used for testing. For each matrix pair, the following 13 tests
-*> will be performed and compared with the threshhold THRESH except
+*> will be performed and compared with the threshold THRESH except
*> the tests (5), (11) and (13).
*>
*>
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2011
+*> \date June 2016
*
*> \ingroup double_eig
*
$ ALPHAI, BETA, WORK, LWORK, RESULT, BWORK,
$ INFO )
*
-* -- LAPACK test routine (version 3.4.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2011
+* June 2016
*
* .. Scalar Arguments ..
INTEGER INFO, LDA, LDQ, LWORK, NOUNIT, NSIZES, NTYPES
*> number of matrix "TYPES" are specified. For each size ("N")
*> and each TYPE of matrix, a pair of matrices (A, B) will be generated
*> and used for testing. For each matrix pair, the following 13 tests
-*> will be performed and compared with the threshhold THRESH except
+*> will be performed and compared with the threshold THRESH except
*> the tests (5), (11) and (13).
*>
*>
$ ALPHAI, BETA, WORK, LWORK, RESULT, BWORK,
$ INFO )
*
-* -- LAPACK test routine (version 3.6.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
* February 2015
*> number of matrix "types" are specified. For each size ("n")
*> and each type of matrix, a pair of matrices (A, B) will be generated
*> and used for testing. For each matrix pair, the following tests
-*> will be performed and compared with the threshhold THRESH.
+*> will be performed and compared with the threshold THRESH.
*>
*> Results from DGGEV:
*>
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2015
+*> \date June 2016
*
*> \ingroup double_eig
*
$ ALPHAR, ALPHAI, BETA, ALPHR1, ALPHI1, BETA1,
$ WORK, LWORK, RESULT, INFO )
*
-* -- LAPACK test routine (version 3.6.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2015
+* June 2016
*
* .. Scalar Arguments ..
INTEGER INFO, LDA, LDQ, LDQE, LWORK, NOUNIT, NSIZES,
*> number of matrix "types" are specified. For each size ("n")
*> and each type of matrix, a pair of matrices (A, B) will be generated
*> and used for testing. For each matrix pair, the following tests
-*> will be performed and compared with the threshhold THRESH.
+*> will be performed and compared with the threshold THRESH.
*>
*> Results from DGGEV3:
*>
$ ALPHAR, ALPHAI, BETA, ALPHR1, ALPHI1, BETA1,
$ WORK, LWORK, RESULT, INFO )
*
-* -- LAPACK test routine (version 3.6.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
* February 2015
*> to test DGGESX.
*>
*> For each matrix pair, the following tests will be performed and
-*> compared with the threshhold THRESH except for the tests (7) and (9):
+*> compared with the threshold THRESH except for the tests (7) and (9):
*>
*> (1) | A - Q S Z' | / ( |A| n ulp )
*>
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2011
+*> \date June 2016
*
*> \ingroup double_eig
*
$ BI, Z, Q, ALPHAR, ALPHAI, BETA, C, LDC, S,
$ WORK, LWORK, IWORK, LIWORK, BWORK, INFO )
*
-* -- LAPACK test routine (version 3.4.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2011
+* June 2016
*
* .. Scalar Arguments ..
INTEGER INFO, LDA, LDC, LIWORK, LWORK, NCMAX, NIN,
*> ``exactly'' (see DLATM6).
*>
*> For each matrix pair, the following tests will be performed and
-*> compared with the threshhold THRESH.
+*> compared with the threshold THRESH.
*>
*> (1) max over all left eigenvalue/-vector pairs (beta/alpha,l) of
*>
$ RSCALE, S, DTRU, DIF, DIFTRU, WORK, LWORK,
$ IWORK, LIWORK, RESULT, BWORK, INFO )
*
-* -- LAPACK test routine (version 3.4.1) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
* April 2012
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2015
+*> \date June 2016
*
*> \ingroup double_eig
*
$ A, LDA, U, LDU, VT, LDVT, ASAV, USAV, VTSAV, S,
$ SSAV, E, WORK, LWORK, IWORK, NOUT, INFO )
*
-* -- LAPACK test routine (version 3.6.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2015
+* June 2016
*
* .. Scalar Arguments ..
INTEGER INFO, LDA, LDU, LDVT, LWORK, NOUT, NSIZES,
* ..
* .. External Subroutines ..
EXTERNAL ALASVM, DBDT01, DGEJSV, DGESDD, DGESVD,
- $ DGESVDX, DGESVJ, DLABAD, DLACPY, DLASET, DLATMS,
- $ DORT01, DORT03, XERBLA
+ $ DGESVDX, DGESVJ, DLABAD, DLACPY, DLASET,
+ $ DLATMS, DORT01, DORT03, XERBLA
* ..
* .. Intrinsic Functions ..
INTRINSIC ABS, DBLE, INT, MAX, MIN
$ / '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 ) ',
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2011
+*> \date June 2016
*
*> \ingroup double_eig
*
$ NOUNIT, A, LDA, H, HT, WR, WI, WRT, WIT, VS,
$ LDVS, RESULT, WORK, NWORK, IWORK, BWORK, INFO )
*
-* -- LAPACK test routine (version 3.4.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2011
+* June 2016
*
* .. Scalar Arguments ..
INTEGER INFO, LDA, LDVS, NOUNIT, NSIZES, NTYPES, NWORK
$ ' 1/ulp otherwise', /
$ ' 12 = 0 if WR, WI same no matter if VS computed (sort),',
$ ' 1/ulp otherwise', /
- $ ' 13 = 0 if sorting succesful, 1/ulp otherwise', / )
+ $ ' 13 = 0 if sorting successful, 1/ulp otherwise', / )
9993 FORMAT( ' N=', I5, ', IWK=', I2, ', seed=', 4( I4, ',' ),
$ ' type ', I2, ', test(', I2, ')=', G10.3 )
9992 FORMAT( ' DDRVES: ', A, ' returned INFO=', I6, '.', / 9X, 'N=',
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2011
+*> \date June 2016
*
*> \ingroup double_eig
*
$ WIT, WRTMP, WITMP, VS, LDVS, VS1, RESULT, WORK,
$ LWORK, IWORK, BWORK, INFO )
*
-* -- LAPACK test routine (version 3.4.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2011
+* June 2016
*
* .. Scalar Arguments ..
INTEGER INFO, LDA, LDVS, LWORK, NIUNIT, NOUNIT, NSIZES,
$ ' 1/ulp otherwise', /
$ ' 12 = 0 if WR, WI same no matter what else computed ',
$ '(sort), 1/ulp otherwise', /
- $ ' 13 = 0 if sorting succesful, 1/ulp otherwise',
+ $ ' 13 = 0 if sorting successful, 1/ulp otherwise',
$ / ' 14 = 0 if RCONDE same no matter what else computed,',
$ ' 1/ulp otherwise', /
$ ' 15 = 0 if RCONDv same no matter what else computed,',
*> \verbatim
*> INFO is INTEGER
*> If 0, then successful exit.
-*> If <0, then input paramter -INFO is incorrect.
+*> If <0, then input parameter -INFO is incorrect.
*> If >0, DLATMR, SLATMS, SLATME or DGET23 returned an error
*> code, and INFO is its absolute value.
*>
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2011
+*> \date June 2016
*
*> \ingroup double_eig
*
$ RCDVIN, RCONDE, RCNDE1, RCDEIN, SCALE, SCALE1,
$ RESULT, WORK, NWORK, IWORK, INFO )
*
-* -- LAPACK test routine (version 3.4.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2011
+* June 2016
*
* .. Scalar Arguments ..
INTEGER INFO, LDA, LDLRE, LDVL, LDVR, NIUNIT, NOUNIT,
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2015
+*> \date June 2016
*
*> \ingroup double_eig
*
* =====================================================================
SUBROUTINE DERRBD( PATH, NUNIT )
*
-* -- LAPACK test routine (version 3.6.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2015
+* June 2016
*
* .. Scalar Arguments ..
CHARACTER*3 PATH
EXTERNAL LSAMEN
* ..
* .. External Subroutines ..
- EXTERNAL CHKXER, DBDSDC, DBDSQR, DBDSVDX, DGEBD2, DGEBRD, DORGBR,
- $ DORMBR
+ EXTERNAL CHKXER, DBDSDC, DBDSQR, DBDSVDX, DGEBD2,
+ $ DGEBRD, DORGBR, DORMBR
* ..
* .. Scalars in Common ..
LOGICAL LERR, OK
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2015
+*> \date June 2016
*
*> \ingroup double_eig
*
* =====================================================================
SUBROUTINE DERRED( PATH, NUNIT )
*
-* -- LAPACK test routine (version 3.6.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2015
+* June 2016
*
* .. Scalar Arguments ..
CHARACTER*3 PATH
CALL DGESVDX( 'V', 'N', 'A', 2, 2, A, 2, ZERO, ZERO,
$ 0, 0, NS, S, U, 1, VT, 1, W, 1, IW, INFO )
CALL CHKXER( 'DGESVDX', INFOT, NOUT, LERR, OK )
- INFOT = 16
+ INFOT = 17
CALL DGESVDX( 'N', 'V', 'A', 2, 2, A, 2, ZERO, ZERO,
$ 0, 0, NS, S, U, 1, VT, 1, W, 1, IW, INFO )
CALL CHKXER( 'DGESVDX', INFOT, NOUT, LERR, OK )
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2015
+*> \date June 2016
*
*> \ingroup double_eig
*
* =====================================================================
SUBROUTINE DERRGG( PATH, NUNIT )
*
-* -- LAPACK test routine (version 3.6.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2015
+* June 2016
*
* .. Scalar Arguments ..
CHARACTER*3 PATH
* ..
* .. 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 ),
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
*
*> .FALSE..
*>
*> It is used by the test routine DDRGES to test whether the driver
-*> routine DGGES succesfully sorts eigenvalues.
+*> routine DGGES successfully sorts eigenvalues.
*> \endverbatim
*
* Arguments:
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2011
+*> \date June 2016
*
*> \ingroup double_eig
*
* =====================================================================
LOGICAL FUNCTION DLCTES( ZR, ZI, D )
*
-* -- LAPACK test routine (version 3.4.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2011
+* June 2016
*
* .. Scalar Arguments ..
DOUBLE PRECISION D, ZI, ZR
*>
*> DSLECT returns .TRUE. if the eigenvalue ZR+sqrt(-1)*ZI is to be
*> selected, and otherwise it returns .FALSE.
-*> It is used by DCHK41 to test if DGEES succesfully sorts eigenvalues,
-*> and by DCHK43 to test if DGEESX succesfully sorts eigenvalues.
+*> It is used by DCHK41 to test if DGEES successfully sorts eigenvalues,
+*> and by DCHK43 to test if DGEESX successfully sorts eigenvalues.
*>
*> The common block /SSLCT/ controls how eigenvalues are selected.
*> If SELOPT = 0, then DSLECT return .TRUE. when ZR is less than zero,
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2011
+*> \date June 2016
*
*> \ingroup double_eig
*
* =====================================================================
LOGICAL FUNCTION DSLECT( ZR, ZI )
*
-* -- LAPACK test routine (version 3.4.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2011
+* June 2016
*
* .. Scalar Arguments ..
DOUBLE PRECISION ZI, ZR
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2015
+*> \date June 2016
*
*> \ingroup single_eig
*
$ Y, Z, Q, LDQ, PT, LDPT, U, VT, WORK, LWORK,
$ IWORK, NOUT, INFO )
*
-* -- LAPACK test routine (version 3.6.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2015
+* June 2016
*
* .. Scalar Arguments ..
INTEGER INFO, LDA, LDPT, LDQ, LDX, LWORK, NOUT, NRHS,
CHARACTER*3 PATH
INTEGER I, IINFO, IL, IMODE, ITEMP, ITYPE, IU, IWBD,
$ IWBE, IWBS, IWBZ, IWWORK, J, JCOL, JSIZE,
- $ JTYPE, LOG2UI, M, MINWRK, MMAX, MNMAX, MNMIN,
+ $ JTYPE, LOG2UI, M, MINWRK, MMAX, MNMAX, MNMIN,
$ MNMIN2, MQ, MTYPES, N, NFAIL, NMAX,
$ NS1, NS2, NTEST
REAL ABSTOL, AMNINV, ANORM, COND, OVFL, RTOVFL,
EXTERNAL SLAMCH, SLARND, SSXT1
* ..
* .. External Subroutines ..
- EXTERNAL ALASUM, SBDSDC, SBDSQR, SBDSVDX, SBDT01, SBDT02,
- $ SBDT03, SBDT04, SCOPY, SGEBRD, SGEMM, SLABAD,
- $ SLACPY, SLAHD2, SLASET, SLATMR, SLATMS,
- $ SORGBR, SORT01, XERBLA
+ EXTERNAL ALASUM, SBDSDC, SBDSQR, SBDSVDX, SBDT01,
+ $ SBDT02, SBDT03, SBDT04, SCOPY, SGEBRD,
+ $ SGEMM, SLABAD, SLACPY, SLAHD2, SLASET,
+ $ SLATMR, SLATMS, SORGBR, SORT01, XERBLA
* ..
* .. Intrinsic Functions ..
INTRINSIC ABS, EXP, INT, LOG, MAX, MIN, SQRT
COMMON / SRNAMC / SRNAMT
* ..
* .. Data statements ..
- DATA KTYPE / 1, 2, 5*4, 5*6, 3*9, 10 /
- DATA KMAGN / 2*1, 3*1, 2, 3, 3*1, 2, 3, 1, 2, 3, 0 /
- DATA KMODE / 2*0, 4, 3, 1, 4, 4, 4, 3, 1, 4, 4, 0,
+ DATA KTYPE / 1, 2, 5*4, 5*6, 3*9, 10 /
+ DATA KMAGN / 2*1, 3*1, 2, 3, 3*1, 2, 3, 1, 2, 3, 0 /
+ DATA KMODE / 2*0, 4, 3, 1, 4, 4, 4, 3, 1, 4, 4, 0,
$ 0, 0, 0 /
* ..
* .. Executable Statements ..
IWBD = IWBS + MNMIN
IWBE = IWBD + MNMIN
IWBZ = IWBE + MNMIN
- IWWORK = IWBZ + MNMIN*(MNMIN*2+1)
+ IWWORK = IWBZ + 2*MNMIN*(MNMIN+1)
MNMIN2 = MAX( 1,MNMIN*2 )
*
CALL SCOPY( MNMIN, BD, 1, WORK( IWBD ), 1 )
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2015
+*> \date June 2016
*
*> \ingroup single_eig
*
* =====================================================================
PROGRAM SCHKEE
*
-* -- LAPACK test routine (version 3.6.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2015
+* June 2016
*
* =====================================================================
*
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 ) )
*
* 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 ),
* 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,
*> and each type of matrix, one matrix will be generated and used
*> to test the nonsymmetric eigenroutines. For each matrix, 15
*> tests will be performed. The first twelve "test ratios" should be
-*> small -- O(1). They will be compared with the threshhold THRESH:
+*> small -- O(1). They will be compared with the threshold THRESH:
*>
*> T
*> (1) | A - U H V | / ( |A| n ulp )
*> |beta(Q,Z computed) - beta(Q,Z not computed)|/|P| ) / ulp
*>
*> In addition, the normalization of L and R are checked, and compared
-*> with the threshhold THRSHN.
+*> with the threshold THRSHN.
*>
*> Test Matrices
*> ---- --------
*> \param[in] THRSHN
*> \verbatim
*> THRSHN is REAL
-*> Threshhold for reporting eigenvector normalization error.
+*> Threshold for reporting eigenvector normalization error.
*> If the normalization of any eigenvector differs from 1 by
*> more than THRSHN*ulp, then a special error message will be
*> printed. (This is handled separately from the other tests,
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2011
+*> \date June 2016
*
*> \ingroup single_eig
*
$ BETA1, ALPHR3, ALPHI3, BETA3, EVECTL, EVECTR,
$ WORK, LWORK, LLWORK, RESULT, INFO )
*
-* -- LAPACK test routine (version 3.4.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2011
+* June 2016
*
* .. Scalar Arguments ..
LOGICAL TSTDIF
*> number of matrix "TYPES" are specified. For each size ("N")
*> and each TYPE of matrix, a pair of matrices (A, B) will be generated
*> and used for testing. For each matrix pair, the following 13 tests
-*> will be performed and compared with the threshhold THRESH except
+*> will be performed and compared with the threshold THRESH except
*> the tests (5), (11) and (13).
*>
*>
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2011
+*> \date June 2016
*
*> \ingroup single_eig
*
$ ALPHAI, BETA, WORK, LWORK, RESULT, BWORK,
$ INFO )
*
-* -- LAPACK test routine (version 3.4.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2011
+* June 2016
*
* .. Scalar Arguments ..
INTEGER INFO, LDA, LDQ, LWORK, NOUNIT, NSIZES, NTYPES
*> number of matrix "TYPES" are specified. For each size ("N")
*> and each TYPE of matrix, a pair of matrices (A, B) will be generated
*> and used for testing. For each matrix pair, the following 13 tests
-*> will be performed and compared with the threshhold THRESH except
+*> will be performed and compared with the threshold THRESH except
*> the tests (5), (11) and (13).
*>
*>
$ ALPHAI, BETA, WORK, LWORK, RESULT, BWORK,
$ INFO )
*
-* -- LAPACK test routine (version 3.6.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
* February 2015
*> number of matrix "types" are specified. For each size ("n")
*> and each type of matrix, a pair of matrices (A, B) will be generated
*> and used for testing. For each matrix pair, the following tests
-*> will be performed and compared with the threshhold THRESH.
+*> will be performed and compared with the threshold THRESH.
*>
*> Results from SGGEV:
*>
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2015
+*> \date June 2016
*
*> \ingroup single_eig
*
$ ALPHAR, ALPHAI, BETA, ALPHR1, ALPHI1, BETA1,
$ WORK, LWORK, RESULT, INFO )
*
-* -- LAPACK test routine (version 3.6.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2015
+* June 2016
*
* .. Scalar Arguments ..
INTEGER INFO, LDA, LDQ, LDQE, LWORK, NOUNIT, NSIZES,
*> number of matrix "types" are specified. For each size ("n")
*> and each type of matrix, a pair of matrices (A, B) will be generated
*> and used for testing. For each matrix pair, the following tests
-*> will be performed and compared with the threshhold THRESH.
+*> will be performed and compared with the threshold THRESH.
*>
*> Results from SGGEV3:
*>
$ ALPHAR, ALPHAI, BETA, ALPHR1, ALPHI1, BETA1,
$ WORK, LWORK, RESULT, INFO )
*
-* -- LAPACK test routine (version 3.6.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
* February 2015
*> to test SGGESX.
*>
*> For each matrix pair, the following tests will be performed and
-*> compared with the threshhold THRESH except for the tests (7) and (9):
+*> compared with the threshold THRESH except for the tests (7) and (9):
*>
*> (1) | A - Q S Z' | / ( |A| n ulp )
*>
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2011
+*> \date June 2016
*
*> \ingroup single_eig
*
$ AI, BI, Z, Q, ALPHAR, ALPHAI, BETA, C, LDC, S,
$ WORK, LWORK, IWORK, LIWORK, BWORK, INFO )
*
-* -- LAPACK test routine (version 3.4.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2011
+* June 2016
*
* .. Scalar Arguments ..
INTEGER INFO, LDA, LDC, LIWORK, LWORK, NCMAX, NIN,
*> ``exactly'' (see SLATM6).
*>
*> For each matrix pair, the following tests will be performed and
-*> compared with the threshhold THRESH.
+*> compared with the threshold THRESH.
*>
*> (1) max over all left eigenvalue/-vector pairs (beta/alpha,l) of
*>
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2011
+*> \date June 2016
*
*> \ingroup single_eig
*
$ RSCALE, S, STRU, DIF, DIFTRU, WORK, LWORK,
$ IWORK, LIWORK, RESULT, BWORK, INFO )
*
-* -- LAPACK test routine (version 3.4.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2011
+* June 2016
*
* .. Scalar Arguments ..
INTEGER IHI, ILO, INFO, LDA, LIWORK, LWORK, NIN, NOUT,
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2015
+*> \date June 2016
*
*> \ingroup single_eig
*
$ A, LDA, U, LDU, VT, LDVT, ASAV, USAV, VTSAV, S,
$ SSAV, E, WORK, LWORK, IWORK, NOUT, INFO )
*
-* -- LAPACK test routine (version 3.6.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2015
+* June 2016
*
* .. Scalar Arguments ..
INTEGER INFO, LDA, LDU, LDVT, LWORK, NOUT, NSIZES,
* ..
* .. External Subroutines ..
EXTERNAL ALASVM, SBDT01, SGEJSV, SGESDD, SGESVD,
- $ SGESVDX, SGESVJ, SLABAD, SLACPY, SLASET, SLATMS,
- $ SORT01, SORT03, XERBLA
+ $ SGESVDX, SGESVJ, SLABAD, SLACPY, SLASET,
+ $ SLATMS, SORT01, SORT03, XERBLA
* ..
* .. Intrinsic Functions ..
INTRINSIC ABS, REAL, INT, MAX, MIN
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2011
+*> \date June 2016
*
*> \ingroup single_eig
*
$ NOUNIT, A, LDA, H, HT, WR, WI, WRT, WIT, VS,
$ LDVS, RESULT, WORK, NWORK, IWORK, BWORK, INFO )
*
-* -- LAPACK test routine (version 3.4.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2011
+* June 2016
*
* .. Scalar Arguments ..
INTEGER INFO, LDA, LDVS, NOUNIT, NSIZES, NTYPES, NWORK
$ ' 1/ulp otherwise', /
$ ' 12 = 0 if WR, WI same no matter if VS computed (sort),',
$ ' 1/ulp otherwise', /
- $ ' 13 = 0 if sorting succesful, 1/ulp otherwise', / )
+ $ ' 13 = 0 if sorting successful, 1/ulp otherwise', / )
9993 FORMAT( ' N=', I5, ', IWK=', I2, ', seed=', 4( I4, ',' ),
$ ' type ', I2, ', test(', I2, ')=', G10.3 )
9992 FORMAT( ' SDRVES: ', A, ' returned INFO=', I6, '.', / 9X, 'N=',
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2011
+*> \date June 2016
*
*> \ingroup single_eig
*
$ WIT, WRTMP, WITMP, VS, LDVS, VS1, RESULT, WORK,
$ LWORK, IWORK, BWORK, INFO )
*
-* -- LAPACK test routine (version 3.4.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2011
+* June 2016
*
* .. Scalar Arguments ..
INTEGER INFO, LDA, LDVS, LWORK, NIUNIT, NOUNIT, NSIZES,
$ ' 1/ulp otherwise', /
$ ' 12 = 0 if WR, WI same no matter what else computed ',
$ '(sort), 1/ulp otherwise', /
- $ ' 13 = 0 if sorting succesful, 1/ulp otherwise',
+ $ ' 13 = 0 if sorting successful, 1/ulp otherwise',
$ / ' 14 = 0 if RCONDE same no matter what else computed,',
$ ' 1/ulp otherwise', /
$ ' 15 = 0 if RCONDv same no matter what else computed,',
*> \verbatim
*> INFO is INTEGER
*> If 0, then successful exit.
-*> If <0, then input paramter -INFO is incorrect.
+*> If <0, then input parameter -INFO is incorrect.
*> If >0, SLATMR, SLATMS, SLATME or SGET23 returned an error
*> code, and INFO is its absolute value.
*>
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2011
+*> \date June 2016
*
*> \ingroup single_eig
*
$ RCDVIN, RCONDE, RCNDE1, RCDEIN, SCALE, SCALE1,
$ RESULT, WORK, NWORK, IWORK, INFO )
*
-* -- LAPACK test routine (version 3.4.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2011
+* June 2016
*
* .. Scalar Arguments ..
INTEGER INFO, LDA, LDLRE, LDVL, LDVR, NIUNIT, NOUNIT,
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2015
+*> \date June 2016
*
*> \ingroup single_eig
*
* =====================================================================
SUBROUTINE SERRBD( PATH, NUNIT )
*
-* -- LAPACK test routine (version 3.6.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2015
+* June 2016
*
* .. Scalar Arguments ..
CHARACTER*3 PATH
EXTERNAL LSAMEN
* ..
* .. External Subroutines ..
- EXTERNAL CHKXER, SBDSDC, SBDSQR, SBDSVDX, SGEBD2, SGEBRD, SORGBR,
- $ SORMBR
+ EXTERNAL CHKXER, SBDSDC, SBDSQR, SBDSVDX, SGEBD2,
+ $ SGEBRD, SORGBR, SORMBR
* ..
* .. Scalars in Common ..
LOGICAL LERR, OK
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2015
+*> \date June 2016
*
*> \ingroup single_eig
*
* =====================================================================
SUBROUTINE SERRED( PATH, NUNIT )
*
-* -- LAPACK test routine (version 3.6.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2015
+* June 2016
*
* .. Scalar Arguments ..
CHARACTER*3 PATH
CALL SGESVDX( 'V', 'N', 'A', 2, 2, A, 2, ZERO, ZERO,
$ 0, 0, NS, S, U, 1, VT, 1, W, 1, IW, INFO )
CALL CHKXER( 'SGESVDX', INFOT, NOUT, LERR, OK )
- INFOT = 16
+ INFOT = 17
CALL SGESVDX( 'N', 'V', 'A', 2, 2, A, 2, ZERO, ZERO,
$ 0, 0, NS, S, U, 1, VT, 1, W, 1, IW, INFO )
CALL CHKXER( 'SGESVDX', INFOT, NOUT, LERR, OK )
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2015
+*> \date June 2016
*
*> \ingroup single_eig
*
* =====================================================================
SUBROUTINE SERRGG( PATH, NUNIT )
*
-* -- LAPACK test routine (version 3.6.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2015
+* June 2016
*
* .. Scalar Arguments ..
CHARACTER*3 PATH
* .. 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 ),
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
*
*> .FALSE..
*>
*> It is used by the test routine SDRGES to test whether the driver
-*> routine SGGES succesfully sorts eigenvalues.
+*> routine SGGES successfully sorts eigenvalues.
*> \endverbatim
*
* Arguments:
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2011
+*> \date June 2016
*
*> \ingroup single_eig
*
* =====================================================================
LOGICAL FUNCTION SLCTES( ZR, ZI, D )
*
-* -- LAPACK test routine (version 3.4.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2011
+* June 2016
*
* .. Scalar Arguments ..
REAL D, ZI, ZR
*>
*> SSLECT returns .TRUE. if the eigenvalue ZR+sqrt(-1)*ZI is to be
*> selected, and otherwise it returns .FALSE.
-*> It is used by SCHK41 to test if SGEES succesfully sorts eigenvalues,
-*> and by SCHK43 to test if SGEESX succesfully sorts eigenvalues.
+*> It is used by SCHK41 to test if SGEES successfully sorts eigenvalues,
+*> and by SCHK43 to test if SGEESX successfully sorts eigenvalues.
*>
*> The common block /SSLCT/ controls how eigenvalues are selected.
*> If SELOPT = 0, then SSLECT return .TRUE. when ZR is less than zero,
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2011
+*> \date June 2016
*
*> \ingroup single_eig
*
* =====================================================================
LOGICAL FUNCTION SSLECT( ZR, ZI )
*
-* -- LAPACK test routine (version 3.4.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2011
+* June 2016
*
* .. Scalar Arguments ..
REAL ZI, ZR
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2011
+*> \date June 2016
*
*> \ingroup complex16_eig
*
$ Y, Z, Q, LDQ, PT, LDPT, U, VT, WORK, LWORK,
$ RWORK, NOUT, INFO )
*
-* -- LAPACK test routine (version 3.4.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2011
+* June 2016
*
* .. Scalar Arguments ..
INTEGER INFO, LDA, LDPT, LDQ, LDX, LWORK, NOUT, NRHS,
COMMON / SRNAMC / SRNAMT
* ..
* .. Data statements ..
- DATA KTYPE / 1, 2, 5*4, 5*6, 3*9, 10 /
- DATA KMAGN / 2*1, 3*1, 2, 3, 3*1, 2, 3, 1, 2, 3, 0 /
- DATA KMODE / 2*0, 4, 3, 1, 4, 4, 4, 3, 1, 4, 4, 0,
+ DATA KTYPE / 1, 2, 5*4, 5*6, 3*9, 10 /
+ DATA KMAGN / 2*1, 3*1, 2, 3, 3*1, 2, 3, 1, 2, 3, 0 /
+ DATA KMODE / 2*0, 4, 3, 1, 4, 4, 4, 3, 1, 4, 4, 0,
$ 0, 0, 0 /
* ..
* .. Executable Statements ..
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2015
+*> \date June 2016
*
*> \ingroup complex16_eig
*
* =====================================================================
PROGRAM ZCHKEE
*
-* -- LAPACK test routine (version 3.6.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2015
+* June 2016
*
* =====================================================================
*
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
*
* 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 ),
* 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,
*> and each type of matrix, one matrix will be generated and used
*> to test the nonsymmetric eigenroutines. For each matrix, 13
*> tests will be performed. The first twelve "test ratios" should be
-*> small -- O(1). They will be compared with the threshhold THRESH:
+*> small -- O(1). They will be compared with the threshold THRESH:
*>
*> H
*> (1) | A - U H V | / ( |A| n ulp )
*> |beta(Q,Z computed) - beta(Q,Z not computed)|/|P| ) / ulp
*>
*> In addition, the normalization of L and R are checked, and compared
-*> with the threshhold THRSHN.
+*> with the threshold THRSHN.
*>
*> Test Matrices
*> ---- --------
*> \param[in] THRSHN
*> \verbatim
*> THRSHN is DOUBLE PRECISION
-*> Threshhold for reporting eigenvector normalization error.
+*> Threshold for reporting eigenvector normalization error.
*> If the normalization of any eigenvector differs from 1 by
*> more than THRSHN*ulp, then a special error message will be
*> printed. (This is handled separately from the other tests,
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2011
+*> \date June 2016
*
*> \ingroup complex16_eig
*
$ ALPHA3, BETA3, EVECTL, EVECTR, WORK, LWORK,
$ RWORK, LLWORK, RESULT, INFO )
*
-* -- LAPACK test routine (version 3.4.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2011
+* June 2016
*
* .. Scalar Arguments ..
LOGICAL TSTDIF
*> number of matrix "TYPES" are specified. For each size ("N")
*> and each TYPE of matrix, a pair of matrices (A, B) will be generated
*> and used for testing. For each matrix pair, the following 13 tests
-*> will be performed and compared with the threshhold THRESH except
+*> will be performed and compared with the threshold THRESH except
*> the tests (5), (11) and (13).
*>
*>
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2011
+*> \date June 2016
*
*> \ingroup complex16_eig
*
$ NOUNIT, A, LDA, B, S, T, Q, LDQ, Z, ALPHA,
$ BETA, WORK, LWORK, RWORK, RESULT, BWORK, INFO )
*
-* -- LAPACK test routine (version 3.4.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2011
+* June 2016
*
* .. Scalar Arguments ..
INTEGER INFO, LDA, LDQ, LWORK, NOUNIT, NSIZES, NTYPES
*> number of matrix "TYPES" are specified. For each size ("N")
*> and each TYPE of matrix, a pair of matrices (A, B) will be generated
*> and used for testing. For each matrix pair, the following 13 tests
-*> will be performed and compared with the threshhold THRESH except
+*> will be performed and compared with the threshold THRESH except
*> the tests (5), (11) and (13).
*>
*>
$ BETA, WORK, LWORK, RWORK, RESULT, BWORK,
$ INFO )
*
-* -- LAPACK test routine (version 3.6.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
* February 2015
*> number of matrix "types" are specified. For each size ("n")
*> and each type of matrix, a pair of matrices (A, B) will be generated
*> and used for testing. For each matrix pair, the following tests
-*> will be performed and compared with the threshhold THRESH.
+*> will be performed and compared with the threshold THRESH.
*>
*> Results from ZGGEV:
*>
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2015
+*> \date June 2016
*
*> \ingroup complex16_eig
*
$ ALPHA, BETA, ALPHA1, BETA1, WORK, LWORK, RWORK,
$ RESULT, INFO )
*
-* -- LAPACK test routine (version 3.6.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2015
+* June 2016
*
* .. Scalar Arguments ..
INTEGER INFO, LDA, LDQ, LDQE, LWORK, NOUNIT, NSIZES,
*> number of matrix "types" are specified. For each size ("n")
*> and each type of matrix, a pair of matrices (A, B) will be generated
*> and used for testing. For each matrix pair, the following tests
-*> will be performed and compared with the threshhold THRESH.
+*> will be performed and compared with the threshold THRESH.
*>
*> Results from ZGGEV3:
*>
$ ALPHA, BETA, ALPHA1, BETA1, WORK, LWORK,
$ RWORK, RESULT, INFO )
*
-* -- LAPACK test routine (version 3.6.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
* February 2015
*> (need more details on what kind of read-in data are needed).
*>
*> For each matrix pair, the following tests will be performed and
-*> compared with the threshhold THRESH except for the tests (7) and (9):
+*> compared with the threshold THRESH except for the tests (7) and (9):
*>
*> (1) | A - Q S Z' | / ( |A| n ulp )
*>
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2011
+*> \date June 2016
*
*> \ingroup complex16_eig
*
$ BI, Z, Q, ALPHA, BETA, C, LDC, S, WORK, LWORK,
$ RWORK, IWORK, LIWORK, BWORK, INFO )
*
-* -- LAPACK test routine (version 3.4.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2011
+* June 2016
*
* .. Scalar Arguments ..
INTEGER INFO, LDA, LDC, LIWORK, LWORK, NCMAX, NIN,
*> corresponding the first and last eigenvalues are also know
*> ``exactly'' (see ZLATM6).
*> For each matrix pair, the following tests will be performed and
-*> compared with the threshhold THRESH.
+*> compared with the threshold THRESH.
*>
*> (1) max over all left eigenvalue/-vector pairs (beta/alpha,l) of
*>
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2011
+*> \date June 2016
*
*> \ingroup complex16_eig
*
$ S, DTRU, DIF, DIFTRU, WORK, LWORK, RWORK,
$ IWORK, LIWORK, RESULT, BWORK, INFO )
*
-* -- LAPACK test routine (version 3.4.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2011
+* June 2016
*
* .. Scalar Arguments ..
INTEGER IHI, ILO, INFO, LDA, LIWORK, LWORK, NIN, NOUT,
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2015
+*> \date June 2016
*
*> \ingroup complex16_eig
*
$ SSAV, E, WORK, LWORK, RWORK, IWORK, NOUNIT,
$ INFO )
*
-* -- LAPACK test routine (version 3.6.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2015
+* June 2016
*
* .. Scalar Arguments ..
INTEGER INFO, LDA, LDU, LDVT, LWORK, NOUNIT, NSIZES,
* .. Local Scalars ..
LOGICAL BADMM, BADNN
CHARACTER JOBQ, JOBU, JOBVT, RANGE
- INTEGER I, IINFO, IJQ, IJU, IJVT, IL, IU, ITEMP, IWSPC,
- $ IWTMP, J, JSIZE, JTYPE, LSWORK, M, MINWRK,
- $ MMAX, MNMAX, MNMIN, MTYPES, N, NERRS, NFAIL,
- $ NMAX, NS, NSI, NSV, NTEST, NTESTF, NTESTT,
- $ LRWORK
+ INTEGER I, IINFO, IJQ, IJU, IJVT, IL, IU, ITEMP,
+ $ IWSPC, IWTMP, J, JSIZE, JTYPE, LSWORK, M,
+ $ MINWRK, MMAX, MNMAX, MNMIN, MTYPES, N,
+ $ NERRS, NFAIL, NMAX, NS, NSI, NSV, NTEST,
+ $ NTESTF, NTESTT, LRWORK
DOUBLE PRECISION ANORM, DIF, DIV, OVFL, RTUNFL, ULP, ULPINV,
$ UNFL, VL, VU
* ..
EXTERNAL DLAMCH, DLARND
* ..
* .. External Subroutines ..
- EXTERNAL ALASVM, XERBLA, ZBDT01, ZBDT05, ZGESDD, ZGESVD,
- $ ZGESVJ, ZGEJSV, ZGESVDX, ZLACPY, ZLASET, ZLATMS,
- $ ZUNT01, ZUNT03
+ EXTERNAL ALASVM, XERBLA, ZBDT01, ZBDT05, ZGESDD,
+ $ ZGESVD, ZGESVJ, ZGEJSV, ZGESVDX, ZLACPY,
+ $ ZLASET, ZLATMS, ZUNT01, ZUNT03
* ..
* .. Intrinsic Functions ..
INTRINSIC ABS, DBLE, MAX, MIN
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2011
+*> \date June 2016
*
*> \ingroup complex16_eig
*
$ NOUNIT, A, LDA, H, HT, W, WT, VS, LDVS, RESULT,
$ WORK, NWORK, RWORK, IWORK, BWORK, INFO )
*
-* -- LAPACK test routine (version 3.4.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2011
+* June 2016
*
* .. Scalar Arguments ..
INTEGER INFO, LDA, LDVS, NOUNIT, NSIZES, NTYPES, NWORK
$ ' 1/ulp otherwise', /
$ ' 12 = 0 if W same no matter if VS computed (sort),',
$ ' 1/ulp otherwise', /
- $ ' 13 = 0 if sorting succesful, 1/ulp otherwise', / )
+ $ ' 13 = 0 if sorting successful, 1/ulp otherwise', / )
9993 FORMAT( ' N=', I5, ', IWK=', I2, ', seed=', 4( I4, ',' ),
$ ' type ', I2, ', test(', I2, ')=', G10.3 )
9992 FORMAT( ' ZDRVES: ', A, ' returned INFO=', I6, '.', / 9X, 'N=',
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2011
+*> \date June 2016
*
*> \ingroup complex16_eig
*
$ LDVS, VS1, RESULT, WORK, LWORK, RWORK, BWORK,
$ INFO )
*
-* -- LAPACK test routine (version 3.4.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2011
+* June 2016
*
* .. Scalar Arguments ..
INTEGER INFO, LDA, LDVS, LWORK, NIUNIT, NOUNIT, NSIZES,
$ ' 1/ulp otherwise', /
$ ' 12 = 0 if W same no matter what else computed ',
$ '(sort), 1/ulp otherwise', /
- $ ' 13 = 0 if sorting succesful, 1/ulp otherwise',
+ $ ' 13 = 0 if sorting successful, 1/ulp otherwise',
$ / ' 14 = 0 if RCONDE same no matter what else computed,',
$ ' 1/ulp otherwise', /
$ ' 15 = 0 if RCONDv same no matter what else computed,',
*> \verbatim
*> INFO is INTEGER
*> If 0, then successful exit.
-*> If <0, then input paramter -INFO is incorrect.
+*> If <0, then input parameter -INFO is incorrect.
*> If >0, ZLATMR, CLATMS, CLATME or ZGET23 returned an error
*> code, and INFO is its absolute value.
*>
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2011
+*> \date June 2016
*
*> \ingroup complex16_eig
*
$ RCONDE, RCNDE1, RCDEIN, SCALE, SCALE1, RESULT,
$ WORK, NWORK, RWORK, INFO )
*
-* -- LAPACK test routine (version 3.4.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2011
+* June 2016
*
* .. Scalar Arguments ..
INTEGER INFO, LDA, LDLRE, LDVL, LDVR, NIUNIT, NOUNIT,
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2015
+*> \date June 2016
*
*> \ingroup complex16_eig
*
* =====================================================================
SUBROUTINE ZERRED( PATH, NUNIT )
*
-* -- LAPACK test routine (version 3.6.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2015
+* June 2016
*
* .. Scalar Arguments ..
CHARACTER*3 PATH
CALL ZGESVDX( 'V', 'N', 'A', 2, 2, A, 2, ZERO, ZERO,
$ 0, 0, NS, S, U, 1, VT, 1, W, 1, RW, IW, INFO )
CALL CHKXER( 'ZGESVDX', INFOT, NOUT, LERR, OK )
- INFOT = 16
+ INFOT = 17
CALL ZGESVDX( 'N', 'V', 'A', 2, 2, A, 2, ZERO, ZERO,
$ 0, 0, NS, S, U, 1, VT, 1, W, 1, RW, IW, INFO )
CALL CHKXER( 'ZGESVDX', INFOT, NOUT, LERR, OK )
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2015
+*> \date June 2016
*
*> \ingroup complex16_eig
*
* =====================================================================
SUBROUTINE ZERRGG( PATH, NUNIT )
*
-* -- LAPACK test routine (version 3.6.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2015
+* June 2016
*
* .. Scalar Arguments ..
CHARACTER*3 PATH
* ..
* .. 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 ),
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
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 )
*> eigenvalue is negative), and otherwise it returns .FALSE..
*>
*> It is used by the test routine ZDRGES to test whether the driver
-*> routine ZGGES succesfully sorts eigenvalues.
+*> routine ZGGES successfully sorts eigenvalues.
*> \endverbatim
*
* Arguments:
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2011
+*> \date June 2016
*
*> \ingroup complex16_eig
*
* =====================================================================
LOGICAL FUNCTION ZLCTES( Z, D )
*
-* -- LAPACK test routine (version 3.4.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2011
+* June 2016
*
* .. Scalar Arguments ..
COMPLEX*16 D, Z
*>
*> ZSLECT returns .TRUE. if the eigenvalue Z is to be selected,
*> otherwise it returns .FALSE.
-*> It is used by ZCHK41 to test if ZGEES succesfully sorts eigenvalues,
-*> and by ZCHK43 to test if ZGEESX succesfully sorts eigenvalues.
+*> It is used by ZCHK41 to test if ZGEES successfully sorts eigenvalues,
+*> and by ZCHK43 to test if ZGEESX successfully sorts eigenvalues.
*>
*> The common block /SSLCT/ controls how eigenvalues are selected.
*> If SELOPT = 0, then ZSLECT return .TRUE. when real(Z) is less than
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2011
+*> \date June 2016
*
*> \ingroup complex16_eig
*
* =====================================================================
LOGICAL FUNCTION ZSLECT( Z )
*
-* -- LAPACK test routine (version 3.4.0) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2011
+* June 2016
*
* .. Scalar Arguments ..
COMPLEX*16 Z
.f.o:
$(FORTRAN) $(OPTS) -c $< -o $@
+
+.NOTPARALLEL:
SUBROUTINE CQRT04(M,N,NB,RESULT)
IMPLICIT NONE
*
-* -- LAPACK test routine (version 3.4.1) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
* April 2012
* ..
* .. Local allocatable arrays
COMPLEX, ALLOCATABLE :: AF(:,:), Q(:,:),
- $ R(:,:), RWORK(:), WORK( : ), T(:,:),
+ $ R(:,:), WORK( : ), T(:,:),
$ CF(:,:), DF(:,:), A(:,:), C(:,:), D(:,:)
+ REAL, ALLOCATABLE :: RWORK(:)
*
* .. Parameters ..
REAL ZERO
SUBROUTINE CQRT05(M,N,L,NB,RESULT)
IMPLICIT NONE
*
-* -- LAPACK test routine (version 3.4.1) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
* April 2012
* ..
* .. Local allocatable arrays
COMPLEX, ALLOCATABLE :: AF(:,:), Q(:,:),
- $ R(:,:), RWORK(:), WORK( : ), T(:,:),
+ $ R(:,:), WORK( : ), T(:,:),
$ CF(:,:), DF(:,:), A(:,:), C(:,:), D(:,:)
+ REAL, ALLOCATABLE :: RWORK(:)
*
* .. Parameters ..
REAL ZERO
SUBROUTINE ZQRT04(M,N,NB,RESULT)
IMPLICIT NONE
*
-* -- LAPACK test routine (version 3.4.1) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
* April 2012
* ..
* .. Local allocatable arrays
COMPLEX*16, ALLOCATABLE :: AF(:,:), Q(:,:),
- $ R(:,:), RWORK(:), WORK( : ), T(:,:),
+ $ R(:,:), WORK( : ), T(:,:),
$ CF(:,:), DF(:,:), A(:,:), C(:,:), D(:,:)
+ DOUBLE PRECISION, ALLOCATABLE :: RWORK(:)
*
* .. Parameters ..
DOUBLE PRECISION ZERO
SUBROUTINE ZQRT05(M,N,L,NB,RESULT)
IMPLICIT NONE
*
-* -- LAPACK test routine (version 3.4.1) --
+* -- LAPACK test routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
* April 2012
* ..
* .. Local allocatable arrays
COMPLEX*16, ALLOCATABLE :: AF(:,:), Q(:,:),
- $ R(:,:), RWORK(:), WORK( : ), T(:,:),
+ $ R(:,:), WORK( : ), T(:,:),
$ CF(:,:), DF(:,:), A(:,:), C(:,:), D(:,:)
+ DOUBLE PRECISION, ALLOCATABLE :: RWORK(:)
*
* .. Parameters ..
DOUBLE PRECISION ZERO
*> \verbatim
*>
*> CLATM2 returns the (I,J) entry of a random matrix of dimension
-*> (M, N) described by the other paramters. It is called by the
+*> (M, N) described by the other parameters. It is called by the
*> CLATMR routine in order to build random test matrices. No error
*> checking on parameters is done, because this routine is called in
*> a tight loop by CLATMR which has already checked the parameters.
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2011
+*> \date June 2016
*
*> \ingroup complex_matgen
*
COMPLEX FUNCTION CLATM2( M, N, I, J, KL, KU, IDIST, ISEED, D,
$ IGRADE, DL, DR, IPVTNG, IWORK, SPARSE )
*
-* -- LAPACK auxiliary routine (version 3.4.0) --
+* -- LAPACK auxiliary routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2011
+* June 2016
*
* .. Scalar Arguments ..
*
*> \verbatim
*>
*> CLATM3 returns the (ISUB,JSUB) entry of a random matrix of
-*> dimension (M, N) described by the other paramters. (ISUB,JSUB)
+*> dimension (M, N) described by the other parameters. (ISUB,JSUB)
*> is the final position of the (I,J) entry after pivoting
*> according to IPVTNG and IWORK. CLATM3 is called by the
*> CLATMR routine in order to build random test matrices. No error
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2011
+*> \date June 2016
*
*> \ingroup complex_matgen
*
$ ISEED, D, IGRADE, DL, DR, IPVTNG, IWORK,
$ SPARSE )
*
-* -- LAPACK auxiliary routine (version 3.4.0) --
+* -- LAPACK auxiliary routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2011
+* June 2016
*
* .. Scalar Arguments ..
*
*> \param[in] PRTYPE
*> \verbatim
*> PRTYPE is INTEGER
-*> "Points" to a certian type of the matrices to generate
+*> "Points" to a certain type of the matrices to generate
*> (see futher details).
*> \endverbatim
*>
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2011
+*> \date June 2016
*
*> \ingroup complex_matgen
*
$ E, LDE, F, LDF, R, LDR, L, LDL, ALPHA, QBLCKA,
$ QBLCKB )
*
-* -- LAPACK computational routine (version 3.4.0) --
+* -- LAPACK computational routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2011
+* June 2016
*
* .. Scalar Arguments ..
INTEGER LDA, LDB, LDC, LDD, LDE, LDF, LDL, LDR, M, N,
*> \verbatim
*>
*> DLATM2 returns the (I,J) entry of a random matrix of dimension
-*> (M, N) described by the other paramters. It is called by the
+*> (M, N) described by the other parameters. It is called by the
*> DLATMR routine in order to build random test matrices. No error
*> checking on parameters is done, because this routine is called in
*> a tight loop by DLATMR which has already checked the parameters.
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2011
+*> \date June 2016
*
*> \ingroup double_matgen
*
DOUBLE PRECISION FUNCTION DLATM2( M, N, I, J, KL, KU, IDIST,
$ ISEED, D, IGRADE, DL, DR, IPVTNG, IWORK, SPARSE )
*
-* -- LAPACK auxiliary routine (version 3.4.0) --
+* -- LAPACK auxiliary routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2011
+* June 2016
*
* .. Scalar Arguments ..
*
*> \verbatim
*>
*> DLATM3 returns the (ISUB,JSUB) entry of a random matrix of
-*> dimension (M, N) described by the other paramters. (ISUB,JSUB)
+*> dimension (M, N) described by the other parameters. (ISUB,JSUB)
*> is the final position of the (I,J) entry after pivoting
*> according to IPVTNG and IWORK. DLATM3 is called by the
*> DLATMR routine in order to build random test matrices. No error
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2011
+*> \date June 2016
*
*> \ingroup double_matgen
*
$ IDIST, ISEED, D, IGRADE, DL, DR, IPVTNG, IWORK,
$ SPARSE )
*
-* -- LAPACK auxiliary routine (version 3.4.0) --
+* -- LAPACK auxiliary routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2011
+* June 2016
*
* .. Scalar Arguments ..
*
*> \param[in] PRTYPE
*> \verbatim
*> PRTYPE is INTEGER
-*> "Points" to a certian type of the matrices to generate
+*> "Points" to a certain type of the matrices to generate
*> (see futher details).
*> \endverbatim
*>
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2011
+*> \date June 2016
*
*> \ingroup double_matgen
*
$ E, LDE, F, LDF, R, LDR, L, LDL, ALPHA, QBLCKA,
$ QBLCKB )
*
-* -- LAPACK computational routine (version 3.4.0) --
+* -- LAPACK computational routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2011
+* June 2016
*
* .. Scalar Arguments ..
INTEGER LDA, LDB, LDC, LDD, LDE, LDF, LDL, LDR, M, N,
*> \verbatim
*>
*> SLATM2 returns the (I,J) entry of a random matrix of dimension
-*> (M, N) described by the other paramters. It is called by the
+*> (M, N) described by the other parameters. It is called by the
*> SLATMR routine in order to build random test matrices. No error
*> checking on parameters is done, because this routine is called in
*> a tight loop by SLATMR which has already checked the parameters.
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2011
+*> \date June 2016
*
*> \ingroup real_matgen
*
REAL FUNCTION SLATM2( M, N, I, J, KL, KU, IDIST,
$ ISEED, D, IGRADE, DL, DR, IPVTNG, IWORK, SPARSE )
*
-* -- LAPACK auxiliary routine (version 3.4.0) --
+* -- LAPACK auxiliary routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2011
+* June 2016
*
* .. Scalar Arguments ..
*
*> \verbatim
*>
*> SLATM3 returns the (ISUB,JSUB) entry of a random matrix of
-*> dimension (M, N) described by the other paramters. (ISUB,JSUB)
+*> dimension (M, N) described by the other parameters. (ISUB,JSUB)
*> is the final position of the (I,J) entry after pivoting
*> according to IPVTNG and IWORK. SLATM3 is called by the
*> SLATMR routine in order to build random test matrices. No error
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2011
+*> \date June 2016
*
*> \ingroup real_matgen
*
$ IDIST, ISEED, D, IGRADE, DL, DR, IPVTNG, IWORK,
$ SPARSE )
*
-* -- LAPACK auxiliary routine (version 3.4.0) --
+* -- LAPACK auxiliary routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2011
+* June 2016
*
* .. Scalar Arguments ..
*
*> \param[in] PRTYPE
*> \verbatim
*> PRTYPE is INTEGER
-*> "Points" to a certian type of the matrices to generate
+*> "Points" to a certain type of the matrices to generate
*> (see futher details).
*> \endverbatim
*>
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2011
+*> \date June 2016
*
*> \ingroup real_matgen
*
$ E, LDE, F, LDF, R, LDR, L, LDL, ALPHA, QBLCKA,
$ QBLCKB )
*
-* -- LAPACK computational routine (version 3.4.0) --
+* -- LAPACK computational routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2011
+* June 2016
*
* .. Scalar Arguments ..
INTEGER LDA, LDB, LDC, LDD, LDE, LDF, LDL, LDR, M, N,
*> \verbatim
*>
*> ZLATM2 returns the (I,J) entry of a random matrix of dimension
-*> (M, N) described by the other paramters. It is called by the
+*> (M, N) described by the other parameters. It is called by the
*> ZLATMR routine in order to build random test matrices. No error
*> checking on parameters is done, because this routine is called in
*> a tight loop by ZLATMR which has already checked the parameters.
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2011
+*> \date June 2016
*
*> \ingroup complex16_matgen
*
COMPLEX*16 FUNCTION ZLATM2( M, N, I, J, KL, KU, IDIST,
$ ISEED, D, IGRADE, DL, DR, IPVTNG, IWORK, SPARSE )
*
-* -- LAPACK auxiliary routine (version 3.4.0) --
+* -- LAPACK auxiliary routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2011
+* June 2016
*
* .. Scalar Arguments ..
*
*> \verbatim
*>
*> ZLATM3 returns the (ISUB,JSUB) entry of a random matrix of
-*> dimension (M, N) described by the other paramters. (ISUB,JSUB)
+*> dimension (M, N) described by the other parameters. (ISUB,JSUB)
*> is the final position of the (I,J) entry after pivoting
*> according to IPVTNG and IWORK. ZLATM3 is called by the
*> ZLATMR routine in order to build random test matrices. No error
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2011
+*> \date June 2016
*
*> \ingroup complex16_matgen
*
$ IDIST, ISEED, D, IGRADE, DL, DR, IPVTNG, IWORK,
$ SPARSE )
*
-* -- LAPACK auxiliary routine (version 3.4.0) --
+* -- LAPACK auxiliary routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2011
+* June 2016
*
* .. Scalar Arguments ..
*
*> \param[in] PRTYPE
*> \verbatim
*> PRTYPE is INTEGER
-*> "Points" to a certian type of the matrices to generate
+*> "Points" to a certain type of the matrices to generate
*> (see futher details).
*> \endverbatim
*>
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
-*> \date November 2011
+*> \date June 2016
*
*> \ingroup complex16_matgen
*
$ E, LDE, F, LDF, R, LDR, L, LDL, ALPHA, QBLCKA,
$ QBLCKB )
*
-* -- LAPACK computational routine (version 3.4.0) --
+* -- LAPACK computational routine (version 3.6.1) --
* -- LAPACK is a software package provided by Univ. of Tennessee, --
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-* November 2011
+* June 2016
*
* .. Scalar Arguments ..
INTEGER LDA, LDB, LDC, LDD, LDE, LDF, LDL, LDR, M, N,
FRC:
@FRC=$(FRC)
+
+.NOTPARALLEL:
1 2 15 Values of NRHS (number of right hand sides)
9 Number of matrix types (list types on next line if 0 < NTYPES < 9)
1 2 3 4 5 6 7 8 9 Matrix Types
-60.0 Threshold value of test ratio
+30.0 Threshold value of test ratio
T Put T to test the error exits
0 5 7 3 200 Values of INIBL (nibble crossover point)
1 2 4 2 1 Values of ISHFTS (number of simultaneous shifts)
0 1 2 0 1 Values of IACC22 (select structured matrix multiply: 0, 1 or 2)
-40.0 Threshold value
+20.0 Threshold value
T Put T to test the error exits
1 Code to interpret the seed
NEP 21
SEP: Data file for testing Symmetric Eigenvalue Problem routines
-8 Number of values of N
-0 1 2 3 5 16 19 20 Values of N (dimension)
+6 Number of values of N
+0 1 2 3 5 20 Values of N (dimension)
5 Number of values of NB
1 3 3 3 10 Values of NB (blocksize)
2 2 2 2 2 Values of NBMIN (minimum blocksize)
1 0 5 9 1 Values of NX (crossover point)
-600.0 Threshold value
+50.0 Threshold value
T Put T to test the LAPACK routines
T Put T to test the driver routines
T Put T to test the error exits
2 2 2 2 2 Values of NBMIN (minimum blocksize)
1 0 5 9 1 Values of NX (crossover point)
2 0 2 2 2 Values of NRHS
-100.0 Threshold value
+50.0 Threshold value
T Put T to test the LAPACK routines
T Put T to test the driver routines
T Put T to test the error exits