diff options
Diffstat (limited to 'TESTING/EIG')
108 files changed, 715 insertions, 1430 deletions
diff --git a/TESTING/EIG/cchkbb.f b/TESTING/EIG/cchkbb.f index a930fc10..0f7cfbec 100644 --- a/TESTING/EIG/cchkbb.f +++ b/TESTING/EIG/cchkbb.f @@ -316,11 +316,9 @@ *> \verbatim *> INFO is INTEGER *> If 0, then everything ran OK. -*> \endverbatim -*> \verbatim +*> *>----------------------------------------------------------------------- -*> \endverbatim -*> \verbatim +*> *> Some Local Variables and Parameters: *> ---- ----- --------- --- ---------- *> ZERO, ONE Real 0 and 1. @@ -334,8 +332,7 @@ *> so far. *> COND, IMODE Values to be passed to the matrix generators. *> ANORM Norm of A; passed to matrix generators. -*> \endverbatim -*> \verbatim +*> *> OVFL, UNFL Overflow and underflow thresholds. *> ULP, ULPINV Finest relative precision and its inverse. *> RTOVFL, RTUNFL Square roots of the previous 2 values. diff --git a/TESTING/EIG/cchkbd.f b/TESTING/EIG/cchkbd.f index 659bbb98..604a43a9 100644 --- a/TESTING/EIG/cchkbd.f +++ b/TESTING/EIG/cchkbd.f @@ -367,15 +367,12 @@ *> If CLATMR, CLATMS, CGEBRD, CUNGBR, or CBDSQR, *> returns an error code, the *> absolute value of it is returned. -*> \endverbatim -*> \verbatim +*> *>----------------------------------------------------------------------- -*> \endverbatim -*> \verbatim +*> *> Some Local Variables and Parameters: *> ---- ----- --------- --- ---------- -*> \endverbatim -*> \verbatim +*> *> ZERO, ONE Real 0 and 1. *> MAXTYP The number of types defined. *> NTEST The number of tests performed, or which can @@ -388,13 +385,11 @@ *> NFAIL The number of tests which have exceeded THRESH *> COND, IMODE Values to be passed to the matrix generators. *> ANORM Norm of A; passed to matrix generators. -*> \endverbatim -*> \verbatim +*> *> OVFL, UNFL Overflow and underflow thresholds. *> RTOVFL, RTUNFL Square roots of the previous 2 values. *> ULP, ULPINV Finest relative precision and its inverse. -*> \endverbatim -*> \verbatim +*> *> The following four arrays decode JTYPE: *> KTYPE(j) The general type (1-10) for type "j". *> KMODE(j) The MODE value to be passed to the matrix diff --git a/TESTING/EIG/cchkhb.f b/TESTING/EIG/cchkhb.f index 52d8c981..80e8b0af 100644 --- a/TESTING/EIG/cchkhb.f +++ b/TESTING/EIG/cchkhb.f @@ -254,11 +254,9 @@ *> \verbatim *> INFO is INTEGER *> If 0, then everything ran OK. -*> \endverbatim -*> \verbatim +*> *>----------------------------------------------------------------------- -*> \endverbatim -*> \verbatim +*> *> Some Local Variables and Parameters: *> ---- ----- --------- --- ---------- *> ZERO, ONE Real 0 and 1. @@ -272,8 +270,7 @@ *> so far. *> COND, IMODE Values to be passed to the matrix generators. *> ANORM Norm of A; passed to matrix generators. -*> \endverbatim -*> \verbatim +*> *> OVFL, UNFL Overflow and underflow thresholds. *> ULP, ULPINV Finest relative precision and its inverse. *> RTOVFL, RTUNFL Square roots of the previous 2 values. diff --git a/TESTING/EIG/cchkhs.f b/TESTING/EIG/cchkhs.f index 66663338..f931e7fe 100644 --- a/TESTING/EIG/cchkhs.f +++ b/TESTING/EIG/cchkhs.f @@ -176,15 +176,13 @@ *> The number of sizes of matrices to use. If it is zero, *> CCHKHS does nothing. It must be at least zero. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> NN - INTEGER array, dimension (NSIZES) *> An array containing the sizes to be used for the matrices. *> Zero values will be skipped. The values must be at least *> zero. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> NTYPES - INTEGER *> The number of elements in DOTYPE. If it is zero, CCHKHS *> does nothing. It must be at least zero. If it is MAXTYP+1 @@ -193,8 +191,7 @@ *> is only useful if DOTYPE(1:MAXTYP) is .FALSE. and *> DOTYPE(MAXTYP+1) is .TRUE. . *> Not modified. -*> \endverbatim -*> \verbatim +*> *> DOTYPE - LOGICAL array, dimension (NTYPES) *> If DOTYPE(j) is .TRUE., then for each size in NN a *> matrix of that size and of type j will be generated. @@ -204,8 +201,7 @@ *> than MAXTYP, DOTYPE(MAXTYP+1) through DOTYPE(NTYPES) *> will be ignored. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> ISEED - INTEGER array, dimension (4) *> On entry ISEED specifies the seed of the random number *> generator. The array elements should be between 0 and 4095; @@ -217,8 +213,7 @@ *> next call to CCHKHS to continue the same random number *> sequence. *> Modified. -*> \endverbatim -*> \verbatim +*> *> THRESH - REAL *> A test will count as "failed" if the "error", computed as *> described above, exceeds THRESH. Note that the error @@ -227,74 +222,62 @@ *> it should not depend on the precision (single vs. double) *> or the size of the matrix. It must be at least zero. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> NOUNIT - INTEGER *> The FORTRAN unit number for printing out error messages *> (e.g., if a routine returns IINFO not equal to 0.) *> Not modified. -*> \endverbatim -*> \verbatim +*> *> A - COMPLEX array, dimension (LDA,max(NN)) *> Used to hold the matrix whose eigenvalues are to be *> computed. On exit, A contains the last matrix actually *> used. *> Modified. -*> \endverbatim -*> \verbatim +*> *> LDA - INTEGER *> The leading dimension of A, H, T1 and T2. It must be at *> least 1 and at least max( NN ). *> Not modified. -*> \endverbatim -*> \verbatim +*> *> H - COMPLEX array, dimension (LDA,max(NN)) *> The upper hessenberg matrix computed by CGEHRD. On exit, *> H contains the Hessenberg form of the matrix in A. *> Modified. -*> \endverbatim -*> \verbatim +*> *> T1 - COMPLEX array, dimension (LDA,max(NN)) *> The Schur (="quasi-triangular") matrix computed by CHSEQR *> if Z is computed. On exit, T1 contains the Schur form of *> the matrix in A. *> Modified. -*> \endverbatim -*> \verbatim +*> *> T2 - COMPLEX array, dimension (LDA,max(NN)) *> The Schur matrix computed by CHSEQR when Z is not computed. *> This should be identical to T1. *> Modified. -*> \endverbatim -*> \verbatim +*> *> LDU - INTEGER *> The leading dimension of U, Z, UZ and UU. It must be at *> least 1 and at least max( NN ). *> Not modified. -*> \endverbatim -*> \verbatim +*> *> U - COMPLEX array, dimension (LDU,max(NN)) *> The unitary matrix computed by CGEHRD. *> Modified. -*> \endverbatim -*> \verbatim +*> *> Z - COMPLEX array, dimension (LDU,max(NN)) *> The unitary matrix computed by CHSEQR. *> Modified. -*> \endverbatim -*> \verbatim +*> *> UZ - COMPLEX array, dimension (LDU,max(NN)) *> The product of U times Z. *> Modified. -*> \endverbatim -*> \verbatim +*> *> W1 - COMPLEX array, dimension (max(NN)) *> The eigenvalues of A, as computed by a full Schur *> decomposition H = Z T Z'. On exit, W1 contains the *> eigenvalues of the matrix in A. *> Modified. -*> \endverbatim -*> \verbatim +*> *> W3 - COMPLEX array, dimension (max(NN)) *> The eigenvalues of A, as computed by a partial Schur *> decomposition (Z not computed, T only computed as much @@ -302,72 +285,59 @@ *> W3 contains the eigenvalues of the matrix in A, possibly *> perturbed by CHSEIN. *> Modified. -*> \endverbatim -*> \verbatim +*> *> EVECTL - COMPLEX array, dimension (LDU,max(NN)) *> The conjugate transpose of the (upper triangular) left *> eigenvector matrix for the matrix in T1. *> Modified. -*> \endverbatim -*> \verbatim +*> *> EVECTR - COMPLEX array, dimension (LDU,max(NN)) *> The (upper triangular) right eigenvector matrix for the *> matrix in T1. *> Modified. -*> \endverbatim -*> \verbatim +*> *> EVECTY - COMPLEX array, dimension (LDU,max(NN)) *> The conjugate transpose of the left eigenvector matrix *> for the matrix in H. *> Modified. -*> \endverbatim -*> \verbatim +*> *> EVECTX - COMPLEX array, dimension (LDU,max(NN)) *> The right eigenvector matrix for the matrix in H. *> Modified. -*> \endverbatim -*> \verbatim +*> *> UU - COMPLEX array, dimension (LDU,max(NN)) *> Details of the unitary matrix computed by CGEHRD. *> Modified. -*> \endverbatim -*> \verbatim +*> *> TAU - COMPLEX array, dimension (max(NN)) *> Further details of the unitary matrix computed by CGEHRD. *> Modified. -*> \endverbatim -*> \verbatim +*> *> WORK - COMPLEX array, dimension (NWORK) *> Workspace. *> Modified. -*> \endverbatim -*> \verbatim +*> *> NWORK - INTEGER *> The number of entries in WORK. NWORK >= 4*NN(j)*NN(j) + 2. -*> \endverbatim -*> \verbatim +*> *> RWORK - REAL array, dimension (max(NN)) *> Workspace. Could be equivalenced to IWORK, but not SELECT. *> Modified. -*> \endverbatim -*> \verbatim +*> *> IWORK - INTEGER array, dimension (max(NN)) *> Workspace. *> Modified. -*> \endverbatim -*> \verbatim +*> *> SELECT - LOGICAL array, dimension (max(NN)) *> Workspace. Could be equivalenced to IWORK, but not RWORK. *> Modified. -*> \endverbatim -*> \verbatim +*> *> RESULT - REAL array, dimension (14) *> The values computed by the fourteen tests described above. *> The values are currently limited to 1/ulp, to avoid *> overflow. *> Modified. -*> \endverbatim -*> \verbatim +*> *> INFO - INTEGER *> If 0, then everything ran OK. *> -1: NSIZES < 0 @@ -384,15 +354,12 @@ *> If >2, then 30*N iterations were not enough to find an *> eigenvalue or to decompose the problem. *> Modified. -*> \endverbatim -*> \verbatim +*> *>----------------------------------------------------------------------- -*> \endverbatim -*> \verbatim +*> *> Some Local Variables and Parameters: *> ---- ----- --------- --- ---------- -*> \endverbatim -*> \verbatim +*> *> ZERO, ONE Real 0 and 1. *> MAXTYP The number of types defined. *> MTEST The number of tests defined: care must be taken @@ -411,14 +378,12 @@ *> COND, CONDS, *> IMODE Values to be passed to the matrix generators. *> ANORM Norm of A; passed to matrix generators. -*> \endverbatim -*> \verbatim +*> *> OVFL, UNFL Overflow and underflow thresholds. *> ULP, ULPINV Finest relative precision and its inverse. *> RTOVFL, RTUNFL, *> RTULP, RTULPI Square roots of the previous 4 values. -*> \endverbatim -*> \verbatim +*> *> The following four arrays decode JTYPE: *> KTYPE(j) The general type (1-10) for type "j". *> KMODE(j) The MODE value to be passed to the matrix diff --git a/TESTING/EIG/cchkst.f b/TESTING/EIG/cchkst.f index d72c02b5..90caf9d6 100644 --- a/TESTING/EIG/cchkst.f +++ b/TESTING/EIG/cchkst.f @@ -557,11 +557,9 @@ *> If CLATMR, CLATMS, CHETRD, CUNGTR, CSTEQR, SSTERF, *> or CUNMC2 returns an error code, the *> absolute value of it is returned. -*> \endverbatim -*> \verbatim +*> *>----------------------------------------------------------------------- -*> \endverbatim -*> \verbatim +*> *> Some Local Variables and Parameters: *> ---- ----- --------- --- ---------- *> ZERO, ONE Real 0 and 1. @@ -576,8 +574,7 @@ *> so far. *> COND, IMODE Values to be passed to the matrix generators. *> ANORM Norm of A; passed to matrix generators. -*> \endverbatim -*> \verbatim +*> *> OVFL, UNFL Overflow and underflow thresholds. *> ULP, ULPINV Finest relative precision and its inverse. *> RTOVFL, RTUNFL Square roots of the previous 2 values. diff --git a/TESTING/EIG/cdrges.f b/TESTING/EIG/cdrges.f index 1d76f265..0f0a265a 100644 --- a/TESTING/EIG/cdrges.f +++ b/TESTING/EIG/cdrges.f @@ -321,8 +321,7 @@ *> \param[out] BETA *> \verbatim *> BETA is COMPLEX array, dimension (max(NN)) -*> \endverbatim -*> \verbatim +*> *> The generalized eigenvalues of (A,B) computed by CGGES. *> ALPHA(k) / BETA(k) is the k-th generalized eigenvalue of A *> and B. diff --git a/TESTING/EIG/cdrgev.f b/TESTING/EIG/cdrgev.f index 21cb4bf0..15e872b7 100644 --- a/TESTING/EIG/cdrgev.f +++ b/TESTING/EIG/cdrgev.f @@ -328,8 +328,7 @@ *> \param[out] BETA *> \verbatim *> BETA is COMPLEX array, dimension (max(NN)) -*> \endverbatim -*> \verbatim +*> *> The generalized eigenvalues of (A,B) computed by CGGEV. *> ( ALPHAR(k)+ALPHAI(k)*i ) / BETA(k) is the k-th *> generalized eigenvalue of A and B. @@ -343,8 +342,7 @@ *> \param[out] BETA1 *> \verbatim *> BETA1 is COMPLEX array, dimension (max(NN)) -*> \endverbatim -*> \verbatim +*> *> Like ALPHAR, ALPHAI, BETA, these arrays contain the *> eigenvalues of A and B, but those computed when CGGEV only *> computes a partial eigendecomposition, i.e. not the diff --git a/TESTING/EIG/cdrgsx.f b/TESTING/EIG/cdrgsx.f index d73f6683..c285195c 100644 --- a/TESTING/EIG/cdrgsx.f +++ b/TESTING/EIG/cdrgsx.f @@ -268,8 +268,7 @@ *> \param[out] BETA *> \verbatim *> BETA is COMPLEX array, dimension (NSIZE) -*> \endverbatim -*> \verbatim +*> *> On exit, ALPHA/BETA are the eigenvalues. *> \endverbatim *> diff --git a/TESTING/EIG/cdrgvx.f b/TESTING/EIG/cdrgvx.f index 264b4cd1..77bb5aa5 100644 --- a/TESTING/EIG/cdrgvx.f +++ b/TESTING/EIG/cdrgvx.f @@ -180,8 +180,7 @@ *> \param[out] BETA *> \verbatim *> BETA is COMPLEX array, dimension (NSIZE) -*> \endverbatim -*> \verbatim +*> *> On exit, ALPHA/BETA are the eigenvalues. *> \endverbatim *> diff --git a/TESTING/EIG/cdrves.f b/TESTING/EIG/cdrves.f index a79cfe14..11bf31da 100644 --- a/TESTING/EIG/cdrves.f +++ b/TESTING/EIG/cdrves.f @@ -336,11 +336,9 @@ *> -18: NWORK too small. *> If CLATMR, CLATMS, CLATME or CGEES returns an error code, *> the absolute value of it is returned. -*> \endverbatim -*> \verbatim +*> *>----------------------------------------------------------------------- -*> \endverbatim -*> \verbatim +*> *> Some Local Variables and Parameters: *> ---- ----- --------- --- ---------- *> ZERO, ONE Real 0 and 1. @@ -350,8 +348,7 @@ *> COND, CONDS, *> IMODE Values to be passed to the matrix generators. *> ANORM Norm of A; passed to matrix generators. -*> \endverbatim -*> \verbatim +*> *> OVFL, UNFL Overflow and underflow thresholds. *> ULP, ULPINV Finest relative precision and its inverse. *> RTULP, RTULPI Square roots of the previous 4 values. diff --git a/TESTING/EIG/cdrvev.f b/TESTING/EIG/cdrvev.f index 412554c1..7ac5cade 100644 --- a/TESTING/EIG/cdrvev.f +++ b/TESTING/EIG/cdrvev.f @@ -346,15 +346,12 @@ *> -21: NWORK too small. *> If CLATMR, CLATMS, CLATME or CGEEV returns an error code, *> the absolute value of it is returned. -*> \endverbatim -*> \verbatim +*> *>----------------------------------------------------------------------- -*> \endverbatim -*> \verbatim +*> *> Some Local Variables and Parameters: *> ---- ----- --------- --- ---------- -*> \endverbatim -*> \verbatim +*> *> ZERO, ONE Real 0 and 1. *> MAXTYP The number of types defined. *> NMAX Largest value in NN. @@ -362,13 +359,11 @@ *> COND, CONDS, *> IMODE Values to be passed to the matrix generators. *> ANORM Norm of A; passed to matrix generators. -*> \endverbatim -*> \verbatim +*> *> OVFL, UNFL Overflow and underflow thresholds. *> ULP, ULPINV Finest relative precision and its inverse. *> RTULP, RTULPI Square roots of the previous 4 values. -*> \endverbatim -*> \verbatim +*> *> The following four arrays decode JTYPE: *> KTYPE(j) The general type (1-10) for type "j". *> KMODE(j) The MODE value to be passed to the matrix diff --git a/TESTING/EIG/cdrvgg.f b/TESTING/EIG/cdrvgg.f index 20dba2de..ebe3f9a3 100644 --- a/TESTING/EIG/cdrvgg.f +++ b/TESTING/EIG/cdrvgg.f @@ -333,8 +333,7 @@ *> \param[out] BETA1 *> \verbatim *> BETA1 is COMPLEX array, dimension (max(NN)) -*> \endverbatim -*> \verbatim +*> *> The generalized eigenvalues of (A,B) computed by CGEGS. *> ALPHA1(k) / BETA1(k) is the k-th generalized eigenvalue of *> the matrices in A and B. @@ -348,8 +347,7 @@ *> \param[out] BETA2 *> \verbatim *> BETA2 is COMPLEX array, dimension (max(NN)) -*> \endverbatim -*> \verbatim +*> *> The generalized eigenvalues of (A,B) computed by CGEGV. *> ALPHA2(k) / BETA2(k) is the k-th generalized eigenvalue of *> the matrices in A and B. diff --git a/TESTING/EIG/cdrvsg.f b/TESTING/EIG/cdrvsg.f index 31d6fb67..daac3882 100644 --- a/TESTING/EIG/cdrvsg.f +++ b/TESTING/EIG/cdrvsg.f @@ -172,15 +172,13 @@ *> The number of sizes of matrices to use. If it is zero, *> CDRVSG does nothing. It must be at least zero. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> NN INTEGER array, dimension (NSIZES) *> An array containing the sizes to be used for the matrices. *> Zero values will be skipped. The values must be at least *> zero. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> NTYPES INTEGER *> The number of elements in DOTYPE. If it is zero, CDRVSG *> does nothing. It must be at least zero. If it is MAXTYP+1 @@ -189,8 +187,7 @@ *> is only useful if DOTYPE(1:MAXTYP) is .FALSE. and *> DOTYPE(MAXTYP+1) is .TRUE. . *> Not modified. -*> \endverbatim -*> \verbatim +*> *> DOTYPE LOGICAL array, dimension (NTYPES) *> If DOTYPE(j) is .TRUE., then for each size in NN a *> matrix of that size and of type j will be generated. @@ -200,8 +197,7 @@ *> than MAXTYP, DOTYPE(MAXTYP+1) through DOTYPE(NTYPES) *> will be ignored. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> ISEED INTEGER array, dimension (4) *> On entry ISEED specifies the seed of the random number *> generator. The array elements should be between 0 and 4095; @@ -213,8 +209,7 @@ *> next call to CDRVSG to continue the same random number *> sequence. *> Modified. -*> \endverbatim -*> \verbatim +*> *> THRESH REAL *> A test will count as "failed" if the "error", computed as *> described above, exceeds THRESH. Note that the error @@ -223,118 +218,98 @@ *> it should not depend on the precision (single vs. double) *> or the size of the matrix. It must be at least zero. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> NOUNIT INTEGER *> The FORTRAN unit number for printing out error messages *> (e.g., if a routine returns IINFO not equal to 0.) *> Not modified. -*> \endverbatim -*> \verbatim +*> *> A COMPLEX array, dimension (LDA , max(NN)) *> Used to hold the matrix whose eigenvalues are to be *> computed. On exit, A contains the last matrix actually *> used. *> Modified. -*> \endverbatim -*> \verbatim +*> *> LDA INTEGER *> The leading dimension of A. It must be at *> least 1 and at least max( NN ). *> Not modified. -*> \endverbatim -*> \verbatim +*> *> B COMPLEX array, dimension (LDB , max(NN)) *> Used to hold the Hermitian positive definite matrix for *> the generailzed problem. *> On exit, B contains the last matrix actually *> used. *> Modified. -*> \endverbatim -*> \verbatim +*> *> LDB INTEGER *> The leading dimension of B. It must be at *> least 1 and at least max( NN ). *> Not modified. -*> \endverbatim -*> \verbatim +*> *> D REAL array, dimension (max(NN)) *> The eigenvalues of A. On exit, the eigenvalues in D *> correspond with the matrix in A. *> Modified. -*> \endverbatim -*> \verbatim +*> *> Z COMPLEX array, dimension (LDZ, max(NN)) *> The matrix of eigenvectors. *> Modified. -*> \endverbatim -*> \verbatim +*> *> LDZ INTEGER *> The leading dimension of ZZ. It must be at least 1 and *> at least max( NN ). *> Not modified. -*> \endverbatim -*> \verbatim +*> *> AB COMPLEX array, dimension (LDA, max(NN)) *> Workspace. *> Modified. -*> \endverbatim -*> \verbatim +*> *> BB COMPLEX array, dimension (LDB, max(NN)) *> Workspace. *> Modified. -*> \endverbatim -*> \verbatim +*> *> AP COMPLEX array, dimension (max(NN)**2) *> Workspace. *> Modified. -*> \endverbatim -*> \verbatim +*> *> BP COMPLEX array, dimension (max(NN)**2) *> Workspace. *> Modified. -*> \endverbatim -*> \verbatim +*> *> WORK COMPLEX array, dimension (NWORK) *> Workspace. *> Modified. -*> \endverbatim -*> \verbatim +*> *> NWORK INTEGER *> The number of entries in WORK. This must be at least *> 2*N + N**2 where N = max( NN(j), 2 ). *> Not modified. -*> \endverbatim -*> \verbatim +*> *> RWORK REAL array, dimension (LRWORK) *> Workspace. *> Modified. -*> \endverbatim -*> \verbatim +*> *> LRWORK INTEGER *> The number of entries in RWORK. This must be at least *> max( 7*N, 1 + 4*N + 2*N*lg(N) + 3*N**2 ) where *> N = max( NN(j) ) and lg( N ) = smallest integer k such *> that 2**k >= N . *> Not modified. -*> \endverbatim -*> \verbatim +*> *> IWORK INTEGER array, dimension (LIWORK)) *> Workspace. *> Modified. -*> \endverbatim -*> \verbatim +*> *> LIWORK INTEGER *> The number of entries in IWORK. This must be at least *> 2 + 5*max( NN(j) ). *> Not modified. -*> \endverbatim -*> \verbatim +*> *> RESULT REAL array, dimension (70) *> The values computed by the 70 tests described above. *> Modified. -*> \endverbatim -*> \verbatim +*> *> INFO INTEGER *> If 0, then everything ran OK. *> -1: NSIZES < 0 @@ -350,11 +325,9 @@ *> CHPGVD, CHEGVX, CHPGVX, CHBGVX returns an error code, *> the absolute value of it is returned. *> Modified. -*> \endverbatim -*> \verbatim +*> *>----------------------------------------------------------------------- -*> \endverbatim -*> \verbatim +*> *> Some Local Variables and Parameters: *> ---- ----- --------- --- ---------- *> ZERO, ONE Real 0 and 1. @@ -368,8 +341,7 @@ *> so far (computed by SLAFTS). *> COND, IMODE Values to be passed to the matrix generators. *> ANORM Norm of A; passed to matrix generators. -*> \endverbatim -*> \verbatim +*> *> OVFL, UNFL Overflow and underflow thresholds. *> ULP, ULPINV Finest relative precision and its inverse. *> RTOVFL, RTUNFL Square roots of the previous 2 values. diff --git a/TESTING/EIG/cdrvst.f b/TESTING/EIG/cdrvst.f index 837341fe..94949b38 100644 --- a/TESTING/EIG/cdrvst.f +++ b/TESTING/EIG/cdrvst.f @@ -141,15 +141,13 @@ *> The number of sizes of matrices to use. If it is zero, *> CDRVST does nothing. It must be at least zero. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> NN INTEGER array, dimension (NSIZES) *> An array containing the sizes to be used for the matrices. *> Zero values will be skipped. The values must be at least *> zero. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> NTYPES INTEGER *> The number of elements in DOTYPE. If it is zero, CDRVST *> does nothing. It must be at least zero. If it is MAXTYP+1 @@ -158,8 +156,7 @@ *> is only useful if DOTYPE(1:MAXTYP) is .FALSE. and *> DOTYPE(MAXTYP+1) is .TRUE. . *> Not modified. -*> \endverbatim -*> \verbatim +*> *> DOTYPE LOGICAL array, dimension (NTYPES) *> If DOTYPE(j) is .TRUE., then for each size in NN a *> matrix of that size and of type j will be generated. @@ -169,8 +166,7 @@ *> than MAXTYP, DOTYPE(MAXTYP+1) through DOTYPE(NTYPES) *> will be ignored. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> ISEED INTEGER array, dimension (4) *> On entry ISEED specifies the seed of the random number *> generator. The array elements should be between 0 and 4095; @@ -182,8 +178,7 @@ *> next call to CDRVST to continue the same random number *> sequence. *> Modified. -*> \endverbatim -*> \verbatim +*> *> THRESH REAL *> A test will count as "failed" if the "error", computed as *> described above, exceeds THRESH. Note that the error @@ -192,121 +187,99 @@ *> it should not depend on the precision (single vs. double) *> or the size of the matrix. It must be at least zero. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> NOUNIT INTEGER *> The FORTRAN unit number for printing out error messages *> (e.g., if a routine returns IINFO not equal to 0.) *> Not modified. -*> \endverbatim -*> \verbatim +*> *> A COMPLEX array, dimension (LDA , max(NN)) *> Used to hold the matrix whose eigenvalues are to be *> computed. On exit, A contains the last matrix actually *> used. *> Modified. -*> \endverbatim -*> \verbatim +*> *> LDA INTEGER *> The leading dimension of A. It must be at *> least 1 and at least max( NN ). *> Not modified. -*> \endverbatim -*> \verbatim +*> *> D1 REAL array, dimension (max(NN)) *> The eigenvalues of A, as computed by CSTEQR simlutaneously *> with Z. On exit, the eigenvalues in D1 correspond with the *> matrix in A. *> Modified. -*> \endverbatim -*> \verbatim +*> *> D2 REAL array, dimension (max(NN)) *> The eigenvalues of A, as computed by CSTEQR if Z is not *> computed. On exit, the eigenvalues in D2 correspond with *> the matrix in A. *> Modified. -*> \endverbatim -*> \verbatim +*> *> D3 REAL array, dimension (max(NN)) *> The eigenvalues of A, as computed by SSTERF. On exit, the *> eigenvalues in D3 correspond with the matrix in A. *> Modified. -*> \endverbatim -*> \verbatim +*> *> WA1 REAL array, dimension -*> \endverbatim -*> \verbatim +*> *> WA2 REAL array, dimension -*> \endverbatim -*> \verbatim +*> *> WA3 REAL array, dimension -*> \endverbatim -*> \verbatim +*> *> U COMPLEX array, dimension (LDU, max(NN)) *> The unitary matrix computed by CHETRD + CUNGC3. *> Modified. -*> \endverbatim -*> \verbatim +*> *> LDU INTEGER *> The leading dimension of U, Z, and V. It must be at *> least 1 and at least max( NN ). *> Not modified. -*> \endverbatim -*> \verbatim +*> *> V COMPLEX array, dimension (LDU, max(NN)) *> The Housholder vectors computed by CHETRD in reducing A to *> tridiagonal form. *> Modified. -*> \endverbatim -*> \verbatim +*> *> TAU COMPLEX array, dimension (max(NN)) *> The Householder factors computed by CHETRD in reducing A *> to tridiagonal form. *> Modified. -*> \endverbatim -*> \verbatim +*> *> Z COMPLEX array, dimension (LDU, max(NN)) *> The unitary matrix of eigenvectors computed by CHEEVD, *> CHEEVX, CHPEVD, CHPEVX, CHBEVD, and CHBEVX. *> Modified. -*> \endverbatim -*> \verbatim +*> *> WORK - COMPLEX array of dimension ( LWORK ) *> Workspace. *> Modified. -*> \endverbatim -*> \verbatim +*> *> LWORK - INTEGER *> The number of entries in WORK. This must be at least *> 2*max( NN(j), 2 )**2. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> RWORK REAL array, dimension (3*max(NN)) *> Workspace. *> Modified. -*> \endverbatim -*> \verbatim +*> *> LRWORK - INTEGER *> The number of entries in RWORK. -*> \endverbatim -*> \verbatim +*> *> IWORK INTEGER array, dimension (6*max(NN)) *> Workspace. *> Modified. -*> \endverbatim -*> \verbatim +*> *> LIWORK - INTEGER *> The number of entries in IWORK. -*> \endverbatim -*> \verbatim +*> *> RESULT REAL array, dimension (??) *> The values computed by the tests described above. *> The values are currently limited to 1/ulp, to avoid *> overflow. *> Modified. -*> \endverbatim -*> \verbatim +*> *> INFO INTEGER *> If 0, then everything ran OK. *> -1: NSIZES < 0 @@ -320,11 +293,9 @@ *> or SORMC2 returns an error code, the *> absolute value of it is returned. *> Modified. -*> \endverbatim -*> \verbatim +*> *>----------------------------------------------------------------------- -*> \endverbatim -*> \verbatim +*> *> Some Local Variables and Parameters: *> ---- ----- --------- --- ---------- *> ZERO, ONE Real 0 and 1. @@ -338,8 +309,7 @@ *> so far (computed by SLAFTS). *> COND, IMODE Values to be passed to the matrix generators. *> ANORM Norm of A; passed to matrix generators. -*> \endverbatim -*> \verbatim +*> *> OVFL, UNFL Overflow and underflow thresholds. *> ULP, ULPINV Finest relative precision and its inverse. *> RTOVFL, RTUNFL Square roots of the previous 2 values. diff --git a/TESTING/EIG/cdrvsx.f b/TESTING/EIG/cdrvsx.f index 9aafeed7..659c26f9 100644 --- a/TESTING/EIG/cdrvsx.f +++ b/TESTING/EIG/cdrvsx.f @@ -392,11 +392,9 @@ *> <0, input parameter -INFO is incorrect *> >0, CLATMR, CLATMS, CLATME or CGET24 returned an error *> code and INFO is its absolute value -*> \endverbatim -*> \verbatim +*> *>----------------------------------------------------------------------- -*> \endverbatim -*> \verbatim +*> *> Some Local Variables and Parameters: *> ---- ----- --------- --- ---------- *> ZERO, ONE Real 0 and 1. @@ -406,8 +404,7 @@ *> COND, CONDS, *> IMODE Values to be passed to the matrix generators. *> ANORM Norm of A; passed to matrix generators. -*> \endverbatim -*> \verbatim +*> *> OVFL, UNFL Overflow and underflow thresholds. *> ULP, ULPINV Finest relative precision and its inverse. *> RTULP, RTULPI Square roots of the previous 4 values. diff --git a/TESTING/EIG/cdrvvx.f b/TESTING/EIG/cdrvvx.f index 3a15d8f6..7aa0ce03 100644 --- a/TESTING/EIG/cdrvvx.f +++ b/TESTING/EIG/cdrvvx.f @@ -450,15 +450,12 @@ *> If <0, then input paramter -INFO is incorrect. *> If >0, CLATMR, CLATMS, CLATME or CGET23 returned an error *> code, and INFO is its absolute value. -*> \endverbatim -*> \verbatim +*> *>----------------------------------------------------------------------- -*> \endverbatim -*> \verbatim +*> *> Some Local Variables and Parameters: *> ---- ----- --------- --- ---------- -*> \endverbatim -*> \verbatim +*> *> ZERO, ONE Real 0 and 1. *> MAXTYP The number of types defined. *> NMAX Largest value in NN or 12. @@ -466,13 +463,11 @@ *> COND, CONDS, *> IMODE Values to be passed to the matrix generators. *> ANORM Norm of A; passed to matrix generators. -*> \endverbatim -*> \verbatim +*> *> OVFL, UNFL Overflow and underflow thresholds. *> ULP, ULPINV Finest relative precision and its inverse. *> RTULP, RTULPI Square roots of the previous 4 values. -*> \endverbatim -*> \verbatim +*> *> The following four arrays decode JTYPE: *> KTYPE(j) The general type (1-10) for type "j". *> KMODE(j) The MODE value to be passed to the matrix diff --git a/TESTING/EIG/cgsvts.f b/TESTING/EIG/cgsvts.f index 89bcd753..b981bad7 100644 --- a/TESTING/EIG/cgsvts.f +++ b/TESTING/EIG/cgsvts.f @@ -140,8 +140,7 @@ *> \param[out] BETA *> \verbatim *> BETA is REAL array, dimension (N) -*> \endverbatim -*> \verbatim +*> *> The generalized singular value pairs of A and B, the *> ``diagonal'' matrices D1 and D2 are constructed from *> ALPHA and BETA, see subroutine CGGSVD for details. diff --git a/TESTING/EIG/chet21.f b/TESTING/EIG/chet21.f index 11f3fbab..c5b0f054 100644 --- a/TESTING/EIG/chet21.f +++ b/TESTING/EIG/chet21.f @@ -68,12 +68,10 @@ *> Specifies the type of tests to be performed. *> 1: U expressed as a dense unitary matrix: *> RESULT(1) = | A - U S U* | / ( |A| n ulp ) *andC> RESULT(2) = | I - UU* | / ( n ulp ) -*> \endverbatim -*> \verbatim +*> *> 2: U expressed as a product V of Housholder transformations: *> RESULT(1) = | A - V S V* | / ( |A| n ulp ) -*> \endverbatim -*> \verbatim +*> *> 3: U expressed both as a dense unitary matrix and *> as a product of Housholder transformations: *> RESULT(1) = | I - UV* | / ( n ulp ) diff --git a/TESTING/EIG/chet22.f b/TESTING/EIG/chet22.f index db8a768b..b6e5b495 100644 --- a/TESTING/EIG/chet22.f +++ b/TESTING/EIG/chet22.f @@ -54,104 +54,88 @@ *> Specifies the type of tests to be performed. *> 1: U expressed as a dense orthogonal matrix: *> RESULT(1) = | A - U S U' | / ( |A| n ulp ) *andC> RESULT(2) = | I - UU' | / ( n ulp ) -*> \endverbatim -*> \verbatim +*> *> UPLO CHARACTER *> If UPLO='U', the upper triangle of A will be used and the *> (strictly) lower triangle will not be referenced. If *> UPLO='L', the lower triangle of A will be used and the *> (strictly) upper triangle will not be referenced. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> N INTEGER *> The size of the matrix. If it is zero, CHET22 does nothing. *> It must be at least zero. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> M INTEGER *> The number of columns of U. If it is zero, CHET22 does *> nothing. It must be at least zero. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> KBAND INTEGER *> The bandwidth of the matrix. It may only be zero or one. *> If zero, then S is diagonal, and E is not referenced. If *> one, then S is symmetric tri-diagonal. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> A COMPLEX array, dimension (LDA , N) *> The original (unfactored) matrix. It is assumed to be *> symmetric, and only the upper (UPLO='U') or only the lower *> (UPLO='L') will be referenced. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> LDA INTEGER *> The leading dimension of A. It must be at least 1 *> and at least N. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> D REAL array, dimension (N) *> The diagonal of the (symmetric tri-) diagonal matrix. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> E REAL array, dimension (N) *> The off-diagonal of the (symmetric tri-) diagonal matrix. *> E(1) is ignored, E(2) is the (1,2) and (2,1) element, etc. *> Not referenced if KBAND=0. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> U COMPLEX array, dimension (LDU, N) *> If ITYPE=1, this contains the orthogonal matrix in *> the decomposition, expressed as a dense matrix. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> LDU INTEGER *> The leading dimension of U. LDU must be at least N and *> at least 1. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> V COMPLEX array, dimension (LDV, N) *> If ITYPE=2 or 3, the lower triangle of this array contains *> the Householder vectors used to describe the orthogonal *> matrix in the decomposition. If ITYPE=1, then it is not *> referenced. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> LDV INTEGER *> The leading dimension of V. LDV must be at least N and *> at least 1. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> TAU COMPLEX array, dimension (N) *> If ITYPE >= 2, then TAU(j) is the scalar factor of *> v(j) v(j)' in the Householder transformation H(j) of *> the product U = H(1)...H(n-2) *> If ITYPE < 2, then TAU is not referenced. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> WORK COMPLEX array, dimension (2*N**2) *> Workspace. *> Modified. -*> \endverbatim -*> \verbatim +*> *> RWORK REAL array, dimension (N) *> Workspace. *> Modified. -*> \endverbatim -*> \verbatim +*> *> RESULT REAL array, dimension (2) *> The values computed by the two tests described above. The *> values are currently limited to 1/ulp, to avoid overflow. diff --git a/TESTING/EIG/chpt21.f b/TESTING/EIG/chpt21.f index 337af63f..9ad6d7c6 100644 --- a/TESTING/EIG/chpt21.f +++ b/TESTING/EIG/chpt21.f @@ -93,12 +93,10 @@ *> Specifies the type of tests to be performed. *> 1: U expressed as a dense unitary matrix: *> RESULT(1) = | A - U S U* | / ( |A| n ulp ) *andC> RESULT(2) = | I - UU* | / ( n ulp ) -*> \endverbatim -*> \verbatim +*> *> 2: U expressed as a product V of Housholder transformations: *> RESULT(1) = | A - V S V* | / ( |A| n ulp ) -*> \endverbatim -*> \verbatim +*> *> 3: U expressed both as a dense unitary matrix and *> as a product of Housholder transformations: *> RESULT(1) = | I - UV* | / ( n ulp ) diff --git a/TESTING/EIG/chst01.f b/TESTING/EIG/chst01.f index ad373917..a007da90 100644 --- a/TESTING/EIG/chst01.f +++ b/TESTING/EIG/chst01.f @@ -57,8 +57,7 @@ *> \param[in] IHI *> \verbatim *> IHI is INTEGER -*> \endverbatim -*> \verbatim +*> *> A is assumed to be upper triangular in rows and columns *> 1:ILO-1 and IHI+1:N, so Q differs from the identity only in *> rows and columns ILO+1:IHI. diff --git a/TESTING/EIG/clarhs.f b/TESTING/EIG/clarhs.f index f861dc4e..d7fb11bf 100644 --- a/TESTING/EIG/clarhs.f +++ b/TESTING/EIG/clarhs.f @@ -118,12 +118,10 @@ *> KU is INTEGER *> Used only if A is a general band matrix or if A is *> triangular. -*> \endverbatim -*> \verbatim +*> *> If PATH = xGB, specifies the number of superdiagonals of A, *> and 0 <= KU <= N-1. -*> \endverbatim -*> \verbatim +*> *> If PATH = xTR, xTP, or xTB, specifies whether or not the *> matrix has unit diagonal: *> = 1: matrix has non-unit diagonal (default) diff --git a/TESTING/EIG/clatm4.f b/TESTING/EIG/clatm4.f index 17ee2dc9..c4a39940 100644 --- a/TESTING/EIG/clatm4.f +++ b/TESTING/EIG/clatm4.f @@ -49,8 +49,7 @@ *> If ITYPE < 0, then type abs(ITYPE) is generated and then *> swapped end for end (A(I,J) := A'(N-J,N-I).) See also *> the description of AMAGN and RSIGN. -*> \endverbatim -*> \verbatim +*> *> Special types: *> = 0: the zero matrix. *> = 1: the identity. @@ -59,8 +58,7 @@ *> followed by a k x k identity block, where k=(N-1)/2. *> If N is even, then k=(N-2)/2, and a zero diagonal entry *> is tacked onto the end. -*> \endverbatim -*> \verbatim +*> *> Diagonal types. The diagonal consists of NZ1 zeros, then *> k=N-NZ1-NZ2 nonzeros. The subdiagonal is zero. ITYPE *> specifies the nonzero diagonal entries as follows: diff --git a/TESTING/EIG/clctsx.f b/TESTING/EIG/clctsx.f index c2208012..01ed84f7 100644 --- a/TESTING/EIG/clctsx.f +++ b/TESTING/EIG/clctsx.f @@ -38,8 +38,7 @@ *> \param[in] BETA *> \verbatim *> BETA is COMPLEX -*> \endverbatim -*> \verbatim +*> *> parameters to decide whether the pair (ALPHA, BETA) is *> selected. *> \endverbatim diff --git a/TESTING/EIG/csbmv.f b/TESTING/EIG/csbmv.f index a0b6558d..5cfff04f 100644 --- a/TESTING/EIG/csbmv.f +++ b/TESTING/EIG/csbmv.f @@ -43,36 +43,29 @@ *> On entry, UPLO specifies whether the upper or lower *> triangular part of the band matrix A is being supplied as *> follows: -*> \endverbatim -*> \verbatim +*> *> UPLO = 'U' or 'u' The upper triangular part of A is *> being supplied. -*> \endverbatim -*> \verbatim +*> *> UPLO = 'L' or 'l' The lower triangular part of A is *> being supplied. -*> \endverbatim -*> \verbatim +*> *> Unchanged on exit. -*> \endverbatim -*> \verbatim +*> *> N - INTEGER *> On entry, N specifies the order of the matrix A. *> N must be at least zero. *> Unchanged on exit. -*> \endverbatim -*> \verbatim +*> *> K - INTEGER *> On entry, K specifies the number of super-diagonals of the *> matrix A. K must satisfy 0 .le. K. *> Unchanged on exit. -*> \endverbatim -*> \verbatim +*> *> ALPHA - COMPLEX *> On entry, ALPHA specifies the scalar alpha. *> Unchanged on exit. -*> \endverbatim -*> \verbatim +*> *> A - COMPLEX array, dimension( LDA, N ) *> Before entry with UPLO = 'U' or 'u', the leading ( k + 1 ) *> by n part of the array A must contain the upper triangular @@ -84,16 +77,14 @@ *> The following program segment will transfer the upper *> triangular part of a symmetric band matrix from conventional *> full matrix storage to band storage: -*> \endverbatim -*> \verbatim +*> *> DO 20, J = 1, N *> M = K + 1 - J *> DO 10, I = MAX( 1, J - K ), J *> A( M + I, J ) = matrix( I, J ) *> 10 CONTINUE *> 20 CONTINUE -*> \endverbatim -*> \verbatim +*> *> Before entry with UPLO = 'L' or 'l', the leading ( k + 1 ) *> by n part of the array A must contain the lower triangular *> band part of the symmetric matrix, supplied column by @@ -104,50 +95,42 @@ *> The following program segment will transfer the lower *> triangular part of a symmetric band matrix from conventional *> full matrix storage to band storage: -*> \endverbatim -*> \verbatim +*> *> DO 20, J = 1, N *> M = 1 - J *> DO 10, I = J, MIN( N, J + K ) *> A( M + I, J ) = matrix( I, J ) *> 10 CONTINUE *> 20 CONTINUE -*> \endverbatim -*> \verbatim +*> *> Unchanged on exit. -*> \endverbatim -*> \verbatim +*> *> LDA - INTEGER *> On entry, LDA specifies the first dimension of A as declared *> in the calling (sub) program. LDA must be at least *> ( k + 1 ). *> Unchanged on exit. -*> \endverbatim -*> \verbatim +*> *> X - COMPLEX array, dimension at least *> ( 1 + ( N - 1 )*abs( INCX ) ). *> Before entry, the incremented array X must contain the *> vector x. *> Unchanged on exit. -*> \endverbatim -*> \verbatim +*> *> INCX - INTEGER *> On entry, INCX specifies the increment for the elements of *> X. INCX must not be zero. *> Unchanged on exit. -*> \endverbatim -*> \verbatim +*> *> BETA - COMPLEX *> On entry, BETA specifies the scalar beta. *> Unchanged on exit. -*> \endverbatim -*> \verbatim +*> *> Y - COMPLEX array, dimension at least *> ( 1 + ( N - 1 )*abs( INCY ) ). *> Before entry, the incremented array Y must contain the *> vector y. On exit, Y is overwritten by the updated vector y. -*> \endverbatim -*> \verbatim +*> *> INCY - INTEGER *> On entry, INCY specifies the increment for the elements of *> Y. INCY must not be zero. diff --git a/TESTING/EIG/dchkbb.f b/TESTING/EIG/dchkbb.f index e040dc7b..ca9b5375 100644 --- a/TESTING/EIG/dchkbb.f +++ b/TESTING/EIG/dchkbb.f @@ -310,11 +310,9 @@ *> \verbatim *> INFO is INTEGER *> If 0, then everything ran OK. -*> \endverbatim -*> \verbatim +*> *>----------------------------------------------------------------------- -*> \endverbatim -*> \verbatim +*> *> Some Local Variables and Parameters: *> ---- ----- --------- --- ---------- *> ZERO, ONE Real 0 and 1. @@ -328,8 +326,7 @@ *> so far. *> COND, IMODE Values to be passed to the matrix generators. *> ANORM Norm of A; passed to matrix generators. -*> \endverbatim -*> \verbatim +*> *> OVFL, UNFL Overflow and underflow thresholds. *> ULP, ULPINV Finest relative precision and its inverse. *> RTOVFL, RTUNFL Square roots of the previous 2 values. diff --git a/TESTING/EIG/dchkbd.f b/TESTING/EIG/dchkbd.f index 760b6bf9..32c8dbd3 100644 --- a/TESTING/EIG/dchkbd.f +++ b/TESTING/EIG/dchkbd.f @@ -388,15 +388,12 @@ *> If DLATMR, SLATMS, DGEBRD, DORGBR, or DBDSQR, *> returns an error code, the *> absolute value of it is returned. -*> \endverbatim -*> \verbatim +*> *>----------------------------------------------------------------------- -*> \endverbatim -*> \verbatim +*> *> Some Local Variables and Parameters: *> ---- ----- --------- --- ---------- -*> \endverbatim -*> \verbatim +*> *> ZERO, ONE Real 0 and 1. *> MAXTYP The number of types defined. *> NTEST The number of tests performed, or which can @@ -409,13 +406,11 @@ *> NFAIL The number of tests which have exceeded THRESH *> COND, IMODE Values to be passed to the matrix generators. *> ANORM Norm of A; passed to matrix generators. -*> \endverbatim -*> \verbatim +*> *> OVFL, UNFL Overflow and underflow thresholds. *> RTOVFL, RTUNFL Square roots of the previous 2 values. *> ULP, ULPINV Finest relative precision and its inverse. -*> \endverbatim -*> \verbatim +*> *> The following four arrays decode JTYPE: *> KTYPE(j) The general type (1-10) for type "j". *> KMODE(j) The MODE value to be passed to the matrix diff --git a/TESTING/EIG/dchkgg.f b/TESTING/EIG/dchkgg.f index 9c9202d7..15a375ab 100644 --- a/TESTING/EIG/dchkgg.f +++ b/TESTING/EIG/dchkgg.f @@ -423,8 +423,7 @@ *> \param[out] BETA1 *> \verbatim *> BETA1 is DOUBLE PRECISION array, dimension (max(NN)) -*> \endverbatim -*> \verbatim +*> *> The generalized eigenvalues of (A,B) computed by DHGEQZ *> when Q, Z, and the full Schur matrices are computed. *> On exit, ( ALPHR1(k)+ALPHI1(k)*i ) / BETA1(k) is the k-th diff --git a/TESTING/EIG/dchkhs.f b/TESTING/EIG/dchkhs.f index d1cc39ab..0657f682 100644 --- a/TESTING/EIG/dchkhs.f +++ b/TESTING/EIG/dchkhs.f @@ -166,15 +166,13 @@ *> The number of sizes of matrices to use. If it is zero, *> DCHKHS does nothing. It must be at least zero. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> NN - INTEGER array, dimension (NSIZES) *> An array containing the sizes to be used for the matrices. *> Zero values will be skipped. The values must be at least *> zero. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> NTYPES - INTEGER *> The number of elements in DOTYPE. If it is zero, DCHKHS *> does nothing. It must be at least zero. If it is MAXTYP+1 @@ -183,8 +181,7 @@ *> is only useful if DOTYPE(1:MAXTYP) is .FALSE. and *> DOTYPE(MAXTYP+1) is .TRUE. . *> Not modified. -*> \endverbatim -*> \verbatim +*> *> DOTYPE - LOGICAL array, dimension (NTYPES) *> If DOTYPE(j) is .TRUE., then for each size in NN a *> matrix of that size and of type j will be generated. @@ -194,8 +191,7 @@ *> than MAXTYP, DOTYPE(MAXTYP+1) through DOTYPE(NTYPES) *> will be ignored. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> ISEED - INTEGER array, dimension (4) *> On entry ISEED specifies the seed of the random number *> generator. The array elements should be between 0 and 4095; @@ -207,8 +203,7 @@ *> next call to DCHKHS to continue the same random number *> sequence. *> Modified. -*> \endverbatim -*> \verbatim +*> *> THRESH - DOUBLE PRECISION *> A test will count as "failed" if the "error", computed as *> described above, exceeds THRESH. Note that the error @@ -217,75 +212,63 @@ *> it should not depend on the precision (single vs. double) *> or the size of the matrix. It must be at least zero. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> NOUNIT - INTEGER *> The FORTRAN unit number for printing out error messages *> (e.g., if a routine returns IINFO not equal to 0.) *> Not modified. -*> \endverbatim -*> \verbatim +*> *> A - DOUBLE PRECISION array, dimension (LDA,max(NN)) *> Used to hold the matrix whose eigenvalues are to be *> computed. On exit, A contains the last matrix actually *> used. *> Modified. -*> \endverbatim -*> \verbatim +*> *> LDA - INTEGER *> The leading dimension of A, H, T1 and T2. It must be at *> least 1 and at least max( NN ). *> Not modified. -*> \endverbatim -*> \verbatim +*> *> H - DOUBLE PRECISION array, dimension (LDA,max(NN)) *> The upper hessenberg matrix computed by DGEHRD. On exit, *> H contains the Hessenberg form of the matrix in A. *> Modified. -*> \endverbatim -*> \verbatim +*> *> T1 - DOUBLE PRECISION array, dimension (LDA,max(NN)) *> The Schur (="quasi-triangular") matrix computed by DHSEQR *> if Z is computed. On exit, T1 contains the Schur form of *> the matrix in A. *> Modified. -*> \endverbatim -*> \verbatim +*> *> T2 - DOUBLE PRECISION array, dimension (LDA,max(NN)) *> The Schur matrix computed by DHSEQR when Z is not computed. *> This should be identical to T1. *> Modified. -*> \endverbatim -*> \verbatim +*> *> LDU - INTEGER *> The leading dimension of U, Z, UZ and UU. It must be at *> least 1 and at least max( NN ). *> Not modified. -*> \endverbatim -*> \verbatim +*> *> U - DOUBLE PRECISION array, dimension (LDU,max(NN)) *> The orthogonal matrix computed by DGEHRD. *> Modified. -*> \endverbatim -*> \verbatim +*> *> Z - DOUBLE PRECISION array, dimension (LDU,max(NN)) *> The orthogonal matrix computed by DHSEQR. *> Modified. -*> \endverbatim -*> \verbatim +*> *> UZ - DOUBLE PRECISION array, dimension (LDU,max(NN)) *> The product of U times Z. *> Modified. -*> \endverbatim -*> \verbatim +*> *> WR1 - DOUBLE PRECISION array, dimension (max(NN)) *> WI1 - DOUBLE PRECISION array, dimension (max(NN)) *> The real and imaginary parts of the eigenvalues of A, *> as computed when Z is computed. *> On exit, WR1 + WI1*i are the eigenvalues of the matrix in A. *> Modified. -*> \endverbatim -*> \verbatim +*> *> WR3 - DOUBLE PRECISION array, dimension (max(NN)) *> WI3 - DOUBLE PRECISION array, dimension (max(NN)) *> Like WR1, WI1, these arrays contain the eigenvalues of A, @@ -294,72 +277,60 @@ *> Schur form than is necessary for computing the *> eigenvalues. *> Modified. -*> \endverbatim -*> \verbatim +*> *> EVECTL - DOUBLE PRECISION array, dimension (LDU,max(NN)) *> The (upper triangular) left eigenvector matrix for the *> matrix in T1. For complex conjugate pairs, the real part *> is stored in one row and the imaginary part in the next. *> Modified. -*> \endverbatim -*> \verbatim +*> *> EVEZTR - DOUBLE PRECISION array, dimension (LDU,max(NN)) *> The (upper triangular) right eigenvector matrix for the *> matrix in T1. For complex conjugate pairs, the real part *> is stored in one column and the imaginary part in the next. *> Modified. -*> \endverbatim -*> \verbatim +*> *> EVECTY - DOUBLE PRECISION array, dimension (LDU,max(NN)) *> The left eigenvector matrix for the *> matrix in H. For complex conjugate pairs, the real part *> is stored in one row and the imaginary part in the next. *> Modified. -*> \endverbatim -*> \verbatim +*> *> EVECTX - DOUBLE PRECISION array, dimension (LDU,max(NN)) *> The right eigenvector matrix for the *> matrix in H. For complex conjugate pairs, the real part *> is stored in one column and the imaginary part in the next. *> Modified. -*> \endverbatim -*> \verbatim +*> *> UU - DOUBLE PRECISION array, dimension (LDU,max(NN)) *> Details of the orthogonal matrix computed by DGEHRD. *> Modified. -*> \endverbatim -*> \verbatim +*> *> TAU - DOUBLE PRECISION array, dimension(max(NN)) *> Further details of the orthogonal matrix computed by DGEHRD. *> Modified. -*> \endverbatim -*> \verbatim +*> *> WORK - DOUBLE PRECISION array, dimension (NWORK) *> Workspace. *> Modified. -*> \endverbatim -*> \verbatim +*> *> NWORK - INTEGER *> The number of entries in WORK. NWORK >= 4*NN(j)*NN(j) + 2. -*> \endverbatim -*> \verbatim +*> *> IWORK - INTEGER array, dimension (max(NN)) *> Workspace. *> Modified. -*> \endverbatim -*> \verbatim +*> *> SELECT - LOGICAL array, dimension (max(NN)) *> Workspace. *> Modified. -*> \endverbatim -*> \verbatim +*> *> RESULT - DOUBLE PRECISION array, dimension (14) *> The values computed by the fourteen tests described above. *> The values are currently limited to 1/ulp, to avoid *> overflow. *> Modified. -*> \endverbatim -*> \verbatim +*> *> INFO - INTEGER *> If 0, then everything ran OK. *> -1: NSIZES < 0 @@ -376,15 +347,12 @@ *> If >2, then 30*N iterations were not enough to find an *> eigenvalue or to decompose the problem. *> Modified. -*> \endverbatim -*> \verbatim +*> *>----------------------------------------------------------------------- -*> \endverbatim -*> \verbatim +*> *> Some Local Variables and Parameters: *> ---- ----- --------- --- ---------- -*> \endverbatim -*> \verbatim +*> *> ZERO, ONE Real 0 and 1. *> MAXTYP The number of types defined. *> MTEST The number of tests defined: care must be taken @@ -403,14 +371,12 @@ *> COND, CONDS, *> IMODE Values to be passed to the matrix generators. *> ANORM Norm of A; passed to matrix generators. -*> \endverbatim -*> \verbatim +*> *> OVFL, UNFL Overflow and underflow thresholds. *> ULP, ULPINV Finest relative precision and its inverse. *> RTOVFL, RTUNFL, *> RTULP, RTULPI Square roots of the previous 4 values. -*> \endverbatim -*> \verbatim +*> *> The following four arrays decode JTYPE: *> KTYPE(j) The general type (1-10) for type "j". *> KMODE(j) The MODE value to be passed to the matrix diff --git a/TESTING/EIG/dchksb.f b/TESTING/EIG/dchksb.f index d3b12bd8..cda3fae6 100644 --- a/TESTING/EIG/dchksb.f +++ b/TESTING/EIG/dchksb.f @@ -249,11 +249,9 @@ *> \verbatim *> INFO is INTEGER *> If 0, then everything ran OK. -*> \endverbatim -*> \verbatim +*> *>----------------------------------------------------------------------- -*> \endverbatim -*> \verbatim +*> *> Some Local Variables and Parameters: *> ---- ----- --------- --- ---------- *> ZERO, ONE Real 0 and 1. @@ -267,8 +265,7 @@ *> so far. *> COND, IMODE Values to be passed to the matrix generators. *> ANORM Norm of A; passed to matrix generators. -*> \endverbatim -*> \verbatim +*> *> OVFL, UNFL Overflow and underflow thresholds. *> ULP, ULPINV Finest relative precision and its inverse. *> RTOVFL, RTUNFL Square roots of the previous 2 values. diff --git a/TESTING/EIG/dchkst.f b/TESTING/EIG/dchkst.f index c585fe3d..76b5752d 100644 --- a/TESTING/EIG/dchkst.f +++ b/TESTING/EIG/dchkst.f @@ -545,11 +545,9 @@ *> If DLATMR, SLATMS, DSYTRD, DORGTR, DSTEQR, SSTERF, *> or DORMC2 returns an error code, the *> absolute value of it is returned. -*> \endverbatim -*> \verbatim +*> *>----------------------------------------------------------------------- -*> \endverbatim -*> \verbatim +*> *> Some Local Variables and Parameters: *> ---- ----- --------- --- ---------- *> ZERO, ONE Real 0 and 1. @@ -564,8 +562,7 @@ *> so far. *> COND, IMODE Values to be passed to the matrix generators. *> ANORM Norm of A; passed to matrix generators. -*> \endverbatim -*> \verbatim +*> *> OVFL, UNFL Overflow and underflow thresholds. *> ULP, ULPINV Finest relative precision and its inverse. *> RTOVFL, RTUNFL Square roots of the previous 2 values. diff --git a/TESTING/EIG/ddrges.f b/TESTING/EIG/ddrges.f index 8eaec08e..dd96e7c2 100644 --- a/TESTING/EIG/ddrges.f +++ b/TESTING/EIG/ddrges.f @@ -346,8 +346,7 @@ *> \param[out] BETA *> \verbatim *> BETA is DOUBLE PRECISION array, dimension (max(NN)) -*> \endverbatim -*> \verbatim +*> *> The generalized eigenvalues of (A,B) computed by DGGES. *> ( ALPHAR(k)+ALPHAI(k)*i ) / BETA(k) is the k-th *> generalized eigenvalue of A and B. diff --git a/TESTING/EIG/ddrgev.f b/TESTING/EIG/ddrgev.f index 45ba288f..bf138b72 100644 --- a/TESTING/EIG/ddrgev.f +++ b/TESTING/EIG/ddrgev.f @@ -338,8 +338,7 @@ *> \param[out] BETA *> \verbatim *> BETA is DOUBLE PRECISION array, dimension (max(NN)) -*> \endverbatim -*> \verbatim +*> *> The generalized eigenvalues of (A,B) computed by DGGEV. *> ( ALPHAR(k)+ALPHAI(k)*i ) / BETA(k) is the k-th *> generalized eigenvalue of A and B. @@ -358,8 +357,7 @@ *> \param[out] BETA1 *> \verbatim *> BETA1 is DOUBLE PRECISION array, dimension (max(NN)) -*> \endverbatim -*> \verbatim +*> *> Like ALPHAR, ALPHAI, BETA, these arrays contain the *> eigenvalues of A and B, but those computed when DGGEV only *> computes a partial eigendecomposition, i.e. not the diff --git a/TESTING/EIG/ddrgsx.f b/TESTING/EIG/ddrgsx.f index 99e77861..6f2f36a5 100644 --- a/TESTING/EIG/ddrgsx.f +++ b/TESTING/EIG/ddrgsx.f @@ -283,8 +283,7 @@ *> \param[out] BETA *> \verbatim *> BETA is DOUBLE PRECISION array, dimension (NSIZE) -*> \endverbatim -*> \verbatim +*> *> On exit, (ALPHAR + ALPHAI*i)/BETA are the eigenvalues. *> \endverbatim *> diff --git a/TESTING/EIG/ddrgvx.f b/TESTING/EIG/ddrgvx.f index f038d75d..e1486261 100644 --- a/TESTING/EIG/ddrgvx.f +++ b/TESTING/EIG/ddrgvx.f @@ -187,8 +187,7 @@ *> \param[out] BETA *> \verbatim *> BETA is DOUBLE PRECISION array, dimension (NSIZE) -*> \endverbatim -*> \verbatim +*> *> On exit, (ALPHAR + ALPHAI*i)/BETA are the eigenvalues. *> \endverbatim *> diff --git a/TESTING/EIG/ddrves.f b/TESTING/EIG/ddrves.f index 93fa2b6d..6cf863b3 100644 --- a/TESTING/EIG/ddrves.f +++ b/TESTING/EIG/ddrves.f @@ -271,8 +271,7 @@ *> \param[out] WI *> \verbatim *> WI is DOUBLE PRECISION array, dimension (max(NN)) -*> \endverbatim -*> \verbatim +*> *> The real and imaginary parts of the eigenvalues of A. *> On exit, WR + WI*i are the eigenvalues of the matrix in A. *> \endverbatim @@ -285,8 +284,7 @@ *> \param[out] WIT *> \verbatim *> WIT is DOUBLE PRECISION array, dimension (max(NN)) -*> \endverbatim -*> \verbatim +*> *> Like WR, WI, these arrays contain the eigenvalues of A, *> but those computed when DGEES only computes a partial *> eigendecomposition, i.e. not Schur vectors @@ -346,15 +344,12 @@ *> -20: NWORK too small. *> If DLATMR, SLATMS, SLATME or DGEES returns an error code, *> the absolute value of it is returned. -*> \endverbatim -*> \verbatim +*> *>----------------------------------------------------------------------- -*> \endverbatim -*> \verbatim +*> *> Some Local Variables and Parameters: *> ---- ----- --------- --- ---------- -*> \endverbatim -*> \verbatim +*> *> ZERO, ONE Real 0 and 1. *> MAXTYP The number of types defined. *> NMAX Largest value in NN. @@ -362,13 +357,11 @@ *> COND, CONDS, *> IMODE Values to be passed to the matrix generators. *> ANORM Norm of A; passed to matrix generators. -*> \endverbatim -*> \verbatim +*> *> OVFL, UNFL Overflow and underflow thresholds. *> ULP, ULPINV Finest relative precision and its inverse. *> RTULP, RTULPI Square roots of the previous 4 values. -*> \endverbatim -*> \verbatim +*> *> The following four arrays decode JTYPE: *> KTYPE(j) The general type (1-10) for type "j". *> KMODE(j) The MODE value to be passed to the matrix diff --git a/TESTING/EIG/ddrvev.f b/TESTING/EIG/ddrvev.f index b79b1fa0..66054207 100644 --- a/TESTING/EIG/ddrvev.f +++ b/TESTING/EIG/ddrvev.f @@ -266,8 +266,7 @@ *> \param[out] WI *> \verbatim *> WI is DOUBLE PRECISION array, dimension (max(NN)) -*> \endverbatim -*> \verbatim +*> *> The real and imaginary parts of the eigenvalues of A. *> On exit, WR + WI*i are the eigenvalues of the matrix in A. *> \endverbatim @@ -280,8 +279,7 @@ *> \param[out] WI1 *> \verbatim *> WI1 is DOUBLE PRECISION array, dimension (max(NN)) -*> \endverbatim -*> \verbatim +*> *> Like WR, WI, these arrays contain the eigenvalues of A, *> but those computed when DGEEV only computes a partial *> eigendecomposition, i.e. not the eigenvalues and left @@ -363,15 +361,12 @@ *> -23: NWORK too small. *> If DLATMR, SLATMS, SLATME or DGEEV returns an error code, *> the absolute value of it is returned. -*> \endverbatim -*> \verbatim +*> *>----------------------------------------------------------------------- -*> \endverbatim -*> \verbatim +*> *> Some Local Variables and Parameters: *> ---- ----- --------- --- ---------- -*> \endverbatim -*> \verbatim +*> *> ZERO, ONE Real 0 and 1. *> MAXTYP The number of types defined. *> NMAX Largest value in NN. @@ -379,13 +374,11 @@ *> COND, CONDS, *> IMODE Values to be passed to the matrix generators. *> ANORM Norm of A; passed to matrix generators. -*> \endverbatim -*> \verbatim +*> *> OVFL, UNFL Overflow and underflow thresholds. *> ULP, ULPINV Finest relative precision and its inverse. *> RTULP, RTULPI Square roots of the previous 4 values. -*> \endverbatim -*> \verbatim +*> *> The following four arrays decode JTYPE: *> KTYPE(j) The general type (1-10) for type "j". *> KMODE(j) The MODE value to be passed to the matrix diff --git a/TESTING/EIG/ddrvgg.f b/TESTING/EIG/ddrvgg.f index e94d6b95..9ab41fb6 100644 --- a/TESTING/EIG/ddrvgg.f +++ b/TESTING/EIG/ddrvgg.f @@ -363,8 +363,7 @@ *> \param[out] BETA1 *> \verbatim *> BETA1 is DOUBLE PRECISION array, dimension (max(NN)) -*> \endverbatim -*> \verbatim +*> *> The generalized eigenvalues of (A,B) computed by DGEGS. *> ( ALPHR1(k)+ALPHI1(k)*i ) / BETA1(k) is the k-th *> generalized eigenvalue of the matrices in A and B. @@ -383,8 +382,7 @@ *> \param[out] BETA2 *> \verbatim *> BETA2 is DOUBLE PRECISION array, dimension (max(NN)) -*> \endverbatim -*> \verbatim +*> *> The generalized eigenvalues of (A,B) computed by DGEGV. *> ( ALPHR2(k)+ALPHI2(k)*i ) / BETA2(k) is the k-th *> generalized eigenvalue of the matrices in A and B. diff --git a/TESTING/EIG/ddrvsg.f b/TESTING/EIG/ddrvsg.f index 6afe70dc..cfe511ee 100644 --- a/TESTING/EIG/ddrvsg.f +++ b/TESTING/EIG/ddrvsg.f @@ -170,15 +170,13 @@ *> The number of sizes of matrices to use. If it is zero, *> DDRVSG does nothing. It must be at least zero. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> NN INTEGER array, dimension (NSIZES) *> An array containing the sizes to be used for the matrices. *> Zero values will be skipped. The values must be at least *> zero. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> NTYPES INTEGER *> The number of elements in DOTYPE. If it is zero, DDRVSG *> does nothing. It must be at least zero. If it is MAXTYP+1 @@ -187,8 +185,7 @@ *> is only useful if DOTYPE(1:MAXTYP) is .FALSE. and *> DOTYPE(MAXTYP+1) is .TRUE. . *> Not modified. -*> \endverbatim -*> \verbatim +*> *> DOTYPE LOGICAL array, dimension (NTYPES) *> If DOTYPE(j) is .TRUE., then for each size in NN a *> matrix of that size and of type j will be generated. @@ -198,8 +195,7 @@ *> than MAXTYP, DOTYPE(MAXTYP+1) through DOTYPE(NTYPES) *> will be ignored. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> ISEED INTEGER array, dimension (4) *> On entry ISEED specifies the seed of the random number *> generator. The array elements should be between 0 and 4095; @@ -211,8 +207,7 @@ *> next call to DDRVSG to continue the same random number *> sequence. *> Modified. -*> \endverbatim -*> \verbatim +*> *> THRESH DOUBLE PRECISION *> A test will count as "failed" if the "error", computed as *> described above, exceeds THRESH. Note that the error @@ -221,105 +216,87 @@ *> it should not depend on the precision (single vs. double) *> or the size of the matrix. It must be at least zero. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> NOUNIT INTEGER *> The FORTRAN unit number for printing out error messages *> (e.g., if a routine returns IINFO not equal to 0.) *> Not modified. -*> \endverbatim -*> \verbatim +*> *> A DOUBLE PRECISION array, dimension (LDA , max(NN)) *> Used to hold the matrix whose eigenvalues are to be *> computed. On exit, A contains the last matrix actually *> used. *> Modified. -*> \endverbatim -*> \verbatim +*> *> LDA INTEGER *> The leading dimension of A and AB. It must be at *> least 1 and at least max( NN ). *> Not modified. -*> \endverbatim -*> \verbatim +*> *> B DOUBLE PRECISION array, dimension (LDB , max(NN)) *> Used to hold the symmetric positive definite matrix for *> the generailzed problem. *> On exit, B contains the last matrix actually *> used. *> Modified. -*> \endverbatim -*> \verbatim +*> *> LDB INTEGER *> The leading dimension of B and BB. It must be at *> least 1 and at least max( NN ). *> Not modified. -*> \endverbatim -*> \verbatim +*> *> D DOUBLE PRECISION array, dimension (max(NN)) *> The eigenvalues of A. On exit, the eigenvalues in D *> correspond with the matrix in A. *> Modified. -*> \endverbatim -*> \verbatim +*> *> Z DOUBLE PRECISION array, dimension (LDZ, max(NN)) *> The matrix of eigenvectors. *> Modified. -*> \endverbatim -*> \verbatim +*> *> LDZ INTEGER *> The leading dimension of Z. It must be at least 1 and *> at least max( NN ). *> Not modified. -*> \endverbatim -*> \verbatim +*> *> AB DOUBLE PRECISION array, dimension (LDA, max(NN)) *> Workspace. *> Modified. -*> \endverbatim -*> \verbatim +*> *> BB DOUBLE PRECISION array, dimension (LDB, max(NN)) *> Workspace. *> Modified. -*> \endverbatim -*> \verbatim +*> *> AP DOUBLE PRECISION array, dimension (max(NN)**2) *> Workspace. *> Modified. -*> \endverbatim -*> \verbatim +*> *> BP DOUBLE PRECISION array, dimension (max(NN)**2) *> Workspace. *> Modified. -*> \endverbatim -*> \verbatim +*> *> WORK DOUBLE PRECISION array, dimension (NWORK) *> Workspace. *> Modified. -*> \endverbatim -*> \verbatim +*> *> NWORK INTEGER *> The number of entries in WORK. This must be at least *> 1+5*N+2*N*lg(N)+3*N**2 where N = max( NN(j) ) and *> lg( N ) = smallest integer k such that 2**k >= N. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> IWORK INTEGER array, dimension (LIWORK) *> Workspace. *> Modified. -*> \endverbatim -*> \verbatim +*> *> LIWORK INTEGER *> The number of entries in WORK. This must be at least 6*N. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> RESULT DOUBLE PRECISION array, dimension (70) *> The values computed by the 70 tests described above. *> Modified. -*> \endverbatim -*> \verbatim +*> *> INFO INTEGER *> If 0, then everything ran OK. *> -1: NSIZES < 0 @@ -334,11 +311,9 @@ *> DSBGVD, DSYGVX, DSPGVX or SSBGVX returns an error code, *> the absolute value of it is returned. *> Modified. -*> \endverbatim -*> \verbatim +*> *> ---------------------------------------------------------------------- -*> \endverbatim -*> \verbatim +*> *> Some Local Variables and Parameters: *> ---- ----- --------- --- ---------- *> ZERO, ONE Real 0 and 1. @@ -352,8 +327,7 @@ *> so far (computed by DLAFTS). *> COND, IMODE Values to be passed to the matrix generators. *> ANORM Norm of A; passed to matrix generators. -*> \endverbatim -*> \verbatim +*> *> OVFL, UNFL Overflow and underflow thresholds. *> ULP, ULPINV Finest relative precision and its inverse. *> RTOVFL, RTUNFL Square roots of the previous 2 values. diff --git a/TESTING/EIG/ddrvst.f b/TESTING/EIG/ddrvst.f index 623ea3d8..e1319caf 100644 --- a/TESTING/EIG/ddrvst.f +++ b/TESTING/EIG/ddrvst.f @@ -151,15 +151,13 @@ *> The number of sizes of matrices to use. If it is zero, *> DDRVST does nothing. It must be at least zero. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> NN INTEGER array, dimension (NSIZES) *> An array containing the sizes to be used for the matrices. *> Zero values will be skipped. The values must be at least *> zero. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> NTYPES INTEGER *> The number of elements in DOTYPE. If it is zero, DDRVST *> does nothing. It must be at least zero. If it is MAXTYP+1 @@ -168,8 +166,7 @@ *> is only useful if DOTYPE(1:MAXTYP) is .FALSE. and *> DOTYPE(MAXTYP+1) is .TRUE. . *> Not modified. -*> \endverbatim -*> \verbatim +*> *> DOTYPE LOGICAL array, dimension (NTYPES) *> If DOTYPE(j) is .TRUE., then for each size in NN a *> matrix of that size and of type j will be generated. @@ -179,8 +176,7 @@ *> than MAXTYP, DOTYPE(MAXTYP+1) through DOTYPE(NTYPES) *> will be ignored. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> ISEED INTEGER array, dimension (4) *> On entry ISEED specifies the seed of the random number *> generator. The array elements should be between 0 and 4095; @@ -192,8 +188,7 @@ *> next call to DDRVST to continue the same random number *> sequence. *> Modified. -*> \endverbatim -*> \verbatim +*> *> THRESH DOUBLE PRECISION *> A test will count as "failed" if the "error", computed as *> described above, exceeds THRESH. Note that the error @@ -202,120 +197,99 @@ *> it should not depend on the precision (single vs. double) *> or the size of the matrix. It must be at least zero. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> NOUNIT INTEGER *> The FORTRAN unit number for printing out error messages *> (e.g., if a routine returns IINFO not equal to 0.) *> Not modified. -*> \endverbatim -*> \verbatim +*> *> A DOUBLE PRECISION array, dimension (LDA , max(NN)) *> Used to hold the matrix whose eigenvalues are to be *> computed. On exit, A contains the last matrix actually *> used. *> Modified. -*> \endverbatim -*> \verbatim +*> *> LDA INTEGER *> The leading dimension of A. It must be at *> least 1 and at least max( NN ). *> Not modified. -*> \endverbatim -*> \verbatim +*> *> D1 DOUBLE PRECISION array, dimension (max(NN)) *> The eigenvalues of A, as computed by DSTEQR simlutaneously *> with Z. On exit, the eigenvalues in D1 correspond with the *> matrix in A. *> Modified. -*> \endverbatim -*> \verbatim +*> *> D2 DOUBLE PRECISION array, dimension (max(NN)) *> The eigenvalues of A, as computed by DSTEQR if Z is not *> computed. On exit, the eigenvalues in D2 correspond with *> the matrix in A. *> Modified. -*> \endverbatim -*> \verbatim +*> *> D3 DOUBLE PRECISION array, dimension (max(NN)) *> The eigenvalues of A, as computed by DSTERF. On exit, the *> eigenvalues in D3 correspond with the matrix in A. *> Modified. -*> \endverbatim -*> \verbatim +*> *> D4 DOUBLE PRECISION array, dimension -*> \endverbatim -*> \verbatim +*> *> EVEIGS DOUBLE PRECISION array, dimension (max(NN)) *> The eigenvalues as computed by DSTEV('N', ... ) *> (I reserve the right to change this to the output of *> whichever algorithm computes the most accurate eigenvalues). -*> \endverbatim -*> \verbatim +*> *> WA1 DOUBLE PRECISION array, dimension -*> \endverbatim -*> \verbatim +*> *> WA2 DOUBLE PRECISION array, dimension -*> \endverbatim -*> \verbatim +*> *> WA3 DOUBLE PRECISION array, dimension -*> \endverbatim -*> \verbatim +*> *> U DOUBLE PRECISION array, dimension (LDU, max(NN)) *> The orthogonal matrix computed by DSYTRD + DORGTR. *> Modified. -*> \endverbatim -*> \verbatim +*> *> LDU INTEGER *> The leading dimension of U, Z, and V. It must be at *> least 1 and at least max( NN ). *> Not modified. -*> \endverbatim -*> \verbatim +*> *> V DOUBLE PRECISION array, dimension (LDU, max(NN)) *> The Housholder vectors computed by DSYTRD in reducing A to *> tridiagonal form. *> Modified. -*> \endverbatim -*> \verbatim +*> *> TAU DOUBLE PRECISION array, dimension (max(NN)) *> The Householder factors computed by DSYTRD in reducing A *> to tridiagonal form. *> Modified. -*> \endverbatim -*> \verbatim +*> *> Z DOUBLE PRECISION array, dimension (LDU, max(NN)) *> The orthogonal matrix of eigenvectors computed by DSTEQR, *> DPTEQR, and DSTEIN. *> Modified. -*> \endverbatim -*> \verbatim +*> *> WORK DOUBLE PRECISION array, dimension (LWORK) *> Workspace. *> Modified. -*> \endverbatim -*> \verbatim +*> *> LWORK INTEGER *> The number of entries in WORK. This must be at least *> 1 + 4 * Nmax + 2 * Nmax * lg Nmax + 4 * Nmax**2 *> where Nmax = max( NN(j), 2 ) and lg = log base 2. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> IWORK INTEGER array, *> dimension (6 + 6*Nmax + 5 * Nmax * lg Nmax ) *> where Nmax = max( NN(j), 2 ) and lg = log base 2. *> Workspace. *> Modified. -*> \endverbatim -*> \verbatim +*> *> RESULT DOUBLE PRECISION array, dimension (105) *> The values computed by the tests described above. *> The values are currently limited to 1/ulp, to avoid *> overflow. *> Modified. -*> \endverbatim -*> \verbatim +*> *> INFO INTEGER *> If 0, then everything ran OK. *> -1: NSIZES < 0 @@ -329,11 +303,9 @@ *> or DORMTR returns an error code, the *> absolute value of it is returned. *> Modified. -*> \endverbatim -*> \verbatim +*> *>----------------------------------------------------------------------- -*> \endverbatim -*> \verbatim +*> *> Some Local Variables and Parameters: *> ---- ----- --------- --- ---------- *> ZERO, ONE Real 0 and 1. @@ -347,8 +319,7 @@ *> so far (computed by DLAFTS). *> COND, IMODE Values to be passed to the matrix generators. *> ANORM Norm of A; passed to matrix generators. -*> \endverbatim -*> \verbatim +*> *> OVFL, UNFL Overflow and underflow thresholds. *> ULP, ULPINV Finest relative precision and its inverse. *> RTOVFL, RTUNFL Square roots of the previous 2 values. @@ -358,8 +329,7 @@ *> generator for type "j". *> KMAGN(j) The order of magnitude ( O(1), *> O(overflow^(1/2) ), O(underflow^(1/2) ) -*> \endverbatim -*> \verbatim +*> *> The tests performed are: Routine tested *> 1= | A - U S U' | / ( |A| n ulp ) DSTEV('V', ... ) *> 2= | I - U U' | / ( n ulp ) DSTEV('V', ... ) @@ -385,8 +355,7 @@ *> 22= | A - U S U' | / ( |A| n ulp ) DSTEVR('V','V', ... ) *> 23= | I - U U' | / ( n ulp ) DSTEVR('V','V', ... ) *> 24= |D(with Z) - D(w/o Z)| / (|D| ulp) DSTEVR('N','V', ... ) -*> \endverbatim -*> \verbatim +*> *> 25= | A - U S U' | / ( |A| n ulp ) DSYEV('L','V', ... ) *> 26= | I - U U' | / ( n ulp ) DSYEV('L','V', ... ) *> 27= |D(with Z) - D(w/o Z)| / (|D| ulp) DSYEV('L','N', ... ) @@ -441,15 +410,12 @@ *> 76= | A - U S U' | / ( |A| n ulp ) DSYEVR('L','V','V', ... ) *> 77= | I - U U' | / ( n ulp ) DSYEVR('L','V','V', ... ) *> 78= |D(with Z) - D(w/o Z)| / (|D| ulp) DSYEVR('L','N','V', ... ) -*> \endverbatim -*> \verbatim +*> *> Tests 25 through 78 are repeated (as tests 79 through 132) *> with UPLO='U' -*> \endverbatim -*> \verbatim +*> *> To be added in 1999 -*> \endverbatim -*> \verbatim +*> *> 79= | A - U S U' | / ( |A| n ulp ) DSPEVR('L','V','A', ... ) *> 80= | I - U U' | / ( n ulp ) DSPEVR('L','V','A', ... ) *> 81= |D(with Z) - D(w/o Z)| / (|D| ulp) DSPEVR('L','N','A', ... ) diff --git a/TESTING/EIG/ddrvsx.f b/TESTING/EIG/ddrvsx.f index 2644eba1..11688d6f 100644 --- a/TESTING/EIG/ddrvsx.f +++ b/TESTING/EIG/ddrvsx.f @@ -326,8 +326,7 @@ *> \param[out] WI *> \verbatim *> WI is DOUBLE PRECISION array, dimension (max(NN)) -*> \endverbatim -*> \verbatim +*> *> The real and imaginary parts of the eigenvalues of A. *> On exit, WR + WI*i are the eigenvalues of the matrix in A. *> \endverbatim @@ -340,8 +339,7 @@ *> \param[out] WIT *> \verbatim *> WIT is DOUBLE PRECISION array, dimension (max(NN)) -*> \endverbatim -*> \verbatim +*> *> Like WR, WI, these arrays contain the eigenvalues of A, *> but those computed when DGEESX only computes a partial *> eigendecomposition, i.e. not Schur vectors @@ -355,8 +353,7 @@ *> \param[out] WITMP *> \verbatim *> WITMP is DOUBLE PRECISION array, dimension (max(NN)) -*> \endverbatim -*> \verbatim +*> *> More temporary storage for eigenvalues. *> \endverbatim *> @@ -414,11 +411,9 @@ *> <0, input parameter -INFO is incorrect *> >0, DLATMR, SLATMS, SLATME or DGET24 returned an error *> code and INFO is its absolute value -*> \endverbatim -*> \verbatim +*> *>----------------------------------------------------------------------- -*> \endverbatim -*> \verbatim +*> *> Some Local Variables and Parameters: *> ---- ----- --------- --- ---------- *> ZERO, ONE Real 0 and 1. @@ -428,8 +423,7 @@ *> COND, CONDS, *> IMODE Values to be passed to the matrix generators. *> ANORM Norm of A; passed to matrix generators. -*> \endverbatim -*> \verbatim +*> *> OVFL, UNFL Overflow and underflow thresholds. *> ULP, ULPINV Finest relative precision and its inverse. *> RTULP, RTULPI Square roots of the previous 4 values. diff --git a/TESTING/EIG/ddrvvx.f b/TESTING/EIG/ddrvvx.f index 230912b3..4b754313 100644 --- a/TESTING/EIG/ddrvvx.f +++ b/TESTING/EIG/ddrvvx.f @@ -328,8 +328,7 @@ *> \param[out] WI *> \verbatim *> WI is DOUBLE PRECISION array, dimension (max(NN)) -*> \endverbatim -*> \verbatim +*> *> The real and imaginary parts of the eigenvalues of A. *> On exit, WR + WI*i are the eigenvalues of the matrix in A. *> \endverbatim @@ -342,8 +341,7 @@ *> \param[out] WI1 *> \verbatim *> WI1 is DOUBLE PRECISION array, dimension (max(NN,12)) -*> \endverbatim -*> \verbatim +*> *> Like WR, WI, these arrays contain the eigenvalues of A, *> but those computed when DGEEVX only computes a partial *> eigendecomposition, i.e. not the eigenvalues and left @@ -477,15 +475,12 @@ *> If <0, then input paramter -INFO is incorrect. *> If >0, DLATMR, SLATMS, SLATME or DGET23 returned an error *> code, and INFO is its absolute value. -*> \endverbatim -*> \verbatim +*> *>----------------------------------------------------------------------- -*> \endverbatim -*> \verbatim +*> *> Some Local Variables and Parameters: *> ---- ----- --------- --- ---------- -*> \endverbatim -*> \verbatim +*> *> ZERO, ONE Real 0 and 1. *> MAXTYP The number of types defined. *> NMAX Largest value in NN or 12. @@ -493,13 +488,11 @@ *> COND, CONDS, *> IMODE Values to be passed to the matrix generators. *> ANORM Norm of A; passed to matrix generators. -*> \endverbatim -*> \verbatim +*> *> OVFL, UNFL Overflow and underflow thresholds. *> ULP, ULPINV Finest relative precision and its inverse. *> RTULP, RTULPI Square roots of the previous 4 values. -*> \endverbatim -*> \verbatim +*> *> The following four arrays decode JTYPE: *> KTYPE(j) The general type (1-10) for type "j". *> KMODE(j) The MODE value to be passed to the matrix diff --git a/TESTING/EIG/dget22.f b/TESTING/EIG/dget22.f index 8749c6f5..39aa95d0 100644 --- a/TESTING/EIG/dget22.f +++ b/TESTING/EIG/dget22.f @@ -131,8 +131,7 @@ *> \param[in] WI *> \verbatim *> WI is DOUBLE PRECISION array, dimension (N) -*> \endverbatim -*> \verbatim +*> *> The real and imaginary parts of the eigenvalues of A. *> Purely real eigenvalues are indicated by WI(j) = 0. *> Complex conjugate pairs are indicated by WR(j)=WR(j+1) and diff --git a/TESTING/EIG/dget23.f b/TESTING/EIG/dget23.f index bd1471dc..c055d620 100644 --- a/TESTING/EIG/dget23.f +++ b/TESTING/EIG/dget23.f @@ -214,8 +214,7 @@ *> \param[out] WI *> \verbatim *> WI is DOUBLE PRECISION array, dimension (N) -*> \endverbatim -*> \verbatim +*> *> The real and imaginary parts of the eigenvalues of A. *> On exit, WR + WI*i are the eigenvalues of the matrix in A. *> \endverbatim @@ -228,8 +227,7 @@ *> \param[out] WI1 *> \verbatim *> WI1 is DOUBLE PRECISION array, dimension (N) -*> \endverbatim -*> \verbatim +*> *> Like WR, WI, these arrays contain the eigenvalues of A, *> but those computed when DGEEVX only computes a partial *> eigendecomposition, i.e. not the eigenvalues and left diff --git a/TESTING/EIG/dget24.f b/TESTING/EIG/dget24.f index d5901759..5020f56c 100644 --- a/TESTING/EIG/dget24.f +++ b/TESTING/EIG/dget24.f @@ -205,8 +205,7 @@ *> \param[out] WI *> \verbatim *> WI is DOUBLE PRECISION array, dimension (N) -*> \endverbatim -*> \verbatim +*> *> The real and imaginary parts of the eigenvalues of A. *> On exit, WR + WI*i are the eigenvalues of the matrix in A. *> \endverbatim @@ -219,8 +218,7 @@ *> \param[out] WIT *> \verbatim *> WIT is DOUBLE PRECISION array, dimension (N) -*> \endverbatim -*> \verbatim +*> *> Like WR, WI, these arrays contain the eigenvalues of A, *> but those computed when DGEESX only computes a partial *> eigendecomposition, i.e. not Schur vectors @@ -234,8 +232,7 @@ *> \param[out] WITMP *> \verbatim *> WITMP is DOUBLE PRECISION array, dimension (N) -*> \endverbatim -*> \verbatim +*> *> Like WR, WI, these arrays contain the eigenvalues of A, *> but sorted by increasing real part. *> \endverbatim diff --git a/TESTING/EIG/dgsvts.f b/TESTING/EIG/dgsvts.f index 8e410520..7f989f25 100644 --- a/TESTING/EIG/dgsvts.f +++ b/TESTING/EIG/dgsvts.f @@ -141,8 +141,7 @@ *> \param[out] BETA *> \verbatim *> BETA is DOUBLE PRECISION array, dimension (N) -*> \endverbatim -*> \verbatim +*> *> The generalized singular value pairs of A and B, the *> ``diagonal'' matrices D1 and D2 are constructed from *> ALPHA and BETA, see subroutine DGGSVD for details. diff --git a/TESTING/EIG/dhst01.f b/TESTING/EIG/dhst01.f index 5ddffc1a..9af8fd9b 100644 --- a/TESTING/EIG/dhst01.f +++ b/TESTING/EIG/dhst01.f @@ -56,8 +56,7 @@ *> \param[in] IHI *> \verbatim *> IHI is INTEGER -*> \endverbatim -*> \verbatim +*> *> A is assumed to be upper triangular in rows and columns *> 1:ILO-1 and IHI+1:N, so Q differs from the identity only in *> rows and columns ILO+1:IHI. diff --git a/TESTING/EIG/dlafts.f b/TESTING/EIG/dlafts.f index a71a832a..29cf46f5 100644 --- a/TESTING/EIG/dlafts.f +++ b/TESTING/EIG/dlafts.f @@ -41,51 +41,43 @@ *> On entry, TYPE specifies the matrix type to be used in the *> printed messages. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> N - INTEGER *> On entry, N specifies the order of the test matrix. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> IMAT - INTEGER *> On entry, IMAT specifies the type of the test matrix. *> A listing of the different types is printed by DLAHD2 *> to the output file if a test fails to pass the threshold. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> NTESTS - INTEGER *> On entry, NTESTS is the number of tests performed on the *> subroutines in the path given by TYPE. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> RESULT - DOUBLE PRECISION array of dimension( NTESTS ) *> On entry, RESULT contains the test ratios from the tests *> performed in the calling program. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> ISEED - INTEGER array of dimension( 4 ) *> Contains the random seed that generated the matrix used *> for the tests whose ratios are in RESULT. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> THRESH - DOUBLE PRECISION *> On entry, THRESH specifies the acceptable threshold of the *> test ratios. If RESULT( K ) > THRESH, then the K-th test *> did not pass the threshold and a message will be printed. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> IOUNIT - INTEGER *> On entry, IOUNIT specifies the unit number of the file *> to which the messages are printed. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> IE - INTEGER *> On entry, IE contains the number of tests which have *> failed to pass the threshold so far. diff --git a/TESTING/EIG/dlahd2.f b/TESTING/EIG/dlahd2.f index 33844fdf..9d2cac9a 100644 --- a/TESTING/EIG/dlahd2.f +++ b/TESTING/EIG/dlahd2.f @@ -40,15 +40,13 @@ *> PATH is CHARACTER*3. *> On entry, PATH contains the name of the path for which the *> header information is to be printed. Current paths are -*> \endverbatim -*> \verbatim +*> *> DHS, ZHS: Non-symmetric eigenproblem. *> DST, ZST: Symmetric eigenproblem. *> DSG, ZSG: Symmetric Generalized eigenproblem. *> DBD, ZBD: Singular Value Decomposition (SVD) *> DBB, ZBB: General Banded reduction to bidiagonal form -*> \endverbatim -*> \verbatim +*> *> These paths also are supplied in double precision (replace *> leading S by D and leading C by Z in path names). *> \endverbatim diff --git a/TESTING/EIG/dlarhs.f b/TESTING/EIG/dlarhs.f index c1331dad..b071ccd9 100644 --- a/TESTING/EIG/dlarhs.f +++ b/TESTING/EIG/dlarhs.f @@ -113,12 +113,10 @@ *> KU is INTEGER *> Used only if A is a general band matrix or if A is *> triangular. -*> \endverbatim -*> \verbatim +*> *> If PATH = xGB, specifies the number of superdiagonals of A, *> and 0 <= KU <= N-1. -*> \endverbatim -*> \verbatim +*> *> If PATH = xTR, xTP, or xTB, specifies whether or not the *> matrix has unit diagonal: *> = 1: matrix has non-unit diagonal (default) diff --git a/TESTING/EIG/dlatm4.f b/TESTING/EIG/dlatm4.f index d692f067..dcb3d9c1 100644 --- a/TESTING/EIG/dlatm4.f +++ b/TESTING/EIG/dlatm4.f @@ -48,8 +48,7 @@ *> If ITYPE < 0, then type abs(ITYPE) is generated and then *> swapped end for end (A(I,J) := A'(N-J,N-I).) See also *> the description of AMAGN and ISIGN. -*> \endverbatim -*> \verbatim +*> *> Special types: *> = 0: the zero matrix. *> = 1: the identity. @@ -58,8 +57,7 @@ *> followed by a k x k identity block, where k=(N-1)/2. *> If N is even, then k=(N-2)/2, and a zero diagonal entry *> is tacked onto the end. -*> \endverbatim -*> \verbatim +*> *> Diagonal types. The diagonal consists of NZ1 zeros, then *> k=N-NZ1-NZ2 nonzeros. The subdiagonal is zero. ITYPE *> specifies the nonzero diagonal entries as follows: diff --git a/TESTING/EIG/dspt21.f b/TESTING/EIG/dspt21.f index f310bb26..862b26e2 100644 --- a/TESTING/EIG/dspt21.f +++ b/TESTING/EIG/dspt21.f @@ -95,12 +95,10 @@ *> Specifies the type of tests to be performed. *> 1: U expressed as a dense orthogonal matrix: *> RESULT(1) = | A - U S U' | / ( |A| n ulp ) *andC> RESULT(2) = | I - UU' | / ( n ulp ) -*> \endverbatim -*> \verbatim +*> *> 2: U expressed as a product V of Housholder transformations: *> RESULT(1) = | A - V S V' | / ( |A| n ulp ) -*> \endverbatim -*> \verbatim +*> *> 3: U expressed both as a dense orthogonal matrix and *> as a product of Housholder transformations: *> RESULT(1) = | I - VU' | / ( n ulp ) diff --git a/TESTING/EIG/dsyt21.f b/TESTING/EIG/dsyt21.f index c908b3cd..c6ee194a 100644 --- a/TESTING/EIG/dsyt21.f +++ b/TESTING/EIG/dsyt21.f @@ -67,12 +67,10 @@ *> Specifies the type of tests to be performed. *> 1: U expressed as a dense orthogonal matrix: *> RESULT(1) = | A - U S U' | / ( |A| n ulp ) *andC> RESULT(2) = | I - UU' | / ( n ulp ) -*> \endverbatim -*> \verbatim +*> *> 2: U expressed as a product V of Housholder transformations: *> RESULT(1) = | A - V S V' | / ( |A| n ulp ) -*> \endverbatim -*> \verbatim +*> *> 3: U expressed both as a dense orthogonal matrix and *> as a product of Housholder transformations: *> RESULT(1) = | I - VU' | / ( n ulp ) diff --git a/TESTING/EIG/dsyt22.f b/TESTING/EIG/dsyt22.f index 23ecdbf5..dc0e5ba4 100644 --- a/TESTING/EIG/dsyt22.f +++ b/TESTING/EIG/dsyt22.f @@ -53,100 +53,85 @@ *> Specifies the type of tests to be performed. *> 1: U expressed as a dense orthogonal matrix: *> RESULT(1) = | A - U S U' | / ( |A| n ulp ) *andC> RESULT(2) = | I - UU' | / ( n ulp ) -*> \endverbatim -*> \verbatim +*> *> UPLO CHARACTER *> If UPLO='U', the upper triangle of A will be used and the *> (strictly) lower triangle will not be referenced. If *> UPLO='L', the lower triangle of A will be used and the *> (strictly) upper triangle will not be referenced. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> N INTEGER *> The size of the matrix. If it is zero, DSYT22 does nothing. *> It must be at least zero. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> M INTEGER *> The number of columns of U. If it is zero, DSYT22 does *> nothing. It must be at least zero. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> KBAND INTEGER *> The bandwidth of the matrix. It may only be zero or one. *> If zero, then S is diagonal, and E is not referenced. If *> one, then S is symmetric tri-diagonal. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> A DOUBLE PRECISION array, dimension (LDA , N) *> The original (unfactored) matrix. It is assumed to be *> symmetric, and only the upper (UPLO='U') or only the lower *> (UPLO='L') will be referenced. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> LDA INTEGER *> The leading dimension of A. It must be at least 1 *> and at least N. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> D DOUBLE PRECISION array, dimension (N) *> The diagonal of the (symmetric tri-) diagonal matrix. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> E DOUBLE PRECISION array, dimension (N) *> The off-diagonal of the (symmetric tri-) diagonal matrix. *> E(1) is ignored, E(2) is the (1,2) and (2,1) element, etc. *> Not referenced if KBAND=0. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> U DOUBLE PRECISION array, dimension (LDU, N) *> If ITYPE=1 or 3, this contains the orthogonal matrix in *> the decomposition, expressed as a dense matrix. If ITYPE=2, *> then it is not referenced. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> LDU INTEGER *> The leading dimension of U. LDU must be at least N and *> at least 1. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> V DOUBLE PRECISION array, dimension (LDV, N) *> If ITYPE=2 or 3, the lower triangle of this array contains *> the Householder vectors used to describe the orthogonal *> matrix in the decomposition. If ITYPE=1, then it is not *> referenced. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> LDV INTEGER *> The leading dimension of V. LDV must be at least N and *> at least 1. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> TAU DOUBLE PRECISION array, dimension (N) *> If ITYPE >= 2, then TAU(j) is the scalar factor of *> v(j) v(j)' in the Householder transformation H(j) of *> the product U = H(1)...H(n-2) *> If ITYPE < 2, then TAU is not referenced. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> WORK DOUBLE PRECISION array, dimension (2*N**2) *> Workspace. *> Modified. -*> \endverbatim -*> \verbatim +*> *> RESULT DOUBLE PRECISION array, dimension (2) *> The values computed by the two tests described above. The *> values are currently limited to 1/ulp, to avoid overflow. diff --git a/TESTING/EIG/ilaenv.f b/TESTING/EIG/ilaenv.f index 442a0e36..09b45a99 100644 --- a/TESTING/EIG/ilaenv.f +++ b/TESTING/EIG/ilaenv.f @@ -66,8 +66,7 @@ *> 12 <= ISPEC <= 16: *> xHSEQR or one of its subroutines, *> see IPARMQ for detailed explanation -*> \endverbatim -*> \verbatim +*> *> Other specifications (up to 100) can be added later. *> \endverbatim *> @@ -104,8 +103,7 @@ *> \param[in] N4 *> \verbatim *> N4 is INTEGER -*> \endverbatim -*> \verbatim +*> *> Problem dimensions for the subroutine NAME; these may not all *> be required. *> \endverbatim diff --git a/TESTING/EIG/schkbb.f b/TESTING/EIG/schkbb.f index a13ff601..59292703 100644 --- a/TESTING/EIG/schkbb.f +++ b/TESTING/EIG/schkbb.f @@ -310,11 +310,9 @@ *> \verbatim *> INFO is INTEGER *> If 0, then everything ran OK. -*> \endverbatim -*> \verbatim +*> *>----------------------------------------------------------------------- -*> \endverbatim -*> \verbatim +*> *> Some Local Variables and Parameters: *> ---- ----- --------- --- ---------- *> ZERO, ONE Real 0 and 1. @@ -328,8 +326,7 @@ *> so far. *> COND, IMODE Values to be passed to the matrix generators. *> ANORM Norm of A; passed to matrix generators. -*> \endverbatim -*> \verbatim +*> *> OVFL, UNFL Overflow and underflow thresholds. *> ULP, ULPINV Finest relative precision and its inverse. *> RTOVFL, RTUNFL Square roots of the previous 2 values. diff --git a/TESTING/EIG/schkbd.f b/TESTING/EIG/schkbd.f index 59d6edeb..831c6679 100644 --- a/TESTING/EIG/schkbd.f +++ b/TESTING/EIG/schkbd.f @@ -388,15 +388,12 @@ *> If SLATMR, SLATMS, SGEBRD, SORGBR, or SBDSQR, *> returns an error code, the *> absolute value of it is returned. -*> \endverbatim -*> \verbatim +*> *>----------------------------------------------------------------------- -*> \endverbatim -*> \verbatim +*> *> Some Local Variables and Parameters: *> ---- ----- --------- --- ---------- -*> \endverbatim -*> \verbatim +*> *> ZERO, ONE Real 0 and 1. *> MAXTYP The number of types defined. *> NTEST The number of tests performed, or which can @@ -409,13 +406,11 @@ *> NFAIL The number of tests which have exceeded THRESH *> COND, IMODE Values to be passed to the matrix generators. *> ANORM Norm of A; passed to matrix generators. -*> \endverbatim -*> \verbatim +*> *> OVFL, UNFL Overflow and underflow thresholds. *> RTOVFL, RTUNFL Square roots of the previous 2 values. *> ULP, ULPINV Finest relative precision and its inverse. -*> \endverbatim -*> \verbatim +*> *> The following four arrays decode JTYPE: *> KTYPE(j) The general type (1-10) for type "j". *> KMODE(j) The MODE value to be passed to the matrix diff --git a/TESTING/EIG/schkgg.f b/TESTING/EIG/schkgg.f index ef04fd8d..e458510f 100644 --- a/TESTING/EIG/schkgg.f +++ b/TESTING/EIG/schkgg.f @@ -423,8 +423,7 @@ *> \param[out] BETA1 *> \verbatim *> BETA1 is REAL array, dimension (max(NN)) -*> \endverbatim -*> \verbatim +*> *> The generalized eigenvalues of (A,B) computed by SHGEQZ *> when Q, Z, and the full Schur matrices are computed. *> On exit, ( ALPHR1(k)+ALPHI1(k)*i ) / BETA1(k) is the k-th diff --git a/TESTING/EIG/schkhs.f b/TESTING/EIG/schkhs.f index 60043c2d..47b6a079 100644 --- a/TESTING/EIG/schkhs.f +++ b/TESTING/EIG/schkhs.f @@ -166,15 +166,13 @@ *> The number of sizes of matrices to use. If it is zero, *> SCHKHS does nothing. It must be at least zero. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> NN - INTEGER array, dimension (NSIZES) *> An array containing the sizes to be used for the matrices. *> Zero values will be skipped. The values must be at least *> zero. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> NTYPES - INTEGER *> The number of elements in DOTYPE. If it is zero, SCHKHS *> does nothing. It must be at least zero. If it is MAXTYP+1 @@ -183,8 +181,7 @@ *> is only useful if DOTYPE(1:MAXTYP) is .FALSE. and *> DOTYPE(MAXTYP+1) is .TRUE. . *> Not modified. -*> \endverbatim -*> \verbatim +*> *> DOTYPE - LOGICAL array, dimension (NTYPES) *> If DOTYPE(j) is .TRUE., then for each size in NN a *> matrix of that size and of type j will be generated. @@ -194,8 +191,7 @@ *> than MAXTYP, DOTYPE(MAXTYP+1) through DOTYPE(NTYPES) *> will be ignored. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> ISEED - INTEGER array, dimension (4) *> On entry ISEED specifies the seed of the random number *> generator. The array elements should be between 0 and 4095; @@ -207,8 +203,7 @@ *> next call to SCHKHS to continue the same random number *> sequence. *> Modified. -*> \endverbatim -*> \verbatim +*> *> THRESH - REAL *> A test will count as "failed" if the "error", computed as *> described above, exceeds THRESH. Note that the error @@ -217,75 +212,63 @@ *> it should not depend on the precision (single vs. double) *> or the size of the matrix. It must be at least zero. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> NOUNIT - INTEGER *> The FORTRAN unit number for printing out error messages *> (e.g., if a routine returns IINFO not equal to 0.) *> Not modified. -*> \endverbatim -*> \verbatim +*> *> A - REAL array, dimension (LDA,max(NN)) *> Used to hold the matrix whose eigenvalues are to be *> computed. On exit, A contains the last matrix actually *> used. *> Modified. -*> \endverbatim -*> \verbatim +*> *> LDA - INTEGER *> The leading dimension of A, H, T1 and T2. It must be at *> least 1 and at least max( NN ). *> Not modified. -*> \endverbatim -*> \verbatim +*> *> H - REAL array, dimension (LDA,max(NN)) *> The upper hessenberg matrix computed by SGEHRD. On exit, *> H contains the Hessenberg form of the matrix in A. *> Modified. -*> \endverbatim -*> \verbatim +*> *> T1 - REAL array, dimension (LDA,max(NN)) *> The Schur (="quasi-triangular") matrix computed by SHSEQR *> if Z is computed. On exit, T1 contains the Schur form of *> the matrix in A. *> Modified. -*> \endverbatim -*> \verbatim +*> *> T2 - REAL array, dimension (LDA,max(NN)) *> The Schur matrix computed by SHSEQR when Z is not computed. *> This should be identical to T1. *> Modified. -*> \endverbatim -*> \verbatim +*> *> LDU - INTEGER *> The leading dimension of U, Z, UZ and UU. It must be at *> least 1 and at least max( NN ). *> Not modified. -*> \endverbatim -*> \verbatim +*> *> U - REAL array, dimension (LDU,max(NN)) *> The orthogonal matrix computed by SGEHRD. *> Modified. -*> \endverbatim -*> \verbatim +*> *> Z - REAL array, dimension (LDU,max(NN)) *> The orthogonal matrix computed by SHSEQR. *> Modified. -*> \endverbatim -*> \verbatim +*> *> UZ - REAL array, dimension (LDU,max(NN)) *> The product of U times Z. *> Modified. -*> \endverbatim -*> \verbatim +*> *> WR1 - REAL array, dimension (max(NN)) *> WI1 - REAL array, dimension (max(NN)) *> The real and imaginary parts of the eigenvalues of A, *> as computed when Z is computed. *> On exit, WR1 + WI1*i are the eigenvalues of the matrix in A. *> Modified. -*> \endverbatim -*> \verbatim +*> *> WR3 - REAL array, dimension (max(NN)) *> WI3 - REAL array, dimension (max(NN)) *> Like WR1, WI1, these arrays contain the eigenvalues of A, @@ -294,72 +277,60 @@ *> Schur form than is necessary for computing the *> eigenvalues. *> Modified. -*> \endverbatim -*> \verbatim +*> *> EVECTL - REAL array, dimension (LDU,max(NN)) *> The (upper triangular) left eigenvector matrix for the *> matrix in T1. For complex conjugate pairs, the real part *> is stored in one row and the imaginary part in the next. *> Modified. -*> \endverbatim -*> \verbatim +*> *> EVECTR - REAL array, dimension (LDU,max(NN)) *> The (upper triangular) right eigenvector matrix for the *> matrix in T1. For complex conjugate pairs, the real part *> is stored in one column and the imaginary part in the next. *> Modified. -*> \endverbatim -*> \verbatim +*> *> EVECTY - REAL array, dimension (LDU,max(NN)) *> The left eigenvector matrix for the *> matrix in H. For complex conjugate pairs, the real part *> is stored in one row and the imaginary part in the next. *> Modified. -*> \endverbatim -*> \verbatim +*> *> EVECTX - REAL array, dimension (LDU,max(NN)) *> The right eigenvector matrix for the *> matrix in H. For complex conjugate pairs, the real part *> is stored in one column and the imaginary part in the next. *> Modified. -*> \endverbatim -*> \verbatim +*> *> UU - REAL array, dimension (LDU,max(NN)) *> Details of the orthogonal matrix computed by SGEHRD. *> Modified. -*> \endverbatim -*> \verbatim +*> *> TAU - REAL array, dimension(max(NN)) *> Further details of the orthogonal matrix computed by SGEHRD. *> Modified. -*> \endverbatim -*> \verbatim +*> *> WORK - REAL array, dimension (NWORK) *> Workspace. *> Modified. -*> \endverbatim -*> \verbatim +*> *> NWORK - INTEGER *> The number of entries in WORK. NWORK >= 4*NN(j)*NN(j) + 2. -*> \endverbatim -*> \verbatim +*> *> IWORK - INTEGER array, dimension (max(NN)) *> Workspace. *> Modified. -*> \endverbatim -*> \verbatim +*> *> SELECT - LOGICAL array, dimension (max(NN)) *> Workspace. *> Modified. -*> \endverbatim -*> \verbatim +*> *> RESULT - REAL array, dimension (14) *> The values computed by the fourteen tests described above. *> The values are currently limited to 1/ulp, to avoid *> overflow. *> Modified. -*> \endverbatim -*> \verbatim +*> *> INFO - INTEGER *> If 0, then everything ran OK. *> -1: NSIZES < 0 @@ -376,15 +347,12 @@ *> If >2, then 30*N iterations were not enough to find an *> eigenvalue or to decompose the problem. *> Modified. -*> \endverbatim -*> \verbatim +*> *>----------------------------------------------------------------------- -*> \endverbatim -*> \verbatim +*> *> Some Local Variables and Parameters: *> ---- ----- --------- --- ---------- -*> \endverbatim -*> \verbatim +*> *> ZERO, ONE Real 0 and 1. *> MAXTYP The number of types defined. *> MTEST The number of tests defined: care must be taken @@ -403,14 +371,12 @@ *> COND, CONDS, *> IMODE Values to be passed to the matrix generators. *> ANORM Norm of A; passed to matrix generators. -*> \endverbatim -*> \verbatim +*> *> OVFL, UNFL Overflow and underflow thresholds. *> ULP, ULPINV Finest relative precision and its inverse. *> RTOVFL, RTUNFL, *> RTULP, RTULPI Square roots of the previous 4 values. -*> \endverbatim -*> \verbatim +*> *> The following four arrays decode JTYPE: *> KTYPE(j) The general type (1-10) for type "j". *> KMODE(j) The MODE value to be passed to the matrix diff --git a/TESTING/EIG/schksb.f b/TESTING/EIG/schksb.f index 128428fd..2dafac47 100644 --- a/TESTING/EIG/schksb.f +++ b/TESTING/EIG/schksb.f @@ -249,11 +249,9 @@ *> \verbatim *> INFO is INTEGER *> If 0, then everything ran OK. -*> \endverbatim -*> \verbatim +*> *>----------------------------------------------------------------------- -*> \endverbatim -*> \verbatim +*> *> Some Local Variables and Parameters: *> ---- ----- --------- --- ---------- *> ZERO, ONE Real 0 and 1. @@ -267,8 +265,7 @@ *> so far. *> COND, IMODE Values to be passed to the matrix generators. *> ANORM Norm of A; passed to matrix generators. -*> \endverbatim -*> \verbatim +*> *> OVFL, UNFL Overflow and underflow thresholds. *> ULP, ULPINV Finest relative precision and its inverse. *> RTOVFL, RTUNFL Square roots of the previous 2 values. diff --git a/TESTING/EIG/schkst.f b/TESTING/EIG/schkst.f index c937cbae..1f516284 100644 --- a/TESTING/EIG/schkst.f +++ b/TESTING/EIG/schkst.f @@ -545,11 +545,9 @@ *> If SLATMR, SLATMS, SSYTRD, SORGTR, SSTEQR, SSTERF, *> or SORMC2 returns an error code, the *> absolute value of it is returned. -*> \endverbatim -*> \verbatim +*> *>----------------------------------------------------------------------- -*> \endverbatim -*> \verbatim +*> *> Some Local Variables and Parameters: *> ---- ----- --------- --- ---------- *> ZERO, ONE Real 0 and 1. @@ -564,8 +562,7 @@ *> so far. *> COND, IMODE Values to be passed to the matrix generators. *> ANORM Norm of A; passed to matrix generators. -*> \endverbatim -*> \verbatim +*> *> OVFL, UNFL Overflow and underflow thresholds. *> ULP, ULPINV Finest relative precision and its inverse. *> RTOVFL, RTUNFL Square roots of the previous 2 values. diff --git a/TESTING/EIG/sdrges.f b/TESTING/EIG/sdrges.f index f4823667..7fd779e1 100644 --- a/TESTING/EIG/sdrges.f +++ b/TESTING/EIG/sdrges.f @@ -346,8 +346,7 @@ *> \param[out] BETA *> \verbatim *> BETA is REAL array, dimension (max(NN)) -*> \endverbatim -*> \verbatim +*> *> The generalized eigenvalues of (A,B) computed by SGGES. *> ( ALPHAR(k)+ALPHAI(k)*i ) / BETA(k) is the k-th *> generalized eigenvalue of A and B. diff --git a/TESTING/EIG/sdrgev.f b/TESTING/EIG/sdrgev.f index c926dcc3..2adcd115 100644 --- a/TESTING/EIG/sdrgev.f +++ b/TESTING/EIG/sdrgev.f @@ -357,8 +357,7 @@ *> \param[out] BETA1 *> \verbatim *> BETA1 is REAL array, dimension (max(NN)) -*> \endverbatim -*> \verbatim +*> *> Like ALPHAR, ALPHAI, BETA, these arrays contain the *> eigenvalues of A and B, but those computed when SGGEV only *> computes a partial eigendecomposition, i.e. not the diff --git a/TESTING/EIG/sdrgvx.f b/TESTING/EIG/sdrgvx.f index d214e4ce..f1aad543 100644 --- a/TESTING/EIG/sdrgvx.f +++ b/TESTING/EIG/sdrgvx.f @@ -188,8 +188,7 @@ *> \param[out] BETA *> \verbatim *> BETA is REAL array, dimension (NSIZE) -*> \endverbatim -*> \verbatim +*> *> On exit, (ALPHAR + ALPHAI*i)/BETA are the eigenvalues. *> \endverbatim *> diff --git a/TESTING/EIG/sdrves.f b/TESTING/EIG/sdrves.f index ad3d1784..f4cdd51c 100644 --- a/TESTING/EIG/sdrves.f +++ b/TESTING/EIG/sdrves.f @@ -271,8 +271,7 @@ *> \param[out] WI *> \verbatim *> WI is REAL array, dimension (max(NN)) -*> \endverbatim -*> \verbatim +*> *> The real and imaginary parts of the eigenvalues of A. *> On exit, WR + WI*i are the eigenvalues of the matrix in A. *> \endverbatim @@ -285,8 +284,7 @@ *> \param[out] WIT *> \verbatim *> WIT is REAL array, dimension (max(NN)) -*> \endverbatim -*> \verbatim +*> *> Like WR, WI, these arrays contain the eigenvalues of A, *> but those computed when SGEES only computes a partial *> eigendecomposition, i.e. not Schur vectors @@ -346,15 +344,12 @@ *> -20: NWORK too small. *> If SLATMR, SLATMS, SLATME or SGEES returns an error code, *> the absolute value of it is returned. -*> \endverbatim -*> \verbatim +*> *>----------------------------------------------------------------------- -*> \endverbatim -*> \verbatim +*> *> Some Local Variables and Parameters: *> ---- ----- --------- --- ---------- -*> \endverbatim -*> \verbatim +*> *> ZERO, ONE Real 0 and 1. *> MAXTYP The number of types defined. *> NMAX Largest value in NN. @@ -362,13 +357,11 @@ *> COND, CONDS, *> IMODE Values to be passed to the matrix generators. *> ANORM Norm of A; passed to matrix generators. -*> \endverbatim -*> \verbatim +*> *> OVFL, UNFL Overflow and underflow thresholds. *> ULP, ULPINV Finest relative precision and its inverse. *> RTULP, RTULPI Square roots of the previous 4 values. -*> \endverbatim -*> \verbatim +*> *> The following four arrays decode JTYPE: *> KTYPE(j) The general type (1-10) for type "j". *> KMODE(j) The MODE value to be passed to the matrix diff --git a/TESTING/EIG/sdrvev.f b/TESTING/EIG/sdrvev.f index aa97cc53..f75d43bc 100644 --- a/TESTING/EIG/sdrvev.f +++ b/TESTING/EIG/sdrvev.f @@ -266,8 +266,7 @@ *> \param[out] WI *> \verbatim *> WI is REAL array, dimension (max(NN)) -*> \endverbatim -*> \verbatim +*> *> The real and imaginary parts of the eigenvalues of A. *> On exit, WR + WI*i are the eigenvalues of the matrix in A. *> \endverbatim @@ -280,8 +279,7 @@ *> \param[out] WI1 *> \verbatim *> WI1 is REAL array, dimension (max(NN)) -*> \endverbatim -*> \verbatim +*> *> Like WR, WI, these arrays contain the eigenvalues of A, *> but those computed when SGEEV only computes a partial *> eigendecomposition, i.e. not the eigenvalues and left @@ -363,15 +361,12 @@ *> -23: NWORK too small. *> If SLATMR, SLATMS, SLATME or SGEEV returns an error code, *> the absolute value of it is returned. -*> \endverbatim -*> \verbatim +*> *>----------------------------------------------------------------------- -*> \endverbatim -*> \verbatim +*> *> Some Local Variables and Parameters: *> ---- ----- --------- --- ---------- -*> \endverbatim -*> \verbatim +*> *> ZERO, ONE Real 0 and 1. *> MAXTYP The number of types defined. *> NMAX Largest value in NN. @@ -379,13 +374,11 @@ *> COND, CONDS, *> IMODE Values to be passed to the matrix generators. *> ANORM Norm of A; passed to matrix generators. -*> \endverbatim -*> \verbatim +*> *> OVFL, UNFL Overflow and underflow thresholds. *> ULP, ULPINV Finest relative precision and its inverse. *> RTULP, RTULPI Square roots of the previous 4 values. -*> \endverbatim -*> \verbatim +*> *> The following four arrays decode JTYPE: *> KTYPE(j) The general type (1-10) for type "j". *> KMODE(j) The MODE value to be passed to the matrix diff --git a/TESTING/EIG/sdrvgg.f b/TESTING/EIG/sdrvgg.f index 88d33d8e..87cfaca6 100644 --- a/TESTING/EIG/sdrvgg.f +++ b/TESTING/EIG/sdrvgg.f @@ -363,8 +363,7 @@ *> \param[out] BETA1 *> \verbatim *> BETA1 is REAL array, dimension (max(NN)) -*> \endverbatim -*> \verbatim +*> *> The generalized eigenvalues of (A,B) computed by SGEGS. *> ( ALPHR1(k)+ALPHI1(k)*i ) / BETA1(k) is the k-th *> generalized eigenvalue of the matrices in A and B. @@ -383,8 +382,7 @@ *> \param[out] BETA2 *> \verbatim *> BETA2 is REAL array, dimension (max(NN)) -*> \endverbatim -*> \verbatim +*> *> The generalized eigenvalues of (A,B) computed by SGEGV. *> ( ALPHR2(k)+ALPHI2(k)*i ) / BETA2(k) is the k-th *> generalized eigenvalue of the matrices in A and B. diff --git a/TESTING/EIG/sdrvsg.f b/TESTING/EIG/sdrvsg.f index 5b25720b..57703cdc 100644 --- a/TESTING/EIG/sdrvsg.f +++ b/TESTING/EIG/sdrvsg.f @@ -170,15 +170,13 @@ *> The number of sizes of matrices to use. If it is zero, *> SDRVSG does nothing. It must be at least zero. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> NN INTEGER array, dimension (NSIZES) *> An array containing the sizes to be used for the matrices. *> Zero values will be skipped. The values must be at least *> zero. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> NTYPES INTEGER *> The number of elements in DOTYPE. If it is zero, SDRVSG *> does nothing. It must be at least zero. If it is MAXTYP+1 @@ -187,8 +185,7 @@ *> is only useful if DOTYPE(1:MAXTYP) is .FALSE. and *> DOTYPE(MAXTYP+1) is .TRUE. . *> Not modified. -*> \endverbatim -*> \verbatim +*> *> DOTYPE LOGICAL array, dimension (NTYPES) *> If DOTYPE(j) is .TRUE., then for each size in NN a *> matrix of that size and of type j will be generated. @@ -198,8 +195,7 @@ *> than MAXTYP, DOTYPE(MAXTYP+1) through DOTYPE(NTYPES) *> will be ignored. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> ISEED INTEGER array, dimension (4) *> On entry ISEED specifies the seed of the random number *> generator. The array elements should be between 0 and 4095; @@ -211,8 +207,7 @@ *> next call to SDRVSG to continue the same random number *> sequence. *> Modified. -*> \endverbatim -*> \verbatim +*> *> THRESH REAL *> A test will count as "failed" if the "error", computed as *> described above, exceeds THRESH. Note that the error @@ -221,105 +216,87 @@ *> it should not depend on the precision (single vs. double) *> or the size of the matrix. It must be at least zero. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> NOUNIT INTEGER *> The FORTRAN unit number for printing out error messages *> (e.g., if a routine returns IINFO not equal to 0.) *> Not modified. -*> \endverbatim -*> \verbatim +*> *> A REAL array, dimension (LDA , max(NN)) *> Used to hold the matrix whose eigenvalues are to be *> computed. On exit, A contains the last matrix actually *> used. *> Modified. -*> \endverbatim -*> \verbatim +*> *> LDA INTEGER *> The leading dimension of A and AB. It must be at *> least 1 and at least max( NN ). *> Not modified. -*> \endverbatim -*> \verbatim +*> *> B REAL array, dimension (LDB , max(NN)) *> Used to hold the symmetric positive definite matrix for *> the generailzed problem. *> On exit, B contains the last matrix actually *> used. *> Modified. -*> \endverbatim -*> \verbatim +*> *> LDB INTEGER *> The leading dimension of B and BB. It must be at *> least 1 and at least max( NN ). *> Not modified. -*> \endverbatim -*> \verbatim +*> *> D REAL array, dimension (max(NN)) *> The eigenvalues of A. On exit, the eigenvalues in D *> correspond with the matrix in A. *> Modified. -*> \endverbatim -*> \verbatim +*> *> Z REAL array, dimension (LDZ, max(NN)) *> The matrix of eigenvectors. *> Modified. -*> \endverbatim -*> \verbatim +*> *> LDZ INTEGER *> The leading dimension of Z. It must be at least 1 and *> at least max( NN ). *> Not modified. -*> \endverbatim -*> \verbatim +*> *> AB REAL array, dimension (LDA, max(NN)) *> Workspace. *> Modified. -*> \endverbatim -*> \verbatim +*> *> BB REAL array, dimension (LDB, max(NN)) *> Workspace. *> Modified. -*> \endverbatim -*> \verbatim +*> *> AP REAL array, dimension (max(NN)**2) *> Workspace. *> Modified. -*> \endverbatim -*> \verbatim +*> *> BP REAL array, dimension (max(NN)**2) *> Workspace. *> Modified. -*> \endverbatim -*> \verbatim +*> *> WORK REAL array, dimension (NWORK) *> Workspace. *> Modified. -*> \endverbatim -*> \verbatim +*> *> NWORK INTEGER *> The number of entries in WORK. This must be at least *> 1+5*N+2*N*lg(N)+3*N**2 where N = max( NN(j) ) and *> lg( N ) = smallest integer k such that 2**k >= N. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> IWORK INTEGER array, dimension (LIWORK) *> Workspace. *> Modified. -*> \endverbatim -*> \verbatim +*> *> LIWORK INTEGER *> The number of entries in WORK. This must be at least 6*N. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> RESULT REAL array, dimension (70) *> The values computed by the 70 tests described above. *> Modified. -*> \endverbatim -*> \verbatim +*> *> INFO INTEGER *> If 0, then everything ran OK. *> -1: NSIZES < 0 @@ -334,11 +311,9 @@ *> SSBGVD, SSYGVX, SSPGVX or SSBGVX returns an error code, *> the absolute value of it is returned. *> Modified. -*> \endverbatim -*> \verbatim +*> *> ---------------------------------------------------------------------- -*> \endverbatim -*> \verbatim +*> *> Some Local Variables and Parameters: *> ---- ----- --------- --- ---------- *> ZERO, ONE Real 0 and 1. @@ -352,8 +327,7 @@ *> so far (computed by SLAFTS). *> COND, IMODE Values to be passed to the matrix generators. *> ANORM Norm of A; passed to matrix generators. -*> \endverbatim -*> \verbatim +*> *> OVFL, UNFL Overflow and underflow thresholds. *> ULP, ULPINV Finest relative precision and its inverse. *> RTOVFL, RTUNFL Square roots of the previous 2 values. diff --git a/TESTING/EIG/sdrvst.f b/TESTING/EIG/sdrvst.f index 7f603c21..cd554aff 100644 --- a/TESTING/EIG/sdrvst.f +++ b/TESTING/EIG/sdrvst.f @@ -151,15 +151,13 @@ *> The number of sizes of matrices to use. If it is zero, *> SDRVST does nothing. It must be at least zero. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> NN INTEGER array, dimension (NSIZES) *> An array containing the sizes to be used for the matrices. *> Zero values will be skipped. The values must be at least *> zero. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> NTYPES INTEGER *> The number of elements in DOTYPE. If it is zero, SDRVST *> does nothing. It must be at least zero. If it is MAXTYP+1 @@ -168,8 +166,7 @@ *> is only useful if DOTYPE(1:MAXTYP) is .FALSE. and *> DOTYPE(MAXTYP+1) is .TRUE. . *> Not modified. -*> \endverbatim -*> \verbatim +*> *> DOTYPE LOGICAL array, dimension (NTYPES) *> If DOTYPE(j) is .TRUE., then for each size in NN a *> matrix of that size and of type j will be generated. @@ -179,8 +176,7 @@ *> than MAXTYP, DOTYPE(MAXTYP+1) through DOTYPE(NTYPES) *> will be ignored. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> ISEED INTEGER array, dimension (4) *> On entry ISEED specifies the seed of the random number *> generator. The array elements should be between 0 and 4095; @@ -192,8 +188,7 @@ *> next call to SDRVST to continue the same random number *> sequence. *> Modified. -*> \endverbatim -*> \verbatim +*> *> THRESH REAL *> A test will count as "failed" if the "error", computed as *> described above, exceeds THRESH. Note that the error @@ -202,120 +197,99 @@ *> it should not depend on the precision (single vs. double) *> or the size of the matrix. It must be at least zero. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> NOUNIT INTEGER *> The FORTRAN unit number for printing out error messages *> (e.g., if a routine returns IINFO not equal to 0.) *> Not modified. -*> \endverbatim -*> \verbatim +*> *> A REAL array, dimension (LDA , max(NN)) *> Used to hold the matrix whose eigenvalues are to be *> computed. On exit, A contains the last matrix actually *> used. *> Modified. -*> \endverbatim -*> \verbatim +*> *> LDA INTEGER *> The leading dimension of A. It must be at *> least 1 and at least max( NN ). *> Not modified. -*> \endverbatim -*> \verbatim +*> *> D1 REAL array, dimension (max(NN)) *> The eigenvalues of A, as computed by SSTEQR simlutaneously *> with Z. On exit, the eigenvalues in D1 correspond with the *> matrix in A. *> Modified. -*> \endverbatim -*> \verbatim +*> *> D2 REAL array, dimension (max(NN)) *> The eigenvalues of A, as computed by SSTEQR if Z is not *> computed. On exit, the eigenvalues in D2 correspond with *> the matrix in A. *> Modified. -*> \endverbatim -*> \verbatim +*> *> D3 REAL array, dimension (max(NN)) *> The eigenvalues of A, as computed by SSTERF. On exit, the *> eigenvalues in D3 correspond with the matrix in A. *> Modified. -*> \endverbatim -*> \verbatim +*> *> D4 REAL array, dimension -*> \endverbatim -*> \verbatim +*> *> EVEIGS REAL array, dimension (max(NN)) *> The eigenvalues as computed by SSTEV('N', ... ) *> (I reserve the right to change this to the output of *> whichever algorithm computes the most accurate eigenvalues). -*> \endverbatim -*> \verbatim +*> *> WA1 REAL array, dimension -*> \endverbatim -*> \verbatim +*> *> WA2 REAL array, dimension -*> \endverbatim -*> \verbatim +*> *> WA3 REAL array, dimension -*> \endverbatim -*> \verbatim +*> *> U REAL array, dimension (LDU, max(NN)) *> The orthogonal matrix computed by SSYTRD + SORGTR. *> Modified. -*> \endverbatim -*> \verbatim +*> *> LDU INTEGER *> The leading dimension of U, Z, and V. It must be at *> least 1 and at least max( NN ). *> Not modified. -*> \endverbatim -*> \verbatim +*> *> V REAL array, dimension (LDU, max(NN)) *> The Housholder vectors computed by SSYTRD in reducing A to *> tridiagonal form. *> Modified. -*> \endverbatim -*> \verbatim +*> *> TAU REAL array, dimension (max(NN)) *> The Householder factors computed by SSYTRD in reducing A *> to tridiagonal form. *> Modified. -*> \endverbatim -*> \verbatim +*> *> Z REAL array, dimension (LDU, max(NN)) *> The orthogonal matrix of eigenvectors computed by SSTEQR, *> SPTEQR, and SSTEIN. *> Modified. -*> \endverbatim -*> \verbatim +*> *> WORK REAL array, dimension (LWORK) *> Workspace. *> Modified. -*> \endverbatim -*> \verbatim +*> *> LWORK INTEGER *> The number of entries in WORK. This must be at least *> 1 + 4 * Nmax + 2 * Nmax * lg Nmax + 4 * Nmax**2 *> where Nmax = max( NN(j), 2 ) and lg = log base 2. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> IWORK INTEGER array, *> dimension (6 + 6*Nmax + 5 * Nmax * lg Nmax ) *> where Nmax = max( NN(j), 2 ) and lg = log base 2. *> Workspace. *> Modified. -*> \endverbatim -*> \verbatim +*> *> RESULT REAL array, dimension (105) *> The values computed by the tests described above. *> The values are currently limited to 1/ulp, to avoid *> overflow. *> Modified. -*> \endverbatim -*> \verbatim +*> *> INFO INTEGER *> If 0, then everything ran OK. *> -1: NSIZES < 0 @@ -329,11 +303,9 @@ *> or SORMTR returns an error code, the *> absolute value of it is returned. *> Modified. -*> \endverbatim -*> \verbatim +*> *>----------------------------------------------------------------------- -*> \endverbatim -*> \verbatim +*> *> Some Local Variables and Parameters: *> ---- ----- --------- --- ---------- *> ZERO, ONE Real 0 and 1. @@ -347,8 +319,7 @@ *> so far (computed by SLAFTS). *> COND, IMODE Values to be passed to the matrix generators. *> ANORM Norm of A; passed to matrix generators. -*> \endverbatim -*> \verbatim +*> *> OVFL, UNFL Overflow and underflow thresholds. *> ULP, ULPINV Finest relative precision and its inverse. *> RTOVFL, RTUNFL Square roots of the previous 2 values. @@ -358,8 +329,7 @@ *> generator for type "j". *> KMAGN(j) The order of magnitude ( O(1), *> O(overflow^(1/2) ), O(underflow^(1/2) ) -*> \endverbatim -*> \verbatim +*> *> The tests performed are: Routine tested *> 1= | A - U S U' | / ( |A| n ulp ) SSTEV('V', ... ) *> 2= | I - U U' | / ( n ulp ) SSTEV('V', ... ) @@ -385,8 +355,7 @@ *> 22= | A - U S U' | / ( |A| n ulp ) SSTEVR('V','V', ... ) *> 23= | I - U U' | / ( n ulp ) SSTEVR('V','V', ... ) *> 24= |D(with Z) - D(w/o Z)| / (|D| ulp) SSTEVR('N','V', ... ) -*> \endverbatim -*> \verbatim +*> *> 25= | A - U S U' | / ( |A| n ulp ) SSYEV('L','V', ... ) *> 26= | I - U U' | / ( n ulp ) SSYEV('L','V', ... ) *> 27= |D(with Z) - D(w/o Z)| / (|D| ulp) SSYEV('L','N', ... ) @@ -441,15 +410,12 @@ *> 76= | A - U S U' | / ( |A| n ulp ) SSYEVR('L','V','V', ... ) *> 77= | I - U U' | / ( n ulp ) SSYEVR('L','V','V', ... ) *> 78= |D(with Z) - D(w/o Z)| / (|D| ulp) SSYEVR('L','N','V', ... ) -*> \endverbatim -*> \verbatim +*> *> Tests 25 through 78 are repeated (as tests 79 through 132) *> with UPLO='U' -*> \endverbatim -*> \verbatim +*> *> To be added in 1999 -*> \endverbatim -*> \verbatim +*> *> 79= | A - U S U' | / ( |A| n ulp ) SSPEVR('L','V','A', ... ) *> 80= | I - U U' | / ( n ulp ) SSPEVR('L','V','A', ... ) *> 81= |D(with Z) - D(w/o Z)| / (|D| ulp) SSPEVR('L','N','A', ... ) diff --git a/TESTING/EIG/sdrvsx.f b/TESTING/EIG/sdrvsx.f index ad9355a3..a1b98e18 100644 --- a/TESTING/EIG/sdrvsx.f +++ b/TESTING/EIG/sdrvsx.f @@ -326,8 +326,7 @@ *> \param[out] WI *> \verbatim *> WI is REAL array, dimension (max(NN)) -*> \endverbatim -*> \verbatim +*> *> The real and imaginary parts of the eigenvalues of A. *> On exit, WR + WI*i are the eigenvalues of the matrix in A. *> \endverbatim @@ -340,8 +339,7 @@ *> \param[out] WIT *> \verbatim *> WIT is REAL array, dimension (max(NN)) -*> \endverbatim -*> \verbatim +*> *> Like WR, WI, these arrays contain the eigenvalues of A, *> but those computed when SGEESX only computes a partial *> eigendecomposition, i.e. not Schur vectors @@ -355,8 +353,7 @@ *> \param[out] WITMP *> \verbatim *> WITMP is REAL array, dimension (max(NN)) -*> \endverbatim -*> \verbatim +*> *> More temporary storage for eigenvalues. *> \endverbatim *> @@ -414,11 +411,9 @@ *> <0, input parameter -INFO is incorrect *> >0, SLATMR, SLATMS, SLATME or SGET24 returned an error *> code and INFO is its absolute value -*> \endverbatim -*> \verbatim +*> *>----------------------------------------------------------------------- -*> \endverbatim -*> \verbatim +*> *> Some Local Variables and Parameters: *> ---- ----- --------- --- ---------- *> ZERO, ONE Real 0 and 1. @@ -428,8 +423,7 @@ *> COND, CONDS, *> IMODE Values to be passed to the matrix generators. *> ANORM Norm of A; passed to matrix generators. -*> \endverbatim -*> \verbatim +*> *> OVFL, UNFL Overflow and underflow thresholds. *> ULP, ULPINV Finest relative precision and its inverse. *> RTULP, RTULPI Square roots of the previous 4 values. diff --git a/TESTING/EIG/sdrvvx.f b/TESTING/EIG/sdrvvx.f index 28b8b7bd..9b32d066 100644 --- a/TESTING/EIG/sdrvvx.f +++ b/TESTING/EIG/sdrvvx.f @@ -340,8 +340,7 @@ *> \param[out] WI1 *> \verbatim *> WI1 is REAL array, dimension (max(NN,12)) -*> \endverbatim -*> \verbatim +*> *> Like WR, WI, these arrays contain the eigenvalues of A, *> but those computed when SGEEVX only computes a partial *> eigendecomposition, i.e. not the eigenvalues and left @@ -475,15 +474,12 @@ *> If <0, then input paramter -INFO is incorrect. *> If >0, SLATMR, SLATMS, SLATME or SGET23 returned an error *> code, and INFO is its absolute value. -*> \endverbatim -*> \verbatim +*> *>----------------------------------------------------------------------- -*> \endverbatim -*> \verbatim +*> *> Some Local Variables and Parameters: *> ---- ----- --------- --- ---------- -*> \endverbatim -*> \verbatim +*> *> ZERO, ONE Real 0 and 1. *> MAXTYP The number of types defined. *> NMAX Largest value in NN or 12. @@ -491,13 +487,11 @@ *> COND, CONDS, *> IMODE Values to be passed to the matrix generators. *> ANORM Norm of A; passed to matrix generators. -*> \endverbatim -*> \verbatim +*> *> OVFL, UNFL Overflow and underflow thresholds. *> ULP, ULPINV Finest relative precision and its inverse. *> RTULP, RTULPI Square roots of the previous 4 values. -*> \endverbatim -*> \verbatim +*> *> The following four arrays decode JTYPE: *> KTYPE(j) The general type (1-10) for type "j". *> KMODE(j) The MODE value to be passed to the matrix diff --git a/TESTING/EIG/sget22.f b/TESTING/EIG/sget22.f index 96359a13..7db0f4ac 100644 --- a/TESTING/EIG/sget22.f +++ b/TESTING/EIG/sget22.f @@ -131,8 +131,7 @@ *> \param[in] WI *> \verbatim *> WI is REAL array, dimension (N) -*> \endverbatim -*> \verbatim +*> *> The real and imaginary parts of the eigenvalues of A. *> Purely real eigenvalues are indicated by WI(j) = 0. *> Complex conjugate pairs are indicated by WR(j)=WR(j+1) and diff --git a/TESTING/EIG/sget23.f b/TESTING/EIG/sget23.f index 5311eca9..f614b5f2 100644 --- a/TESTING/EIG/sget23.f +++ b/TESTING/EIG/sget23.f @@ -214,8 +214,7 @@ *> \param[out] WI *> \verbatim *> WI is REAL array, dimension (N) -*> \endverbatim -*> \verbatim +*> *> The real and imaginary parts of the eigenvalues of A. *> On exit, WR + WI*i are the eigenvalues of the matrix in A. *> \endverbatim @@ -228,8 +227,7 @@ *> \param[out] WI1 *> \verbatim *> WI1 is REAL array, dimension (N) -*> \endverbatim -*> \verbatim +*> *> Like WR, WI, these arrays contain the eigenvalues of A, *> but those computed when SGEEVX only computes a partial *> eigendecomposition, i.e. not the eigenvalues and left diff --git a/TESTING/EIG/sget24.f b/TESTING/EIG/sget24.f index 5bc08f99..48c35b99 100644 --- a/TESTING/EIG/sget24.f +++ b/TESTING/EIG/sget24.f @@ -205,8 +205,7 @@ *> \param[out] WI *> \verbatim *> WI is REAL array, dimension (N) -*> \endverbatim -*> \verbatim +*> *> The real and imaginary parts of the eigenvalues of A. *> On exit, WR + WI*i are the eigenvalues of the matrix in A. *> \endverbatim @@ -219,8 +218,7 @@ *> \param[out] WIT *> \verbatim *> WIT is REAL array, dimension (N) -*> \endverbatim -*> \verbatim +*> *> Like WR, WI, these arrays contain the eigenvalues of A, *> but those computed when SGEESX only computes a partial *> eigendecomposition, i.e. not Schur vectors @@ -234,8 +232,7 @@ *> \param[out] WITMP *> \verbatim *> WITMP is REAL array, dimension (N) -*> \endverbatim -*> \verbatim +*> *> Like WR, WI, these arrays contain the eigenvalues of A, *> but sorted by increasing real part. *> \endverbatim diff --git a/TESTING/EIG/sgsvts.f b/TESTING/EIG/sgsvts.f index 14418f54..7a74ca9b 100644 --- a/TESTING/EIG/sgsvts.f +++ b/TESTING/EIG/sgsvts.f @@ -141,8 +141,7 @@ *> \param[out] BETA *> \verbatim *> BETA is REAL array, dimension (N) -*> \endverbatim -*> \verbatim +*> *> The generalized singular value pairs of A and B, the *> ``diagonal'' matrices D1 and D2 are constructed from *> ALPHA and BETA, see subroutine SGGSVD for details. diff --git a/TESTING/EIG/shst01.f b/TESTING/EIG/shst01.f index b6937948..a36ac4dd 100644 --- a/TESTING/EIG/shst01.f +++ b/TESTING/EIG/shst01.f @@ -56,8 +56,7 @@ *> \param[in] IHI *> \verbatim *> IHI is INTEGER -*> \endverbatim -*> \verbatim +*> *> A is assumed to be upper triangular in rows and columns *> 1:ILO-1 and IHI+1:N, so Q differs from the identity only in *> rows and columns ILO+1:IHI. diff --git a/TESTING/EIG/slafts.f b/TESTING/EIG/slafts.f index 1566b3dd..7e49023a 100644 --- a/TESTING/EIG/slafts.f +++ b/TESTING/EIG/slafts.f @@ -41,51 +41,43 @@ *> On entry, TYPE specifies the matrix type to be used in the *> printed messages. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> N - INTEGER *> On entry, N specifies the order of the test matrix. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> IMAT - INTEGER *> On entry, IMAT specifies the type of the test matrix. *> A listing of the different types is printed by SLAHD2 *> to the output file if a test fails to pass the threshold. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> NTESTS - INTEGER *> On entry, NTESTS is the number of tests performed on the *> subroutines in the path given by TYPE. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> RESULT - REAL array of dimension( NTESTS ) *> On entry, RESULT contains the test ratios from the tests *> performed in the calling program. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> ISEED - INTEGER array of dimension( 4 ) *> Contains the random seed that generated the matrix used *> for the tests whose ratios are in RESULT. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> THRESH - REAL *> On entry, THRESH specifies the acceptable threshold of the *> test ratios. If RESULT( K ) > THRESH, then the K-th test *> did not pass the threshold and a message will be printed. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> IOUNIT - INTEGER *> On entry, IOUNIT specifies the unit number of the file *> to which the messages are printed. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> IE - INTEGER *> On entry, IE contains the number of tests which have *> failed to pass the threshold so far. diff --git a/TESTING/EIG/slahd2.f b/TESTING/EIG/slahd2.f index 325a9731..2a09afc5 100644 --- a/TESTING/EIG/slahd2.f +++ b/TESTING/EIG/slahd2.f @@ -40,15 +40,13 @@ *> PATH is CHARACTER*3. *> On entry, PATH contains the name of the path for which the *> header information is to be printed. Current paths are -*> \endverbatim -*> \verbatim +*> *> SHS, CHS: Non-symmetric eigenproblem. *> SST, CST: Symmetric eigenproblem. *> SSG, CSG: Symmetric Generalized eigenproblem. *> SBD, CBD: Singular Value Decomposition (SVD) *> SBB, CBB: General Banded reduction to bidiagonal form -*> \endverbatim -*> \verbatim +*> *> These paths also are supplied in double precision (replace *> leading S by D and leading C by Z in path names). *> \endverbatim diff --git a/TESTING/EIG/slarhs.f b/TESTING/EIG/slarhs.f index d6df8fd2..f6b1473d 100644 --- a/TESTING/EIG/slarhs.f +++ b/TESTING/EIG/slarhs.f @@ -113,12 +113,10 @@ *> KU is INTEGER *> Used only if A is a general band matrix or if A is *> triangular. -*> \endverbatim -*> \verbatim +*> *> If PATH = xGB, specifies the number of superdiagonals of A, *> and 0 <= KU <= N-1. -*> \endverbatim -*> \verbatim +*> *> If PATH = xTR, xTP, or xTB, specifies whether or not the *> matrix has unit diagonal: *> = 1: matrix has non-unit diagonal (default) diff --git a/TESTING/EIG/slatm4.f b/TESTING/EIG/slatm4.f index d3ea61e3..757231a6 100644 --- a/TESTING/EIG/slatm4.f +++ b/TESTING/EIG/slatm4.f @@ -48,8 +48,7 @@ *> If ITYPE < 0, then type abs(ITYPE) is generated and then *> swapped end for end (A(I,J) := A'(N-J,N-I).) See also *> the description of AMAGN and ISIGN. -*> \endverbatim -*> \verbatim +*> *> Special types: *> = 0: the zero matrix. *> = 1: the identity. @@ -58,8 +57,7 @@ *> followed by a k x k identity block, where k=(N-1)/2. *> If N is even, then k=(N-2)/2, and a zero diagonal entry *> is tacked onto the end. -*> \endverbatim -*> \verbatim +*> *> Diagonal types. The diagonal consists of NZ1 zeros, then *> k=N-NZ1-NZ2 nonzeros. The subdiagonal is zero. ITYPE *> specifies the nonzero diagonal entries as follows: diff --git a/TESTING/EIG/sspt21.f b/TESTING/EIG/sspt21.f index 3387ada7..dcbc4194 100644 --- a/TESTING/EIG/sspt21.f +++ b/TESTING/EIG/sspt21.f @@ -95,12 +95,10 @@ *> Specifies the type of tests to be performed. *> 1: U expressed as a dense orthogonal matrix: *> RESULT(1) = | A - U S U' | / ( |A| n ulp ) *andC> RESULT(2) = | I - UU' | / ( n ulp ) -*> \endverbatim -*> \verbatim +*> *> 2: U expressed as a product V of Housholder transformations: *> RESULT(1) = | A - V S V' | / ( |A| n ulp ) -*> \endverbatim -*> \verbatim +*> *> 3: U expressed both as a dense orthogonal matrix and *> as a product of Housholder transformations: *> RESULT(1) = | I - VU' | / ( n ulp ) diff --git a/TESTING/EIG/ssyt21.f b/TESTING/EIG/ssyt21.f index 3b38f33a..5fbde6a1 100644 --- a/TESTING/EIG/ssyt21.f +++ b/TESTING/EIG/ssyt21.f @@ -67,12 +67,10 @@ *> Specifies the type of tests to be performed. *> 1: U expressed as a dense orthogonal matrix: *> RESULT(1) = | A - U S U' | / ( |A| n ulp ) *andC> RESULT(2) = | I - UU' | / ( n ulp ) -*> \endverbatim -*> \verbatim +*> *> 2: U expressed as a product V of Housholder transformations: *> RESULT(1) = | A - V S V' | / ( |A| n ulp ) -*> \endverbatim -*> \verbatim +*> *> 3: U expressed both as a dense orthogonal matrix and *> as a product of Housholder transformations: *> RESULT(1) = | I - VU' | / ( n ulp ) diff --git a/TESTING/EIG/ssyt22.f b/TESTING/EIG/ssyt22.f index 33a54a78..7abdea83 100644 --- a/TESTING/EIG/ssyt22.f +++ b/TESTING/EIG/ssyt22.f @@ -53,100 +53,85 @@ *> Specifies the type of tests to be performed. *> 1: U expressed as a dense orthogonal matrix: *> RESULT(1) = | A - U S U' | / ( |A| n ulp ) *andC> RESULT(2) = | I - UU' | / ( n ulp ) -*> \endverbatim -*> \verbatim +*> *> UPLO CHARACTER *> If UPLO='U', the upper triangle of A will be used and the *> (strictly) lower triangle will not be referenced. If *> UPLO='L', the lower triangle of A will be used and the *> (strictly) upper triangle will not be referenced. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> N INTEGER *> The size of the matrix. If it is zero, SSYT22 does nothing. *> It must be at least zero. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> M INTEGER *> The number of columns of U. If it is zero, SSYT22 does *> nothing. It must be at least zero. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> KBAND INTEGER *> The bandwidth of the matrix. It may only be zero or one. *> If zero, then S is diagonal, and E is not referenced. If *> one, then S is symmetric tri-diagonal. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> A REAL array, dimension (LDA , N) *> The original (unfactored) matrix. It is assumed to be *> symmetric, and only the upper (UPLO='U') or only the lower *> (UPLO='L') will be referenced. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> LDA INTEGER *> The leading dimension of A. It must be at least 1 *> and at least N. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> D REAL array, dimension (N) *> The diagonal of the (symmetric tri-) diagonal matrix. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> E REAL array, dimension (N) *> The off-diagonal of the (symmetric tri-) diagonal matrix. *> E(1) is ignored, E(2) is the (1,2) and (2,1) element, etc. *> Not referenced if KBAND=0. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> U REAL array, dimension (LDU, N) *> If ITYPE=1 or 3, this contains the orthogonal matrix in *> the decomposition, expressed as a dense matrix. If ITYPE=2, *> then it is not referenced. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> LDU INTEGER *> The leading dimension of U. LDU must be at least N and *> at least 1. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> V REAL array, dimension (LDV, N) *> If ITYPE=2 or 3, the lower triangle of this array contains *> the Householder vectors used to describe the orthogonal *> matrix in the decomposition. If ITYPE=1, then it is not *> referenced. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> LDV INTEGER *> The leading dimension of V. LDV must be at least N and *> at least 1. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> TAU REAL array, dimension (N) *> If ITYPE >= 2, then TAU(j) is the scalar factor of *> v(j) v(j)' in the Householder transformation H(j) of *> the product U = H(1)...H(n-2) *> If ITYPE < 2, then TAU is not referenced. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> WORK REAL array, dimension (2*N**2) *> Workspace. *> Modified. -*> \endverbatim -*> \verbatim +*> *> RESULT REAL array, dimension (2) *> The values computed by the two tests described above. The *> values are currently limited to 1/ulp, to avoid overflow. diff --git a/TESTING/EIG/zchkbb.f b/TESTING/EIG/zchkbb.f index fd3173e5..51636112 100644 --- a/TESTING/EIG/zchkbb.f +++ b/TESTING/EIG/zchkbb.f @@ -316,11 +316,9 @@ *> \verbatim *> INFO is INTEGER *> If 0, then everything ran OK. -*> \endverbatim -*> \verbatim +*> *>----------------------------------------------------------------------- -*> \endverbatim -*> \verbatim +*> *> Some Local Variables and Parameters: *> ---- ----- --------- --- ---------- *> ZERO, ONE Real 0 and 1. @@ -334,8 +332,7 @@ *> so far. *> COND, IMODE Values to be passed to the matrix generators. *> ANORM Norm of A; passed to matrix generators. -*> \endverbatim -*> \verbatim +*> *> OVFL, UNFL Overflow and underflow thresholds. *> ULP, ULPINV Finest relative precision and its inverse. *> RTOVFL, RTUNFL Square roots of the previous 2 values. diff --git a/TESTING/EIG/zchkbd.f b/TESTING/EIG/zchkbd.f index f38e29e5..376d53f1 100644 --- a/TESTING/EIG/zchkbd.f +++ b/TESTING/EIG/zchkbd.f @@ -367,15 +367,12 @@ *> If ZLATMR, CLATMS, ZGEBRD, ZUNGBR, or ZBDSQR, *> returns an error code, the *> absolute value of it is returned. -*> \endverbatim -*> \verbatim +*> *>----------------------------------------------------------------------- -*> \endverbatim -*> \verbatim +*> *> Some Local Variables and Parameters: *> ---- ----- --------- --- ---------- -*> \endverbatim -*> \verbatim +*> *> ZERO, ONE Real 0 and 1. *> MAXTYP The number of types defined. *> NTEST The number of tests performed, or which can @@ -388,13 +385,11 @@ *> NFAIL The number of tests which have exceeded THRESH *> COND, IMODE Values to be passed to the matrix generators. *> ANORM Norm of A; passed to matrix generators. -*> \endverbatim -*> \verbatim +*> *> OVFL, UNFL Overflow and underflow thresholds. *> RTOVFL, RTUNFL Square roots of the previous 2 values. *> ULP, ULPINV Finest relative precision and its inverse. -*> \endverbatim -*> \verbatim +*> *> The following four arrays decode JTYPE: *> KTYPE(j) The general type (1-10) for type "j". *> KMODE(j) The MODE value to be passed to the matrix diff --git a/TESTING/EIG/zchkhb.f b/TESTING/EIG/zchkhb.f index 72103636..6e879f1f 100644 --- a/TESTING/EIG/zchkhb.f +++ b/TESTING/EIG/zchkhb.f @@ -254,11 +254,9 @@ *> \verbatim *> INFO is INTEGER *> If 0, then everything ran OK. -*> \endverbatim -*> \verbatim +*> *>----------------------------------------------------------------------- -*> \endverbatim -*> \verbatim +*> *> Some Local Variables and Parameters: *> ---- ----- --------- --- ---------- *> ZERO, ONE Real 0 and 1. @@ -272,8 +270,7 @@ *> so far. *> COND, IMODE Values to be passed to the matrix generators. *> ANORM Norm of A; passed to matrix generators. -*> \endverbatim -*> \verbatim +*> *> OVFL, UNFL Overflow and underflow thresholds. *> ULP, ULPINV Finest relative precision and its inverse. *> RTOVFL, RTUNFL Square roots of the previous 2 values. diff --git a/TESTING/EIG/zchkhs.f b/TESTING/EIG/zchkhs.f index cdf951e4..b6afc0a8 100644 --- a/TESTING/EIG/zchkhs.f +++ b/TESTING/EIG/zchkhs.f @@ -176,15 +176,13 @@ *> The number of sizes of matrices to use. If it is zero, *> ZCHKHS does nothing. It must be at least zero. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> NN - INTEGER array, dimension (NSIZES) *> An array containing the sizes to be used for the matrices. *> Zero values will be skipped. The values must be at least *> zero. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> NTYPES - INTEGER *> The number of elements in DOTYPE. If it is zero, ZCHKHS *> does nothing. It must be at least zero. If it is MAXTYP+1 @@ -193,8 +191,7 @@ *> is only useful if DOTYPE(1:MAXTYP) is .FALSE. and *> DOTYPE(MAXTYP+1) is .TRUE. . *> Not modified. -*> \endverbatim -*> \verbatim +*> *> DOTYPE - LOGICAL array, dimension (NTYPES) *> If DOTYPE(j) is .TRUE., then for each size in NN a *> matrix of that size and of type j will be generated. @@ -204,8 +201,7 @@ *> than MAXTYP, DOTYPE(MAXTYP+1) through DOTYPE(NTYPES) *> will be ignored. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> ISEED - INTEGER array, dimension (4) *> On entry ISEED specifies the seed of the random number *> generator. The array elements should be between 0 and 4095; @@ -217,8 +213,7 @@ *> next call to ZCHKHS to continue the same random number *> sequence. *> Modified. -*> \endverbatim -*> \verbatim +*> *> THRESH - DOUBLE PRECISION *> A test will count as "failed" if the "error", computed as *> described above, exceeds THRESH. Note that the error @@ -227,74 +222,62 @@ *> it should not depend on the precision (single vs. double) *> or the size of the matrix. It must be at least zero. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> NOUNIT - INTEGER *> The FORTRAN unit number for printing out error messages *> (e.g., if a routine returns IINFO not equal to 0.) *> Not modified. -*> \endverbatim -*> \verbatim +*> *> A - COMPLEX*16 array, dimension (LDA,max(NN)) *> Used to hold the matrix whose eigenvalues are to be *> computed. On exit, A contains the last matrix actually *> used. *> Modified. -*> \endverbatim -*> \verbatim +*> *> LDA - INTEGER *> The leading dimension of A, H, T1 and T2. It must be at *> least 1 and at least max( NN ). *> Not modified. -*> \endverbatim -*> \verbatim +*> *> H - COMPLEX*16 array, dimension (LDA,max(NN)) *> The upper hessenberg matrix computed by ZGEHRD. On exit, *> H contains the Hessenberg form of the matrix in A. *> Modified. -*> \endverbatim -*> \verbatim +*> *> T1 - COMPLEX*16 array, dimension (LDA,max(NN)) *> The Schur (="quasi-triangular") matrix computed by ZHSEQR *> if Z is computed. On exit, T1 contains the Schur form of *> the matrix in A. *> Modified. -*> \endverbatim -*> \verbatim +*> *> T2 - COMPLEX*16 array, dimension (LDA,max(NN)) *> The Schur matrix computed by ZHSEQR when Z is not computed. *> This should be identical to T1. *> Modified. -*> \endverbatim -*> \verbatim +*> *> LDU - INTEGER *> The leading dimension of U, Z, UZ and UU. It must be at *> least 1 and at least max( NN ). *> Not modified. -*> \endverbatim -*> \verbatim +*> *> U - COMPLEX*16 array, dimension (LDU,max(NN)) *> The unitary matrix computed by ZGEHRD. *> Modified. -*> \endverbatim -*> \verbatim +*> *> Z - COMPLEX*16 array, dimension (LDU,max(NN)) *> The unitary matrix computed by ZHSEQR. *> Modified. -*> \endverbatim -*> \verbatim +*> *> UZ - COMPLEX*16 array, dimension (LDU,max(NN)) *> The product of U times Z. *> Modified. -*> \endverbatim -*> \verbatim +*> *> W1 - COMPLEX*16 array, dimension (max(NN)) *> The eigenvalues of A, as computed by a full Schur *> decomposition H = Z T Z'. On exit, W1 contains the *> eigenvalues of the matrix in A. *> Modified. -*> \endverbatim -*> \verbatim +*> *> W3 - COMPLEX*16 array, dimension (max(NN)) *> The eigenvalues of A, as computed by a partial Schur *> decomposition (Z not computed, T only computed as much @@ -302,72 +285,59 @@ *> W3 contains the eigenvalues of the matrix in A, possibly *> perturbed by ZHSEIN. *> Modified. -*> \endverbatim -*> \verbatim +*> *> EVECTL - COMPLEX*16 array, dimension (LDU,max(NN)) *> The conjugate transpose of the (upper triangular) left *> eigenvector matrix for the matrix in T1. *> Modified. -*> \endverbatim -*> \verbatim +*> *> EVEZTR - COMPLEX*16 array, dimension (LDU,max(NN)) *> The (upper triangular) right eigenvector matrix for the *> matrix in T1. *> Modified. -*> \endverbatim -*> \verbatim +*> *> EVECTY - COMPLEX*16 array, dimension (LDU,max(NN)) *> The conjugate transpose of the left eigenvector matrix *> for the matrix in H. *> Modified. -*> \endverbatim -*> \verbatim +*> *> EVECTX - COMPLEX*16 array, dimension (LDU,max(NN)) *> The right eigenvector matrix for the matrix in H. *> Modified. -*> \endverbatim -*> \verbatim +*> *> UU - COMPLEX*16 array, dimension (LDU,max(NN)) *> Details of the unitary matrix computed by ZGEHRD. *> Modified. -*> \endverbatim -*> \verbatim +*> *> TAU - COMPLEX*16 array, dimension (max(NN)) *> Further details of the unitary matrix computed by ZGEHRD. *> Modified. -*> \endverbatim -*> \verbatim +*> *> WORK - COMPLEX*16 array, dimension (NWORK) *> Workspace. *> Modified. -*> \endverbatim -*> \verbatim +*> *> NWORK - INTEGER *> The number of entries in WORK. NWORK >= 4*NN(j)*NN(j) + 2. -*> \endverbatim -*> \verbatim +*> *> RWORK - DOUBLE PRECISION array, dimension (max(NN)) *> Workspace. Could be equivalenced to IWORK, but not SELECT. *> Modified. -*> \endverbatim -*> \verbatim +*> *> IWORK - INTEGER array, dimension (max(NN)) *> Workspace. *> Modified. -*> \endverbatim -*> \verbatim +*> *> SELECT - LOGICAL array, dimension (max(NN)) *> Workspace. Could be equivalenced to IWORK, but not RWORK. *> Modified. -*> \endverbatim -*> \verbatim +*> *> RESULT - DOUBLE PRECISION array, dimension (14) *> The values computed by the fourteen tests described above. *> The values are currently limited to 1/ulp, to avoid *> overflow. *> Modified. -*> \endverbatim -*> \verbatim +*> *> INFO - INTEGER *> If 0, then everything ran OK. *> -1: NSIZES < 0 @@ -384,15 +354,12 @@ *> If >2, then 30*N iterations were not enough to find an *> eigenvalue or to decompose the problem. *> Modified. -*> \endverbatim -*> \verbatim +*> *>----------------------------------------------------------------------- -*> \endverbatim -*> \verbatim +*> *> Some Local Variables and Parameters: *> ---- ----- --------- --- ---------- -*> \endverbatim -*> \verbatim +*> *> ZERO, ONE Real 0 and 1. *> MAXTYP The number of types defined. *> MTEST The number of tests defined: care must be taken @@ -411,14 +378,12 @@ *> COND, CONDS, *> IMODE Values to be passed to the matrix generators. *> ANORM Norm of A; passed to matrix generators. -*> \endverbatim -*> \verbatim +*> *> OVFL, UNFL Overflow and underflow thresholds. *> ULP, ULPINV Finest relative precision and its inverse. *> RTOVFL, RTUNFL, *> RTULP, RTULPI Square roots of the previous 4 values. -*> \endverbatim -*> \verbatim +*> *> The following four arrays decode JTYPE: *> KTYPE(j) The general type (1-10) for type "j". *> KMODE(j) The MODE value to be passed to the matrix diff --git a/TESTING/EIG/zchkst.f b/TESTING/EIG/zchkst.f index 0e0b8dbb..a0f2a8b8 100644 --- a/TESTING/EIG/zchkst.f +++ b/TESTING/EIG/zchkst.f @@ -557,11 +557,9 @@ *> If ZLATMR, CLATMS, ZHETRD, ZUNGTR, ZSTEQR, DSTERF, *> or ZUNMC2 returns an error code, the *> absolute value of it is returned. -*> \endverbatim -*> \verbatim +*> *>----------------------------------------------------------------------- -*> \endverbatim -*> \verbatim +*> *> Some Local Variables and Parameters: *> ---- ----- --------- --- ---------- *> ZERO, ONE Real 0 and 1. @@ -576,8 +574,7 @@ *> so far. *> COND, IMODE Values to be passed to the matrix generators. *> ANORM Norm of A; passed to matrix generators. -*> \endverbatim -*> \verbatim +*> *> OVFL, UNFL Overflow and underflow thresholds. *> ULP, ULPINV Finest relative precision and its inverse. *> RTOVFL, RTUNFL Square roots of the previous 2 values. diff --git a/TESTING/EIG/zdrges.f b/TESTING/EIG/zdrges.f index db51af7e..dd9c57aa 100644 --- a/TESTING/EIG/zdrges.f +++ b/TESTING/EIG/zdrges.f @@ -321,8 +321,7 @@ *> \param[out] BETA *> \verbatim *> BETA is COMPLEX*16 array, dimension (max(NN)) -*> \endverbatim -*> \verbatim +*> *> The generalized eigenvalues of (A,B) computed by ZGGES. *> ALPHA(k) / BETA(k) is the k-th generalized eigenvalue of A *> and B. diff --git a/TESTING/EIG/zdrgev.f b/TESTING/EIG/zdrgev.f index e02b8c88..166b8489 100644 --- a/TESTING/EIG/zdrgev.f +++ b/TESTING/EIG/zdrgev.f @@ -328,8 +328,7 @@ *> \param[out] BETA *> \verbatim *> BETA is COMPLEX*16 array, dimension (max(NN)) -*> \endverbatim -*> \verbatim +*> *> The generalized eigenvalues of (A,B) computed by ZGGEV. *> ( ALPHAR(k)+ALPHAI(k)*i ) / BETA(k) is the k-th *> generalized eigenvalue of A and B. @@ -343,8 +342,7 @@ *> \param[out] BETA1 *> \verbatim *> BETA1 is COMPLEX*16 array, dimension (max(NN)) -*> \endverbatim -*> \verbatim +*> *> Like ALPHAR, ALPHAI, BETA, these arrays contain the *> eigenvalues of A and B, but those computed when ZGGEV only *> computes a partial eigendecomposition, i.e. not the diff --git a/TESTING/EIG/zdrgsx.f b/TESTING/EIG/zdrgsx.f index bd82889c..60fd6760 100644 --- a/TESTING/EIG/zdrgsx.f +++ b/TESTING/EIG/zdrgsx.f @@ -268,8 +268,7 @@ *> \param[out] BETA *> \verbatim *> BETA is COMPLEX*16 array, dimension (NSIZE) -*> \endverbatim -*> \verbatim +*> *> On exit, ALPHA/BETA are the eigenvalues. *> \endverbatim *> diff --git a/TESTING/EIG/zdrgvx.f b/TESTING/EIG/zdrgvx.f index 99fe8597..4e4d2f46 100644 --- a/TESTING/EIG/zdrgvx.f +++ b/TESTING/EIG/zdrgvx.f @@ -179,8 +179,7 @@ *> \param[out] BETA *> \verbatim *> BETA is COMPLEX*16 array, dimension (NSIZE) -*> \endverbatim -*> \verbatim +*> *> On exit, ALPHA/BETA are the eigenvalues. *> \endverbatim *> diff --git a/TESTING/EIG/zdrves.f b/TESTING/EIG/zdrves.f index 9f6f1423..034f5040 100644 --- a/TESTING/EIG/zdrves.f +++ b/TESTING/EIG/zdrves.f @@ -336,11 +336,9 @@ *> -18: NWORK too small. *> If ZLATMR, CLATMS, CLATME or ZGEES returns an error code, *> the absolute value of it is returned. -*> \endverbatim -*> \verbatim +*> *>----------------------------------------------------------------------- -*> \endverbatim -*> \verbatim +*> *> Some Local Variables and Parameters: *> ---- ----- --------- --- ---------- *> ZERO, ONE Real 0 and 1. @@ -350,8 +348,7 @@ *> COND, CONDS, *> IMODE Values to be passed to the matrix generators. *> ANORM Norm of A; passed to matrix generators. -*> \endverbatim -*> \verbatim +*> *> OVFL, UNFL Overflow and underflow thresholds. *> ULP, ULPINV Finest relative precision and its inverse. *> RTULP, RTULPI Square roots of the previous 4 values. diff --git a/TESTING/EIG/zdrvev.f b/TESTING/EIG/zdrvev.f index 227b4e1a..af692d5c 100644 --- a/TESTING/EIG/zdrvev.f +++ b/TESTING/EIG/zdrvev.f @@ -346,15 +346,12 @@ *> -21: NWORK too small. *> If ZLATMR, CLATMS, CLATME or ZGEEV returns an error code, *> the absolute value of it is returned. -*> \endverbatim -*> \verbatim +*> *>----------------------------------------------------------------------- -*> \endverbatim -*> \verbatim +*> *> Some Local Variables and Parameters: *> ---- ----- --------- --- ---------- -*> \endverbatim -*> \verbatim +*> *> ZERO, ONE Real 0 and 1. *> MAXTYP The number of types defined. *> NMAX Largest value in NN. @@ -362,13 +359,11 @@ *> COND, CONDS, *> IMODE Values to be passed to the matrix generators. *> ANORM Norm of A; passed to matrix generators. -*> \endverbatim -*> \verbatim +*> *> OVFL, UNFL Overflow and underflow thresholds. *> ULP, ULPINV Finest relative precision and its inverse. *> RTULP, RTULPI Square roots of the previous 4 values. -*> \endverbatim -*> \verbatim +*> *> The following four arrays decode JTYPE: *> KTYPE(j) The general type (1-10) for type "j". *> KMODE(j) The MODE value to be passed to the matrix diff --git a/TESTING/EIG/zdrvgg.f b/TESTING/EIG/zdrvgg.f index 26c6b430..88577740 100644 --- a/TESTING/EIG/zdrvgg.f +++ b/TESTING/EIG/zdrvgg.f @@ -333,8 +333,7 @@ *> \param[out] BETA1 *> \verbatim *> BETA1 is COMPLEX*16 array, dimension (max(NN)) -*> \endverbatim -*> \verbatim +*> *> The generalized eigenvalues of (A,B) computed by ZGEGS. *> ALPHA1(k) / BETA1(k) is the k-th generalized eigenvalue of *> the matrices in A and B. @@ -348,8 +347,7 @@ *> \param[out] BETA2 *> \verbatim *> BETA2 is COMPLEX*16 array, dimension (max(NN)) -*> \endverbatim -*> \verbatim +*> *> The generalized eigenvalues of (A,B) computed by ZGEGV. *> ALPHA2(k) / BETA2(k) is the k-th generalized eigenvalue of *> the matrices in A and B. diff --git a/TESTING/EIG/zdrvsg.f b/TESTING/EIG/zdrvsg.f index 10cb0b37..d7734bf9 100644 --- a/TESTING/EIG/zdrvsg.f +++ b/TESTING/EIG/zdrvsg.f @@ -172,15 +172,13 @@ *> The number of sizes of matrices to use. If it is zero, *> ZDRVSG does nothing. It must be at least zero. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> NN INTEGER array, dimension (NSIZES) *> An array containing the sizes to be used for the matrices. *> Zero values will be skipped. The values must be at least *> zero. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> NTYPES INTEGER *> The number of elements in DOTYPE. If it is zero, ZDRVSG *> does nothing. It must be at least zero. If it is MAXTYP+1 @@ -189,8 +187,7 @@ *> is only useful if DOTYPE(1:MAXTYP) is .FALSE. and *> DOTYPE(MAXTYP+1) is .TRUE. . *> Not modified. -*> \endverbatim -*> \verbatim +*> *> DOTYPE LOGICAL array, dimension (NTYPES) *> If DOTYPE(j) is .TRUE., then for each size in NN a *> matrix of that size and of type j will be generated. @@ -200,8 +197,7 @@ *> than MAXTYP, DOTYPE(MAXTYP+1) through DOTYPE(NTYPES) *> will be ignored. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> ISEED INTEGER array, dimension (4) *> On entry ISEED specifies the seed of the random number *> generator. The array elements should be between 0 and 4095; @@ -213,8 +209,7 @@ *> next call to ZDRVSG to continue the same random number *> sequence. *> Modified. -*> \endverbatim -*> \verbatim +*> *> THRESH DOUBLE PRECISION *> A test will count as "failed" if the "error", computed as *> described above, exceeds THRESH. Note that the error @@ -223,118 +218,98 @@ *> it should not depend on the precision (single vs. double) *> or the size of the matrix. It must be at least zero. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> NOUNIT INTEGER *> The FORTRAN unit number for printing out error messages *> (e.g., if a routine returns IINFO not equal to 0.) *> Not modified. -*> \endverbatim -*> \verbatim +*> *> A COMPLEX*16 array, dimension (LDA , max(NN)) *> Used to hold the matrix whose eigenvalues are to be *> computed. On exit, A contains the last matrix actually *> used. *> Modified. -*> \endverbatim -*> \verbatim +*> *> LDA INTEGER *> The leading dimension of A. It must be at *> least 1 and at least max( NN ). *> Not modified. -*> \endverbatim -*> \verbatim +*> *> B COMPLEX*16 array, dimension (LDB , max(NN)) *> Used to hold the Hermitian positive definite matrix for *> the generailzed problem. *> On exit, B contains the last matrix actually *> used. *> Modified. -*> \endverbatim -*> \verbatim +*> *> LDB INTEGER *> The leading dimension of B. It must be at *> least 1 and at least max( NN ). *> Not modified. -*> \endverbatim -*> \verbatim +*> *> D DOUBLE PRECISION array, dimension (max(NN)) *> The eigenvalues of A. On exit, the eigenvalues in D *> correspond with the matrix in A. *> Modified. -*> \endverbatim -*> \verbatim +*> *> Z COMPLEX*16 array, dimension (LDZ, max(NN)) *> The matrix of eigenvectors. *> Modified. -*> \endverbatim -*> \verbatim +*> *> LDZ INTEGER *> The leading dimension of ZZ. It must be at least 1 and *> at least max( NN ). *> Not modified. -*> \endverbatim -*> \verbatim +*> *> AB COMPLEX*16 array, dimension (LDA, max(NN)) *> Workspace. *> Modified. -*> \endverbatim -*> \verbatim +*> *> BB COMPLEX*16 array, dimension (LDB, max(NN)) *> Workspace. *> Modified. -*> \endverbatim -*> \verbatim +*> *> AP COMPLEX*16 array, dimension (max(NN)**2) *> Workspace. *> Modified. -*> \endverbatim -*> \verbatim +*> *> BP COMPLEX*16 array, dimension (max(NN)**2) *> Workspace. *> Modified. -*> \endverbatim -*> \verbatim +*> *> WORK COMPLEX*16 array, dimension (NWORK) *> Workspace. *> Modified. -*> \endverbatim -*> \verbatim +*> *> NWORK INTEGER *> The number of entries in WORK. This must be at least *> 2*N + N**2 where N = max( NN(j), 2 ). *> Not modified. -*> \endverbatim -*> \verbatim +*> *> RWORK DOUBLE PRECISION array, dimension (LRWORK) *> Workspace. *> Modified. -*> \endverbatim -*> \verbatim +*> *> LRWORK INTEGER *> The number of entries in RWORK. This must be at least *> max( 7*N, 1 + 4*N + 2*N*lg(N) + 3*N**2 ) where *> N = max( NN(j) ) and lg( N ) = smallest integer k such *> that 2**k >= N . *> Not modified. -*> \endverbatim -*> \verbatim +*> *> IWORK INTEGER array, dimension (LIWORK)) *> Workspace. *> Modified. -*> \endverbatim -*> \verbatim +*> *> LIWORK INTEGER *> The number of entries in IWORK. This must be at least *> 2 + 5*max( NN(j) ). *> Not modified. -*> \endverbatim -*> \verbatim +*> *> RESULT DOUBLE PRECISION array, dimension (70) *> The values computed by the 70 tests described above. *> Modified. -*> \endverbatim -*> \verbatim +*> *> INFO INTEGER *> If 0, then everything ran OK. *> -1: NSIZES < 0 @@ -350,11 +325,9 @@ *> ZHPGVD, ZHEGVX, CHPGVX, ZHBGVX returns an error code, *> the absolute value of it is returned. *> Modified. -*> \endverbatim -*> \verbatim +*> *>----------------------------------------------------------------------- -*> \endverbatim -*> \verbatim +*> *> Some Local Variables and Parameters: *> ---- ----- --------- --- ---------- *> ZERO, ONE Real 0 and 1. @@ -368,8 +341,7 @@ *> so far (computed by DLAFTS). *> COND, IMODE Values to be passed to the matrix generators. *> ANORM Norm of A; passed to matrix generators. -*> \endverbatim -*> \verbatim +*> *> OVFL, UNFL Overflow and underflow thresholds. *> ULP, ULPINV Finest relative precision and its inverse. *> RTOVFL, RTUNFL Square roots of the previous 2 values. diff --git a/TESTING/EIG/zdrvst.f b/TESTING/EIG/zdrvst.f index ddc7774f..923322b2 100644 --- a/TESTING/EIG/zdrvst.f +++ b/TESTING/EIG/zdrvst.f @@ -141,15 +141,13 @@ *> The number of sizes of matrices to use. If it is zero, *> ZDRVST does nothing. It must be at least zero. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> NN INTEGER array, dimension (NSIZES) *> An array containing the sizes to be used for the matrices. *> Zero values will be skipped. The values must be at least *> zero. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> NTYPES INTEGER *> The number of elements in DOTYPE. If it is zero, ZDRVST *> does nothing. It must be at least zero. If it is MAXTYP+1 @@ -158,8 +156,7 @@ *> is only useful if DOTYPE(1:MAXTYP) is .FALSE. and *> DOTYPE(MAXTYP+1) is .TRUE. . *> Not modified. -*> \endverbatim -*> \verbatim +*> *> DOTYPE LOGICAL array, dimension (NTYPES) *> If DOTYPE(j) is .TRUE., then for each size in NN a *> matrix of that size and of type j will be generated. @@ -169,8 +166,7 @@ *> than MAXTYP, DOTYPE(MAXTYP+1) through DOTYPE(NTYPES) *> will be ignored. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> ISEED INTEGER array, dimension (4) *> On entry ISEED specifies the seed of the random number *> generator. The array elements should be between 0 and 4095; @@ -182,8 +178,7 @@ *> next call to ZDRVST to continue the same random number *> sequence. *> Modified. -*> \endverbatim -*> \verbatim +*> *> THRESH DOUBLE PRECISION *> A test will count as "failed" if the "error", computed as *> described above, exceeds THRESH. Note that the error @@ -192,121 +187,99 @@ *> it should not depend on the precision (single vs. double) *> or the size of the matrix. It must be at least zero. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> NOUNIT INTEGER *> The FORTRAN unit number for printing out error messages *> (e.g., if a routine returns IINFO not equal to 0.) *> Not modified. -*> \endverbatim -*> \verbatim +*> *> A COMPLEX*16 array, dimension (LDA , max(NN)) *> Used to hold the matrix whose eigenvalues are to be *> computed. On exit, A contains the last matrix actually *> used. *> Modified. -*> \endverbatim -*> \verbatim +*> *> LDA INTEGER *> The leading dimension of A. It must be at *> least 1 and at least max( NN ). *> Not modified. -*> \endverbatim -*> \verbatim +*> *> D1 DOUBLE PRECISION array, dimension (max(NN)) *> The eigenvalues of A, as computed by ZSTEQR simlutaneously *> with Z. On exit, the eigenvalues in D1 correspond with the *> matrix in A. *> Modified. -*> \endverbatim -*> \verbatim +*> *> D2 DOUBLE PRECISION array, dimension (max(NN)) *> The eigenvalues of A, as computed by ZSTEQR if Z is not *> computed. On exit, the eigenvalues in D2 correspond with *> the matrix in A. *> Modified. -*> \endverbatim -*> \verbatim +*> *> D3 DOUBLE PRECISION array, dimension (max(NN)) *> The eigenvalues of A, as computed by DSTERF. On exit, the *> eigenvalues in D3 correspond with the matrix in A. *> Modified. -*> \endverbatim -*> \verbatim +*> *> WA1 DOUBLE PRECISION array, dimension -*> \endverbatim -*> \verbatim +*> *> WA2 DOUBLE PRECISION array, dimension -*> \endverbatim -*> \verbatim +*> *> WA3 DOUBLE PRECISION array, dimension -*> \endverbatim -*> \verbatim +*> *> U COMPLEX*16 array, dimension (LDU, max(NN)) *> The unitary matrix computed by ZHETRD + ZUNGC3. *> Modified. -*> \endverbatim -*> \verbatim +*> *> LDU INTEGER *> The leading dimension of U, Z, and V. It must be at *> least 1 and at least max( NN ). *> Not modified. -*> \endverbatim -*> \verbatim +*> *> V COMPLEX*16 array, dimension (LDU, max(NN)) *> The Housholder vectors computed by ZHETRD in reducing A to *> tridiagonal form. *> Modified. -*> \endverbatim -*> \verbatim +*> *> TAU COMPLEX*16 array, dimension (max(NN)) *> The Householder factors computed by ZHETRD in reducing A *> to tridiagonal form. *> Modified. -*> \endverbatim -*> \verbatim +*> *> Z COMPLEX*16 array, dimension (LDU, max(NN)) *> The unitary matrix of eigenvectors computed by ZHEEVD, *> ZHEEVX, ZHPEVD, CHPEVX, ZHBEVD, and CHBEVX. *> Modified. -*> \endverbatim -*> \verbatim +*> *> WORK - COMPLEX*16 array of dimension ( LWORK ) *> Workspace. *> Modified. -*> \endverbatim -*> \verbatim +*> *> LWORK - INTEGER *> The number of entries in WORK. This must be at least *> 2*max( NN(j), 2 )**2. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> RWORK DOUBLE PRECISION array, dimension (3*max(NN)) *> Workspace. *> Modified. -*> \endverbatim -*> \verbatim +*> *> LRWORK - INTEGER *> The number of entries in RWORK. -*> \endverbatim -*> \verbatim +*> *> IWORK INTEGER array, dimension (6*max(NN)) *> Workspace. *> Modified. -*> \endverbatim -*> \verbatim +*> *> LIWORK - INTEGER *> The number of entries in IWORK. -*> \endverbatim -*> \verbatim +*> *> RESULT DOUBLE PRECISION array, dimension (??) *> The values computed by the tests described above. *> The values are currently limited to 1/ulp, to avoid *> overflow. *> Modified. -*> \endverbatim -*> \verbatim +*> *> INFO INTEGER *> If 0, then everything ran OK. *> -1: NSIZES < 0 @@ -320,11 +293,9 @@ *> or DORMC2 returns an error code, the *> absolute value of it is returned. *> Modified. -*> \endverbatim -*> \verbatim +*> *>----------------------------------------------------------------------- -*> \endverbatim -*> \verbatim +*> *> Some Local Variables and Parameters: *> ---- ----- --------- --- ---------- *> ZERO, ONE Real 0 and 1. @@ -338,8 +309,7 @@ *> so far (computed by DLAFTS). *> COND, IMODE Values to be passed to the matrix generators. *> ANORM Norm of A; passed to matrix generators. -*> \endverbatim -*> \verbatim +*> *> OVFL, UNFL Overflow and underflow thresholds. *> ULP, ULPINV Finest relative precision and its inverse. *> RTOVFL, RTUNFL Square roots of the previous 2 values. diff --git a/TESTING/EIG/zdrvsx.f b/TESTING/EIG/zdrvsx.f index 4b8f5578..39e51dd4 100644 --- a/TESTING/EIG/zdrvsx.f +++ b/TESTING/EIG/zdrvsx.f @@ -392,11 +392,9 @@ *> <0, input parameter -INFO is incorrect *> >0, ZLATMR, CLATMS, CLATME or ZGET24 returned an error *> code and INFO is its absolute value -*> \endverbatim -*> \verbatim +*> *>----------------------------------------------------------------------- -*> \endverbatim -*> \verbatim +*> *> Some Local Variables and Parameters: *> ---- ----- --------- --- ---------- *> ZERO, ONE Real 0 and 1. @@ -406,8 +404,7 @@ *> COND, CONDS, *> IMODE Values to be passed to the matrix generators. *> ANORM Norm of A; passed to matrix generators. -*> \endverbatim -*> \verbatim +*> *> OVFL, UNFL Overflow and underflow thresholds. *> ULP, ULPINV Finest relative precision and its inverse. *> RTULP, RTULPI Square roots of the previous 4 values. diff --git a/TESTING/EIG/zdrvvx.f b/TESTING/EIG/zdrvvx.f index e9794cf2..9b3763d4 100644 --- a/TESTING/EIG/zdrvvx.f +++ b/TESTING/EIG/zdrvvx.f @@ -450,15 +450,12 @@ *> If <0, then input paramter -INFO is incorrect. *> If >0, ZLATMR, CLATMS, CLATME or ZGET23 returned an error *> code, and INFO is its absolute value. -*> \endverbatim -*> \verbatim +*> *>----------------------------------------------------------------------- -*> \endverbatim -*> \verbatim +*> *> Some Local Variables and Parameters: *> ---- ----- --------- --- ---------- -*> \endverbatim -*> \verbatim +*> *> ZERO, ONE Real 0 and 1. *> MAXTYP The number of types defined. *> NMAX Largest value in NN or 12. @@ -466,13 +463,11 @@ *> COND, CONDS, *> IMODE Values to be passed to the matrix generators. *> ANORM Norm of A; passed to matrix generators. -*> \endverbatim -*> \verbatim +*> *> OVFL, UNFL Overflow and underflow thresholds. *> ULP, ULPINV Finest relative precision and its inverse. *> RTULP, RTULPI Square roots of the previous 4 values. -*> \endverbatim -*> \verbatim +*> *> The following four arrays decode JTYPE: *> KTYPE(j) The general type (1-10) for type "j". *> KMODE(j) The MODE value to be passed to the matrix diff --git a/TESTING/EIG/zgsvts.f b/TESTING/EIG/zgsvts.f index e4b6ef31..bd099c6b 100644 --- a/TESTING/EIG/zgsvts.f +++ b/TESTING/EIG/zgsvts.f @@ -140,8 +140,7 @@ *> \param[out] BETA *> \verbatim *> BETA is DOUBLE PRECISION array, dimension (N) -*> \endverbatim -*> \verbatim +*> *> The generalized singular value pairs of A and B, the *> ``diagonal'' matrices D1 and D2 are constructed from *> ALPHA and BETA, see subroutine ZGGSVD for details. diff --git a/TESTING/EIG/zhet21.f b/TESTING/EIG/zhet21.f index 0457a014..b86d05df 100644 --- a/TESTING/EIG/zhet21.f +++ b/TESTING/EIG/zhet21.f @@ -68,12 +68,10 @@ *> Specifies the type of tests to be performed. *> 1: U expressed as a dense unitary matrix: *> RESULT(1) = | A - U S U* | / ( |A| n ulp ) *andC> RESULT(2) = | I - UU* | / ( n ulp ) -*> \endverbatim -*> \verbatim +*> *> 2: U expressed as a product V of Housholder transformations: *> RESULT(1) = | A - V S V* | / ( |A| n ulp ) -*> \endverbatim -*> \verbatim +*> *> 3: U expressed both as a dense unitary matrix and *> as a product of Housholder transformations: *> RESULT(1) = | I - UV* | / ( n ulp ) diff --git a/TESTING/EIG/zhet22.f b/TESTING/EIG/zhet22.f index 42d2b68b..36affd02 100644 --- a/TESTING/EIG/zhet22.f +++ b/TESTING/EIG/zhet22.f @@ -54,104 +54,88 @@ *> Specifies the type of tests to be performed. *> 1: U expressed as a dense orthogonal matrix: *> RESULT(1) = | A - U S U' | / ( |A| n ulp ) *andC> RESULT(2) = | I - UU' | / ( n ulp ) -*> \endverbatim -*> \verbatim +*> *> UPLO CHARACTER *> If UPLO='U', the upper triangle of A will be used and the *> (strictly) lower triangle will not be referenced. If *> UPLO='L', the lower triangle of A will be used and the *> (strictly) upper triangle will not be referenced. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> N INTEGER *> The size of the matrix. If it is zero, ZHET22 does nothing. *> It must be at least zero. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> M INTEGER *> The number of columns of U. If it is zero, ZHET22 does *> nothing. It must be at least zero. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> KBAND INTEGER *> The bandwidth of the matrix. It may only be zero or one. *> If zero, then S is diagonal, and E is not referenced. If *> one, then S is symmetric tri-diagonal. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> A COMPLEX*16 array, dimension (LDA , N) *> The original (unfactored) matrix. It is assumed to be *> symmetric, and only the upper (UPLO='U') or only the lower *> (UPLO='L') will be referenced. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> LDA INTEGER *> The leading dimension of A. It must be at least 1 *> and at least N. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> D DOUBLE PRECISION array, dimension (N) *> The diagonal of the (symmetric tri-) diagonal matrix. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> E DOUBLE PRECISION array, dimension (N) *> The off-diagonal of the (symmetric tri-) diagonal matrix. *> E(1) is ignored, E(2) is the (1,2) and (2,1) element, etc. *> Not referenced if KBAND=0. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> U COMPLEX*16 array, dimension (LDU, N) *> If ITYPE=1, this contains the orthogonal matrix in *> the decomposition, expressed as a dense matrix. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> LDU INTEGER *> The leading dimension of U. LDU must be at least N and *> at least 1. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> V COMPLEX*16 array, dimension (LDV, N) *> If ITYPE=2 or 3, the lower triangle of this array contains *> the Householder vectors used to describe the orthogonal *> matrix in the decomposition. If ITYPE=1, then it is not *> referenced. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> LDV INTEGER *> The leading dimension of V. LDV must be at least N and *> at least 1. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> TAU COMPLEX*16 array, dimension (N) *> If ITYPE >= 2, then TAU(j) is the scalar factor of *> v(j) v(j)' in the Householder transformation H(j) of *> the product U = H(1)...H(n-2) *> If ITYPE < 2, then TAU is not referenced. *> Not modified. -*> \endverbatim -*> \verbatim +*> *> WORK COMPLEX*16 array, dimension (2*N**2) *> Workspace. *> Modified. -*> \endverbatim -*> \verbatim +*> *> RWORK DOUBLE PRECISION array, dimension (N) *> Workspace. *> Modified. -*> \endverbatim -*> \verbatim +*> *> RESULT DOUBLE PRECISION array, dimension (2) *> The values computed by the two tests described above. The *> values are currently limited to 1/ulp, to avoid overflow. diff --git a/TESTING/EIG/zhpt21.f b/TESTING/EIG/zhpt21.f index 63f3828c..08bf2594 100644 --- a/TESTING/EIG/zhpt21.f +++ b/TESTING/EIG/zhpt21.f @@ -93,12 +93,10 @@ *> Specifies the type of tests to be performed. *> 1: U expressed as a dense unitary matrix: *> RESULT(1) = | A - U S U* | / ( |A| n ulp ) *andC> RESULT(2) = | I - UU* | / ( n ulp ) -*> \endverbatim -*> \verbatim +*> *> 2: U expressed as a product V of Housholder transformations: *> RESULT(1) = | A - V S V* | / ( |A| n ulp ) -*> \endverbatim -*> \verbatim +*> *> 3: U expressed both as a dense unitary matrix and *> as a product of Housholder transformations: *> RESULT(1) = | I - UV* | / ( n ulp ) diff --git a/TESTING/EIG/zhst01.f b/TESTING/EIG/zhst01.f index 051d4f03..ed356665 100644 --- a/TESTING/EIG/zhst01.f +++ b/TESTING/EIG/zhst01.f @@ -57,8 +57,7 @@ *> \param[in] IHI *> \verbatim *> IHI is INTEGER -*> \endverbatim -*> \verbatim +*> *> A is assumed to be upper triangular in rows and columns *> 1:ILO-1 and IHI+1:N, so Q differs from the identity only in *> rows and columns ILO+1:IHI. diff --git a/TESTING/EIG/zlarhs.f b/TESTING/EIG/zlarhs.f index fdeccbbf..330a9af7 100644 --- a/TESTING/EIG/zlarhs.f +++ b/TESTING/EIG/zlarhs.f @@ -118,12 +118,10 @@ *> KU is INTEGER *> Used only if A is a general band matrix or if A is *> triangular. -*> \endverbatim -*> \verbatim +*> *> If PATH = xGB, specifies the number of superdiagonals of A, *> and 0 <= KU <= N-1. -*> \endverbatim -*> \verbatim +*> *> If PATH = xTR, xTP, or xTB, specifies whether or not the *> matrix has unit diagonal: *> = 1: matrix has non-unit diagonal (default) diff --git a/TESTING/EIG/zlatm4.f b/TESTING/EIG/zlatm4.f index 1153b1c3..3fb07fcf 100644 --- a/TESTING/EIG/zlatm4.f +++ b/TESTING/EIG/zlatm4.f @@ -49,8 +49,7 @@ *> If ITYPE < 0, then type abs(ITYPE) is generated and then *> swapped end for end (A(I,J) := A'(N-J,N-I).) See also *> the description of AMAGN and RSIGN. -*> \endverbatim -*> \verbatim +*> *> Special types: *> = 0: the zero matrix. *> = 1: the identity. @@ -59,8 +58,7 @@ *> followed by a k x k identity block, where k=(N-1)/2. *> If N is even, then k=(N-2)/2, and a zero diagonal entry *> is tacked onto the end. -*> \endverbatim -*> \verbatim +*> *> Diagonal types. The diagonal consists of NZ1 zeros, then *> k=N-NZ1-NZ2 nonzeros. The subdiagonal is zero. ITYPE *> specifies the nonzero diagonal entries as follows: diff --git a/TESTING/EIG/zlctsx.f b/TESTING/EIG/zlctsx.f index 7624ed5f..ea5c268e 100644 --- a/TESTING/EIG/zlctsx.f +++ b/TESTING/EIG/zlctsx.f @@ -38,8 +38,7 @@ *> \param[in] BETA *> \verbatim *> BETA is COMPLEX*16 -*> \endverbatim -*> \verbatim +*> *> parameters to decide whether the pair (ALPHA, BETA) is *> selected. *> \endverbatim diff --git a/TESTING/EIG/zsbmv.f b/TESTING/EIG/zsbmv.f index 03193b61..3f4adcf6 100644 --- a/TESTING/EIG/zsbmv.f +++ b/TESTING/EIG/zsbmv.f @@ -43,36 +43,29 @@ *> On entry, UPLO specifies whether the upper or lower *> triangular part of the band matrix A is being supplied as *> follows: -*> \endverbatim -*> \verbatim +*> *> UPLO = 'U' or 'u' The upper triangular part of A is *> being supplied. -*> \endverbatim -*> \verbatim +*> *> UPLO = 'L' or 'l' The lower triangular part of A is *> being supplied. -*> \endverbatim -*> \verbatim +*> *> Unchanged on exit. -*> \endverbatim -*> \verbatim +*> *> N - INTEGER *> On entry, N specifies the order of the matrix A. *> N must be at least zero. *> Unchanged on exit. -*> \endverbatim -*> \verbatim +*> *> K - INTEGER *> On entry, K specifies the number of super-diagonals of the *> matrix A. K must satisfy 0 .le. K. *> Unchanged on exit. -*> \endverbatim -*> \verbatim +*> *> ALPHA - COMPLEX*16 *> On entry, ALPHA specifies the scalar alpha. *> Unchanged on exit. -*> \endverbatim -*> \verbatim +*> *> A - COMPLEX*16 array, dimension( LDA, N ) *> Before entry with UPLO = 'U' or 'u', the leading ( k + 1 ) *> by n part of the array A must contain the upper triangular @@ -84,16 +77,14 @@ *> The following program segment will transfer the upper *> triangular part of a symmetric band matrix from conventional *> full matrix storage to band storage: -*> \endverbatim -*> \verbatim +*> *> DO 20, J = 1, N *> M = K + 1 - J *> DO 10, I = MAX( 1, J - K ), J *> A( M + I, J ) = matrix( I, J ) *> 10 CONTINUE *> 20 CONTINUE -*> \endverbatim -*> \verbatim +*> *> Before entry with UPLO = 'L' or 'l', the leading ( k + 1 ) *> by n part of the array A must contain the lower triangular *> band part of the symmetric matrix, supplied column by @@ -104,50 +95,42 @@ *> The following program segment will transfer the lower *> triangular part of a symmetric band matrix from conventional *> full matrix storage to band storage: -*> \endverbatim -*> \verbatim +*> *> DO 20, J = 1, N *> M = 1 - J *> DO 10, I = J, MIN( N, J + K ) *> A( M + I, J ) = matrix( I, J ) *> 10 CONTINUE *> 20 CONTINUE -*> \endverbatim -*> \verbatim +*> *> Unchanged on exit. -*> \endverbatim -*> \verbatim +*> *> LDA - INTEGER *> On entry, LDA specifies the first dimension of A as declared *> in the calling (sub) program. LDA must be at least *> ( k + 1 ). *> Unchanged on exit. -*> \endverbatim -*> \verbatim +*> *> X - COMPLEX*16 array, dimension at least *> ( 1 + ( N - 1 )*abs( INCX ) ). *> Before entry, the incremented array X must contain the *> vector x. *> Unchanged on exit. -*> \endverbatim -*> \verbatim +*> *> INCX - INTEGER *> On entry, INCX specifies the increment for the elements of *> X. INCX must not be zero. *> Unchanged on exit. -*> \endverbatim -*> \verbatim +*> *> BETA - COMPLEX*16 *> On entry, BETA specifies the scalar beta. *> Unchanged on exit. -*> \endverbatim -*> \verbatim +*> *> Y - COMPLEX*16 array, dimension at least *> ( 1 + ( N - 1 )*abs( INCY ) ). *> Before entry, the incremented array Y must contain the *> vector y. On exit, Y is overwritten by the updated vector y. -*> \endverbatim -*> \verbatim +*> *> INCY - INTEGER *> On entry, INCY specifies the increment for the elements of *> Y. INCY must not be zero. |