Import LAPACK: TESTING directory
This commit is contained in:
parent
92a858e69e
commit
57eee3fa43
|
@ -405,7 +405,7 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2011
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup complex_eig
|
*> \ingroup complex_eig
|
||||||
*
|
*
|
||||||
|
@ -415,10 +415,10 @@
|
||||||
$ Y, Z, Q, LDQ, PT, LDPT, U, VT, WORK, LWORK,
|
$ Y, Z, Q, LDQ, PT, LDPT, U, VT, WORK, LWORK,
|
||||||
$ RWORK, NOUT, INFO )
|
$ 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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2011
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
INTEGER INFO, LDA, LDPT, LDQ, LDX, LWORK, NOUT, NRHS,
|
INTEGER INFO, LDA, LDPT, LDQ, LDX, LWORK, NOUT, NRHS,
|
||||||
|
@ -466,9 +466,10 @@
|
||||||
EXTERNAL SLAMCH, SLARND
|
EXTERNAL SLAMCH, SLARND
|
||||||
* ..
|
* ..
|
||||||
* .. External Subroutines ..
|
* .. External Subroutines ..
|
||||||
EXTERNAL ALASUM, CBDSQR, CBDT01, CBDT02, CBDT03, CGEBRD,
|
EXTERNAL ALASUM, CBDSQR, CBDT01, CBDT02, CBDT03,
|
||||||
$ CGEMM, CLACPY, CLASET, CLATMR, CLATMS, CUNGBR,
|
$ CGEBRD, CGEMM, CLACPY, CLASET, CLATMR,
|
||||||
$ CUNT01, SCOPY, SLABAD, SLAHD2, SSVDCH, XERBLA
|
$ CLATMS, CUNGBR, CUNT01, SCOPY, SLABAD,
|
||||||
|
$ SLAHD2, SSVDCH, XERBLA
|
||||||
* ..
|
* ..
|
||||||
* .. Intrinsic Functions ..
|
* .. Intrinsic Functions ..
|
||||||
INTRINSIC ABS, EXP, INT, LOG, MAX, MIN, SQRT
|
INTRINSIC ABS, EXP, INT, LOG, MAX, MIN, SQRT
|
||||||
|
|
|
@ -1027,17 +1027,17 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2015
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup complex_eig
|
*> \ingroup complex_eig
|
||||||
*
|
*
|
||||||
* =====================================================================
|
* =====================================================================
|
||||||
PROGRAM CCHKEE
|
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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2015
|
* June 2016
|
||||||
*
|
*
|
||||||
* =====================================================================
|
* =====================================================================
|
||||||
*
|
*
|
||||||
|
@ -2098,6 +2098,7 @@
|
||||||
MAXTYP = 26
|
MAXTYP = 26
|
||||||
NTYPES = MIN( MAXTYP, NTYPES )
|
NTYPES = MIN( MAXTYP, NTYPES )
|
||||||
CALL ALAREQ( C3, NTYPES, DOTYPE, MAXTYP, NIN, NOUT )
|
CALL ALAREQ( C3, NTYPES, DOTYPE, MAXTYP, NIN, NOUT )
|
||||||
|
CALL XLAENV(1,1)
|
||||||
IF( TSTCHK .AND. TSTERR )
|
IF( TSTCHK .AND. TSTERR )
|
||||||
$ CALL CERRGG( C3, NOUT )
|
$ CALL CERRGG( C3, NOUT )
|
||||||
DO 350 I = 1, NPARMS
|
DO 350 I = 1, NPARMS
|
||||||
|
@ -2157,6 +2158,7 @@
|
||||||
*
|
*
|
||||||
* Blocked version
|
* Blocked version
|
||||||
*
|
*
|
||||||
|
CALL XLAENV(16,2)
|
||||||
CALL CDRGES3( NN, NVAL, MAXTYP, DOTYPE, ISEED, THRESH, NOUT,
|
CALL CDRGES3( NN, NVAL, MAXTYP, DOTYPE, ISEED, THRESH, NOUT,
|
||||||
$ A( 1, 1 ), NMAX, A( 1, 2 ), A( 1, 3 ),
|
$ A( 1, 1 ), NMAX, A( 1, 2 ), A( 1, 3 ),
|
||||||
$ A( 1, 4 ), A( 1, 7 ), NMAX, A( 1, 8 ),
|
$ A( 1, 4 ), A( 1, 7 ), NMAX, A( 1, 8 ),
|
||||||
|
@ -2223,6 +2225,7 @@
|
||||||
*
|
*
|
||||||
* Blocked version
|
* Blocked version
|
||||||
*
|
*
|
||||||
|
CALL XLAENV(16,2)
|
||||||
CALL CDRGEV3( NN, NVAL, MAXTYP, DOTYPE, ISEED, THRESH, NOUT,
|
CALL CDRGEV3( NN, NVAL, MAXTYP, DOTYPE, ISEED, THRESH, NOUT,
|
||||||
$ A( 1, 1 ), NMAX, A( 1, 2 ), A( 1, 3 ),
|
$ A( 1, 1 ), NMAX, A( 1, 2 ), A( 1, 3 ),
|
||||||
$ A( 1, 4 ), A( 1, 7 ), NMAX, A( 1, 8 ),
|
$ A( 1, 4 ), A( 1, 7 ), NMAX, A( 1, 8 ),
|
||||||
|
@ -2348,6 +2351,7 @@
|
||||||
* GSV: Generalized Singular Value Decomposition
|
* GSV: Generalized Singular Value Decomposition
|
||||||
* ----------------------------------------------
|
* ----------------------------------------------
|
||||||
*
|
*
|
||||||
|
CALL XLAENV(1,1)
|
||||||
IF( TSTERR )
|
IF( TSTERR )
|
||||||
$ CALL CERRGG( 'GSV', NOUT )
|
$ CALL CERRGG( 'GSV', NOUT )
|
||||||
CALL CCKGSV( NN, MVAL, PVAL, NVAL, NTYPES, ISEED, THRESH, NMAX,
|
CALL CCKGSV( NN, MVAL, PVAL, NVAL, NTYPES, ISEED, THRESH, NMAX,
|
||||||
|
|
|
@ -68,7 +68,7 @@
|
||||||
*> and each type of matrix, one matrix will be generated and used
|
*> and each type of matrix, one matrix will be generated and used
|
||||||
*> to test the nonsymmetric eigenroutines. For each matrix, 13
|
*> to test the nonsymmetric eigenroutines. For each matrix, 13
|
||||||
*> tests will be performed. The first twelve "test ratios" should be
|
*> 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
|
*> H
|
||||||
*> (1) | A - U H V | / ( |A| n ulp )
|
*> (1) | A - U H V | / ( |A| n ulp )
|
||||||
|
@ -128,7 +128,7 @@
|
||||||
*> |beta(Q,Z computed) - beta(Q,Z not computed)|/|P| ) / ulp
|
*> |beta(Q,Z computed) - beta(Q,Z not computed)|/|P| ) / ulp
|
||||||
*>
|
*>
|
||||||
*> In addition, the normalization of L and R are checked, and compared
|
*> In addition, the normalization of L and R are checked, and compared
|
||||||
*> with the threshhold THRSHN.
|
*> with the threshold THRSHN.
|
||||||
*>
|
*>
|
||||||
*> Test Matrices
|
*> Test Matrices
|
||||||
*> ---- --------
|
*> ---- --------
|
||||||
|
@ -298,7 +298,7 @@
|
||||||
*> \param[in] THRSHN
|
*> \param[in] THRSHN
|
||||||
*> \verbatim
|
*> \verbatim
|
||||||
*> THRSHN is REAL
|
*> 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
|
*> If the normalization of any eigenvector differs from 1 by
|
||||||
*> more than THRSHN*ulp, then a special error message will be
|
*> more than THRSHN*ulp, then a special error message will be
|
||||||
*> printed. (This is handled separately from the other tests,
|
*> printed. (This is handled separately from the other tests,
|
||||||
|
@ -492,7 +492,7 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2011
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup complex_eig
|
*> \ingroup complex_eig
|
||||||
*
|
*
|
||||||
|
@ -503,10 +503,10 @@
|
||||||
$ ALPHA3, BETA3, EVECTL, EVECTR, WORK, LWORK,
|
$ ALPHA3, BETA3, EVECTL, EVECTR, WORK, LWORK,
|
||||||
$ RWORK, LLWORK, RESULT, INFO )
|
$ 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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2011
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
LOGICAL TSTDIF
|
LOGICAL TSTDIF
|
||||||
|
|
|
@ -50,7 +50,7 @@
|
||||||
*> number of matrix "TYPES" are specified. For each size ("N")
|
*> 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 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
|
*> 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).
|
*> the tests (5), (11) and (13).
|
||||||
*>
|
*>
|
||||||
*>
|
*>
|
||||||
|
@ -372,7 +372,7 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2011
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup complex_eig
|
*> \ingroup complex_eig
|
||||||
*
|
*
|
||||||
|
@ -381,10 +381,10 @@
|
||||||
$ NOUNIT, A, LDA, B, S, T, Q, LDQ, Z, ALPHA,
|
$ NOUNIT, A, LDA, B, S, T, Q, LDQ, Z, ALPHA,
|
||||||
$ BETA, WORK, LWORK, RWORK, RESULT, BWORK, INFO )
|
$ 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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2011
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
INTEGER INFO, LDA, LDQ, LWORK, NOUNIT, NSIZES, NTYPES
|
INTEGER INFO, LDA, LDQ, LWORK, NOUNIT, NSIZES, NTYPES
|
||||||
|
|
|
@ -50,7 +50,7 @@
|
||||||
*> number of matrix "TYPES" are specified. For each size ("N")
|
*> 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 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
|
*> 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).
|
*> the tests (5), (11) and (13).
|
||||||
*>
|
*>
|
||||||
*>
|
*>
|
||||||
|
@ -382,7 +382,7 @@
|
||||||
$ BETA, WORK, LWORK, RWORK, RESULT, BWORK,
|
$ BETA, WORK, LWORK, RWORK, RESULT, BWORK,
|
||||||
$ INFO )
|
$ 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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* February 2015
|
* February 2015
|
||||||
|
|
|
@ -55,7 +55,7 @@
|
||||||
*> number of matrix "types" are specified. For each size ("n")
|
*> 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 each type of matrix, a pair of matrices (A, B) will be generated
|
||||||
*> and used for testing. For each matrix pair, the following tests
|
*> 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:
|
*> Results from CGGEV:
|
||||||
*>
|
*>
|
||||||
|
@ -389,7 +389,7 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2015
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup complex_eig
|
*> \ingroup complex_eig
|
||||||
*
|
*
|
||||||
|
@ -399,10 +399,10 @@
|
||||||
$ ALPHA, BETA, ALPHA1, BETA1, WORK, LWORK, RWORK,
|
$ ALPHA, BETA, ALPHA1, BETA1, WORK, LWORK, RWORK,
|
||||||
$ RESULT, INFO )
|
$ 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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2015
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
INTEGER INFO, LDA, LDQ, LDQE, LWORK, NOUNIT, NSIZES,
|
INTEGER INFO, LDA, LDQ, LDQE, LWORK, NOUNIT, NSIZES,
|
||||||
|
|
|
@ -55,7 +55,7 @@
|
||||||
*> number of matrix "types" are specified. For each size ("n")
|
*> 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 each type of matrix, a pair of matrices (A, B) will be generated
|
||||||
*> and used for testing. For each matrix pair, the following tests
|
*> 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:
|
*> Results from CGGEV3:
|
||||||
*>
|
*>
|
||||||
|
@ -399,7 +399,7 @@
|
||||||
$ ALPHA, BETA, ALPHA1, BETA1, WORK, LWORK,
|
$ ALPHA, BETA, ALPHA1, BETA1, WORK, LWORK,
|
||||||
$ RWORK, RESULT, INFO )
|
$ 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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* January 2015
|
* January 2015
|
||||||
|
|
|
@ -59,7 +59,7 @@
|
||||||
*> (need more details on what kind of read-in data are needed).
|
*> (need more details on what kind of read-in data are needed).
|
||||||
*>
|
*>
|
||||||
*> For each matrix pair, the following tests will be performed and
|
*> 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 )
|
*> (1) | A - Q S Z' | / ( |A| n ulp )
|
||||||
*>
|
*>
|
||||||
|
@ -340,7 +340,7 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2011
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup complex_eig
|
*> \ingroup complex_eig
|
||||||
*
|
*
|
||||||
|
@ -349,10 +349,10 @@
|
||||||
$ AI, BI, Z, Q, ALPHA, BETA, C, LDC, S, WORK,
|
$ AI, BI, Z, Q, ALPHA, BETA, C, LDC, S, WORK,
|
||||||
$ LWORK, RWORK, IWORK, LIWORK, BWORK, INFO )
|
$ 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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2011
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
INTEGER INFO, LDA, LDC, LIWORK, LWORK, NCMAX, NIN,
|
INTEGER INFO, LDA, LDC, LIWORK, LWORK, NCMAX, NIN,
|
||||||
|
|
|
@ -50,7 +50,7 @@
|
||||||
*> corresponding the first and last eigenvalues are also know
|
*> corresponding the first and last eigenvalues are also know
|
||||||
*> ``exactly'' (see CLATM6).
|
*> ``exactly'' (see CLATM6).
|
||||||
*> For each matrix pair, the following tests will be performed and
|
*> 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
|
*> (1) max over all left eigenvalue/-vector pairs (beta/alpha,l) of
|
||||||
*>
|
*>
|
||||||
|
@ -288,7 +288,7 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2011
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup complex_eig
|
*> \ingroup complex_eig
|
||||||
*
|
*
|
||||||
|
@ -298,10 +298,10 @@
|
||||||
$ S, STRU, DIF, DIFTRU, WORK, LWORK, RWORK,
|
$ S, STRU, DIF, DIFTRU, WORK, LWORK, RWORK,
|
||||||
$ IWORK, LIWORK, RESULT, BWORK, INFO )
|
$ 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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2011
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
INTEGER IHI, ILO, INFO, LDA, LIWORK, LWORK, NIN, NOUT,
|
INTEGER IHI, ILO, INFO, LDA, LIWORK, LWORK, NIN, NOUT,
|
||||||
|
|
|
@ -379,7 +379,7 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2015
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup complex_eig
|
*> \ingroup complex_eig
|
||||||
*
|
*
|
||||||
|
@ -389,10 +389,10 @@
|
||||||
$ SSAV, E, WORK, LWORK, RWORK, IWORK, NOUNIT,
|
$ SSAV, E, WORK, LWORK, RWORK, IWORK, NOUNIT,
|
||||||
$ INFO )
|
$ 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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2015
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
INTEGER INFO, LDA, LDU, LDVT, LWORK, NOUNIT, NSIZES,
|
INTEGER INFO, LDA, LDU, LDVT, LWORK, NOUNIT, NSIZES,
|
||||||
|
@ -423,11 +423,11 @@
|
||||||
* .. Local Scalars ..
|
* .. Local Scalars ..
|
||||||
LOGICAL BADMM, BADNN
|
LOGICAL BADMM, BADNN
|
||||||
CHARACTER JOBQ, JOBU, JOBVT, RANGE
|
CHARACTER JOBQ, JOBU, JOBVT, RANGE
|
||||||
INTEGER I, IINFO, IJQ, IJU, IJVT, IL, IU, ITEMP, IWSPC,
|
INTEGER I, IINFO, IJQ, IJU, IJVT, IL, IU, ITEMP,
|
||||||
$ IWTMP, J, JSIZE, JTYPE, LSWORK, M, MINWRK,
|
$ IWSPC, IWTMP, J, JSIZE, JTYPE, LSWORK, M,
|
||||||
$ MMAX, MNMAX, MNMIN, MTYPES, N, NERRS, NFAIL,
|
$ MINWRK, MMAX, MNMAX, MNMIN, MTYPES, N,
|
||||||
$ NMAX, NS, NSI, NSV, NTEST, NTESTF, NTESTT,
|
$ NERRS, NFAIL, NMAX, NS, NSI, NSV, NTEST,
|
||||||
$ LRWORK
|
$ NTESTF, NTESTT, LRWORK
|
||||||
REAL ANORM, DIF, DIV, OVFL, RTUNFL, ULP, ULPINV,
|
REAL ANORM, DIF, DIV, OVFL, RTUNFL, ULP, ULPINV,
|
||||||
$ UNFL, VL, VU
|
$ UNFL, VL, VU
|
||||||
* ..
|
* ..
|
||||||
|
@ -441,9 +441,9 @@
|
||||||
EXTERNAL SLAMCH, SLARND
|
EXTERNAL SLAMCH, SLARND
|
||||||
* ..
|
* ..
|
||||||
* .. External Subroutines ..
|
* .. External Subroutines ..
|
||||||
EXTERNAL ALASVM, XERBLA, CBDT01, CBDT05, CGESDD, CGESVD,
|
EXTERNAL ALASVM, XERBLA, CBDT01, CBDT05, CGESDD,
|
||||||
$ CGESVJ, CGEJSV, CGESVDX, CLACPY, CLASET, CLATMS,
|
$ CGESVD, CGESVJ, CGEJSV, CGESVDX, CLACPY,
|
||||||
$ CUNT01, CUNT03
|
$ CLASET, CLATMS, CUNT01, CUNT03
|
||||||
* ..
|
* ..
|
||||||
* .. Intrinsic Functions ..
|
* .. Intrinsic Functions ..
|
||||||
INTRINSIC ABS, REAL, MAX, MIN
|
INTRINSIC ABS, REAL, MAX, MIN
|
||||||
|
|
|
@ -369,7 +369,7 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2011
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup complex_eig
|
*> \ingroup complex_eig
|
||||||
*
|
*
|
||||||
|
@ -378,10 +378,10 @@
|
||||||
$ NOUNIT, A, LDA, H, HT, W, WT, VS, LDVS, RESULT,
|
$ NOUNIT, A, LDA, H, HT, W, WT, VS, LDVS, RESULT,
|
||||||
$ WORK, NWORK, RWORK, IWORK, BWORK, INFO )
|
$ 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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2011
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
INTEGER INFO, LDA, LDVS, NOUNIT, NSIZES, NTYPES, NWORK
|
INTEGER INFO, LDA, LDVS, NOUNIT, NSIZES, NTYPES, NWORK
|
||||||
|
@ -901,7 +901,7 @@
|
||||||
$ ' 1/ulp otherwise', /
|
$ ' 1/ulp otherwise', /
|
||||||
$ ' 12 = 0 if W same no matter if VS computed (sort),',
|
$ ' 12 = 0 if W same no matter if VS computed (sort),',
|
||||||
$ ' 1/ulp otherwise', /
|
$ ' 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, ',' ),
|
9993 FORMAT( ' N=', I5, ', IWK=', I2, ', seed=', 4( I4, ',' ),
|
||||||
$ ' type ', I2, ', test(', I2, ')=', G10.3 )
|
$ ' type ', I2, ', test(', I2, ')=', G10.3 )
|
||||||
9992 FORMAT( ' CDRVES: ', A, ' returned INFO=', I6, '.', / 9X, 'N=',
|
9992 FORMAT( ' CDRVES: ', A, ' returned INFO=', I6, '.', / 9X, 'N=',
|
||||||
|
|
|
@ -425,7 +425,7 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2011
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup complex_eig
|
*> \ingroup complex_eig
|
||||||
*
|
*
|
||||||
|
@ -435,10 +435,10 @@
|
||||||
$ LDVS, VS1, RESULT, WORK, LWORK, RWORK, BWORK,
|
$ LDVS, VS1, RESULT, WORK, LWORK, RWORK, BWORK,
|
||||||
$ INFO )
|
$ 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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2011
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
INTEGER INFO, LDA, LDVS, LWORK, NIUNIT, NOUNIT, NSIZES,
|
INTEGER INFO, LDA, LDVS, LWORK, NIUNIT, NOUNIT, NSIZES,
|
||||||
|
@ -914,7 +914,7 @@
|
||||||
$ ' 1/ulp otherwise', /
|
$ ' 1/ulp otherwise', /
|
||||||
$ ' 12 = 0 if W same no matter what else computed ',
|
$ ' 12 = 0 if W same no matter what else computed ',
|
||||||
$ '(sort), 1/ulp otherwise', /
|
$ '(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,',
|
$ / ' 14 = 0 if RCONDE same no matter what else computed,',
|
||||||
$ ' 1/ulp otherwise', /
|
$ ' 1/ulp otherwise', /
|
||||||
$ ' 15 = 0 if RCONDv same no matter what else computed,',
|
$ ' 15 = 0 if RCONDv same no matter what else computed,',
|
||||||
|
|
|
@ -446,7 +446,7 @@
|
||||||
*> \verbatim
|
*> \verbatim
|
||||||
*> INFO is INTEGER
|
*> INFO is INTEGER
|
||||||
*> If 0, then successful exit.
|
*> 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
|
*> If >0, CLATMR, CLATMS, CLATME or CGET23 returned an error
|
||||||
*> code, and INFO is its absolute value.
|
*> code, and INFO is its absolute value.
|
||||||
*>
|
*>
|
||||||
|
@ -485,7 +485,7 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2011
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup complex_eig
|
*> \ingroup complex_eig
|
||||||
*
|
*
|
||||||
|
@ -496,10 +496,10 @@
|
||||||
$ RCONDE, RCNDE1, RCDEIN, SCALE, SCALE1, RESULT,
|
$ RCONDE, RCNDE1, RCDEIN, SCALE, SCALE1, RESULT,
|
||||||
$ WORK, NWORK, RWORK, INFO )
|
$ 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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2011
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
INTEGER INFO, LDA, LDLRE, LDVL, LDVR, NIUNIT, NOUNIT,
|
INTEGER INFO, LDA, LDLRE, LDVL, LDVR, NIUNIT, NOUNIT,
|
||||||
|
|
|
@ -61,17 +61,17 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2015
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup complex_eig
|
*> \ingroup complex_eig
|
||||||
*
|
*
|
||||||
* =====================================================================
|
* =====================================================================
|
||||||
SUBROUTINE CERRED( PATH, NUNIT )
|
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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2015
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
CHARACTER*3 PATH
|
CHARACTER*3 PATH
|
||||||
|
@ -100,7 +100,7 @@
|
||||||
$ VT( NMAX, NMAX ), W( 10*NMAX ), X( NMAX )
|
$ VT( NMAX, NMAX ), W( 10*NMAX ), X( NMAX )
|
||||||
* ..
|
* ..
|
||||||
* .. External Subroutines ..
|
* .. External Subroutines ..
|
||||||
EXTERNAL CHKXER, CGEES, CGEESX, CGEEV, CGEEVX, CGEJSV
|
EXTERNAL CHKXER, CGEES, CGEESX, CGEEV, CGEEVX, CGEJSV,
|
||||||
$ CGESDD, CGESVD
|
$ CGESDD, CGESVD
|
||||||
* ..
|
* ..
|
||||||
* .. External Functions ..
|
* .. External Functions ..
|
||||||
|
@ -484,7 +484,7 @@
|
||||||
CALL CGESVDX( 'V', 'N', 'A', 2, 2, A, 2, ZERO, ZERO,
|
CALL CGESVDX( 'V', 'N', 'A', 2, 2, A, 2, ZERO, ZERO,
|
||||||
$ 0, 0, NS, S, U, 1, VT, 1, W, 1, RW, IW, INFO )
|
$ 0, 0, NS, S, U, 1, VT, 1, W, 1, RW, IW, INFO )
|
||||||
CALL CHKXER( 'CGESVDX', INFOT, NOUT, LERR, OK )
|
CALL CHKXER( 'CGESVDX', INFOT, NOUT, LERR, OK )
|
||||||
INFOT = 16
|
INFOT = 17
|
||||||
CALL CGESVDX( 'N', 'V', 'A', 2, 2, A, 2, ZERO, ZERO,
|
CALL CGESVDX( 'N', 'V', 'A', 2, 2, A, 2, ZERO, ZERO,
|
||||||
$ 0, 0, NS, S, U, 1, VT, 1, W, 1, RW, IW, INFO )
|
$ 0, 0, NS, S, U, 1, VT, 1, W, 1, RW, IW, INFO )
|
||||||
CALL CHKXER( 'CGESVDX', INFOT, NOUT, LERR, OK )
|
CALL CHKXER( 'CGESVDX', INFOT, NOUT, LERR, OK )
|
||||||
|
|
|
@ -50,17 +50,17 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2015
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup complex_eig
|
*> \ingroup complex_eig
|
||||||
*
|
*
|
||||||
* =====================================================================
|
* =====================================================================
|
||||||
SUBROUTINE CERRGG( PATH, NUNIT )
|
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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2015
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
CHARACTER*3 PATH
|
CHARACTER*3 PATH
|
||||||
|
@ -83,7 +83,7 @@
|
||||||
* ..
|
* ..
|
||||||
* .. Local Arrays ..
|
* .. Local Arrays ..
|
||||||
LOGICAL BW( NMAX ), SEL( NMAX )
|
LOGICAL BW( NMAX ), SEL( NMAX )
|
||||||
INTEGER IW( LW )
|
INTEGER IW( LW ), IDUM(NMAX)
|
||||||
REAL LS( NMAX ), R1( NMAX ), R2( NMAX ),
|
REAL LS( NMAX ), R1( NMAX ), R2( NMAX ),
|
||||||
$ RCE( NMAX ), RCV( NMAX ), RS( NMAX ), RW( LW )
|
$ RCE( NMAX ), RCV( NMAX ), RS( NMAX ), RW( LW )
|
||||||
COMPLEX A( NMAX, NMAX ), ALPHA( NMAX ),
|
COMPLEX A( NMAX, NMAX ), ALPHA( NMAX ),
|
||||||
|
@ -306,57 +306,57 @@
|
||||||
SRNAMT = 'CGGSVD3'
|
SRNAMT = 'CGGSVD3'
|
||||||
INFOT = 1
|
INFOT = 1
|
||||||
CALL CGGSVD3( '/', 'N', 'N', 0, 0, 0, DUMMYK, DUMMYL, A, 1, B,
|
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 )
|
$ INFO )
|
||||||
CALL CHKXER( 'CGGSVD3', INFOT, NOUT, LERR, OK )
|
CALL CHKXER( 'CGGSVD3', INFOT, NOUT, LERR, OK )
|
||||||
INFOT = 2
|
INFOT = 2
|
||||||
CALL CGGSVD3( 'N', '/', 'N', 0, 0, 0, DUMMYK, DUMMYL, A, 1, B,
|
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 )
|
$ INFO )
|
||||||
CALL CHKXER( 'CGGSVD3', INFOT, NOUT, LERR, OK )
|
CALL CHKXER( 'CGGSVD3', INFOT, NOUT, LERR, OK )
|
||||||
INFOT = 3
|
INFOT = 3
|
||||||
CALL CGGSVD3( 'N', 'N', '/', 0, 0, 0, DUMMYK, DUMMYL, A, 1, B,
|
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 )
|
$ INFO )
|
||||||
CALL CHKXER( 'CGGSVD3', INFOT, NOUT, LERR, OK )
|
CALL CHKXER( 'CGGSVD3', INFOT, NOUT, LERR, OK )
|
||||||
INFOT = 4
|
INFOT = 4
|
||||||
CALL CGGSVD3( 'N', 'N', 'N', -1, 0, 0, DUMMYK, DUMMYL, A, 1, B,
|
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 )
|
$ INFO )
|
||||||
CALL CHKXER( 'CGGSVD3', INFOT, NOUT, LERR, OK )
|
CALL CHKXER( 'CGGSVD3', INFOT, NOUT, LERR, OK )
|
||||||
INFOT = 5
|
INFOT = 5
|
||||||
CALL CGGSVD3( 'N', 'N', 'N', 0, -1, 0, DUMMYK, DUMMYL, A, 1, B,
|
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 )
|
$ INFO )
|
||||||
CALL CHKXER( 'CGGSVD3', INFOT, NOUT, LERR, OK )
|
CALL CHKXER( 'CGGSVD3', INFOT, NOUT, LERR, OK )
|
||||||
INFOT = 6
|
INFOT = 6
|
||||||
CALL CGGSVD3( 'N', 'N', 'N', 0, 0, -1, DUMMYK, DUMMYL, A, 1, B,
|
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 )
|
$ INFO )
|
||||||
CALL CHKXER( 'CGGSVD3', INFOT, NOUT, LERR, OK )
|
CALL CHKXER( 'CGGSVD3', INFOT, NOUT, LERR, OK )
|
||||||
INFOT = 10
|
INFOT = 10
|
||||||
CALL CGGSVD3( 'N', 'N', 'N', 2, 1, 1, DUMMYK, DUMMYL, A, 1, B,
|
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 )
|
$ INFO )
|
||||||
CALL CHKXER( 'CGGSVD3', INFOT, NOUT, LERR, OK )
|
CALL CHKXER( 'CGGSVD3', INFOT, NOUT, LERR, OK )
|
||||||
INFOT = 12
|
INFOT = 12
|
||||||
CALL CGGSVD3( 'N', 'N', 'N', 1, 1, 2, DUMMYK, DUMMYL, A, 1, B,
|
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 )
|
$ INFO )
|
||||||
CALL CHKXER( 'CGGSVD3', INFOT, NOUT, LERR, OK )
|
CALL CHKXER( 'CGGSVD3', INFOT, NOUT, LERR, OK )
|
||||||
INFOT = 16
|
INFOT = 16
|
||||||
CALL CGGSVD3( 'U', 'N', 'N', 2, 2, 2, DUMMYK, DUMMYL, A, 2, B,
|
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 )
|
$ INFO )
|
||||||
CALL CHKXER( 'CGGSVD3', INFOT, NOUT, LERR, OK )
|
CALL CHKXER( 'CGGSVD3', INFOT, NOUT, LERR, OK )
|
||||||
INFOT = 18
|
INFOT = 18
|
||||||
CALL CGGSVD3( 'N', 'V', 'N', 2, 2, 2, DUMMYK, DUMMYL, A, 2, B,
|
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 )
|
$ INFO )
|
||||||
CALL CHKXER( 'CGGSVD3', INFOT, NOUT, LERR, OK )
|
CALL CHKXER( 'CGGSVD3', INFOT, NOUT, LERR, OK )
|
||||||
INFOT = 20
|
INFOT = 20
|
||||||
CALL CGGSVD3( 'N', 'N', 'Q', 2, 2, 2, DUMMYK, DUMMYL, A, 2, B,
|
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 )
|
$ INFO )
|
||||||
CALL CHKXER( 'CGGSVD3', INFOT, NOUT, LERR, OK )
|
CALL CHKXER( 'CGGSVD3', INFOT, NOUT, LERR, OK )
|
||||||
NT = NT + 11
|
NT = NT + 11
|
||||||
|
@ -573,56 +573,56 @@
|
||||||
INFOT = 7
|
INFOT = 7
|
||||||
CALL CUNCSD( 'Y', 'Y', 'Y', 'Y', 'N', 'N',
|
CALL CUNCSD( 'Y', 'Y', 'Y', 'Y', 'N', 'N',
|
||||||
$ -1, 0, 0, A, 1, A,
|
$ -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,
|
$ A, 1, A, 1, A, 1, A,
|
||||||
$ 1, W, LW, RW, LW, IW, INFO )
|
$ 1, W, LW, RW, LW, IW, INFO )
|
||||||
CALL CHKXER( 'CUNCSD', INFOT, NOUT, LERR, OK )
|
CALL CHKXER( 'CUNCSD', INFOT, NOUT, LERR, OK )
|
||||||
INFOT = 8
|
INFOT = 8
|
||||||
CALL CUNCSD( 'Y', 'Y', 'Y', 'Y', 'N', 'N',
|
CALL CUNCSD( 'Y', 'Y', 'Y', 'Y', 'N', 'N',
|
||||||
$ 1, -1, 0, A, 1, A,
|
$ 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,
|
$ A, 1, A, 1, A, 1, A,
|
||||||
$ 1, W, LW, RW, LW, IW, INFO )
|
$ 1, W, LW, RW, LW, IW, INFO )
|
||||||
CALL CHKXER( 'CUNCSD', INFOT, NOUT, LERR, OK )
|
CALL CHKXER( 'CUNCSD', INFOT, NOUT, LERR, OK )
|
||||||
INFOT = 9
|
INFOT = 9
|
||||||
CALL CUNCSD( 'Y', 'Y', 'Y', 'Y', 'N', 'N',
|
CALL CUNCSD( 'Y', 'Y', 'Y', 'Y', 'N', 'N',
|
||||||
$ 1, 1, -1, A, 1, A,
|
$ 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,
|
$ A, 1, A, 1, A, 1, A,
|
||||||
$ 1, W, LW, RW, LW, IW, INFO )
|
$ 1, W, LW, RW, LW, IW, INFO )
|
||||||
CALL CHKXER( 'CUNCSD', INFOT, NOUT, LERR, OK )
|
CALL CHKXER( 'CUNCSD', INFOT, NOUT, LERR, OK )
|
||||||
INFOT = 11
|
INFOT = 11
|
||||||
CALL CUNCSD( 'Y', 'Y', 'Y', 'Y', 'N', 'N',
|
CALL CUNCSD( 'Y', 'Y', 'Y', 'Y', 'N', 'N',
|
||||||
$ 1, 1, 1, A, -1, A,
|
$ 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,
|
$ A, 1, A, 1, A, 1, A,
|
||||||
$ 1, W, LW, RW, LW, IW, INFO )
|
$ 1, W, LW, RW, LW, IW, INFO )
|
||||||
CALL CHKXER( 'CUNCSD', INFOT, NOUT, LERR, OK )
|
CALL CHKXER( 'CUNCSD', INFOT, NOUT, LERR, OK )
|
||||||
INFOT = 20
|
INFOT = 20
|
||||||
CALL CUNCSD( 'Y', 'Y', 'Y', 'Y', 'N', 'N',
|
CALL CUNCSD( 'Y', 'Y', 'Y', 'Y', 'N', 'N',
|
||||||
$ 1, 1, 1, A, 1, A,
|
$ 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,
|
$ A, -1, A, 1, A, 1, A,
|
||||||
$ 1, W, LW, RW, LW, IW, INFO )
|
$ 1, W, LW, RW, LW, IW, INFO )
|
||||||
CALL CHKXER( 'CUNCSD', INFOT, NOUT, LERR, OK )
|
CALL CHKXER( 'CUNCSD', INFOT, NOUT, LERR, OK )
|
||||||
INFOT = 22
|
INFOT = 22
|
||||||
CALL CUNCSD( 'Y', 'Y', 'Y', 'Y', 'N', 'N',
|
CALL CUNCSD( 'Y', 'Y', 'Y', 'Y', 'N', 'N',
|
||||||
$ 1, 1, 1, A, 1, A,
|
$ 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,
|
$ A, 1, A, -1, A, 1, A,
|
||||||
$ 1, W, LW, RW, LW, IW, INFO )
|
$ 1, W, LW, RW, LW, IW, INFO )
|
||||||
CALL CHKXER( 'CUNCSD', INFOT, NOUT, LERR, OK )
|
CALL CHKXER( 'CUNCSD', INFOT, NOUT, LERR, OK )
|
||||||
INFOT = 24
|
INFOT = 24
|
||||||
CALL CUNCSD( 'Y', 'Y', 'Y', 'Y', 'N', 'N',
|
CALL CUNCSD( 'Y', 'Y', 'Y', 'Y', 'N', 'N',
|
||||||
$ 1, 1, 1, A, 1, A,
|
$ 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,
|
$ A, 1, A, 1, A, -1, A,
|
||||||
$ 1, W, LW, RW, LW, IW, INFO )
|
$ 1, W, LW, RW, LW, IW, INFO )
|
||||||
CALL CHKXER( 'CUNCSD', INFOT, NOUT, LERR, OK )
|
CALL CHKXER( 'CUNCSD', INFOT, NOUT, LERR, OK )
|
||||||
INFOT = 26
|
INFOT = 26
|
||||||
CALL CUNCSD( 'Y', 'Y', 'Y', 'Y', 'N', 'N',
|
CALL CUNCSD( 'Y', 'Y', 'Y', 'Y', 'N', 'N',
|
||||||
$ 1, 1, 1, A, 1, A,
|
$ 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,
|
$ A, 1, A, 1, A, 1, A,
|
||||||
$ -1, W, LW, RW, LW, IW, INFO )
|
$ -1, W, LW, RW, LW, IW, INFO )
|
||||||
CALL CHKXER( 'CUNCSD', INFOT, NOUT, LERR, OK )
|
CALL CHKXER( 'CUNCSD', INFOT, NOUT, LERR, OK )
|
||||||
|
|
|
@ -25,7 +25,7 @@
|
||||||
*> eigenvalue is negative), and otherwise it returns .FALSE..
|
*> eigenvalue is negative), and otherwise it returns .FALSE..
|
||||||
*>
|
*>
|
||||||
*> It is used by the test routine CDRGES to test whether the driver
|
*> It is used by the test routine CDRGES to test whether the driver
|
||||||
*> routine CGGES succesfully sorts eigenvalues.
|
*> routine CGGES successfully sorts eigenvalues.
|
||||||
*> \endverbatim
|
*> \endverbatim
|
||||||
*
|
*
|
||||||
* Arguments:
|
* Arguments:
|
||||||
|
@ -51,17 +51,17 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2011
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup complex_eig
|
*> \ingroup complex_eig
|
||||||
*
|
*
|
||||||
* =====================================================================
|
* =====================================================================
|
||||||
LOGICAL FUNCTION CLCTES( Z, D )
|
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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2011
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
COMPLEX D, Z
|
COMPLEX D, Z
|
||||||
|
|
|
@ -22,8 +22,8 @@
|
||||||
*>
|
*>
|
||||||
*> CSLECT returns .TRUE. if the eigenvalue Z is to be selected,
|
*> CSLECT returns .TRUE. if the eigenvalue Z is to be selected,
|
||||||
*> otherwise it returns .FALSE.
|
*> otherwise it returns .FALSE.
|
||||||
*> It is used by CCHK41 to test if CGEES succesfully sorts eigenvalues,
|
*> It is used by CCHK41 to test if CGEES successfully sorts eigenvalues,
|
||||||
*> and by CCHK43 to test if CGEESX succesfully sorts eigenvalues.
|
*> and by CCHK43 to test if CGEESX successfully sorts eigenvalues.
|
||||||
*>
|
*>
|
||||||
*> The common block /SSLCT/ controls how eigenvalues are selected.
|
*> The common block /SSLCT/ controls how eigenvalues are selected.
|
||||||
*> If SELOPT = 0, then CSLECT return .TRUE. when real(Z) is less than
|
*> If SELOPT = 0, then CSLECT return .TRUE. when real(Z) is less than
|
||||||
|
@ -49,17 +49,17 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2011
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup complex_eig
|
*> \ingroup complex_eig
|
||||||
*
|
*
|
||||||
* =====================================================================
|
* =====================================================================
|
||||||
LOGICAL FUNCTION CSLECT( Z )
|
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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2011
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
COMPLEX Z
|
COMPLEX Z
|
||||||
|
|
|
@ -483,7 +483,7 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2015
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup double_eig
|
*> \ingroup double_eig
|
||||||
*
|
*
|
||||||
|
@ -493,10 +493,10 @@
|
||||||
$ Y, Z, Q, LDQ, PT, LDPT, U, VT, WORK, LWORK,
|
$ Y, Z, Q, LDQ, PT, LDPT, U, VT, WORK, LWORK,
|
||||||
$ IWORK, NOUT, INFO )
|
$ 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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2015
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
INTEGER INFO, LDA, LDPT, LDQ, LDX, LWORK, NOUT, NRHS,
|
INTEGER INFO, LDA, LDPT, LDQ, LDX, LWORK, NOUT, NRHS,
|
||||||
|
@ -545,10 +545,10 @@
|
||||||
EXTERNAL DLAMCH, DLARND, DSXT1
|
EXTERNAL DLAMCH, DLARND, DSXT1
|
||||||
* ..
|
* ..
|
||||||
* .. External Subroutines ..
|
* .. External Subroutines ..
|
||||||
EXTERNAL ALASUM, DBDSDC, DBDSQR, DBDSVDX, DBDT01, DBDT02,
|
EXTERNAL ALASUM, DBDSDC, DBDSQR, DBDSVDX, DBDT01,
|
||||||
$ DBDT03, DBDT04, DCOPY, DGEBRD, DGEMM, DLABAD,
|
$ DBDT02, DBDT03, DBDT04, DCOPY, DGEBRD,
|
||||||
$ DLACPY, DLAHD2, DLASET, DLATMR, DLATMS,
|
$ DGEMM, DLABAD, DLACPY, DLAHD2, DLASET,
|
||||||
$ DORGBR, DORT01, XERBLA
|
$ DLATMR, DLATMS, DORGBR, DORT01, XERBLA
|
||||||
* ..
|
* ..
|
||||||
* .. Intrinsic Functions ..
|
* .. Intrinsic Functions ..
|
||||||
INTRINSIC ABS, EXP, INT, LOG, MAX, MIN, SQRT
|
INTRINSIC ABS, EXP, INT, LOG, MAX, MIN, SQRT
|
||||||
|
@ -1143,7 +1143,7 @@
|
||||||
IWBD = IWBS + MNMIN
|
IWBD = IWBS + MNMIN
|
||||||
IWBE = IWBD + MNMIN
|
IWBE = IWBD + MNMIN
|
||||||
IWBZ = IWBE + MNMIN
|
IWBZ = IWBE + MNMIN
|
||||||
IWWORK = IWBZ + MNMIN*(MNMIN*2+1)
|
IWWORK = IWBZ + 2*MNMIN*(MNMIN+1)
|
||||||
MNMIN2 = MAX( 1,MNMIN*2 )
|
MNMIN2 = MAX( 1,MNMIN*2 )
|
||||||
*
|
*
|
||||||
CALL DCOPY( MNMIN, BD, 1, WORK( IWBD ), 1 )
|
CALL DCOPY( MNMIN, BD, 1, WORK( IWBD ), 1 )
|
||||||
|
|
|
@ -1033,17 +1033,17 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2015
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup double_eig
|
*> \ingroup double_eig
|
||||||
*
|
*
|
||||||
* =====================================================================
|
* =====================================================================
|
||||||
PROGRAM DCHKEE
|
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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2015
|
* June 2016
|
||||||
*
|
*
|
||||||
* =====================================================================
|
* =====================================================================
|
||||||
*
|
*
|
||||||
|
@ -2105,6 +2105,7 @@
|
||||||
MAXTYP = 26
|
MAXTYP = 26
|
||||||
NTYPES = MIN( MAXTYP, NTYPES )
|
NTYPES = MIN( MAXTYP, NTYPES )
|
||||||
CALL ALAREQ( C3, NTYPES, DOTYPE, MAXTYP, NIN, NOUT )
|
CALL ALAREQ( C3, NTYPES, DOTYPE, MAXTYP, NIN, NOUT )
|
||||||
|
CALL XLAENV(1,1)
|
||||||
IF( TSTCHK .AND. TSTERR )
|
IF( TSTCHK .AND. TSTERR )
|
||||||
$ CALL DERRGG( C3, NOUT )
|
$ CALL DERRGG( C3, NOUT )
|
||||||
DO 350 I = 1, NPARMS
|
DO 350 I = 1, NPARMS
|
||||||
|
@ -2164,6 +2165,7 @@
|
||||||
*
|
*
|
||||||
* Blocked version
|
* Blocked version
|
||||||
*
|
*
|
||||||
|
CALL XLAENV(16, 2)
|
||||||
CALL DDRGES3( NN, NVAL, MAXTYP, DOTYPE, ISEED, THRESH, NOUT,
|
CALL DDRGES3( NN, NVAL, MAXTYP, DOTYPE, ISEED, THRESH, NOUT,
|
||||||
$ A( 1, 1 ), NMAX, A( 1, 2 ), A( 1, 3 ),
|
$ A( 1, 1 ), NMAX, A( 1, 2 ), A( 1, 3 ),
|
||||||
$ A( 1, 4 ), A( 1, 7 ), NMAX, A( 1, 8 ),
|
$ A( 1, 4 ), A( 1, 7 ), NMAX, A( 1, 8 ),
|
||||||
|
@ -2351,6 +2353,7 @@
|
||||||
* GSV: Generalized Singular Value Decomposition
|
* GSV: Generalized Singular Value Decomposition
|
||||||
* ----------------------------------------------
|
* ----------------------------------------------
|
||||||
*
|
*
|
||||||
|
CALL XLAENV(1,1)
|
||||||
IF( TSTERR )
|
IF( TSTERR )
|
||||||
$ CALL DERRGG( 'GSV', NOUT )
|
$ CALL DERRGG( 'GSV', NOUT )
|
||||||
CALL DCKGSV( NN, MVAL, PVAL, NVAL, NTYPES, ISEED, THRESH, NMAX,
|
CALL DCKGSV( NN, MVAL, PVAL, NVAL, NTYPES, ISEED, THRESH, NMAX,
|
||||||
|
|
|
@ -72,7 +72,7 @@
|
||||||
*> and each type of matrix, one matrix will be generated and used
|
*> and each type of matrix, one matrix will be generated and used
|
||||||
*> to test the nonsymmetric eigenroutines. For each matrix, 15
|
*> to test the nonsymmetric eigenroutines. For each matrix, 15
|
||||||
*> tests will be performed. The first twelve "test ratios" should be
|
*> 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
|
*> T
|
||||||
*> (1) | A - U H V | / ( |A| n ulp )
|
*> (1) | A - U H V | / ( |A| n ulp )
|
||||||
|
@ -132,7 +132,7 @@
|
||||||
*> |beta(Q,Z computed) - beta(Q,Z not computed)|/|P| ) / ulp
|
*> |beta(Q,Z computed) - beta(Q,Z not computed)|/|P| ) / ulp
|
||||||
*>
|
*>
|
||||||
*> In addition, the normalization of L and R are checked, and compared
|
*> In addition, the normalization of L and R are checked, and compared
|
||||||
*> with the threshhold THRSHN.
|
*> with the threshold THRSHN.
|
||||||
*>
|
*>
|
||||||
*> Test Matrices
|
*> Test Matrices
|
||||||
*> ---- --------
|
*> ---- --------
|
||||||
|
@ -298,7 +298,7 @@
|
||||||
*> \param[in] THRSHN
|
*> \param[in] THRSHN
|
||||||
*> \verbatim
|
*> \verbatim
|
||||||
*> THRSHN is DOUBLE PRECISION
|
*> 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
|
*> If the normalization of any eigenvector differs from 1 by
|
||||||
*> more than THRSHN*ulp, then a special error message will be
|
*> more than THRSHN*ulp, then a special error message will be
|
||||||
*> printed. (This is handled separately from the other tests,
|
*> printed. (This is handled separately from the other tests,
|
||||||
|
@ -500,7 +500,7 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2011
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup double_eig
|
*> \ingroup double_eig
|
||||||
*
|
*
|
||||||
|
@ -511,10 +511,10 @@
|
||||||
$ BETA1, ALPHR3, ALPHI3, BETA3, EVECTL, EVECTR,
|
$ BETA1, ALPHR3, ALPHI3, BETA3, EVECTL, EVECTR,
|
||||||
$ WORK, LWORK, LLWORK, RESULT, INFO )
|
$ 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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2011
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
LOGICAL TSTDIF
|
LOGICAL TSTDIF
|
||||||
|
|
|
@ -52,7 +52,7 @@
|
||||||
*> number of matrix "TYPES" are specified. For each size ("N")
|
*> 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 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
|
*> 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).
|
*> the tests (5), (11) and (13).
|
||||||
*>
|
*>
|
||||||
*>
|
*>
|
||||||
|
@ -393,7 +393,7 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2011
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup double_eig
|
*> \ingroup double_eig
|
||||||
*
|
*
|
||||||
|
@ -403,10 +403,10 @@
|
||||||
$ ALPHAI, BETA, WORK, LWORK, RESULT, BWORK,
|
$ ALPHAI, BETA, WORK, LWORK, RESULT, BWORK,
|
||||||
$ INFO )
|
$ 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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2011
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
INTEGER INFO, LDA, LDQ, LWORK, NOUNIT, NSIZES, NTYPES
|
INTEGER INFO, LDA, LDQ, LWORK, NOUNIT, NSIZES, NTYPES
|
||||||
|
|
|
@ -52,7 +52,7 @@
|
||||||
*> number of matrix "TYPES" are specified. For each size ("N")
|
*> 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 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
|
*> 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).
|
*> the tests (5), (11) and (13).
|
||||||
*>
|
*>
|
||||||
*>
|
*>
|
||||||
|
@ -403,7 +403,7 @@
|
||||||
$ ALPHAI, BETA, WORK, LWORK, RESULT, BWORK,
|
$ ALPHAI, BETA, WORK, LWORK, RESULT, BWORK,
|
||||||
$ INFO )
|
$ 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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* February 2015
|
* February 2015
|
||||||
|
|
|
@ -55,7 +55,7 @@
|
||||||
*> number of matrix "types" are specified. For each size ("n")
|
*> 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 each type of matrix, a pair of matrices (A, B) will be generated
|
||||||
*> and used for testing. For each matrix pair, the following tests
|
*> 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:
|
*> Results from DGGEV:
|
||||||
*>
|
*>
|
||||||
|
@ -398,7 +398,7 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2015
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup double_eig
|
*> \ingroup double_eig
|
||||||
*
|
*
|
||||||
|
@ -408,10 +408,10 @@
|
||||||
$ ALPHAR, ALPHAI, BETA, ALPHR1, ALPHI1, BETA1,
|
$ ALPHAR, ALPHAI, BETA, ALPHR1, ALPHI1, BETA1,
|
||||||
$ WORK, LWORK, RESULT, INFO )
|
$ 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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2015
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
INTEGER INFO, LDA, LDQ, LDQE, LWORK, NOUNIT, NSIZES,
|
INTEGER INFO, LDA, LDQ, LDQE, LWORK, NOUNIT, NSIZES,
|
||||||
|
|
|
@ -55,7 +55,7 @@
|
||||||
*> number of matrix "types" are specified. For each size ("n")
|
*> 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 each type of matrix, a pair of matrices (A, B) will be generated
|
||||||
*> and used for testing. For each matrix pair, the following tests
|
*> 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:
|
*> Results from DGGEV3:
|
||||||
*>
|
*>
|
||||||
|
@ -408,7 +408,7 @@
|
||||||
$ ALPHAR, ALPHAI, BETA, ALPHR1, ALPHI1, BETA1,
|
$ ALPHAR, ALPHAI, BETA, ALPHR1, ALPHI1, BETA1,
|
||||||
$ WORK, LWORK, RESULT, INFO )
|
$ 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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* February 2015
|
* February 2015
|
||||||
|
|
|
@ -60,7 +60,7 @@
|
||||||
*> to test DGGESX.
|
*> to test DGGESX.
|
||||||
*>
|
*>
|
||||||
*> For each matrix pair, the following tests will be performed and
|
*> 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 )
|
*> (1) | A - Q S Z' | / ( |A| n ulp )
|
||||||
*>
|
*>
|
||||||
|
@ -350,7 +350,7 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2011
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup double_eig
|
*> \ingroup double_eig
|
||||||
*
|
*
|
||||||
|
@ -359,10 +359,10 @@
|
||||||
$ BI, Z, Q, ALPHAR, ALPHAI, BETA, C, LDC, S,
|
$ BI, Z, Q, ALPHAR, ALPHAI, BETA, C, LDC, S,
|
||||||
$ WORK, LWORK, IWORK, LIWORK, BWORK, INFO )
|
$ 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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2011
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
INTEGER INFO, LDA, LDC, LIWORK, LWORK, NCMAX, NIN,
|
INTEGER INFO, LDA, LDC, LIWORK, LWORK, NCMAX, NIN,
|
||||||
|
|
|
@ -50,7 +50,7 @@
|
||||||
*> ``exactly'' (see DLATM6).
|
*> ``exactly'' (see DLATM6).
|
||||||
*>
|
*>
|
||||||
*> For each matrix pair, the following tests will be performed and
|
*> 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
|
*> (1) max over all left eigenvalue/-vector pairs (beta/alpha,l) of
|
||||||
*>
|
*>
|
||||||
|
@ -300,7 +300,7 @@
|
||||||
$ RSCALE, S, DTRU, DIF, DIFTRU, WORK, LWORK,
|
$ RSCALE, S, DTRU, DIF, DIFTRU, WORK, LWORK,
|
||||||
$ IWORK, LIWORK, RESULT, BWORK, INFO )
|
$ 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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* April 2012
|
* April 2012
|
||||||
|
|
|
@ -346,7 +346,7 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2015
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup double_eig
|
*> \ingroup double_eig
|
||||||
*
|
*
|
||||||
|
@ -355,10 +355,10 @@
|
||||||
$ A, LDA, U, LDU, VT, LDVT, ASAV, USAV, VTSAV, S,
|
$ A, LDA, U, LDU, VT, LDVT, ASAV, USAV, VTSAV, S,
|
||||||
$ SSAV, E, WORK, LWORK, IWORK, NOUT, INFO )
|
$ 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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2015
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
INTEGER INFO, LDA, LDU, LDVT, LWORK, NOUT, NSIZES,
|
INTEGER INFO, LDA, LDU, LDVT, LWORK, NOUT, NSIZES,
|
||||||
|
@ -404,8 +404,8 @@
|
||||||
* ..
|
* ..
|
||||||
* .. External Subroutines ..
|
* .. External Subroutines ..
|
||||||
EXTERNAL ALASVM, DBDT01, DGEJSV, DGESDD, DGESVD,
|
EXTERNAL ALASVM, DBDT01, DGEJSV, DGESDD, DGESVD,
|
||||||
$ DGESVDX, DGESVJ, DLABAD, DLACPY, DLASET, DLATMS,
|
$ DGESVDX, DGESVJ, DLABAD, DLACPY, DLASET,
|
||||||
$ DORT01, DORT03, XERBLA
|
$ DLATMS, DORT01, DORT03, XERBLA
|
||||||
* ..
|
* ..
|
||||||
* .. Intrinsic Functions ..
|
* .. Intrinsic Functions ..
|
||||||
INTRINSIC ABS, DBLE, INT, MAX, MIN
|
INTRINSIC ABS, DBLE, INT, MAX, MIN
|
||||||
|
@ -1141,7 +1141,7 @@
|
||||||
$ / '21 = | S - Spartial | / ( min(M,N) ulp |S| )',
|
$ / '21 = | S - Spartial | / ( min(M,N) ulp |S| )',
|
||||||
$ / '22 = 0 if S contains min(M,N) nonnegative values in',
|
$ / '22 = 0 if S contains min(M,N) nonnegative values in',
|
||||||
$ ' decreasing order, else 1/ulp',
|
$ ' 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) ',
|
$ ' DGESVDX(V,V,A) ',
|
||||||
$ / '24 = | I - U**T U | / ( M ulp ) ',
|
$ / '24 = | I - U**T U | / ( M ulp ) ',
|
||||||
$ / '25 = | I - VT VT**T | / ( N ulp ) ',
|
$ / '25 = | I - VT VT**T | / ( N ulp ) ',
|
||||||
|
|
|
@ -379,7 +379,7 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2011
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup double_eig
|
*> \ingroup double_eig
|
||||||
*
|
*
|
||||||
|
@ -388,10 +388,10 @@
|
||||||
$ NOUNIT, A, LDA, H, HT, WR, WI, WRT, WIT, VS,
|
$ NOUNIT, A, LDA, H, HT, WR, WI, WRT, WIT, VS,
|
||||||
$ LDVS, RESULT, WORK, NWORK, IWORK, BWORK, INFO )
|
$ 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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2011
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
INTEGER INFO, LDA, LDVS, NOUNIT, NSIZES, NTYPES, NWORK
|
INTEGER INFO, LDA, LDVS, NOUNIT, NSIZES, NTYPES, NWORK
|
||||||
|
@ -948,7 +948,7 @@
|
||||||
$ ' 1/ulp otherwise', /
|
$ ' 1/ulp otherwise', /
|
||||||
$ ' 12 = 0 if WR, WI same no matter if VS computed (sort),',
|
$ ' 12 = 0 if WR, WI same no matter if VS computed (sort),',
|
||||||
$ ' 1/ulp otherwise', /
|
$ ' 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, ',' ),
|
9993 FORMAT( ' N=', I5, ', IWK=', I2, ', seed=', 4( I4, ',' ),
|
||||||
$ ' type ', I2, ', test(', I2, ')=', G10.3 )
|
$ ' type ', I2, ', test(', I2, ')=', G10.3 )
|
||||||
9992 FORMAT( ' DDRVES: ', A, ' returned INFO=', I6, '.', / 9X, 'N=',
|
9992 FORMAT( ' DDRVES: ', A, ' returned INFO=', I6, '.', / 9X, 'N=',
|
||||||
|
|
|
@ -444,7 +444,7 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2011
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup double_eig
|
*> \ingroup double_eig
|
||||||
*
|
*
|
||||||
|
@ -454,10 +454,10 @@
|
||||||
$ WIT, WRTMP, WITMP, VS, LDVS, VS1, RESULT, WORK,
|
$ WIT, WRTMP, WITMP, VS, LDVS, VS1, RESULT, WORK,
|
||||||
$ LWORK, IWORK, BWORK, INFO )
|
$ 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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2011
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
INTEGER INFO, LDA, LDVS, LWORK, NIUNIT, NOUNIT, NSIZES,
|
INTEGER INFO, LDA, LDVS, LWORK, NIUNIT, NOUNIT, NSIZES,
|
||||||
|
@ -933,7 +933,7 @@
|
||||||
$ ' 1/ulp otherwise', /
|
$ ' 1/ulp otherwise', /
|
||||||
$ ' 12 = 0 if WR, WI same no matter what else computed ',
|
$ ' 12 = 0 if WR, WI same no matter what else computed ',
|
||||||
$ '(sort), 1/ulp otherwise', /
|
$ '(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,',
|
$ / ' 14 = 0 if RCONDE same no matter what else computed,',
|
||||||
$ ' 1/ulp otherwise', /
|
$ ' 1/ulp otherwise', /
|
||||||
$ ' 15 = 0 if RCONDv same no matter what else computed,',
|
$ ' 15 = 0 if RCONDv same no matter what else computed,',
|
||||||
|
|
|
@ -471,7 +471,7 @@
|
||||||
*> \verbatim
|
*> \verbatim
|
||||||
*> INFO is INTEGER
|
*> INFO is INTEGER
|
||||||
*> If 0, then successful exit.
|
*> 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
|
*> If >0, DLATMR, SLATMS, SLATME or DGET23 returned an error
|
||||||
*> code, and INFO is its absolute value.
|
*> code, and INFO is its absolute value.
|
||||||
*>
|
*>
|
||||||
|
@ -510,7 +510,7 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2011
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup double_eig
|
*> \ingroup double_eig
|
||||||
*
|
*
|
||||||
|
@ -521,10 +521,10 @@
|
||||||
$ RCDVIN, RCONDE, RCNDE1, RCDEIN, SCALE, SCALE1,
|
$ RCDVIN, RCONDE, RCNDE1, RCDEIN, SCALE, SCALE1,
|
||||||
$ RESULT, WORK, NWORK, IWORK, INFO )
|
$ 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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2011
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
INTEGER INFO, LDA, LDLRE, LDVL, LDVR, NIUNIT, NOUNIT,
|
INTEGER INFO, LDA, LDLRE, LDVL, LDVR, NIUNIT, NOUNIT,
|
||||||
|
|
|
@ -48,17 +48,17 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2015
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup double_eig
|
*> \ingroup double_eig
|
||||||
*
|
*
|
||||||
* =====================================================================
|
* =====================================================================
|
||||||
SUBROUTINE DERRBD( PATH, NUNIT )
|
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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2015
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
CHARACTER*3 PATH
|
CHARACTER*3 PATH
|
||||||
|
@ -89,8 +89,8 @@
|
||||||
EXTERNAL LSAMEN
|
EXTERNAL LSAMEN
|
||||||
* ..
|
* ..
|
||||||
* .. External Subroutines ..
|
* .. External Subroutines ..
|
||||||
EXTERNAL CHKXER, DBDSDC, DBDSQR, DBDSVDX, DGEBD2, DGEBRD, DORGBR,
|
EXTERNAL CHKXER, DBDSDC, DBDSQR, DBDSVDX, DGEBD2,
|
||||||
$ DORMBR
|
$ DGEBRD, DORGBR, DORMBR
|
||||||
* ..
|
* ..
|
||||||
* .. Scalars in Common ..
|
* .. Scalars in Common ..
|
||||||
LOGICAL LERR, OK
|
LOGICAL LERR, OK
|
||||||
|
|
|
@ -61,17 +61,17 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2015
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup double_eig
|
*> \ingroup double_eig
|
||||||
*
|
*
|
||||||
* =====================================================================
|
* =====================================================================
|
||||||
SUBROUTINE DERRED( PATH, NUNIT )
|
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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2015
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
CHARACTER*3 PATH
|
CHARACTER*3 PATH
|
||||||
|
@ -475,7 +475,7 @@
|
||||||
CALL DGESVDX( 'V', 'N', 'A', 2, 2, A, 2, ZERO, ZERO,
|
CALL DGESVDX( 'V', 'N', 'A', 2, 2, A, 2, ZERO, ZERO,
|
||||||
$ 0, 0, NS, S, U, 1, VT, 1, W, 1, IW, INFO )
|
$ 0, 0, NS, S, U, 1, VT, 1, W, 1, IW, INFO )
|
||||||
CALL CHKXER( 'DGESVDX', INFOT, NOUT, LERR, OK )
|
CALL CHKXER( 'DGESVDX', INFOT, NOUT, LERR, OK )
|
||||||
INFOT = 16
|
INFOT = 17
|
||||||
CALL DGESVDX( 'N', 'V', 'A', 2, 2, A, 2, ZERO, ZERO,
|
CALL DGESVDX( 'N', 'V', 'A', 2, 2, A, 2, ZERO, ZERO,
|
||||||
$ 0, 0, NS, S, U, 1, VT, 1, W, 1, IW, INFO )
|
$ 0, 0, NS, S, U, 1, VT, 1, W, 1, IW, INFO )
|
||||||
CALL CHKXER( 'DGESVDX', INFOT, NOUT, LERR, OK )
|
CALL CHKXER( 'DGESVDX', INFOT, NOUT, LERR, OK )
|
||||||
|
|
|
@ -50,17 +50,17 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2015
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup double_eig
|
*> \ingroup double_eig
|
||||||
*
|
*
|
||||||
* =====================================================================
|
* =====================================================================
|
||||||
SUBROUTINE DERRGG( PATH, NUNIT )
|
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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2015
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
CHARACTER*3 PATH
|
CHARACTER*3 PATH
|
||||||
|
@ -83,7 +83,7 @@
|
||||||
* ..
|
* ..
|
||||||
* .. Local Arrays ..
|
* .. Local Arrays ..
|
||||||
LOGICAL BW( NMAX ), SEL( NMAX )
|
LOGICAL BW( NMAX ), SEL( NMAX )
|
||||||
INTEGER IW( NMAX )
|
INTEGER IW( NMAX ), IDUM(NMAX)
|
||||||
DOUBLE PRECISION A( NMAX, NMAX ), B( NMAX, NMAX ), LS( NMAX ),
|
DOUBLE PRECISION A( NMAX, NMAX ), B( NMAX, NMAX ), LS( NMAX ),
|
||||||
$ Q( NMAX, NMAX ), R1( NMAX ), R2( NMAX ),
|
$ Q( NMAX, NMAX ), R1( NMAX ), R2( NMAX ),
|
||||||
$ R3( NMAX ), RCE( 2 ), RCV( 2 ), RS( NMAX ),
|
$ R3( NMAX ), RCE( 2 ), RCV( 2 ), RS( NMAX ),
|
||||||
|
@ -305,47 +305,47 @@
|
||||||
SRNAMT = 'DGGSVD3'
|
SRNAMT = 'DGGSVD3'
|
||||||
INFOT = 1
|
INFOT = 1
|
||||||
CALL DGGSVD3( '/', 'N', 'N', 0, 0, 0, DUMMYK, DUMMYL, A, 1, B,
|
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 )
|
CALL CHKXER( 'DGGSVD3', INFOT, NOUT, LERR, OK )
|
||||||
INFOT = 2
|
INFOT = 2
|
||||||
CALL DGGSVD3( 'N', '/', 'N', 0, 0, 0, DUMMYK, DUMMYL, A, 1, B,
|
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 )
|
CALL CHKXER( 'DGGSVD3', INFOT, NOUT, LERR, OK )
|
||||||
INFOT = 3
|
INFOT = 3
|
||||||
CALL DGGSVD3( 'N', 'N', '/', 0, 0, 0, DUMMYK, DUMMYL, A, 1, B,
|
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 )
|
CALL CHKXER( 'DGGSVD3', INFOT, NOUT, LERR, OK )
|
||||||
INFOT = 4
|
INFOT = 4
|
||||||
CALL DGGSVD3( 'N', 'N', 'N', -1, 0, 0, DUMMYK, DUMMYL, A, 1, B,
|
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 )
|
CALL CHKXER( 'DGGSVD3', INFOT, NOUT, LERR, OK )
|
||||||
INFOT = 5
|
INFOT = 5
|
||||||
CALL DGGSVD3( 'N', 'N', 'N', 0, -1, 0, DUMMYK, DUMMYL, A, 1, B,
|
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 )
|
CALL CHKXER( 'DGGSVD3', INFOT, NOUT, LERR, OK )
|
||||||
INFOT = 6
|
INFOT = 6
|
||||||
CALL DGGSVD3( 'N', 'N', 'N', 0, 0, -1, DUMMYK, DUMMYL, A, 1, B,
|
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 )
|
CALL CHKXER( 'DGGSVD3', INFOT, NOUT, LERR, OK )
|
||||||
INFOT = 10
|
INFOT = 10
|
||||||
CALL DGGSVD3( 'N', 'N', 'N', 2, 1, 1, DUMMYK, DUMMYL, A, 1, B,
|
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 )
|
CALL CHKXER( 'DGGSVD3', INFOT, NOUT, LERR, OK )
|
||||||
INFOT = 12
|
INFOT = 12
|
||||||
CALL DGGSVD3( 'N', 'N', 'N', 1, 1, 2, DUMMYK, DUMMYL, A, 1, B,
|
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 )
|
CALL CHKXER( 'DGGSVD3', INFOT, NOUT, LERR, OK )
|
||||||
INFOT = 16
|
INFOT = 16
|
||||||
CALL DGGSVD3( 'U', 'N', 'N', 2, 2, 2, DUMMYK, DUMMYL, A, 2, B,
|
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 )
|
CALL CHKXER( 'DGGSVD3', INFOT, NOUT, LERR, OK )
|
||||||
INFOT = 18
|
INFOT = 18
|
||||||
CALL DGGSVD3( 'N', 'V', 'N', 1, 1, 2, DUMMYK, DUMMYL, A, 1, B,
|
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 )
|
CALL CHKXER( 'DGGSVD3', INFOT, NOUT, LERR, OK )
|
||||||
INFOT = 20
|
INFOT = 20
|
||||||
CALL DGGSVD3( 'N', 'N', 'Q', 1, 2, 1, DUMMYK, DUMMYL, A, 1, B,
|
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 )
|
CALL CHKXER( 'DGGSVD3', INFOT, NOUT, LERR, OK )
|
||||||
NT = NT + 11
|
NT = NT + 11
|
||||||
*
|
*
|
||||||
|
|
|
@ -26,7 +26,7 @@
|
||||||
*> .FALSE..
|
*> .FALSE..
|
||||||
*>
|
*>
|
||||||
*> It is used by the test routine DDRGES to test whether the driver
|
*> It is used by the test routine DDRGES to test whether the driver
|
||||||
*> routine DGGES succesfully sorts eigenvalues.
|
*> routine DGGES successfully sorts eigenvalues.
|
||||||
*> \endverbatim
|
*> \endverbatim
|
||||||
*
|
*
|
||||||
* Arguments:
|
* Arguments:
|
||||||
|
@ -61,17 +61,17 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2011
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup double_eig
|
*> \ingroup double_eig
|
||||||
*
|
*
|
||||||
* =====================================================================
|
* =====================================================================
|
||||||
LOGICAL FUNCTION DLCTES( ZR, ZI, D )
|
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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2011
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
DOUBLE PRECISION D, ZI, ZR
|
DOUBLE PRECISION D, ZI, ZR
|
||||||
|
|
|
@ -22,8 +22,8 @@
|
||||||
*>
|
*>
|
||||||
*> DSLECT returns .TRUE. if the eigenvalue ZR+sqrt(-1)*ZI is to be
|
*> DSLECT returns .TRUE. if the eigenvalue ZR+sqrt(-1)*ZI is to be
|
||||||
*> selected, and otherwise it returns .FALSE.
|
*> selected, and otherwise it returns .FALSE.
|
||||||
*> It is used by DCHK41 to test if DGEES succesfully sorts eigenvalues,
|
*> It is used by DCHK41 to test if DGEES successfully sorts eigenvalues,
|
||||||
*> and by DCHK43 to test if DGEESX succesfully sorts eigenvalues.
|
*> and by DCHK43 to test if DGEESX successfully sorts eigenvalues.
|
||||||
*>
|
*>
|
||||||
*> The common block /SSLCT/ controls how eigenvalues are selected.
|
*> The common block /SSLCT/ controls how eigenvalues are selected.
|
||||||
*> If SELOPT = 0, then DSLECT return .TRUE. when ZR is less than zero,
|
*> If SELOPT = 0, then DSLECT return .TRUE. when ZR is less than zero,
|
||||||
|
@ -55,17 +55,17 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2011
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup double_eig
|
*> \ingroup double_eig
|
||||||
*
|
*
|
||||||
* =====================================================================
|
* =====================================================================
|
||||||
LOGICAL FUNCTION DSLECT( ZR, ZI )
|
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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2011
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
DOUBLE PRECISION ZI, ZR
|
DOUBLE PRECISION ZI, ZR
|
||||||
|
|
|
@ -483,7 +483,7 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2015
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup single_eig
|
*> \ingroup single_eig
|
||||||
*
|
*
|
||||||
|
@ -493,10 +493,10 @@
|
||||||
$ Y, Z, Q, LDQ, PT, LDPT, U, VT, WORK, LWORK,
|
$ Y, Z, Q, LDQ, PT, LDPT, U, VT, WORK, LWORK,
|
||||||
$ IWORK, NOUT, INFO )
|
$ 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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2015
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
INTEGER INFO, LDA, LDPT, LDQ, LDX, LWORK, NOUT, NRHS,
|
INTEGER INFO, LDA, LDPT, LDQ, LDX, LWORK, NOUT, NRHS,
|
||||||
|
@ -545,10 +545,10 @@
|
||||||
EXTERNAL SLAMCH, SLARND, SSXT1
|
EXTERNAL SLAMCH, SLARND, SSXT1
|
||||||
* ..
|
* ..
|
||||||
* .. External Subroutines ..
|
* .. External Subroutines ..
|
||||||
EXTERNAL ALASUM, SBDSDC, SBDSQR, SBDSVDX, SBDT01, SBDT02,
|
EXTERNAL ALASUM, SBDSDC, SBDSQR, SBDSVDX, SBDT01,
|
||||||
$ SBDT03, SBDT04, SCOPY, SGEBRD, SGEMM, SLABAD,
|
$ SBDT02, SBDT03, SBDT04, SCOPY, SGEBRD,
|
||||||
$ SLACPY, SLAHD2, SLASET, SLATMR, SLATMS,
|
$ SGEMM, SLABAD, SLACPY, SLAHD2, SLASET,
|
||||||
$ SORGBR, SORT01, XERBLA
|
$ SLATMR, SLATMS, SORGBR, SORT01, XERBLA
|
||||||
* ..
|
* ..
|
||||||
* .. Intrinsic Functions ..
|
* .. Intrinsic Functions ..
|
||||||
INTRINSIC ABS, EXP, INT, LOG, MAX, MIN, SQRT
|
INTRINSIC ABS, EXP, INT, LOG, MAX, MIN, SQRT
|
||||||
|
@ -1143,7 +1143,7 @@
|
||||||
IWBD = IWBS + MNMIN
|
IWBD = IWBS + MNMIN
|
||||||
IWBE = IWBD + MNMIN
|
IWBE = IWBD + MNMIN
|
||||||
IWBZ = IWBE + MNMIN
|
IWBZ = IWBE + MNMIN
|
||||||
IWWORK = IWBZ + MNMIN*(MNMIN*2+1)
|
IWWORK = IWBZ + 2*MNMIN*(MNMIN+1)
|
||||||
MNMIN2 = MAX( 1,MNMIN*2 )
|
MNMIN2 = MAX( 1,MNMIN*2 )
|
||||||
*
|
*
|
||||||
CALL SCOPY( MNMIN, BD, 1, WORK( IWBD ), 1 )
|
CALL SCOPY( MNMIN, BD, 1, WORK( IWBD ), 1 )
|
||||||
|
|
|
@ -1033,17 +1033,17 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2015
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup single_eig
|
*> \ingroup single_eig
|
||||||
*
|
*
|
||||||
* =====================================================================
|
* =====================================================================
|
||||||
PROGRAM SCHKEE
|
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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2015
|
* June 2016
|
||||||
*
|
*
|
||||||
* =====================================================================
|
* =====================================================================
|
||||||
*
|
*
|
||||||
|
@ -2105,8 +2105,9 @@
|
||||||
MAXTYP = 26
|
MAXTYP = 26
|
||||||
NTYPES = MIN( MAXTYP, NTYPES )
|
NTYPES = MIN( MAXTYP, NTYPES )
|
||||||
CALL ALAREQ( C3, NTYPES, DOTYPE, MAXTYP, NIN, NOUT )
|
CALL ALAREQ( C3, NTYPES, DOTYPE, MAXTYP, NIN, NOUT )
|
||||||
|
CALL XLAENV(1,1)
|
||||||
IF( TSTCHK .AND. TSTERR )
|
IF( TSTCHK .AND. TSTERR )
|
||||||
$ CALL SERRGG( C3, NOUT )
|
& CALL SERRGG( C3, NOUT )
|
||||||
DO 350 I = 1, NPARMS
|
DO 350 I = 1, NPARMS
|
||||||
CALL XLAENV( 1, NBVAL( I ) )
|
CALL XLAENV( 1, NBVAL( I ) )
|
||||||
CALL XLAENV( 2, NBMIN( I ) )
|
CALL XLAENV( 2, NBMIN( I ) )
|
||||||
|
@ -2165,6 +2166,7 @@
|
||||||
*
|
*
|
||||||
* Blocked version
|
* Blocked version
|
||||||
*
|
*
|
||||||
|
CALL XLAENV(16,1)
|
||||||
CALL SDRGES3( NN, NVAL, MAXTYP, DOTYPE, ISEED, THRESH, NOUT,
|
CALL SDRGES3( NN, NVAL, MAXTYP, DOTYPE, ISEED, THRESH, NOUT,
|
||||||
$ A( 1, 1 ), NMAX, A( 1, 2 ), A( 1, 3 ),
|
$ A( 1, 1 ), NMAX, A( 1, 2 ), A( 1, 3 ),
|
||||||
$ A( 1, 4 ), A( 1, 7 ), NMAX, A( 1, 8 ),
|
$ A( 1, 4 ), A( 1, 7 ), NMAX, A( 1, 8 ),
|
||||||
|
@ -2353,6 +2355,7 @@
|
||||||
* GSV: Generalized Singular Value Decomposition
|
* GSV: Generalized Singular Value Decomposition
|
||||||
* ----------------------------------------------
|
* ----------------------------------------------
|
||||||
*
|
*
|
||||||
|
CALL XLAENV( 1, 1 )
|
||||||
IF( TSTERR )
|
IF( TSTERR )
|
||||||
$ CALL SERRGG( 'GSV', NOUT )
|
$ CALL SERRGG( 'GSV', NOUT )
|
||||||
CALL SCKGSV( NN, MVAL, PVAL, NVAL, NTYPES, ISEED, THRESH, NMAX,
|
CALL SCKGSV( NN, MVAL, PVAL, NVAL, NTYPES, ISEED, THRESH, NMAX,
|
||||||
|
|
|
@ -72,7 +72,7 @@
|
||||||
*> and each type of matrix, one matrix will be generated and used
|
*> and each type of matrix, one matrix will be generated and used
|
||||||
*> to test the nonsymmetric eigenroutines. For each matrix, 15
|
*> to test the nonsymmetric eigenroutines. For each matrix, 15
|
||||||
*> tests will be performed. The first twelve "test ratios" should be
|
*> 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
|
*> T
|
||||||
*> (1) | A - U H V | / ( |A| n ulp )
|
*> (1) | A - U H V | / ( |A| n ulp )
|
||||||
|
@ -132,7 +132,7 @@
|
||||||
*> |beta(Q,Z computed) - beta(Q,Z not computed)|/|P| ) / ulp
|
*> |beta(Q,Z computed) - beta(Q,Z not computed)|/|P| ) / ulp
|
||||||
*>
|
*>
|
||||||
*> In addition, the normalization of L and R are checked, and compared
|
*> In addition, the normalization of L and R are checked, and compared
|
||||||
*> with the threshhold THRSHN.
|
*> with the threshold THRSHN.
|
||||||
*>
|
*>
|
||||||
*> Test Matrices
|
*> Test Matrices
|
||||||
*> ---- --------
|
*> ---- --------
|
||||||
|
@ -298,7 +298,7 @@
|
||||||
*> \param[in] THRSHN
|
*> \param[in] THRSHN
|
||||||
*> \verbatim
|
*> \verbatim
|
||||||
*> THRSHN is REAL
|
*> 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
|
*> If the normalization of any eigenvector differs from 1 by
|
||||||
*> more than THRSHN*ulp, then a special error message will be
|
*> more than THRSHN*ulp, then a special error message will be
|
||||||
*> printed. (This is handled separately from the other tests,
|
*> printed. (This is handled separately from the other tests,
|
||||||
|
@ -500,7 +500,7 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2011
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup single_eig
|
*> \ingroup single_eig
|
||||||
*
|
*
|
||||||
|
@ -511,10 +511,10 @@
|
||||||
$ BETA1, ALPHR3, ALPHI3, BETA3, EVECTL, EVECTR,
|
$ BETA1, ALPHR3, ALPHI3, BETA3, EVECTL, EVECTR,
|
||||||
$ WORK, LWORK, LLWORK, RESULT, INFO )
|
$ 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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2011
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
LOGICAL TSTDIF
|
LOGICAL TSTDIF
|
||||||
|
|
|
@ -52,7 +52,7 @@
|
||||||
*> number of matrix "TYPES" are specified. For each size ("N")
|
*> 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 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
|
*> 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).
|
*> the tests (5), (11) and (13).
|
||||||
*>
|
*>
|
||||||
*>
|
*>
|
||||||
|
@ -393,7 +393,7 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2011
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup single_eig
|
*> \ingroup single_eig
|
||||||
*
|
*
|
||||||
|
@ -403,10 +403,10 @@
|
||||||
$ ALPHAI, BETA, WORK, LWORK, RESULT, BWORK,
|
$ ALPHAI, BETA, WORK, LWORK, RESULT, BWORK,
|
||||||
$ INFO )
|
$ 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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2011
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
INTEGER INFO, LDA, LDQ, LWORK, NOUNIT, NSIZES, NTYPES
|
INTEGER INFO, LDA, LDQ, LWORK, NOUNIT, NSIZES, NTYPES
|
||||||
|
|
|
@ -52,7 +52,7 @@
|
||||||
*> number of matrix "TYPES" are specified. For each size ("N")
|
*> 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 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
|
*> 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).
|
*> the tests (5), (11) and (13).
|
||||||
*>
|
*>
|
||||||
*>
|
*>
|
||||||
|
@ -403,7 +403,7 @@
|
||||||
$ ALPHAI, BETA, WORK, LWORK, RESULT, BWORK,
|
$ ALPHAI, BETA, WORK, LWORK, RESULT, BWORK,
|
||||||
$ INFO )
|
$ 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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* February 2015
|
* February 2015
|
||||||
|
|
|
@ -55,7 +55,7 @@
|
||||||
*> number of matrix "types" are specified. For each size ("n")
|
*> 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 each type of matrix, a pair of matrices (A, B) will be generated
|
||||||
*> and used for testing. For each matrix pair, the following tests
|
*> 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:
|
*> Results from SGGEV:
|
||||||
*>
|
*>
|
||||||
|
@ -398,7 +398,7 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2015
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup single_eig
|
*> \ingroup single_eig
|
||||||
*
|
*
|
||||||
|
@ -408,10 +408,10 @@
|
||||||
$ ALPHAR, ALPHAI, BETA, ALPHR1, ALPHI1, BETA1,
|
$ ALPHAR, ALPHAI, BETA, ALPHR1, ALPHI1, BETA1,
|
||||||
$ WORK, LWORK, RESULT, INFO )
|
$ 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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2015
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
INTEGER INFO, LDA, LDQ, LDQE, LWORK, NOUNIT, NSIZES,
|
INTEGER INFO, LDA, LDQ, LDQE, LWORK, NOUNIT, NSIZES,
|
||||||
|
|
|
@ -55,7 +55,7 @@
|
||||||
*> number of matrix "types" are specified. For each size ("n")
|
*> 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 each type of matrix, a pair of matrices (A, B) will be generated
|
||||||
*> and used for testing. For each matrix pair, the following tests
|
*> 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:
|
*> Results from SGGEV3:
|
||||||
*>
|
*>
|
||||||
|
@ -408,7 +408,7 @@
|
||||||
$ ALPHAR, ALPHAI, BETA, ALPHR1, ALPHI1, BETA1,
|
$ ALPHAR, ALPHAI, BETA, ALPHR1, ALPHI1, BETA1,
|
||||||
$ WORK, LWORK, RESULT, INFO )
|
$ 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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* February 2015
|
* February 2015
|
||||||
|
|
|
@ -60,7 +60,7 @@
|
||||||
*> to test SGGESX.
|
*> to test SGGESX.
|
||||||
*>
|
*>
|
||||||
*> For each matrix pair, the following tests will be performed and
|
*> 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 )
|
*> (1) | A - Q S Z' | / ( |A| n ulp )
|
||||||
*>
|
*>
|
||||||
|
@ -350,7 +350,7 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2011
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup single_eig
|
*> \ingroup single_eig
|
||||||
*
|
*
|
||||||
|
@ -359,10 +359,10 @@
|
||||||
$ AI, BI, Z, Q, ALPHAR, ALPHAI, BETA, C, LDC, S,
|
$ AI, BI, Z, Q, ALPHAR, ALPHAI, BETA, C, LDC, S,
|
||||||
$ WORK, LWORK, IWORK, LIWORK, BWORK, INFO )
|
$ 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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2011
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
INTEGER INFO, LDA, LDC, LIWORK, LWORK, NCMAX, NIN,
|
INTEGER INFO, LDA, LDC, LIWORK, LWORK, NCMAX, NIN,
|
||||||
|
|
|
@ -51,7 +51,7 @@
|
||||||
*> ``exactly'' (see SLATM6).
|
*> ``exactly'' (see SLATM6).
|
||||||
*>
|
*>
|
||||||
*> For each matrix pair, the following tests will be performed and
|
*> 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
|
*> (1) max over all left eigenvalue/-vector pairs (beta/alpha,l) of
|
||||||
*>
|
*>
|
||||||
|
@ -291,7 +291,7 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2011
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup single_eig
|
*> \ingroup single_eig
|
||||||
*
|
*
|
||||||
|
@ -301,10 +301,10 @@
|
||||||
$ RSCALE, S, STRU, DIF, DIFTRU, WORK, LWORK,
|
$ RSCALE, S, STRU, DIF, DIFTRU, WORK, LWORK,
|
||||||
$ IWORK, LIWORK, RESULT, BWORK, INFO )
|
$ 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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2011
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
INTEGER IHI, ILO, INFO, LDA, LIWORK, LWORK, NIN, NOUT,
|
INTEGER IHI, ILO, INFO, LDA, LIWORK, LWORK, NIN, NOUT,
|
||||||
|
|
|
@ -346,7 +346,7 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2015
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup single_eig
|
*> \ingroup single_eig
|
||||||
*
|
*
|
||||||
|
@ -355,10 +355,10 @@
|
||||||
$ A, LDA, U, LDU, VT, LDVT, ASAV, USAV, VTSAV, S,
|
$ A, LDA, U, LDU, VT, LDVT, ASAV, USAV, VTSAV, S,
|
||||||
$ SSAV, E, WORK, LWORK, IWORK, NOUT, INFO )
|
$ 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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2015
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
INTEGER INFO, LDA, LDU, LDVT, LWORK, NOUT, NSIZES,
|
INTEGER INFO, LDA, LDU, LDVT, LWORK, NOUT, NSIZES,
|
||||||
|
@ -404,8 +404,8 @@
|
||||||
* ..
|
* ..
|
||||||
* .. External Subroutines ..
|
* .. External Subroutines ..
|
||||||
EXTERNAL ALASVM, SBDT01, SGEJSV, SGESDD, SGESVD,
|
EXTERNAL ALASVM, SBDT01, SGEJSV, SGESDD, SGESVD,
|
||||||
$ SGESVDX, SGESVJ, SLABAD, SLACPY, SLASET, SLATMS,
|
$ SGESVDX, SGESVJ, SLABAD, SLACPY, SLASET,
|
||||||
$ SORT01, SORT03, XERBLA
|
$ SLATMS, SORT01, SORT03, XERBLA
|
||||||
* ..
|
* ..
|
||||||
* .. Intrinsic Functions ..
|
* .. Intrinsic Functions ..
|
||||||
INTRINSIC ABS, REAL, INT, MAX, MIN
|
INTRINSIC ABS, REAL, INT, MAX, MIN
|
||||||
|
|
|
@ -379,7 +379,7 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2011
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup single_eig
|
*> \ingroup single_eig
|
||||||
*
|
*
|
||||||
|
@ -388,10 +388,10 @@
|
||||||
$ NOUNIT, A, LDA, H, HT, WR, WI, WRT, WIT, VS,
|
$ NOUNIT, A, LDA, H, HT, WR, WI, WRT, WIT, VS,
|
||||||
$ LDVS, RESULT, WORK, NWORK, IWORK, BWORK, INFO )
|
$ 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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2011
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
INTEGER INFO, LDA, LDVS, NOUNIT, NSIZES, NTYPES, NWORK
|
INTEGER INFO, LDA, LDVS, NOUNIT, NSIZES, NTYPES, NWORK
|
||||||
|
@ -948,7 +948,7 @@
|
||||||
$ ' 1/ulp otherwise', /
|
$ ' 1/ulp otherwise', /
|
||||||
$ ' 12 = 0 if WR, WI same no matter if VS computed (sort),',
|
$ ' 12 = 0 if WR, WI same no matter if VS computed (sort),',
|
||||||
$ ' 1/ulp otherwise', /
|
$ ' 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, ',' ),
|
9993 FORMAT( ' N=', I5, ', IWK=', I2, ', seed=', 4( I4, ',' ),
|
||||||
$ ' type ', I2, ', test(', I2, ')=', G10.3 )
|
$ ' type ', I2, ', test(', I2, ')=', G10.3 )
|
||||||
9992 FORMAT( ' SDRVES: ', A, ' returned INFO=', I6, '.', / 9X, 'N=',
|
9992 FORMAT( ' SDRVES: ', A, ' returned INFO=', I6, '.', / 9X, 'N=',
|
||||||
|
|
|
@ -444,7 +444,7 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2011
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup single_eig
|
*> \ingroup single_eig
|
||||||
*
|
*
|
||||||
|
@ -454,10 +454,10 @@
|
||||||
$ WIT, WRTMP, WITMP, VS, LDVS, VS1, RESULT, WORK,
|
$ WIT, WRTMP, WITMP, VS, LDVS, VS1, RESULT, WORK,
|
||||||
$ LWORK, IWORK, BWORK, INFO )
|
$ 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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2011
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
INTEGER INFO, LDA, LDVS, LWORK, NIUNIT, NOUNIT, NSIZES,
|
INTEGER INFO, LDA, LDVS, LWORK, NIUNIT, NOUNIT, NSIZES,
|
||||||
|
@ -933,7 +933,7 @@
|
||||||
$ ' 1/ulp otherwise', /
|
$ ' 1/ulp otherwise', /
|
||||||
$ ' 12 = 0 if WR, WI same no matter what else computed ',
|
$ ' 12 = 0 if WR, WI same no matter what else computed ',
|
||||||
$ '(sort), 1/ulp otherwise', /
|
$ '(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,',
|
$ / ' 14 = 0 if RCONDE same no matter what else computed,',
|
||||||
$ ' 1/ulp otherwise', /
|
$ ' 1/ulp otherwise', /
|
||||||
$ ' 15 = 0 if RCONDv same no matter what else computed,',
|
$ ' 15 = 0 if RCONDv same no matter what else computed,',
|
||||||
|
|
|
@ -470,7 +470,7 @@
|
||||||
*> \verbatim
|
*> \verbatim
|
||||||
*> INFO is INTEGER
|
*> INFO is INTEGER
|
||||||
*> If 0, then successful exit.
|
*> 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
|
*> If >0, SLATMR, SLATMS, SLATME or SGET23 returned an error
|
||||||
*> code, and INFO is its absolute value.
|
*> code, and INFO is its absolute value.
|
||||||
*>
|
*>
|
||||||
|
@ -509,7 +509,7 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2011
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup single_eig
|
*> \ingroup single_eig
|
||||||
*
|
*
|
||||||
|
@ -520,10 +520,10 @@
|
||||||
$ RCDVIN, RCONDE, RCNDE1, RCDEIN, SCALE, SCALE1,
|
$ RCDVIN, RCONDE, RCNDE1, RCDEIN, SCALE, SCALE1,
|
||||||
$ RESULT, WORK, NWORK, IWORK, INFO )
|
$ 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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2011
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
INTEGER INFO, LDA, LDLRE, LDVL, LDVR, NIUNIT, NOUNIT,
|
INTEGER INFO, LDA, LDLRE, LDVL, LDVR, NIUNIT, NOUNIT,
|
||||||
|
|
|
@ -48,17 +48,17 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2015
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup single_eig
|
*> \ingroup single_eig
|
||||||
*
|
*
|
||||||
* =====================================================================
|
* =====================================================================
|
||||||
SUBROUTINE SERRBD( PATH, NUNIT )
|
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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2015
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
CHARACTER*3 PATH
|
CHARACTER*3 PATH
|
||||||
|
@ -89,8 +89,8 @@
|
||||||
EXTERNAL LSAMEN
|
EXTERNAL LSAMEN
|
||||||
* ..
|
* ..
|
||||||
* .. External Subroutines ..
|
* .. External Subroutines ..
|
||||||
EXTERNAL CHKXER, SBDSDC, SBDSQR, SBDSVDX, SGEBD2, SGEBRD, SORGBR,
|
EXTERNAL CHKXER, SBDSDC, SBDSQR, SBDSVDX, SGEBD2,
|
||||||
$ SORMBR
|
$ SGEBRD, SORGBR, SORMBR
|
||||||
* ..
|
* ..
|
||||||
* .. Scalars in Common ..
|
* .. Scalars in Common ..
|
||||||
LOGICAL LERR, OK
|
LOGICAL LERR, OK
|
||||||
|
|
|
@ -61,17 +61,17 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2015
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup single_eig
|
*> \ingroup single_eig
|
||||||
*
|
*
|
||||||
* =====================================================================
|
* =====================================================================
|
||||||
SUBROUTINE SERRED( PATH, NUNIT )
|
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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2015
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
CHARACTER*3 PATH
|
CHARACTER*3 PATH
|
||||||
|
@ -475,7 +475,7 @@
|
||||||
CALL SGESVDX( 'V', 'N', 'A', 2, 2, A, 2, ZERO, ZERO,
|
CALL SGESVDX( 'V', 'N', 'A', 2, 2, A, 2, ZERO, ZERO,
|
||||||
$ 0, 0, NS, S, U, 1, VT, 1, W, 1, IW, INFO )
|
$ 0, 0, NS, S, U, 1, VT, 1, W, 1, IW, INFO )
|
||||||
CALL CHKXER( 'SGESVDX', INFOT, NOUT, LERR, OK )
|
CALL CHKXER( 'SGESVDX', INFOT, NOUT, LERR, OK )
|
||||||
INFOT = 16
|
INFOT = 17
|
||||||
CALL SGESVDX( 'N', 'V', 'A', 2, 2, A, 2, ZERO, ZERO,
|
CALL SGESVDX( 'N', 'V', 'A', 2, 2, A, 2, ZERO, ZERO,
|
||||||
$ 0, 0, NS, S, U, 1, VT, 1, W, 1, IW, INFO )
|
$ 0, 0, NS, S, U, 1, VT, 1, W, 1, IW, INFO )
|
||||||
CALL CHKXER( 'SGESVDX', INFOT, NOUT, LERR, OK )
|
CALL CHKXER( 'SGESVDX', INFOT, NOUT, LERR, OK )
|
||||||
|
|
|
@ -50,17 +50,17 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2015
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup single_eig
|
*> \ingroup single_eig
|
||||||
*
|
*
|
||||||
* =====================================================================
|
* =====================================================================
|
||||||
SUBROUTINE SERRGG( PATH, NUNIT )
|
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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2015
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
CHARACTER*3 PATH
|
CHARACTER*3 PATH
|
||||||
|
@ -78,12 +78,12 @@
|
||||||
* .. Local Scalars ..
|
* .. Local Scalars ..
|
||||||
CHARACTER*2 C2
|
CHARACTER*2 C2
|
||||||
INTEGER DUMMYK, DUMMYL, I, IFST, ILO, IHI, ILST, INFO,
|
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
|
REAL ANRM, BNRM, DIF, SCALE, TOLA, TOLB
|
||||||
* ..
|
* ..
|
||||||
* .. Local Arrays ..
|
* .. Local Arrays ..
|
||||||
LOGICAL BW( NMAX ), SEL( NMAX )
|
LOGICAL BW( NMAX ), SEL( NMAX )
|
||||||
INTEGER IW( NMAX )
|
INTEGER IW( NMAX ), IDUM(NMAX)
|
||||||
REAL A( NMAX, NMAX ), B( NMAX, NMAX ), LS( NMAX ),
|
REAL A( NMAX, NMAX ), B( NMAX, NMAX ), LS( NMAX ),
|
||||||
$ Q( NMAX, NMAX ), R1( NMAX ), R2( NMAX ),
|
$ Q( NMAX, NMAX ), R1( NMAX ), R2( NMAX ),
|
||||||
$ R3( NMAX ), RCE( 2 ), RCV( 2 ), RS( NMAX ),
|
$ R3( NMAX ), RCE( 2 ), RCV( 2 ), RS( NMAX ),
|
||||||
|
@ -305,47 +305,47 @@
|
||||||
SRNAMT = 'SGGSVD3'
|
SRNAMT = 'SGGSVD3'
|
||||||
INFOT = 1
|
INFOT = 1
|
||||||
CALL SGGSVD3( '/', 'N', 'N', 0, 0, 0, DUMMYK, DUMMYL, A, 1, B,
|
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 )
|
CALL CHKXER( 'SGGSVD3', INFOT, NOUT, LERR, OK )
|
||||||
INFOT = 2
|
INFOT = 2
|
||||||
CALL SGGSVD3( 'N', '/', 'N', 0, 0, 0, DUMMYK, DUMMYL, A, 1, B,
|
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 )
|
CALL CHKXER( 'SGGSVD3', INFOT, NOUT, LERR, OK )
|
||||||
INFOT = 3
|
INFOT = 3
|
||||||
CALL SGGSVD3( 'N', 'N', '/', 0, 0, 0, DUMMYK, DUMMYL, A, 1, B,
|
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 )
|
CALL CHKXER( 'SGGSVD3', INFOT, NOUT, LERR, OK )
|
||||||
INFOT = 4
|
INFOT = 4
|
||||||
CALL SGGSVD3( 'N', 'N', 'N', -1, 0, 0, DUMMYK, DUMMYL, A, 1, B,
|
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 )
|
CALL CHKXER( 'SGGSVD3', INFOT, NOUT, LERR, OK )
|
||||||
INFOT = 5
|
INFOT = 5
|
||||||
CALL SGGSVD3( 'N', 'N', 'N', 0, -1, 0, DUMMYK, DUMMYL, A, 1, B,
|
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 )
|
CALL CHKXER( 'SGGSVD3', INFOT, NOUT, LERR, OK )
|
||||||
INFOT = 6
|
INFOT = 6
|
||||||
CALL SGGSVD3( 'N', 'N', 'N', 0, 0, -1, DUMMYK, DUMMYL, A, 1, B,
|
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 )
|
CALL CHKXER( 'SGGSVD3', INFOT, NOUT, LERR, OK )
|
||||||
INFOT = 10
|
INFOT = 10
|
||||||
CALL SGGSVD3( 'N', 'N', 'N', 2, 1, 1, DUMMYK, DUMMYL, A, 1, B,
|
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 )
|
CALL CHKXER( 'SGGSVD3', INFOT, NOUT, LERR, OK )
|
||||||
INFOT = 12
|
INFOT = 12
|
||||||
CALL SGGSVD3( 'N', 'N', 'N', 1, 1, 2, DUMMYK, DUMMYL, A, 1, B,
|
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 )
|
CALL CHKXER( 'SGGSVD3', INFOT, NOUT, LERR, OK )
|
||||||
INFOT = 16
|
INFOT = 16
|
||||||
CALL SGGSVD3( 'U', 'N', 'N', 2, 2, 2, DUMMYK, DUMMYL, A, 2, B,
|
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 )
|
CALL CHKXER( 'SGGSVD3', INFOT, NOUT, LERR, OK )
|
||||||
INFOT = 18
|
INFOT = 18
|
||||||
CALL SGGSVD3( 'N', 'V', 'N', 1, 1, 2, DUMMYK, DUMMYL, A, 1, B,
|
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 )
|
CALL CHKXER( 'SGGSVD3', INFOT, NOUT, LERR, OK )
|
||||||
INFOT = 20
|
INFOT = 20
|
||||||
CALL SGGSVD3( 'N', 'N', 'Q', 1, 2, 1, DUMMYK, DUMMYL, A, 1, B,
|
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 )
|
CALL CHKXER( 'SGGSVD3', INFOT, NOUT, LERR, OK )
|
||||||
NT = NT + 11
|
NT = NT + 11
|
||||||
*
|
*
|
||||||
|
|
|
@ -26,7 +26,7 @@
|
||||||
*> .FALSE..
|
*> .FALSE..
|
||||||
*>
|
*>
|
||||||
*> It is used by the test routine SDRGES to test whether the driver
|
*> It is used by the test routine SDRGES to test whether the driver
|
||||||
*> routine SGGES succesfully sorts eigenvalues.
|
*> routine SGGES successfully sorts eigenvalues.
|
||||||
*> \endverbatim
|
*> \endverbatim
|
||||||
*
|
*
|
||||||
* Arguments:
|
* Arguments:
|
||||||
|
@ -61,17 +61,17 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2011
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup single_eig
|
*> \ingroup single_eig
|
||||||
*
|
*
|
||||||
* =====================================================================
|
* =====================================================================
|
||||||
LOGICAL FUNCTION SLCTES( ZR, ZI, D )
|
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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2011
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
REAL D, ZI, ZR
|
REAL D, ZI, ZR
|
||||||
|
|
|
@ -22,8 +22,8 @@
|
||||||
*>
|
*>
|
||||||
*> SSLECT returns .TRUE. if the eigenvalue ZR+sqrt(-1)*ZI is to be
|
*> SSLECT returns .TRUE. if the eigenvalue ZR+sqrt(-1)*ZI is to be
|
||||||
*> selected, and otherwise it returns .FALSE.
|
*> selected, and otherwise it returns .FALSE.
|
||||||
*> It is used by SCHK41 to test if SGEES succesfully sorts eigenvalues,
|
*> It is used by SCHK41 to test if SGEES successfully sorts eigenvalues,
|
||||||
*> and by SCHK43 to test if SGEESX succesfully sorts eigenvalues.
|
*> and by SCHK43 to test if SGEESX successfully sorts eigenvalues.
|
||||||
*>
|
*>
|
||||||
*> The common block /SSLCT/ controls how eigenvalues are selected.
|
*> The common block /SSLCT/ controls how eigenvalues are selected.
|
||||||
*> If SELOPT = 0, then SSLECT return .TRUE. when ZR is less than zero,
|
*> If SELOPT = 0, then SSLECT return .TRUE. when ZR is less than zero,
|
||||||
|
@ -55,17 +55,17 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2011
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup single_eig
|
*> \ingroup single_eig
|
||||||
*
|
*
|
||||||
* =====================================================================
|
* =====================================================================
|
||||||
LOGICAL FUNCTION SSLECT( ZR, ZI )
|
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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2011
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
REAL ZI, ZR
|
REAL ZI, ZR
|
||||||
|
|
|
@ -405,7 +405,7 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2011
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup complex16_eig
|
*> \ingroup complex16_eig
|
||||||
*
|
*
|
||||||
|
@ -415,10 +415,10 @@
|
||||||
$ Y, Z, Q, LDQ, PT, LDPT, U, VT, WORK, LWORK,
|
$ Y, Z, Q, LDQ, PT, LDPT, U, VT, WORK, LWORK,
|
||||||
$ RWORK, NOUT, INFO )
|
$ 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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2011
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
INTEGER INFO, LDA, LDPT, LDQ, LDX, LWORK, NOUT, NRHS,
|
INTEGER INFO, LDA, LDPT, LDQ, LDX, LWORK, NOUT, NRHS,
|
||||||
|
|
|
@ -1027,17 +1027,17 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2015
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup complex16_eig
|
*> \ingroup complex16_eig
|
||||||
*
|
*
|
||||||
* =====================================================================
|
* =====================================================================
|
||||||
PROGRAM ZCHKEE
|
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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2015
|
* June 2016
|
||||||
*
|
*
|
||||||
* =====================================================================
|
* =====================================================================
|
||||||
*
|
*
|
||||||
|
@ -2098,6 +2098,7 @@
|
||||||
MAXTYP = 26
|
MAXTYP = 26
|
||||||
NTYPES = MIN( MAXTYP, NTYPES )
|
NTYPES = MIN( MAXTYP, NTYPES )
|
||||||
CALL ALAREQ( C3, NTYPES, DOTYPE, MAXTYP, NIN, NOUT )
|
CALL ALAREQ( C3, NTYPES, DOTYPE, MAXTYP, NIN, NOUT )
|
||||||
|
CALL XLAENV(1,1)
|
||||||
IF( TSTCHK .AND. TSTERR )
|
IF( TSTCHK .AND. TSTERR )
|
||||||
$ CALL ZERRGG( C3, NOUT )
|
$ CALL ZERRGG( C3, NOUT )
|
||||||
DO 350 I = 1, NPARMS
|
DO 350 I = 1, NPARMS
|
||||||
|
@ -2222,6 +2223,7 @@
|
||||||
*
|
*
|
||||||
* Blocked version
|
* Blocked version
|
||||||
*
|
*
|
||||||
|
CALL XLAENV(16,2)
|
||||||
CALL ZDRGEV3( NN, NVAL, MAXTYP, DOTYPE, ISEED, THRESH, NOUT,
|
CALL ZDRGEV3( NN, NVAL, MAXTYP, DOTYPE, ISEED, THRESH, NOUT,
|
||||||
$ A( 1, 1 ), NMAX, A( 1, 2 ), A( 1, 3 ),
|
$ A( 1, 1 ), NMAX, A( 1, 2 ), A( 1, 3 ),
|
||||||
$ A( 1, 4 ), A( 1, 7 ), NMAX, A( 1, 8 ),
|
$ A( 1, 4 ), A( 1, 7 ), NMAX, A( 1, 8 ),
|
||||||
|
@ -2347,6 +2349,7 @@
|
||||||
* GSV: Generalized Singular Value Decomposition
|
* GSV: Generalized Singular Value Decomposition
|
||||||
* ----------------------------------------------
|
* ----------------------------------------------
|
||||||
*
|
*
|
||||||
|
CALL XLAENV(1,1)
|
||||||
IF( TSTERR )
|
IF( TSTERR )
|
||||||
$ CALL ZERRGG( 'GSV', NOUT )
|
$ CALL ZERRGG( 'GSV', NOUT )
|
||||||
CALL ZCKGSV( NN, MVAL, PVAL, NVAL, NTYPES, ISEED, THRESH, NMAX,
|
CALL ZCKGSV( NN, MVAL, PVAL, NVAL, NTYPES, ISEED, THRESH, NMAX,
|
||||||
|
|
|
@ -68,7 +68,7 @@
|
||||||
*> and each type of matrix, one matrix will be generated and used
|
*> and each type of matrix, one matrix will be generated and used
|
||||||
*> to test the nonsymmetric eigenroutines. For each matrix, 13
|
*> to test the nonsymmetric eigenroutines. For each matrix, 13
|
||||||
*> tests will be performed. The first twelve "test ratios" should be
|
*> 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
|
*> H
|
||||||
*> (1) | A - U H V | / ( |A| n ulp )
|
*> (1) | A - U H V | / ( |A| n ulp )
|
||||||
|
@ -128,7 +128,7 @@
|
||||||
*> |beta(Q,Z computed) - beta(Q,Z not computed)|/|P| ) / ulp
|
*> |beta(Q,Z computed) - beta(Q,Z not computed)|/|P| ) / ulp
|
||||||
*>
|
*>
|
||||||
*> In addition, the normalization of L and R are checked, and compared
|
*> In addition, the normalization of L and R are checked, and compared
|
||||||
*> with the threshhold THRSHN.
|
*> with the threshold THRSHN.
|
||||||
*>
|
*>
|
||||||
*> Test Matrices
|
*> Test Matrices
|
||||||
*> ---- --------
|
*> ---- --------
|
||||||
|
@ -298,7 +298,7 @@
|
||||||
*> \param[in] THRSHN
|
*> \param[in] THRSHN
|
||||||
*> \verbatim
|
*> \verbatim
|
||||||
*> THRSHN is DOUBLE PRECISION
|
*> 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
|
*> If the normalization of any eigenvector differs from 1 by
|
||||||
*> more than THRSHN*ulp, then a special error message will be
|
*> more than THRSHN*ulp, then a special error message will be
|
||||||
*> printed. (This is handled separately from the other tests,
|
*> printed. (This is handled separately from the other tests,
|
||||||
|
@ -492,7 +492,7 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2011
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup complex16_eig
|
*> \ingroup complex16_eig
|
||||||
*
|
*
|
||||||
|
@ -503,10 +503,10 @@
|
||||||
$ ALPHA3, BETA3, EVECTL, EVECTR, WORK, LWORK,
|
$ ALPHA3, BETA3, EVECTL, EVECTR, WORK, LWORK,
|
||||||
$ RWORK, LLWORK, RESULT, INFO )
|
$ 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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2011
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
LOGICAL TSTDIF
|
LOGICAL TSTDIF
|
||||||
|
|
|
@ -50,7 +50,7 @@
|
||||||
*> number of matrix "TYPES" are specified. For each size ("N")
|
*> 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 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
|
*> 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).
|
*> the tests (5), (11) and (13).
|
||||||
*>
|
*>
|
||||||
*>
|
*>
|
||||||
|
@ -372,7 +372,7 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2011
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup complex16_eig
|
*> \ingroup complex16_eig
|
||||||
*
|
*
|
||||||
|
@ -381,10 +381,10 @@
|
||||||
$ NOUNIT, A, LDA, B, S, T, Q, LDQ, Z, ALPHA,
|
$ NOUNIT, A, LDA, B, S, T, Q, LDQ, Z, ALPHA,
|
||||||
$ BETA, WORK, LWORK, RWORK, RESULT, BWORK, INFO )
|
$ 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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2011
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
INTEGER INFO, LDA, LDQ, LWORK, NOUNIT, NSIZES, NTYPES
|
INTEGER INFO, LDA, LDQ, LWORK, NOUNIT, NSIZES, NTYPES
|
||||||
|
|
|
@ -50,7 +50,7 @@
|
||||||
*> number of matrix "TYPES" are specified. For each size ("N")
|
*> 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 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
|
*> 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).
|
*> the tests (5), (11) and (13).
|
||||||
*>
|
*>
|
||||||
*>
|
*>
|
||||||
|
@ -382,7 +382,7 @@
|
||||||
$ BETA, WORK, LWORK, RWORK, RESULT, BWORK,
|
$ BETA, WORK, LWORK, RWORK, RESULT, BWORK,
|
||||||
$ INFO )
|
$ 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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* February 2015
|
* February 2015
|
||||||
|
|
|
@ -55,7 +55,7 @@
|
||||||
*> number of matrix "types" are specified. For each size ("n")
|
*> 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 each type of matrix, a pair of matrices (A, B) will be generated
|
||||||
*> and used for testing. For each matrix pair, the following tests
|
*> 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:
|
*> Results from ZGGEV:
|
||||||
*>
|
*>
|
||||||
|
@ -389,7 +389,7 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2015
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup complex16_eig
|
*> \ingroup complex16_eig
|
||||||
*
|
*
|
||||||
|
@ -399,10 +399,10 @@
|
||||||
$ ALPHA, BETA, ALPHA1, BETA1, WORK, LWORK, RWORK,
|
$ ALPHA, BETA, ALPHA1, BETA1, WORK, LWORK, RWORK,
|
||||||
$ RESULT, INFO )
|
$ 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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2015
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
INTEGER INFO, LDA, LDQ, LDQE, LWORK, NOUNIT, NSIZES,
|
INTEGER INFO, LDA, LDQ, LDQE, LWORK, NOUNIT, NSIZES,
|
||||||
|
|
|
@ -55,7 +55,7 @@
|
||||||
*> number of matrix "types" are specified. For each size ("n")
|
*> 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 each type of matrix, a pair of matrices (A, B) will be generated
|
||||||
*> and used for testing. For each matrix pair, the following tests
|
*> 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:
|
*> Results from ZGGEV3:
|
||||||
*>
|
*>
|
||||||
|
@ -399,7 +399,7 @@
|
||||||
$ ALPHA, BETA, ALPHA1, BETA1, WORK, LWORK,
|
$ ALPHA, BETA, ALPHA1, BETA1, WORK, LWORK,
|
||||||
$ RWORK, RESULT, INFO )
|
$ 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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* February 2015
|
* February 2015
|
||||||
|
|
|
@ -59,7 +59,7 @@
|
||||||
*> (need more details on what kind of read-in data are needed).
|
*> (need more details on what kind of read-in data are needed).
|
||||||
*>
|
*>
|
||||||
*> For each matrix pair, the following tests will be performed and
|
*> 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 )
|
*> (1) | A - Q S Z' | / ( |A| n ulp )
|
||||||
*>
|
*>
|
||||||
|
@ -340,7 +340,7 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2011
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup complex16_eig
|
*> \ingroup complex16_eig
|
||||||
*
|
*
|
||||||
|
@ -349,10 +349,10 @@
|
||||||
$ BI, Z, Q, ALPHA, BETA, C, LDC, S, WORK, LWORK,
|
$ BI, Z, Q, ALPHA, BETA, C, LDC, S, WORK, LWORK,
|
||||||
$ RWORK, IWORK, LIWORK, BWORK, INFO )
|
$ 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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2011
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
INTEGER INFO, LDA, LDC, LIWORK, LWORK, NCMAX, NIN,
|
INTEGER INFO, LDA, LDC, LIWORK, LWORK, NCMAX, NIN,
|
||||||
|
|
|
@ -49,7 +49,7 @@
|
||||||
*> corresponding the first and last eigenvalues are also know
|
*> corresponding the first and last eigenvalues are also know
|
||||||
*> ``exactly'' (see ZLATM6).
|
*> ``exactly'' (see ZLATM6).
|
||||||
*> For each matrix pair, the following tests will be performed and
|
*> 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
|
*> (1) max over all left eigenvalue/-vector pairs (beta/alpha,l) of
|
||||||
*>
|
*>
|
||||||
|
@ -287,7 +287,7 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2011
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup complex16_eig
|
*> \ingroup complex16_eig
|
||||||
*
|
*
|
||||||
|
@ -297,10 +297,10 @@
|
||||||
$ S, DTRU, DIF, DIFTRU, WORK, LWORK, RWORK,
|
$ S, DTRU, DIF, DIFTRU, WORK, LWORK, RWORK,
|
||||||
$ IWORK, LIWORK, RESULT, BWORK, INFO )
|
$ 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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2011
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
INTEGER IHI, ILO, INFO, LDA, LIWORK, LWORK, NIN, NOUT,
|
INTEGER IHI, ILO, INFO, LDA, LIWORK, LWORK, NIN, NOUT,
|
||||||
|
|
|
@ -379,7 +379,7 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2015
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup complex16_eig
|
*> \ingroup complex16_eig
|
||||||
*
|
*
|
||||||
|
@ -389,10 +389,10 @@
|
||||||
$ SSAV, E, WORK, LWORK, RWORK, IWORK, NOUNIT,
|
$ SSAV, E, WORK, LWORK, RWORK, IWORK, NOUNIT,
|
||||||
$ INFO )
|
$ 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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2015
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
INTEGER INFO, LDA, LDU, LDVT, LWORK, NOUNIT, NSIZES,
|
INTEGER INFO, LDA, LDU, LDVT, LWORK, NOUNIT, NSIZES,
|
||||||
|
@ -423,11 +423,11 @@
|
||||||
* .. Local Scalars ..
|
* .. Local Scalars ..
|
||||||
LOGICAL BADMM, BADNN
|
LOGICAL BADMM, BADNN
|
||||||
CHARACTER JOBQ, JOBU, JOBVT, RANGE
|
CHARACTER JOBQ, JOBU, JOBVT, RANGE
|
||||||
INTEGER I, IINFO, IJQ, IJU, IJVT, IL, IU, ITEMP, IWSPC,
|
INTEGER I, IINFO, IJQ, IJU, IJVT, IL, IU, ITEMP,
|
||||||
$ IWTMP, J, JSIZE, JTYPE, LSWORK, M, MINWRK,
|
$ IWSPC, IWTMP, J, JSIZE, JTYPE, LSWORK, M,
|
||||||
$ MMAX, MNMAX, MNMIN, MTYPES, N, NERRS, NFAIL,
|
$ MINWRK, MMAX, MNMAX, MNMIN, MTYPES, N,
|
||||||
$ NMAX, NS, NSI, NSV, NTEST, NTESTF, NTESTT,
|
$ NERRS, NFAIL, NMAX, NS, NSI, NSV, NTEST,
|
||||||
$ LRWORK
|
$ NTESTF, NTESTT, LRWORK
|
||||||
DOUBLE PRECISION ANORM, DIF, DIV, OVFL, RTUNFL, ULP, ULPINV,
|
DOUBLE PRECISION ANORM, DIF, DIV, OVFL, RTUNFL, ULP, ULPINV,
|
||||||
$ UNFL, VL, VU
|
$ UNFL, VL, VU
|
||||||
* ..
|
* ..
|
||||||
|
@ -441,9 +441,9 @@
|
||||||
EXTERNAL DLAMCH, DLARND
|
EXTERNAL DLAMCH, DLARND
|
||||||
* ..
|
* ..
|
||||||
* .. External Subroutines ..
|
* .. External Subroutines ..
|
||||||
EXTERNAL ALASVM, XERBLA, ZBDT01, ZBDT05, ZGESDD, ZGESVD,
|
EXTERNAL ALASVM, XERBLA, ZBDT01, ZBDT05, ZGESDD,
|
||||||
$ ZGESVJ, ZGEJSV, ZGESVDX, ZLACPY, ZLASET, ZLATMS,
|
$ ZGESVD, ZGESVJ, ZGEJSV, ZGESVDX, ZLACPY,
|
||||||
$ ZUNT01, ZUNT03
|
$ ZLASET, ZLATMS, ZUNT01, ZUNT03
|
||||||
* ..
|
* ..
|
||||||
* .. Intrinsic Functions ..
|
* .. Intrinsic Functions ..
|
||||||
INTRINSIC ABS, DBLE, MAX, MIN
|
INTRINSIC ABS, DBLE, MAX, MIN
|
||||||
|
|
|
@ -369,7 +369,7 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2011
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup complex16_eig
|
*> \ingroup complex16_eig
|
||||||
*
|
*
|
||||||
|
@ -378,10 +378,10 @@
|
||||||
$ NOUNIT, A, LDA, H, HT, W, WT, VS, LDVS, RESULT,
|
$ NOUNIT, A, LDA, H, HT, W, WT, VS, LDVS, RESULT,
|
||||||
$ WORK, NWORK, RWORK, IWORK, BWORK, INFO )
|
$ 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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2011
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
INTEGER INFO, LDA, LDVS, NOUNIT, NSIZES, NTYPES, NWORK
|
INTEGER INFO, LDA, LDVS, NOUNIT, NSIZES, NTYPES, NWORK
|
||||||
|
@ -901,7 +901,7 @@
|
||||||
$ ' 1/ulp otherwise', /
|
$ ' 1/ulp otherwise', /
|
||||||
$ ' 12 = 0 if W same no matter if VS computed (sort),',
|
$ ' 12 = 0 if W same no matter if VS computed (sort),',
|
||||||
$ ' 1/ulp otherwise', /
|
$ ' 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, ',' ),
|
9993 FORMAT( ' N=', I5, ', IWK=', I2, ', seed=', 4( I4, ',' ),
|
||||||
$ ' type ', I2, ', test(', I2, ')=', G10.3 )
|
$ ' type ', I2, ', test(', I2, ')=', G10.3 )
|
||||||
9992 FORMAT( ' ZDRVES: ', A, ' returned INFO=', I6, '.', / 9X, 'N=',
|
9992 FORMAT( ' ZDRVES: ', A, ' returned INFO=', I6, '.', / 9X, 'N=',
|
||||||
|
|
|
@ -425,7 +425,7 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2011
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup complex16_eig
|
*> \ingroup complex16_eig
|
||||||
*
|
*
|
||||||
|
@ -435,10 +435,10 @@
|
||||||
$ LDVS, VS1, RESULT, WORK, LWORK, RWORK, BWORK,
|
$ LDVS, VS1, RESULT, WORK, LWORK, RWORK, BWORK,
|
||||||
$ INFO )
|
$ 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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2011
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
INTEGER INFO, LDA, LDVS, LWORK, NIUNIT, NOUNIT, NSIZES,
|
INTEGER INFO, LDA, LDVS, LWORK, NIUNIT, NOUNIT, NSIZES,
|
||||||
|
@ -914,7 +914,7 @@
|
||||||
$ ' 1/ulp otherwise', /
|
$ ' 1/ulp otherwise', /
|
||||||
$ ' 12 = 0 if W same no matter what else computed ',
|
$ ' 12 = 0 if W same no matter what else computed ',
|
||||||
$ '(sort), 1/ulp otherwise', /
|
$ '(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,',
|
$ / ' 14 = 0 if RCONDE same no matter what else computed,',
|
||||||
$ ' 1/ulp otherwise', /
|
$ ' 1/ulp otherwise', /
|
||||||
$ ' 15 = 0 if RCONDv same no matter what else computed,',
|
$ ' 15 = 0 if RCONDv same no matter what else computed,',
|
||||||
|
|
|
@ -446,7 +446,7 @@
|
||||||
*> \verbatim
|
*> \verbatim
|
||||||
*> INFO is INTEGER
|
*> INFO is INTEGER
|
||||||
*> If 0, then successful exit.
|
*> 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
|
*> If >0, ZLATMR, CLATMS, CLATME or ZGET23 returned an error
|
||||||
*> code, and INFO is its absolute value.
|
*> code, and INFO is its absolute value.
|
||||||
*>
|
*>
|
||||||
|
@ -485,7 +485,7 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2011
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup complex16_eig
|
*> \ingroup complex16_eig
|
||||||
*
|
*
|
||||||
|
@ -496,10 +496,10 @@
|
||||||
$ RCONDE, RCNDE1, RCDEIN, SCALE, SCALE1, RESULT,
|
$ RCONDE, RCNDE1, RCDEIN, SCALE, SCALE1, RESULT,
|
||||||
$ WORK, NWORK, RWORK, INFO )
|
$ 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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2011
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
INTEGER INFO, LDA, LDLRE, LDVL, LDVR, NIUNIT, NOUNIT,
|
INTEGER INFO, LDA, LDLRE, LDVL, LDVR, NIUNIT, NOUNIT,
|
||||||
|
|
|
@ -61,17 +61,17 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2015
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup complex16_eig
|
*> \ingroup complex16_eig
|
||||||
*
|
*
|
||||||
* =====================================================================
|
* =====================================================================
|
||||||
SUBROUTINE ZERRED( PATH, NUNIT )
|
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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2015
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
CHARACTER*3 PATH
|
CHARACTER*3 PATH
|
||||||
|
@ -484,7 +484,7 @@
|
||||||
CALL ZGESVDX( 'V', 'N', 'A', 2, 2, A, 2, ZERO, ZERO,
|
CALL ZGESVDX( 'V', 'N', 'A', 2, 2, A, 2, ZERO, ZERO,
|
||||||
$ 0, 0, NS, S, U, 1, VT, 1, W, 1, RW, IW, INFO )
|
$ 0, 0, NS, S, U, 1, VT, 1, W, 1, RW, IW, INFO )
|
||||||
CALL CHKXER( 'ZGESVDX', INFOT, NOUT, LERR, OK )
|
CALL CHKXER( 'ZGESVDX', INFOT, NOUT, LERR, OK )
|
||||||
INFOT = 16
|
INFOT = 17
|
||||||
CALL ZGESVDX( 'N', 'V', 'A', 2, 2, A, 2, ZERO, ZERO,
|
CALL ZGESVDX( 'N', 'V', 'A', 2, 2, A, 2, ZERO, ZERO,
|
||||||
$ 0, 0, NS, S, U, 1, VT, 1, W, 1, RW, IW, INFO )
|
$ 0, 0, NS, S, U, 1, VT, 1, W, 1, RW, IW, INFO )
|
||||||
CALL CHKXER( 'ZGESVDX', INFOT, NOUT, LERR, OK )
|
CALL CHKXER( 'ZGESVDX', INFOT, NOUT, LERR, OK )
|
||||||
|
|
|
@ -50,17 +50,17 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2015
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup complex16_eig
|
*> \ingroup complex16_eig
|
||||||
*
|
*
|
||||||
* =====================================================================
|
* =====================================================================
|
||||||
SUBROUTINE ZERRGG( PATH, NUNIT )
|
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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2015
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
CHARACTER*3 PATH
|
CHARACTER*3 PATH
|
||||||
|
@ -83,7 +83,7 @@
|
||||||
* ..
|
* ..
|
||||||
* .. Local Arrays ..
|
* .. Local Arrays ..
|
||||||
LOGICAL BW( NMAX ), SEL( NMAX )
|
LOGICAL BW( NMAX ), SEL( NMAX )
|
||||||
INTEGER IW( LW )
|
INTEGER IW( LW ), IDUM(NMAX)
|
||||||
DOUBLE PRECISION LS( NMAX ), R1( NMAX ), R2( NMAX ),
|
DOUBLE PRECISION LS( NMAX ), R1( NMAX ), R2( NMAX ),
|
||||||
$ RCE( NMAX ), RCV( NMAX ), RS( NMAX ), RW( LW )
|
$ RCE( NMAX ), RCV( NMAX ), RS( NMAX ), RW( LW )
|
||||||
COMPLEX*16 A( NMAX, NMAX ), ALPHA( NMAX ),
|
COMPLEX*16 A( NMAX, NMAX ), ALPHA( NMAX ),
|
||||||
|
@ -306,57 +306,57 @@
|
||||||
SRNAMT = 'ZGGSVD3'
|
SRNAMT = 'ZGGSVD3'
|
||||||
INFOT = 1
|
INFOT = 1
|
||||||
CALL ZGGSVD3( '/', 'N', 'N', 0, 0, 0, DUMMYK, DUMMYL, A, 1, B,
|
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 )
|
$ INFO )
|
||||||
CALL CHKXER( 'ZGGSVD3', INFOT, NOUT, LERR, OK )
|
CALL CHKXER( 'ZGGSVD3', INFOT, NOUT, LERR, OK )
|
||||||
INFOT = 2
|
INFOT = 2
|
||||||
CALL ZGGSVD3( 'N', '/', 'N', 0, 0, 0, DUMMYK, DUMMYL, A, 1, B,
|
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 )
|
$ INFO )
|
||||||
CALL CHKXER( 'ZGGSVD3', INFOT, NOUT, LERR, OK )
|
CALL CHKXER( 'ZGGSVD3', INFOT, NOUT, LERR, OK )
|
||||||
INFOT = 3
|
INFOT = 3
|
||||||
CALL ZGGSVD3( 'N', 'N', '/', 0, 0, 0, DUMMYK, DUMMYL, A, 1, B,
|
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 )
|
$ INFO )
|
||||||
CALL CHKXER( 'ZGGSVD3', INFOT, NOUT, LERR, OK )
|
CALL CHKXER( 'ZGGSVD3', INFOT, NOUT, LERR, OK )
|
||||||
INFOT = 4
|
INFOT = 4
|
||||||
CALL ZGGSVD3( 'N', 'N', 'N', -1, 0, 0, DUMMYK, DUMMYL, A, 1, B,
|
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 )
|
$ INFO )
|
||||||
CALL CHKXER( 'ZGGSVD3', INFOT, NOUT, LERR, OK )
|
CALL CHKXER( 'ZGGSVD3', INFOT, NOUT, LERR, OK )
|
||||||
INFOT = 5
|
INFOT = 5
|
||||||
CALL ZGGSVD3( 'N', 'N', 'N', 0, -1, 0, DUMMYK, DUMMYL, A, 1, B,
|
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 )
|
$ INFO )
|
||||||
CALL CHKXER( 'ZGGSVD3', INFOT, NOUT, LERR, OK )
|
CALL CHKXER( 'ZGGSVD3', INFOT, NOUT, LERR, OK )
|
||||||
INFOT = 6
|
INFOT = 6
|
||||||
CALL ZGGSVD3( 'N', 'N', 'N', 0, 0, -1, DUMMYK, DUMMYL, A, 1, B,
|
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 )
|
$ INFO )
|
||||||
CALL CHKXER( 'ZGGSVD3', INFOT, NOUT, LERR, OK )
|
CALL CHKXER( 'ZGGSVD3', INFOT, NOUT, LERR, OK )
|
||||||
INFOT = 10
|
INFOT = 10
|
||||||
CALL ZGGSVD3( 'N', 'N', 'N', 2, 1, 1, DUMMYK, DUMMYL, A, 1, B,
|
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 )
|
$ INFO )
|
||||||
CALL CHKXER( 'ZGGSVD3', INFOT, NOUT, LERR, OK )
|
CALL CHKXER( 'ZGGSVD3', INFOT, NOUT, LERR, OK )
|
||||||
INFOT = 12
|
INFOT = 12
|
||||||
CALL ZGGSVD3( 'N', 'N', 'N', 1, 1, 2, DUMMYK, DUMMYL, A, 1, B,
|
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 )
|
$ INFO )
|
||||||
CALL CHKXER( 'ZGGSVD3', INFOT, NOUT, LERR, OK )
|
CALL CHKXER( 'ZGGSVD3', INFOT, NOUT, LERR, OK )
|
||||||
INFOT = 16
|
INFOT = 16
|
||||||
CALL ZGGSVD3( 'U', 'N', 'N', 2, 2, 2, DUMMYK, DUMMYL, A, 2, B,
|
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 )
|
$ INFO )
|
||||||
CALL CHKXER( 'ZGGSVD3', INFOT, NOUT, LERR, OK )
|
CALL CHKXER( 'ZGGSVD3', INFOT, NOUT, LERR, OK )
|
||||||
INFOT = 18
|
INFOT = 18
|
||||||
CALL ZGGSVD3( 'N', 'V', 'N', 2, 2, 2, DUMMYK, DUMMYL, A, 2, B,
|
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 )
|
$ INFO )
|
||||||
CALL CHKXER( 'ZGGSVD3', INFOT, NOUT, LERR, OK )
|
CALL CHKXER( 'ZGGSVD3', INFOT, NOUT, LERR, OK )
|
||||||
INFOT = 20
|
INFOT = 20
|
||||||
CALL ZGGSVD3( 'N', 'N', 'Q', 2, 2, 2, DUMMYK, DUMMYL, A, 2, B,
|
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 )
|
$ INFO )
|
||||||
CALL CHKXER( 'ZGGSVD3', INFOT, NOUT, LERR, OK )
|
CALL CHKXER( 'ZGGSVD3', INFOT, NOUT, LERR, OK )
|
||||||
NT = NT + 11
|
NT = NT + 11
|
||||||
|
@ -573,56 +573,56 @@
|
||||||
INFOT = 7
|
INFOT = 7
|
||||||
CALL ZUNCSD( 'Y', 'Y', 'Y', 'Y', 'N', 'N',
|
CALL ZUNCSD( 'Y', 'Y', 'Y', 'Y', 'N', 'N',
|
||||||
$ -1, 0, 0, A, 1, A,
|
$ -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,
|
$ A, 1, A, 1, A, 1, A,
|
||||||
$ 1, W, LW, RW, LW, IW, INFO )
|
$ 1, W, LW, RW, LW, IW, INFO )
|
||||||
CALL CHKXER( 'ZUNCSD', INFOT, NOUT, LERR, OK )
|
CALL CHKXER( 'ZUNCSD', INFOT, NOUT, LERR, OK )
|
||||||
INFOT = 8
|
INFOT = 8
|
||||||
CALL ZUNCSD( 'Y', 'Y', 'Y', 'Y', 'N', 'N',
|
CALL ZUNCSD( 'Y', 'Y', 'Y', 'Y', 'N', 'N',
|
||||||
$ 1, -1, 0, A, 1, A,
|
$ 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,
|
$ A, 1, A, 1, A, 1, A,
|
||||||
$ 1, W, LW, RW, LW, IW, INFO )
|
$ 1, W, LW, RW, LW, IW, INFO )
|
||||||
CALL CHKXER( 'ZUNCSD', INFOT, NOUT, LERR, OK )
|
CALL CHKXER( 'ZUNCSD', INFOT, NOUT, LERR, OK )
|
||||||
INFOT = 9
|
INFOT = 9
|
||||||
CALL ZUNCSD( 'Y', 'Y', 'Y', 'Y', 'N', 'N',
|
CALL ZUNCSD( 'Y', 'Y', 'Y', 'Y', 'N', 'N',
|
||||||
$ 1, 1, -1, A, 1, A,
|
$ 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,
|
$ A, 1, A, 1, A, 1, A,
|
||||||
$ 1, W, LW, RW, LW, IW, INFO )
|
$ 1, W, LW, RW, LW, IW, INFO )
|
||||||
CALL CHKXER( 'ZUNCSD', INFOT, NOUT, LERR, OK )
|
CALL CHKXER( 'ZUNCSD', INFOT, NOUT, LERR, OK )
|
||||||
INFOT = 11
|
INFOT = 11
|
||||||
CALL ZUNCSD( 'Y', 'Y', 'Y', 'Y', 'N', 'N',
|
CALL ZUNCSD( 'Y', 'Y', 'Y', 'Y', 'N', 'N',
|
||||||
$ 1, 1, 1, A, -1, A,
|
$ 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,
|
$ A, 1, A, 1, A, 1, A,
|
||||||
$ 1, W, LW, RW, LW, IW, INFO )
|
$ 1, W, LW, RW, LW, IW, INFO )
|
||||||
CALL CHKXER( 'ZUNCSD', INFOT, NOUT, LERR, OK )
|
CALL CHKXER( 'ZUNCSD', INFOT, NOUT, LERR, OK )
|
||||||
INFOT = 20
|
INFOT = 20
|
||||||
CALL ZUNCSD( 'Y', 'Y', 'Y', 'Y', 'N', 'N',
|
CALL ZUNCSD( 'Y', 'Y', 'Y', 'Y', 'N', 'N',
|
||||||
$ 1, 1, 1, A, 1, A,
|
$ 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,
|
$ A, -1, A, 1, A, 1, A,
|
||||||
$ 1, W, LW, RW, LW, IW, INFO )
|
$ 1, W, LW, RW, LW, IW, INFO )
|
||||||
CALL CHKXER( 'ZUNCSD', INFOT, NOUT, LERR, OK )
|
CALL CHKXER( 'ZUNCSD', INFOT, NOUT, LERR, OK )
|
||||||
INFOT = 22
|
INFOT = 22
|
||||||
CALL ZUNCSD( 'Y', 'Y', 'Y', 'Y', 'N', 'N',
|
CALL ZUNCSD( 'Y', 'Y', 'Y', 'Y', 'N', 'N',
|
||||||
$ 1, 1, 1, A, 1, A,
|
$ 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,
|
$ A, 1, A, -1, A, 1, A,
|
||||||
$ 1, W, LW, RW, LW, IW, INFO )
|
$ 1, W, LW, RW, LW, IW, INFO )
|
||||||
CALL CHKXER( 'ZUNCSD', INFOT, NOUT, LERR, OK )
|
CALL CHKXER( 'ZUNCSD', INFOT, NOUT, LERR, OK )
|
||||||
INFOT = 24
|
INFOT = 24
|
||||||
CALL ZUNCSD( 'Y', 'Y', 'Y', 'Y', 'N', 'N',
|
CALL ZUNCSD( 'Y', 'Y', 'Y', 'Y', 'N', 'N',
|
||||||
$ 1, 1, 1, A, 1, A,
|
$ 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,
|
$ A, 1, A, 1, A, -1, A,
|
||||||
$ 1, W, LW, RW, LW, IW, INFO )
|
$ 1, W, LW, RW, LW, IW, INFO )
|
||||||
CALL CHKXER( 'ZUNCSD', INFOT, NOUT, LERR, OK )
|
CALL CHKXER( 'ZUNCSD', INFOT, NOUT, LERR, OK )
|
||||||
INFOT = 26
|
INFOT = 26
|
||||||
CALL ZUNCSD( 'Y', 'Y', 'Y', 'Y', 'N', 'N',
|
CALL ZUNCSD( 'Y', 'Y', 'Y', 'Y', 'N', 'N',
|
||||||
$ 1, 1, 1, A, 1, A,
|
$ 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,
|
$ A, 1, A, 1, A, 1, A,
|
||||||
$ -1, W, LW, RW, LW, IW, INFO )
|
$ -1, W, LW, RW, LW, IW, INFO )
|
||||||
CALL CHKXER( 'ZUNCSD', INFOT, NOUT, LERR, OK )
|
CALL CHKXER( 'ZUNCSD', INFOT, NOUT, LERR, OK )
|
||||||
|
|
|
@ -25,7 +25,7 @@
|
||||||
*> eigenvalue is negative), and otherwise it returns .FALSE..
|
*> eigenvalue is negative), and otherwise it returns .FALSE..
|
||||||
*>
|
*>
|
||||||
*> It is used by the test routine ZDRGES to test whether the driver
|
*> It is used by the test routine ZDRGES to test whether the driver
|
||||||
*> routine ZGGES succesfully sorts eigenvalues.
|
*> routine ZGGES successfully sorts eigenvalues.
|
||||||
*> \endverbatim
|
*> \endverbatim
|
||||||
*
|
*
|
||||||
* Arguments:
|
* Arguments:
|
||||||
|
@ -51,17 +51,17 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2011
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup complex16_eig
|
*> \ingroup complex16_eig
|
||||||
*
|
*
|
||||||
* =====================================================================
|
* =====================================================================
|
||||||
LOGICAL FUNCTION ZLCTES( Z, D )
|
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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2011
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
COMPLEX*16 D, Z
|
COMPLEX*16 D, Z
|
||||||
|
|
|
@ -22,8 +22,8 @@
|
||||||
*>
|
*>
|
||||||
*> ZSLECT returns .TRUE. if the eigenvalue Z is to be selected,
|
*> ZSLECT returns .TRUE. if the eigenvalue Z is to be selected,
|
||||||
*> otherwise it returns .FALSE.
|
*> otherwise it returns .FALSE.
|
||||||
*> It is used by ZCHK41 to test if ZGEES succesfully sorts eigenvalues,
|
*> It is used by ZCHK41 to test if ZGEES successfully sorts eigenvalues,
|
||||||
*> and by ZCHK43 to test if ZGEESX succesfully sorts eigenvalues.
|
*> and by ZCHK43 to test if ZGEESX successfully sorts eigenvalues.
|
||||||
*>
|
*>
|
||||||
*> The common block /SSLCT/ controls how eigenvalues are selected.
|
*> The common block /SSLCT/ controls how eigenvalues are selected.
|
||||||
*> If SELOPT = 0, then ZSLECT return .TRUE. when real(Z) is less than
|
*> If SELOPT = 0, then ZSLECT return .TRUE. when real(Z) is less than
|
||||||
|
@ -49,17 +49,17 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2011
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup complex16_eig
|
*> \ingroup complex16_eig
|
||||||
*
|
*
|
||||||
* =====================================================================
|
* =====================================================================
|
||||||
LOGICAL FUNCTION ZSLECT( Z )
|
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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2011
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
COMPLEX*16 Z
|
COMPLEX*16 Z
|
||||||
|
|
|
@ -339,3 +339,5 @@ zchkaa.o: zchkaa.f
|
||||||
|
|
||||||
.f.o:
|
.f.o:
|
||||||
$(FORTRAN) $(OPTS) -c $< -o $@
|
$(FORTRAN) $(OPTS) -c $< -o $@
|
||||||
|
|
||||||
|
.NOTPARALLEL:
|
||||||
|
|
|
@ -74,7 +74,7 @@
|
||||||
SUBROUTINE CQRT04(M,N,NB,RESULT)
|
SUBROUTINE CQRT04(M,N,NB,RESULT)
|
||||||
IMPLICIT NONE
|
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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* April 2012
|
* April 2012
|
||||||
|
@ -89,8 +89,9 @@
|
||||||
* ..
|
* ..
|
||||||
* .. Local allocatable arrays
|
* .. Local allocatable arrays
|
||||||
COMPLEX, ALLOCATABLE :: AF(:,:), Q(:,:),
|
COMPLEX, ALLOCATABLE :: AF(:,:), Q(:,:),
|
||||||
$ R(:,:), RWORK(:), WORK( : ), T(:,:),
|
$ R(:,:), WORK( : ), T(:,:),
|
||||||
$ CF(:,:), DF(:,:), A(:,:), C(:,:), D(:,:)
|
$ CF(:,:), DF(:,:), A(:,:), C(:,:), D(:,:)
|
||||||
|
REAL, ALLOCATABLE :: RWORK(:)
|
||||||
*
|
*
|
||||||
* .. Parameters ..
|
* .. Parameters ..
|
||||||
REAL ZERO
|
REAL ZERO
|
||||||
|
|
|
@ -81,7 +81,7 @@
|
||||||
SUBROUTINE CQRT05(M,N,L,NB,RESULT)
|
SUBROUTINE CQRT05(M,N,L,NB,RESULT)
|
||||||
IMPLICIT NONE
|
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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* April 2012
|
* April 2012
|
||||||
|
@ -96,8 +96,9 @@
|
||||||
* ..
|
* ..
|
||||||
* .. Local allocatable arrays
|
* .. Local allocatable arrays
|
||||||
COMPLEX, ALLOCATABLE :: AF(:,:), Q(:,:),
|
COMPLEX, ALLOCATABLE :: AF(:,:), Q(:,:),
|
||||||
$ R(:,:), RWORK(:), WORK( : ), T(:,:),
|
$ R(:,:), WORK( : ), T(:,:),
|
||||||
$ CF(:,:), DF(:,:), A(:,:), C(:,:), D(:,:)
|
$ CF(:,:), DF(:,:), A(:,:), C(:,:), D(:,:)
|
||||||
|
REAL, ALLOCATABLE :: RWORK(:)
|
||||||
*
|
*
|
||||||
* .. Parameters ..
|
* .. Parameters ..
|
||||||
REAL ZERO
|
REAL ZERO
|
||||||
|
|
|
@ -74,7 +74,7 @@
|
||||||
SUBROUTINE ZQRT04(M,N,NB,RESULT)
|
SUBROUTINE ZQRT04(M,N,NB,RESULT)
|
||||||
IMPLICIT NONE
|
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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* April 2012
|
* April 2012
|
||||||
|
@ -89,8 +89,9 @@
|
||||||
* ..
|
* ..
|
||||||
* .. Local allocatable arrays
|
* .. Local allocatable arrays
|
||||||
COMPLEX*16, ALLOCATABLE :: AF(:,:), Q(:,:),
|
COMPLEX*16, ALLOCATABLE :: AF(:,:), Q(:,:),
|
||||||
$ R(:,:), RWORK(:), WORK( : ), T(:,:),
|
$ R(:,:), WORK( : ), T(:,:),
|
||||||
$ CF(:,:), DF(:,:), A(:,:), C(:,:), D(:,:)
|
$ CF(:,:), DF(:,:), A(:,:), C(:,:), D(:,:)
|
||||||
|
DOUBLE PRECISION, ALLOCATABLE :: RWORK(:)
|
||||||
*
|
*
|
||||||
* .. Parameters ..
|
* .. Parameters ..
|
||||||
DOUBLE PRECISION ZERO
|
DOUBLE PRECISION ZERO
|
||||||
|
|
|
@ -81,7 +81,7 @@
|
||||||
SUBROUTINE ZQRT05(M,N,L,NB,RESULT)
|
SUBROUTINE ZQRT05(M,N,L,NB,RESULT)
|
||||||
IMPLICIT NONE
|
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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* April 2012
|
* April 2012
|
||||||
|
@ -96,8 +96,9 @@
|
||||||
* ..
|
* ..
|
||||||
* .. Local allocatable arrays
|
* .. Local allocatable arrays
|
||||||
COMPLEX*16, ALLOCATABLE :: AF(:,:), Q(:,:),
|
COMPLEX*16, ALLOCATABLE :: AF(:,:), Q(:,:),
|
||||||
$ R(:,:), RWORK(:), WORK( : ), T(:,:),
|
$ R(:,:), WORK( : ), T(:,:),
|
||||||
$ CF(:,:), DF(:,:), A(:,:), C(:,:), D(:,:)
|
$ CF(:,:), DF(:,:), A(:,:), C(:,:), D(:,:)
|
||||||
|
DOUBLE PRECISION, ALLOCATABLE :: RWORK(:)
|
||||||
*
|
*
|
||||||
* .. Parameters ..
|
* .. Parameters ..
|
||||||
DOUBLE PRECISION ZERO
|
DOUBLE PRECISION ZERO
|
||||||
|
|
|
@ -30,7 +30,7 @@
|
||||||
*> \verbatim
|
*> \verbatim
|
||||||
*>
|
*>
|
||||||
*> CLATM2 returns the (I,J) entry of a random matrix of dimension
|
*> 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
|
*> CLATMR routine in order to build random test matrices. No error
|
||||||
*> checking on parameters is done, because this routine is called in
|
*> checking on parameters is done, because this routine is called in
|
||||||
*> a tight loop by CLATMR which has already checked the parameters.
|
*> a tight loop by CLATMR which has already checked the parameters.
|
||||||
|
@ -204,7 +204,7 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2011
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup complex_matgen
|
*> \ingroup complex_matgen
|
||||||
*
|
*
|
||||||
|
@ -212,10 +212,10 @@
|
||||||
COMPLEX FUNCTION CLATM2( M, N, I, J, KL, KU, IDIST, ISEED, D,
|
COMPLEX FUNCTION CLATM2( M, N, I, J, KL, KU, IDIST, ISEED, D,
|
||||||
$ IGRADE, DL, DR, IPVTNG, IWORK, SPARSE )
|
$ 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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2011
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
*
|
*
|
||||||
|
|
|
@ -32,7 +32,7 @@
|
||||||
*> \verbatim
|
*> \verbatim
|
||||||
*>
|
*>
|
||||||
*> CLATM3 returns the (ISUB,JSUB) entry of a random matrix of
|
*> 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
|
*> is the final position of the (I,J) entry after pivoting
|
||||||
*> according to IPVTNG and IWORK. CLATM3 is called by the
|
*> according to IPVTNG and IWORK. CLATM3 is called by the
|
||||||
*> CLATMR routine in order to build random test matrices. No error
|
*> CLATMR routine in order to build random test matrices. No error
|
||||||
|
@ -220,7 +220,7 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2011
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup complex_matgen
|
*> \ingroup complex_matgen
|
||||||
*
|
*
|
||||||
|
@ -229,10 +229,10 @@
|
||||||
$ ISEED, D, IGRADE, DL, DR, IPVTNG, IWORK,
|
$ ISEED, D, IGRADE, DL, DR, IPVTNG, IWORK,
|
||||||
$ SPARSE )
|
$ 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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2011
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
*
|
*
|
||||||
|
|
|
@ -48,7 +48,7 @@
|
||||||
*> \param[in] PRTYPE
|
*> \param[in] PRTYPE
|
||||||
*> \verbatim
|
*> \verbatim
|
||||||
*> PRTYPE is INTEGER
|
*> 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).
|
*> (see futher details).
|
||||||
*> \endverbatim
|
*> \endverbatim
|
||||||
*>
|
*>
|
||||||
|
@ -192,7 +192,7 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2011
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup complex_matgen
|
*> \ingroup complex_matgen
|
||||||
*
|
*
|
||||||
|
@ -268,10 +268,10 @@
|
||||||
$ E, LDE, F, LDF, R, LDR, L, LDL, ALPHA, QBLCKA,
|
$ E, LDE, F, LDF, R, LDR, L, LDL, ALPHA, QBLCKA,
|
||||||
$ QBLCKB )
|
$ 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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2011
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
INTEGER LDA, LDB, LDC, LDD, LDE, LDF, LDL, LDR, M, N,
|
INTEGER LDA, LDB, LDC, LDD, LDE, LDF, LDL, LDR, M, N,
|
||||||
|
|
|
@ -30,7 +30,7 @@
|
||||||
*> \verbatim
|
*> \verbatim
|
||||||
*>
|
*>
|
||||||
*> DLATM2 returns the (I,J) entry of a random matrix of dimension
|
*> 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
|
*> DLATMR routine in order to build random test matrices. No error
|
||||||
*> checking on parameters is done, because this routine is called in
|
*> checking on parameters is done, because this routine is called in
|
||||||
*> a tight loop by DLATMR which has already checked the parameters.
|
*> a tight loop by DLATMR which has already checked the parameters.
|
||||||
|
@ -200,7 +200,7 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2011
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup double_matgen
|
*> \ingroup double_matgen
|
||||||
*
|
*
|
||||||
|
@ -208,10 +208,10 @@
|
||||||
DOUBLE PRECISION FUNCTION DLATM2( M, N, I, J, KL, KU, IDIST,
|
DOUBLE PRECISION FUNCTION DLATM2( M, N, I, J, KL, KU, IDIST,
|
||||||
$ ISEED, D, IGRADE, DL, DR, IPVTNG, IWORK, SPARSE )
|
$ 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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2011
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
*
|
*
|
||||||
|
|
|
@ -32,7 +32,7 @@
|
||||||
*> \verbatim
|
*> \verbatim
|
||||||
*>
|
*>
|
||||||
*> DLATM3 returns the (ISUB,JSUB) entry of a random matrix of
|
*> 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
|
*> is the final position of the (I,J) entry after pivoting
|
||||||
*> according to IPVTNG and IWORK. DLATM3 is called by the
|
*> according to IPVTNG and IWORK. DLATM3 is called by the
|
||||||
*> DLATMR routine in order to build random test matrices. No error
|
*> DLATMR routine in order to build random test matrices. No error
|
||||||
|
@ -217,7 +217,7 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2011
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup double_matgen
|
*> \ingroup double_matgen
|
||||||
*
|
*
|
||||||
|
@ -226,10 +226,10 @@
|
||||||
$ IDIST, ISEED, D, IGRADE, DL, DR, IPVTNG, IWORK,
|
$ IDIST, ISEED, D, IGRADE, DL, DR, IPVTNG, IWORK,
|
||||||
$ SPARSE )
|
$ 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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2011
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
*
|
*
|
||||||
|
|
|
@ -48,7 +48,7 @@
|
||||||
*> \param[in] PRTYPE
|
*> \param[in] PRTYPE
|
||||||
*> \verbatim
|
*> \verbatim
|
||||||
*> PRTYPE is INTEGER
|
*> 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).
|
*> (see futher details).
|
||||||
*> \endverbatim
|
*> \endverbatim
|
||||||
*>
|
*>
|
||||||
|
@ -192,7 +192,7 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2011
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup double_matgen
|
*> \ingroup double_matgen
|
||||||
*
|
*
|
||||||
|
@ -268,10 +268,10 @@
|
||||||
$ E, LDE, F, LDF, R, LDR, L, LDL, ALPHA, QBLCKA,
|
$ E, LDE, F, LDF, R, LDR, L, LDL, ALPHA, QBLCKA,
|
||||||
$ QBLCKB )
|
$ 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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2011
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
INTEGER LDA, LDB, LDC, LDD, LDE, LDF, LDL, LDR, M, N,
|
INTEGER LDA, LDB, LDC, LDD, LDE, LDF, LDL, LDR, M, N,
|
||||||
|
|
|
@ -30,7 +30,7 @@
|
||||||
*> \verbatim
|
*> \verbatim
|
||||||
*>
|
*>
|
||||||
*> SLATM2 returns the (I,J) entry of a random matrix of dimension
|
*> 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
|
*> SLATMR routine in order to build random test matrices. No error
|
||||||
*> checking on parameters is done, because this routine is called in
|
*> checking on parameters is done, because this routine is called in
|
||||||
*> a tight loop by SLATMR which has already checked the parameters.
|
*> a tight loop by SLATMR which has already checked the parameters.
|
||||||
|
@ -200,7 +200,7 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2011
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup real_matgen
|
*> \ingroup real_matgen
|
||||||
*
|
*
|
||||||
|
@ -208,10 +208,10 @@
|
||||||
REAL FUNCTION SLATM2( M, N, I, J, KL, KU, IDIST,
|
REAL FUNCTION SLATM2( M, N, I, J, KL, KU, IDIST,
|
||||||
$ ISEED, D, IGRADE, DL, DR, IPVTNG, IWORK, SPARSE )
|
$ 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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2011
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
*
|
*
|
||||||
|
|
|
@ -32,7 +32,7 @@
|
||||||
*> \verbatim
|
*> \verbatim
|
||||||
*>
|
*>
|
||||||
*> SLATM3 returns the (ISUB,JSUB) entry of a random matrix of
|
*> 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
|
*> is the final position of the (I,J) entry after pivoting
|
||||||
*> according to IPVTNG and IWORK. SLATM3 is called by the
|
*> according to IPVTNG and IWORK. SLATM3 is called by the
|
||||||
*> SLATMR routine in order to build random test matrices. No error
|
*> SLATMR routine in order to build random test matrices. No error
|
||||||
|
@ -217,7 +217,7 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2011
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup real_matgen
|
*> \ingroup real_matgen
|
||||||
*
|
*
|
||||||
|
@ -226,10 +226,10 @@
|
||||||
$ IDIST, ISEED, D, IGRADE, DL, DR, IPVTNG, IWORK,
|
$ IDIST, ISEED, D, IGRADE, DL, DR, IPVTNG, IWORK,
|
||||||
$ SPARSE )
|
$ 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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2011
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
*
|
*
|
||||||
|
|
|
@ -48,7 +48,7 @@
|
||||||
*> \param[in] PRTYPE
|
*> \param[in] PRTYPE
|
||||||
*> \verbatim
|
*> \verbatim
|
||||||
*> PRTYPE is INTEGER
|
*> 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).
|
*> (see futher details).
|
||||||
*> \endverbatim
|
*> \endverbatim
|
||||||
*>
|
*>
|
||||||
|
@ -192,7 +192,7 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2011
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup real_matgen
|
*> \ingroup real_matgen
|
||||||
*
|
*
|
||||||
|
@ -268,10 +268,10 @@
|
||||||
$ E, LDE, F, LDF, R, LDR, L, LDL, ALPHA, QBLCKA,
|
$ E, LDE, F, LDF, R, LDR, L, LDL, ALPHA, QBLCKA,
|
||||||
$ QBLCKB )
|
$ 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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2011
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
INTEGER LDA, LDB, LDC, LDD, LDE, LDF, LDL, LDR, M, N,
|
INTEGER LDA, LDB, LDC, LDD, LDE, LDF, LDL, LDR, M, N,
|
||||||
|
|
|
@ -30,7 +30,7 @@
|
||||||
*> \verbatim
|
*> \verbatim
|
||||||
*>
|
*>
|
||||||
*> ZLATM2 returns the (I,J) entry of a random matrix of dimension
|
*> 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
|
*> ZLATMR routine in order to build random test matrices. No error
|
||||||
*> checking on parameters is done, because this routine is called in
|
*> checking on parameters is done, because this routine is called in
|
||||||
*> a tight loop by ZLATMR which has already checked the parameters.
|
*> a tight loop by ZLATMR which has already checked the parameters.
|
||||||
|
@ -203,7 +203,7 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2011
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup complex16_matgen
|
*> \ingroup complex16_matgen
|
||||||
*
|
*
|
||||||
|
@ -211,10 +211,10 @@
|
||||||
COMPLEX*16 FUNCTION ZLATM2( M, N, I, J, KL, KU, IDIST,
|
COMPLEX*16 FUNCTION ZLATM2( M, N, I, J, KL, KU, IDIST,
|
||||||
$ ISEED, D, IGRADE, DL, DR, IPVTNG, IWORK, SPARSE )
|
$ 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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2011
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
*
|
*
|
||||||
|
|
|
@ -32,7 +32,7 @@
|
||||||
*> \verbatim
|
*> \verbatim
|
||||||
*>
|
*>
|
||||||
*> ZLATM3 returns the (ISUB,JSUB) entry of a random matrix of
|
*> 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
|
*> is the final position of the (I,J) entry after pivoting
|
||||||
*> according to IPVTNG and IWORK. ZLATM3 is called by the
|
*> according to IPVTNG and IWORK. ZLATM3 is called by the
|
||||||
*> ZLATMR routine in order to build random test matrices. No error
|
*> ZLATMR routine in order to build random test matrices. No error
|
||||||
|
@ -220,7 +220,7 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2011
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup complex16_matgen
|
*> \ingroup complex16_matgen
|
||||||
*
|
*
|
||||||
|
@ -229,10 +229,10 @@
|
||||||
$ IDIST, ISEED, D, IGRADE, DL, DR, IPVTNG, IWORK,
|
$ IDIST, ISEED, D, IGRADE, DL, DR, IPVTNG, IWORK,
|
||||||
$ SPARSE )
|
$ 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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2011
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
*
|
*
|
||||||
|
|
|
@ -48,7 +48,7 @@
|
||||||
*> \param[in] PRTYPE
|
*> \param[in] PRTYPE
|
||||||
*> \verbatim
|
*> \verbatim
|
||||||
*> PRTYPE is INTEGER
|
*> 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).
|
*> (see futher details).
|
||||||
*> \endverbatim
|
*> \endverbatim
|
||||||
*>
|
*>
|
||||||
|
@ -192,7 +192,7 @@
|
||||||
*> \author Univ. of Colorado Denver
|
*> \author Univ. of Colorado Denver
|
||||||
*> \author NAG Ltd.
|
*> \author NAG Ltd.
|
||||||
*
|
*
|
||||||
*> \date November 2011
|
*> \date June 2016
|
||||||
*
|
*
|
||||||
*> \ingroup complex16_matgen
|
*> \ingroup complex16_matgen
|
||||||
*
|
*
|
||||||
|
@ -268,10 +268,10 @@
|
||||||
$ E, LDE, F, LDF, R, LDR, L, LDL, ALPHA, QBLCKA,
|
$ E, LDE, F, LDF, R, LDR, L, LDL, ALPHA, QBLCKA,
|
||||||
$ QBLCKB )
|
$ 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, --
|
* -- LAPACK is a software package provided by Univ. of Tennessee, --
|
||||||
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
|
||||||
* November 2011
|
* June 2016
|
||||||
*
|
*
|
||||||
* .. Scalar Arguments ..
|
* .. Scalar Arguments ..
|
||||||
INTEGER LDA, LDB, LDC, LDD, LDE, LDF, LDL, LDR, M, N,
|
INTEGER LDA, LDB, LDC, LDD, LDE, LDF, LDL, LDR, M, N,
|
||||||
|
|
|
@ -582,3 +582,5 @@ FRCEIG:
|
||||||
|
|
||||||
FRC:
|
FRC:
|
||||||
@FRC=$(FRC)
|
@FRC=$(FRC)
|
||||||
|
|
||||||
|
.NOTPARALLEL:
|
||||||
|
|
|
@ -5,5 +5,5 @@ Data file for testing COMPLEX LAPACK linear equation routines RFP format
|
||||||
1 2 15 Values of NRHS (number of right hand sides)
|
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)
|
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
|
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
|
T Put T to test the error exits
|
||||||
|
|
|
@ -10,7 +10,7 @@ NEP: Data file for testing Nonsymmetric Eigenvalue Problem routines
|
||||||
0 5 7 3 200 Values of INIBL (nibble crossover point)
|
0 5 7 3 200 Values of INIBL (nibble crossover point)
|
||||||
1 2 4 2 1 Values of ISHFTS (number of simultaneous shifts)
|
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)
|
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
|
T Put T to test the error exits
|
||||||
1 Code to interpret the seed
|
1 Code to interpret the seed
|
||||||
NEP 21
|
NEP 21
|
||||||
|
|
|
@ -1,11 +1,11 @@
|
||||||
SEP: Data file for testing Symmetric Eigenvalue Problem routines
|
SEP: Data file for testing Symmetric Eigenvalue Problem routines
|
||||||
8 Number of values of N
|
6 Number of values of N
|
||||||
0 1 2 3 5 16 19 20 Values of N (dimension)
|
0 1 2 3 5 20 Values of N (dimension)
|
||||||
5 Number of values of NB
|
5 Number of values of NB
|
||||||
1 3 3 3 10 Values of NB (blocksize)
|
1 3 3 3 10 Values of NB (blocksize)
|
||||||
2 2 2 2 2 Values of NBMIN (minimum blocksize)
|
2 2 2 2 2 Values of NBMIN (minimum blocksize)
|
||||||
1 0 5 9 1 Values of NX (crossover point)
|
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 LAPACK routines
|
||||||
T Put T to test the driver routines
|
T Put T to test the driver routines
|
||||||
T Put T to test the error exits
|
T Put T to test the error exits
|
||||||
|
|
|
@ -7,7 +7,7 @@ SVD: Data file for testing Singular Value Decomposition routines
|
||||||
2 2 2 2 2 Values of NBMIN (minimum blocksize)
|
2 2 2 2 2 Values of NBMIN (minimum blocksize)
|
||||||
1 0 5 9 1 Values of NX (crossover point)
|
1 0 5 9 1 Values of NX (crossover point)
|
||||||
2 0 2 2 2 Values of NRHS
|
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 LAPACK routines
|
||||||
T Put T to test the driver routines
|
T Put T to test the driver routines
|
||||||
T Put T to test the error exits
|
T Put T to test the error exits
|
||||||
|
|
Loading…
Reference in New Issue