/thirdPartyCode/CLAPACK/TESTING/EIG/dchkee.c
C | 3127 lines | 1825 code | 172 blank | 1130 comment | 390 complexity | 28d96679865987373a96538864497978 MD5 | raw file
Possible License(s): BSD-3-Clause, MPL-2.0-no-copyleft-exception, GPL-3.0
Large files files are truncated, but you can click here to view the full file
- #include "blaswrap.h"
- /* -- translated by f2c (version 19990503).
- You must link the resulting object file with the libraries:
- -lf2c -lm (in that order)
- */
- #include "f2c.h"
- /* Common Block Declarations */
- struct {
- integer nproc, nshift, maxb;
- } cenvir_;
- #define cenvir_1 cenvir_
- struct {
- integer infot, nunit;
- logical ok, lerr;
- } infoc_;
- #define infoc_1 infoc_
- struct {
- char srnamt[6];
- } srnamc_;
- #define srnamc_1 srnamc_
- struct {
- integer selopt, seldim;
- logical selval[20];
- doublereal selwr[20], selwi[20];
- } sslct_;
- #define sslct_1 sslct_
- struct {
- integer iparms[100];
- } zlaenv_;
- #define zlaenv_1 zlaenv_
- /* Table of constant values */
- static integer c__1 = 1;
- static integer c__3 = 3;
- static integer c__5 = 5;
- static integer c__6 = 6;
- static integer c__20 = 20;
- static integer c__0 = 0;
- static integer c__132 = 132;
- static integer c__2 = 2;
- static integer c__4 = 4;
- static integer c__8 = 8;
- static integer c__87781 = 87781;
- static integer c__9 = 9;
- static integer c__25 = 25;
- static integer c__89760 = 89760;
- static integer c__18 = 18;
- static integer c__400 = 400;
- static integer c__89758 = 89758;
- static integer c__264 = 264;
- /* Main program */ MAIN__(void)
- {
- /* Initialized data */
- static char intstr[10] = "0123456789";
- static integer ioldsd[4] = { 0,0,0,1 };
- /* Format strings */
- static char fmt_9987[] = "(\002 Tests of the Nonsymmetric Eigenvalue Pro"
- "blem routines\002)";
- static char fmt_9986[] = "(\002 Tests of the Symmetric Eigenvalue Proble"
- "m routines\002)";
- static char fmt_9985[] = "(\002 Tests of the Singular Value Decompositio"
- "n routines\002)";
- static char fmt_9979[] = "(/\002 Tests of the Nonsymmetric Eigenvalue Pr"
- "oblem Driver\002,/\002 DGEEV (eigenvalues and eigevectors)"
- "\002)";
- static char fmt_9978[] = "(/\002 Tests of the Nonsymmetric Eigenvalue Pr"
- "oblem Driver\002,/\002 DGEES (Schur form)\002)";
- static char fmt_9977[] = "(/\002 Tests of the Nonsymmetric Eigenvalue Pr"
- "oblem Expert\002,\002 Driver\002,/\002 DGEEVX (eigenvalues, e"
- "igenvectors and\002,\002 condition numbers)\002)";
- static char fmt_9976[] = "(/\002 Tests of the Nonsymmetric Eigenvalue Pr"
- "oblem Expert\002,\002 Driver\002,/\002 DGEESX (Schur form and"
- " condition\002,\002 numbers)\002)";
- static char fmt_9975[] = "(/\002 Tests of the Generalized Nonsymmetric E"
- "igenvalue \002,\002Problem routines\002)";
- static char fmt_9964[] = "(/\002 Tests of the Generalized Nonsymmetric E"
- "igenvalue \002,\002Problem Driver DGGES\002)";
- static char fmt_9965[] = "(/\002 Tests of the Generalized Nonsymmetric E"
- "igenvalue \002,\002Problem Expert Driver DGGESX\002)";
- static char fmt_9963[] = "(/\002 Tests of the Generalized Nonsymmetric E"
- "igenvalue \002,\002Problem Driver DGGEV\002)";
- static char fmt_9962[] = "(/\002 Tests of the Generalized Nonsymmetric E"
- "igenvalue \002,\002Problem Expert Driver DGGEVX\002)";
- static char fmt_9974[] = "(\002 Tests of DSBTRD\002,/\002 (reduction of "
- "a symmetric band \002,\002matrix to tridiagonal form)\002)";
- static char fmt_9967[] = "(\002 Tests of DGBBRD\002,/\002 (reduction of "
- "a general band \002,\002matrix to real bidiagonal form)\002)";
- static char fmt_9971[] = "(/\002 Tests of the Generalized Linear Regress"
- "ion Model \002,\002routines\002)";
- static char fmt_9970[] = "(/\002 Tests of the Generalized QR and RQ rout"
- "ines\002)";
- static char fmt_9969[] = "(/\002 Tests of the Generalized Singular Valu"
- "e\002,\002 Decomposition routines\002)";
- static char fmt_9968[] = "(/\002 Tests of the Linear Least Squares routi"
- "nes\002)";
- static char fmt_9992[] = "(1x,a3,\002: Unrecognized path name\002)";
- static char fmt_9972[] = "(/\002 LAPACK VERSION 3.0, released June 30, 1"
- "999 \002)";
- static char fmt_9984[] = "(/\002 The following parameter values will be "
- "used:\002)";
- static char fmt_9989[] = "(\002 Invalid input value: \002,a6,\002=\002,i"
- "6,\002; must be >=\002,i6)";
- static char fmt_9988[] = "(\002 Invalid input value: \002,a6,\002=\002,i"
- "6,\002; must be <=\002,i6)";
- static char fmt_9983[] = "(4x,a6,10i6,/10x,10i6)";
- static char fmt_9981[] = "(\002 Relative machine \002,a,\002 is taken to"
- " be\002,d16.6)";
- static char fmt_9982[] = "(/\002 Routines pass computational tests if te"
- "st ratio is \002,\002less than\002,f8.2,/)";
- static char fmt_9999[] = "(/\002 Execution not attempted due to input er"
- "rors\002)";
- static char fmt_9991[] = "(//\002 *** Invalid integer value in column"
- " \002,i2,\002 of input\002,\002 line:\002,/a79)";
- static char fmt_9990[] = "(//1x,a3,\002 routines were not tested\002)";
- static char fmt_9998[] = "(//1x,a3,\002: NB =\002,i4,\002, NBMIN =\002,"
- "i4,\002, NX =\002,i4,\002, NS =\002,i4,\002, MAXB =\002,i4)";
- static char fmt_9980[] = "(\002 *** Error code from \002,a6,\002 = \002,"
- "i4)";
- static char fmt_9997[] = "(//1x,a3,\002: NB =\002,i4,\002, NBMIN =\002,"
- "i4,\002, NX =\002,i4)";
- static char fmt_9995[] = "(//1x,a3,\002: NB =\002,i4,\002, NBMIN =\002,"
- "i4,\002, NX =\002,i4,\002, NRHS =\002,i4)";
- static char fmt_9973[] = "(/1x,71(\002-\002))";
- static char fmt_9996[] = "(//1x,a3,\002: NB =\002,i4,\002, NBMIN =\002,"
- "i4,\002, NS =\002,i4,\002, MAXB =\002,i4,\002, NBCOL =\002,i4)";
- static char fmt_9966[] = "(//1x,a3,\002: NRHS =\002,i4)";
- static char fmt_9994[] = "(//\002 End of tests\002)";
- static char fmt_9993[] = "(\002 Total time used = \002,f12.2,\002 seco"
- "nds\002,/)";
- /* System generated locals */
- integer i__1;
- doublereal d__1;
- /* Builtin functions */
- integer s_rsfe(cilist *), do_fio(integer *, char *, ftnlen), e_rsfe(void);
- /* Subroutine */ int s_copy(char *, char *, ftnlen, ftnlen);
- integer s_cmp(char *, char *, ftnlen, ftnlen), s_wsfe(cilist *), e_wsfe(
- void), s_rsle(cilist *), do_lio(integer *, integer *, char *,
- ftnlen), e_rsle(void), s_wsle(cilist *), e_wsle(void);
- /* Subroutine */ int s_stop(char *, ftnlen);
- integer i_len(char *, ftnlen);
- /* Local variables */
- static char line[80];
- static doublereal taua[132];
- static integer info;
- static char path[3];
- static integer kval[20], lenp, mval[20], nval[20];
- static doublereal taub[132];
- static integer pval[20], itmp, nrhs;
- static doublereal work[87781], a[243936] /* was [17424][14] */, b[
- 87120] /* was [17424][5] */, c__[160000] /* was [400][
- 400] */, d__[1584] /* was [132][12] */;
- static integer i__, k;
- static logical fatal;
- static doublereal x[660];
- static integer iseed[4], nbcol[20], nbval[20], nbmin[20];
- static char vname[6];
- static integer newsd, nsval[20];
- static char c1[1], c3[3];
- static integer nxval[20], i1, iwork[89760];
- static doublereal s1, s2;
- static integer ic;
- extern /* Subroutine */ int dchkbb_(integer *, integer *, integer *,
- integer *, integer *, integer *, logical *, integer *, integer *,
- doublereal *, integer *, doublereal *, integer *, doublereal *,
- integer *, doublereal *, doublereal *, doublereal *, integer *,
- doublereal *, integer *, doublereal *, integer *, doublereal *,
- doublereal *, integer *, doublereal *, integer *), dchkbd_(
- integer *, integer *, integer *, integer *, logical *, integer *,
- integer *, doublereal *, doublereal *, integer *, doublereal *,
- doublereal *, doublereal *, doublereal *, doublereal *, integer *,
- doublereal *, doublereal *, doublereal *, integer *, doublereal *
- , integer *, doublereal *, doublereal *, doublereal *, integer *,
- integer *, integer *, integer *), dchkec_(doublereal *, logical *,
- integer *, integer *), dchkbk_(integer *, integer *), dchkbl_(
- integer *, integer *);
- extern doublereal dlamch_(char *);
- static integer nk;
- extern /* Subroutine */ int dchkgg_(integer *, integer *, integer *,
- logical *, integer *, doublereal *, logical *, doublereal *,
- integer *, doublereal *, integer *, doublereal *, doublereal *,
- doublereal *, doublereal *, doublereal *, doublereal *,
- doublereal *, doublereal *, integer *, doublereal *, doublereal *,
- doublereal *, doublereal *, doublereal *, doublereal *,
- doublereal *, doublereal *, doublereal *, doublereal *,
- doublereal *, doublereal *, integer *, logical *, doublereal *,
- integer *), dchkgk_(integer *, integer *);
- static integer nn;
- extern /* Subroutine */ int dchkgl_(integer *, integer *), dchksb_(
- integer *, integer *, integer *, integer *, integer *, logical *,
- integer *, doublereal *, integer *, doublereal *, integer *,
- doublereal *, doublereal *, doublereal *, integer *, doublereal *,
- integer *, doublereal *, integer *);
- extern doublereal dsecnd_(void);
- extern /* Subroutine */ int dckglm_(integer *, integer *, integer *,
- integer *, integer *, integer *, doublereal *, integer *,
- doublereal *, doublereal *, doublereal *, doublereal *,
- doublereal *, doublereal *, doublereal *, integer *, integer *,
- integer *), derrbd_(char *, integer *), dchkhs_(integer *,
- integer *, integer *, logical *, integer *, doublereal *,
- integer *, doublereal *, integer *, doublereal *, doublereal *,
- doublereal *, doublereal *, integer *, doublereal *, doublereal *,
- doublereal *, doublereal *, doublereal *, doublereal *,
- doublereal *, doublereal *, doublereal *, doublereal *,
- doublereal *, doublereal *, doublereal *, integer *, integer *,
- logical *, doublereal *, integer *), alareq_(char *, integer *,
- logical *, integer *, integer *, integer *), dcklse_(
- integer *, integer *, integer *, integer *, integer *, integer *,
- doublereal *, integer *, doublereal *, doublereal *, doublereal *,
- doublereal *, doublereal *, doublereal *, doublereal *, integer *
- , integer *, integer *), ddrvbd_(integer *, integer *, integer *,
- integer *, logical *, integer *, doublereal *, doublereal *,
- integer *, doublereal *, integer *, doublereal *, integer *,
- doublereal *, doublereal *, doublereal *, doublereal *,
- doublereal *, doublereal *, doublereal *, integer *, integer *,
- integer *, integer *), ddrges_(integer *, integer *, integer *,
- logical *, integer *, doublereal *, integer *, doublereal *,
- integer *, doublereal *, doublereal *, doublereal *, doublereal *,
- integer *, doublereal *, doublereal *, doublereal *, doublereal *
- , doublereal *, integer *, doublereal *, logical *, integer *),
- derred_(char *, integer *), derrgg_(char *, integer *), dckgqr_(integer *, integer *, integer *, integer *,
- integer *, integer *, integer *, integer *, doublereal *, integer
- *, doublereal *, doublereal *, doublereal *, doublereal *,
- doublereal *, doublereal *, doublereal *, doublereal *,
- doublereal *, doublereal *, doublereal *, doublereal *,
- doublereal *, integer *, integer *, integer *), ddrgev_(integer *,
- integer *, integer *, logical *, integer *, doublereal *,
- integer *, doublereal *, integer *, doublereal *, doublereal *,
- doublereal *, doublereal *, integer *, doublereal *, doublereal *,
- integer *, doublereal *, doublereal *, doublereal *, doublereal *
- , doublereal *, doublereal *, doublereal *, integer *, doublereal
- *, integer *), ddrvgg_(integer *, integer *, integer *, logical *,
- integer *, doublereal *, doublereal *, integer *, doublereal *,
- integer *, doublereal *, doublereal *, doublereal *, doublereal *,
- doublereal *, doublereal *, integer *, doublereal *, doublereal *
- , doublereal *, doublereal *, doublereal *, doublereal *,
- doublereal *, doublereal *, doublereal *, doublereal *, integer *,
- doublereal *, integer *);
- extern logical lsamen_(integer *, char *, char *);
- extern /* Subroutine */ int dchkst_(integer *, integer *, integer *,
- logical *, integer *, doublereal *, integer *, doublereal *,
- integer *, doublereal *, doublereal *, doublereal *, doublereal *,
- doublereal *, doublereal *, doublereal *, doublereal *,
- doublereal *, doublereal *, doublereal *, doublereal *,
- doublereal *, integer *, doublereal *, doublereal *, doublereal *,
- doublereal *, doublereal *, integer *, integer *, integer *,
- doublereal *, integer *), dckgsv_(integer *, integer *, integer *,
- integer *, integer *, integer *, doublereal *, integer *,
- doublereal *, doublereal *, doublereal *, doublereal *,
- doublereal *, doublereal *, doublereal *, doublereal *,
- doublereal *, doublereal *, integer *, doublereal *, doublereal *,
- integer *, integer *, integer *), ddrves_(integer *, integer *,
- integer *, logical *, integer *, doublereal *, integer *,
- doublereal *, integer *, doublereal *, doublereal *, doublereal *,
- doublereal *, doublereal *, doublereal *, doublereal *, integer *
- , doublereal *, doublereal *, integer *, integer *, logical *,
- integer *), derrhs_(char *, integer *);
- static integer mxbval[20];
- extern /* Subroutine */ int ddrvev_(integer *, integer *, integer *,
- logical *, integer *, doublereal *, integer *, doublereal *,
- integer *, doublereal *, doublereal *, doublereal *, doublereal *,
- doublereal *, doublereal *, integer *, doublereal *, integer *,
- doublereal *, integer *, doublereal *, doublereal *, integer *,
- integer *, integer *), ddrgsx_(integer *, integer *, doublereal *,
- integer *, integer *, doublereal *, integer *, doublereal *,
- doublereal *, doublereal *, doublereal *, doublereal *,
- doublereal *, doublereal *, doublereal *, doublereal *, integer *,
- doublereal *, doublereal *, integer *, integer *, integer *,
- logical *, integer *), ddrvsg_(integer *, integer *, integer *,
- logical *, integer *, doublereal *, integer *, doublereal *,
- integer *, doublereal *, integer *, doublereal *, doublereal *,
- integer *, doublereal *, doublereal *, doublereal *, doublereal *,
- doublereal *, integer *, integer *, integer *, doublereal *,
- integer *);
- static logical tstdif;
- static doublereal thresh;
- extern /* Subroutine */ int ddrgvx_(integer *, doublereal *, integer *,
- integer *, doublereal *, integer *, doublereal *, doublereal *,
- doublereal *, doublereal *, doublereal *, doublereal *,
- doublereal *, doublereal *, integer *, integer *, doublereal *,
- doublereal *, doublereal *, doublereal *, doublereal *,
- doublereal *, doublereal *, integer *, integer *, integer *,
- doublereal *, logical *, integer *);
- static logical tstchk;
- static integer nparms;
- extern /* Subroutine */ int xlaenv_(integer *, integer *), derrst_(char *,
- integer *);
- static logical dotype[30], logwrk[132];
- static doublereal thrshn;
- extern /* Subroutine */ int ddrvst_(integer *, integer *, integer *,
- logical *, integer *, doublereal *, integer *, doublereal *,
- integer *, doublereal *, doublereal *, doublereal *, doublereal *,
- doublereal *, doublereal *, doublereal *, doublereal *,
- doublereal *, integer *, doublereal *, doublereal *, doublereal *,
- doublereal *, integer *, integer *, integer *, doublereal *,
- integer *), ddrvsx_(integer *, integer *, integer *, logical *,
- integer *, doublereal *, integer *, integer *, doublereal *,
- integer *, doublereal *, doublereal *, doublereal *, doublereal *,
- doublereal *, doublereal *, doublereal *, doublereal *,
- doublereal *, integer *, doublereal *, doublereal *, doublereal *,
- integer *, integer *, logical *, integer *), ddrvvx_(integer *,
- integer *, integer *, logical *, integer *, doublereal *, integer
- *, integer *, doublereal *, integer *, doublereal *, doublereal *,
- doublereal *, doublereal *, doublereal *, doublereal *, integer *
- , doublereal *, integer *, doublereal *, integer *, doublereal *,
- doublereal *, doublereal *, doublereal *, doublereal *,
- doublereal *, doublereal *, doublereal *, doublereal *,
- doublereal *, integer *, integer *, integer *);
- static doublereal result[500];
- static integer maxtyp;
- static logical tsterr;
- static integer ntypes;
- static logical tstdrv, dbb, dbk, dgg, dbl, dgk, dgl, dsb, des, dgs, dev,
- glm, dgv, nep, lse, dgx, sep;
- static doublereal eps;
- static logical gqr, svd, dsx, gsv, dvx, dxv;
- /* Fortran I/O blocks */
- static cilist io___5 = { 0, 5, 1, "(A80)", 0 };
- static cilist io___30 = { 0, 6, 0, fmt_9987, 0 };
- static cilist io___31 = { 0, 6, 0, fmt_9986, 0 };
- static cilist io___32 = { 0, 6, 0, fmt_9985, 0 };
- static cilist io___33 = { 0, 6, 0, fmt_9979, 0 };
- static cilist io___34 = { 0, 6, 0, fmt_9978, 0 };
- static cilist io___35 = { 0, 6, 0, fmt_9977, 0 };
- static cilist io___36 = { 0, 6, 0, fmt_9976, 0 };
- static cilist io___37 = { 0, 6, 0, fmt_9975, 0 };
- static cilist io___38 = { 0, 6, 0, fmt_9964, 0 };
- static cilist io___39 = { 0, 6, 0, fmt_9965, 0 };
- static cilist io___40 = { 0, 6, 0, fmt_9963, 0 };
- static cilist io___41 = { 0, 6, 0, fmt_9962, 0 };
- static cilist io___42 = { 0, 6, 0, fmt_9974, 0 };
- static cilist io___43 = { 0, 6, 0, fmt_9967, 0 };
- static cilist io___44 = { 0, 6, 0, fmt_9971, 0 };
- static cilist io___45 = { 0, 6, 0, fmt_9970, 0 };
- static cilist io___46 = { 0, 6, 0, fmt_9969, 0 };
- static cilist io___47 = { 0, 6, 0, fmt_9968, 0 };
- static cilist io___48 = { 0, 5, 0, 0, 0 };
- static cilist io___51 = { 0, 6, 0, fmt_9992, 0 };
- static cilist io___52 = { 0, 6, 0, fmt_9972, 0 };
- static cilist io___53 = { 0, 6, 0, fmt_9984, 0 };
- static cilist io___54 = { 0, 5, 0, 0, 0 };
- static cilist io___56 = { 0, 6, 0, fmt_9989, 0 };
- static cilist io___57 = { 0, 6, 0, fmt_9988, 0 };
- static cilist io___58 = { 0, 5, 0, 0, 0 };
- static cilist io___62 = { 0, 6, 0, fmt_9989, 0 };
- static cilist io___63 = { 0, 6, 0, fmt_9988, 0 };
- static cilist io___64 = { 0, 6, 0, fmt_9983, 0 };
- static cilist io___65 = { 0, 5, 0, 0, 0 };
- static cilist io___67 = { 0, 6, 0, fmt_9989, 0 };
- static cilist io___68 = { 0, 6, 0, fmt_9988, 0 };
- static cilist io___69 = { 0, 6, 0, fmt_9983, 0 };
- static cilist io___70 = { 0, 5, 0, 0, 0 };
- static cilist io___72 = { 0, 6, 0, fmt_9989, 0 };
- static cilist io___73 = { 0, 6, 0, fmt_9988, 0 };
- static cilist io___74 = { 0, 6, 0, fmt_9983, 0 };
- static cilist io___75 = { 0, 6, 0, fmt_9983, 0 };
- static cilist io___76 = { 0, 5, 0, 0, 0 };
- static cilist io___78 = { 0, 5, 0, 0, 0 };
- static cilist io___80 = { 0, 6, 0, fmt_9989, 0 };
- static cilist io___81 = { 0, 6, 0, fmt_9988, 0 };
- static cilist io___82 = { 0, 6, 0, fmt_9983, 0 };
- static cilist io___83 = { 0, 5, 0, 0, 0 };
- static cilist io___89 = { 0, 6, 0, fmt_9989, 0 };
- static cilist io___90 = { 0, 6, 0, fmt_9989, 0 };
- static cilist io___91 = { 0, 6, 0, fmt_9989, 0 };
- static cilist io___92 = { 0, 6, 0, fmt_9989, 0 };
- static cilist io___93 = { 0, 6, 0, fmt_9989, 0 };
- static cilist io___94 = { 0, 6, 0, fmt_9983, 0 };
- static cilist io___95 = { 0, 6, 0, fmt_9983, 0 };
- static cilist io___96 = { 0, 6, 0, fmt_9983, 0 };
- static cilist io___97 = { 0, 6, 0, fmt_9983, 0 };
- static cilist io___98 = { 0, 6, 0, fmt_9983, 0 };
- static cilist io___99 = { 0, 5, 0, 0, 0 };
- static cilist io___101 = { 0, 6, 0, fmt_9989, 0 };
- static cilist io___102 = { 0, 6, 0, fmt_9988, 0 };
- static cilist io___103 = { 0, 5, 0, 0, 0 };
- static cilist io___104 = { 0, 6, 0, fmt_9989, 0 };
- static cilist io___105 = { 0, 6, 0, fmt_9988, 0 };
- static cilist io___106 = { 0, 6, 0, fmt_9983, 0 };
- static cilist io___107 = { 0, 5, 0, 0, 0 };
- static cilist io___108 = { 0, 6, 0, fmt_9989, 0 };
- static cilist io___109 = { 0, 6, 0, fmt_9988, 0 };
- static cilist io___110 = { 0, 6, 0, fmt_9983, 0 };
- static cilist io___111 = { 0, 5, 0, 0, 0 };
- static cilist io___112 = { 0, 6, 0, fmt_9989, 0 };
- static cilist io___113 = { 0, 6, 0, fmt_9988, 0 };
- static cilist io___114 = { 0, 6, 0, fmt_9983, 0 };
- static cilist io___115 = { 0, 5, 0, 0, 0 };
- static cilist io___116 = { 0, 6, 0, fmt_9989, 0 };
- static cilist io___117 = { 0, 6, 0, fmt_9988, 0 };
- static cilist io___118 = { 0, 6, 0, fmt_9983, 0 };
- static cilist io___119 = { 0, 5, 0, 0, 0 };
- static cilist io___120 = { 0, 6, 0, fmt_9989, 0 };
- static cilist io___121 = { 0, 6, 0, fmt_9988, 0 };
- static cilist io___122 = { 0, 6, 0, fmt_9983, 0 };
- static cilist io___123 = { 0, 5, 0, 0, 0 };
- static cilist io___125 = { 0, 6, 0, fmt_9989, 0 };
- static cilist io___126 = { 0, 6, 0, fmt_9988, 0 };
- static cilist io___127 = { 0, 6, 0, fmt_9983, 0 };
- static cilist io___128 = { 0, 6, 0, 0, 0 };
- static cilist io___130 = { 0, 6, 0, fmt_9981, 0 };
- static cilist io___131 = { 0, 6, 0, fmt_9981, 0 };
- static cilist io___132 = { 0, 6, 0, fmt_9981, 0 };
- static cilist io___133 = { 0, 5, 0, 0, 0 };
- static cilist io___134 = { 0, 6, 0, fmt_9982, 0 };
- static cilist io___135 = { 0, 5, 0, 0, 0 };
- static cilist io___137 = { 0, 5, 0, 0, 0 };
- static cilist io___139 = { 0, 5, 0, 0, 0 };
- static cilist io___140 = { 0, 5, 0, 0, 0 };
- static cilist io___142 = { 0, 5, 0, 0, 0 };
- static cilist io___144 = { 0, 6, 0, fmt_9999, 0 };
- static cilist io___145 = { 0, 5, 1, "(A80)", 0 };
- static cilist io___154 = { 0, 6, 0, fmt_9991, 0 };
- static cilist io___155 = { 0, 6, 0, fmt_9990, 0 };
- static cilist io___158 = { 0, 6, 0, fmt_9998, 0 };
- static cilist io___166 = { 0, 6, 0, fmt_9980, 0 };
- static cilist io___167 = { 0, 6, 0, fmt_9997, 0 };
- static cilist io___168 = { 0, 6, 0, fmt_9980, 0 };
- static cilist io___169 = { 0, 6, 0, fmt_9980, 0 };
- static cilist io___170 = { 0, 6, 0, fmt_9997, 0 };
- static cilist io___171 = { 0, 6, 0, fmt_9980, 0 };
- static cilist io___173 = { 0, 6, 0, fmt_9995, 0 };
- static cilist io___174 = { 0, 6, 0, fmt_9980, 0 };
- static cilist io___175 = { 0, 6, 0, fmt_9990, 0 };
- static cilist io___176 = { 0, 6, 0, fmt_9980, 0 };
- static cilist io___177 = { 0, 6, 0, fmt_9973, 0 };
- static cilist io___178 = { 0, 6, 0, fmt_9990, 0 };
- static cilist io___179 = { 0, 6, 0, fmt_9980, 0 };
- static cilist io___180 = { 0, 6, 0, fmt_9973, 0 };
- static cilist io___181 = { 0, 6, 0, fmt_9990, 0 };
- static cilist io___182 = { 0, 6, 0, fmt_9980, 0 };
- static cilist io___183 = { 0, 6, 0, fmt_9973, 0 };
- static cilist io___184 = { 0, 6, 0, fmt_9990, 0 };
- static cilist io___185 = { 0, 6, 0, fmt_9980, 0 };
- static cilist io___186 = { 0, 6, 0, fmt_9973, 0 };
- static cilist io___187 = { 0, 6, 0, fmt_9996, 0 };
- static cilist io___190 = { 0, 6, 0, fmt_9980, 0 };
- static cilist io___191 = { 0, 6, 0, fmt_9980, 0 };
- static cilist io___192 = { 0, 6, 0, fmt_9990, 0 };
- static cilist io___193 = { 0, 6, 0, fmt_9980, 0 };
- static cilist io___194 = { 0, 6, 0, fmt_9973, 0 };
- static cilist io___195 = { 0, 6, 0, fmt_9990, 0 };
- static cilist io___197 = { 0, 6, 0, fmt_9980, 0 };
- static cilist io___198 = { 0, 6, 0, fmt_9973, 0 };
- static cilist io___199 = { 0, 6, 0, fmt_9990, 0 };
- static cilist io___200 = { 0, 6, 0, fmt_9980, 0 };
- static cilist io___201 = { 0, 6, 0, fmt_9973, 0 };
- static cilist io___202 = { 0, 6, 0, fmt_9990, 0 };
- static cilist io___203 = { 0, 6, 0, fmt_9980, 0 };
- static cilist io___204 = { 0, 6, 0, fmt_9973, 0 };
- static cilist io___205 = { 0, 6, 0, fmt_9980, 0 };
- static cilist io___206 = { 0, 6, 0, fmt_9966, 0 };
- static cilist io___207 = { 0, 6, 0, fmt_9980, 0 };
- static cilist io___210 = { 0, 6, 0, fmt_9980, 0 };
- static cilist io___213 = { 0, 6, 0, fmt_9980, 0 };
- static cilist io___214 = { 0, 6, 0, fmt_9980, 0 };
- static cilist io___215 = { 0, 6, 0, fmt_9980, 0 };
- static cilist io___216 = { 0, 6, 0, 0, 0 };
- static cilist io___217 = { 0, 6, 0, 0, 0 };
- static cilist io___218 = { 0, 6, 0, fmt_9992, 0 };
- static cilist io___219 = { 0, 6, 0, fmt_9994, 0 };
- static cilist io___221 = { 0, 6, 0, fmt_9993, 0 };
- #define a_ref(a_1,a_2) a[(a_2)*17424 + a_1 - 17425]
- #define b_ref(a_1,a_2) b[(a_2)*17424 + a_1 - 17425]
- #define c___ref(a_1,a_2) c__[(a_2)*400 + a_1 - 401]
- #define d___ref(a_1,a_2) d__[(a_2)*132 + a_1 - 133]
- /* -- LAPACK test routine (version 3.0) --
- Univ. of Tennessee, Univ. of California Berkeley, NAG Ltd.,
- Courant Institute, Argonne National Lab, and Rice University
- June 30, 1999
- Purpose
- =======
- DCHKEE tests the DOUBLE PRECISION LAPACK subroutines for the matrix
- eigenvalue problem. The test paths in this version are
- NEP (Nonsymmetric Eigenvalue Problem):
- Test DGEHRD, DORGHR, DHSEQR, DTREVC, DHSEIN, and DORMHR
- SEP (Symmetric Eigenvalue Problem):
- Test DSYTRD, DORGTR, DSTEQR, DSTERF, DSTEIN, DSTEDC,
- and drivers DSYEV(X), DSBEV(X), DSPEV(X), DSTEV(X),
- DSYEVD, DSBEVD, DSPEVD, DSTEVD
- SVD (Singular Value Decomposition):
- Test DGEBRD, DORGBR, DBDSQR, DBDSDC
- and the drivers DGESVD, DGESDD
- DEV (Nonsymmetric Eigenvalue/eigenvector Driver):
- Test DGEEV
- DES (Nonsymmetric Schur form Driver):
- Test DGEES
- DVX (Nonsymmetric Eigenvalue/eigenvector Expert Driver):
- Test DGEEVX
- DSX (Nonsymmetric Schur form Expert Driver):
- Test DGEESX
- DGG (Generalized Nonsymmetric Eigenvalue Problem):
- Test DGGHRD, DGGBAL, DGGBAK, DHGEQZ, and DTGEVC
- and the driver routines DGEGS and DGEGV
- DGS (Generalized Nonsymmetric Schur form Driver):
- Test DGGES
- DGV (Generalized Nonsymmetric Eigenvalue/eigenvector Driver):
- Test DGGEV
- DGX (Generalized Nonsymmetric Schur form Expert Driver):
- Test DGGESX
- DXV (Generalized Nonsymmetric Eigenvalue/eigenvector Expert Driver):
- Test DGGEVX
- DSG (Symmetric Generalized Eigenvalue Problem):
- Test DSYGST, DSYGV, DSYGVD, DSYGVX, DSPGST, DSPGV, DSPGVD,
- DSPGVX, DSBGST, DSBGV, DSBGVD, and DSBGVX
- DSB (Symmetric Band Eigenvalue Problem):
- Test DSBTRD
- DBB (Band Singular Value Decomposition):
- Test DGBBRD
- DEC (Eigencondition estimation):
- Test DLALN2, DLASY2, DLAEQU, DLAEXC, DTRSYL, DTREXC, DTRSNA,
- DTRSEN, and DLAQTR
- DBL (Balancing a general matrix)
- Test DGEBAL
- DBK (Back transformation on a balanced matrix)
- Test DGEBAK
- DGL (Balancing a matrix pair)
- Test DGGBAL
- DGK (Back transformation on a matrix pair)
- Test DGGBAK
- GLM (Generalized Linear Regression Model):
- Tests DGGGLM
- GQR (Generalized QR and RQ factorizations):
- Tests DGGQRF and DGGRQF
- GSV (Generalized Singular Value Decomposition):
- Tests DGGSVD, DGGSVP, DTGSJA, DLAGS2, DLAPLL, and DLAPMT
- LSE (Constrained Linear Least Squares):
- Tests DGGLSE
- Each test path has a different set of inputs, but the data sets for
- the driver routines xEV, xES, xVX, and xSX can be concatenated in a
- single input file. The first line of input should contain one of the
- 3-character path names in columns 1-3. The number of remaining lines
- depends on what is found on the first line.
- The number of matrix types used in testing is often controllable from
- the input file. The number of matrix types for each path, and the
- test routine that describes them, is as follows:
- Path name(s) Types Test routine
- DHS or NEP 21 DCHKHS
- DST or SEP 21 DCHKST (routines)
- 18 DDRVST (drivers)
- DBD or SVD 16 DCHKBD (routines)
- 5 DDRVBD (drivers)
- DEV 21 DDRVEV
- DES 21 DDRVES
- DVX 21 DDRVVX
- DSX 21 DDRVSX
- DGG 26 DCHKGG (routines)
- 26 DDRVGG (drivers)
- DGS 26 DDRGES
- DGX 5 DDRGSX
- DGV 26 DDRGEV
- DXV 2 DDRGVX
- DSG 21 DDRVSG
- DSB 15 DCHKSB
- DBB 15 DCHKBB
- DEC - DCHKEC
- DBL - DCHKBL
- DBK - DCHKBK
- DGL - DCHKGL
- DGK - DCHKGK
- GLM 8 DCKGLM
- GQR 8 DCKGQR
- GSV 8 DCKGSV
- LSE 8 DCKLSE
- -----------------------------------------------------------------------
- NEP input file:
- line 2: NN, INTEGER
- Number of values of N.
- line 3: NVAL, INTEGER array, dimension (NN)
- The values for the matrix dimension N.
- line 4: NPARMS, INTEGER
- Number of values of the parameters NB, NBMIN, NX, NS, and
- MAXB.
- line 5: NBVAL, INTEGER array, dimension (NPARMS)
- The values for the blocksize NB.
- line 6: NBMIN, INTEGER array, dimension (NPARMS)
- The values for the minimum blocksize NBMIN.
- line 7: NXVAL, INTEGER array, dimension (NPARMS)
- The values for the crossover point NX.
- line 8: NSVAL, INTEGER array, dimension (NPARMS)
- The values for the number of shifts.
- line 9: MXBVAL, INTEGER array, dimension (NPARMS)
- The values for MAXB, used in determining minimum blocksize.
- line 10: THRESH
- Threshold value for the test ratios. Information will be
- printed about each test for which the test ratio is greater
- than or equal to the threshold. To have all of the test
- ratios printed, use THRESH = 0.0 .
- line 11: NEWSD, INTEGER
- A code indicating how to set the random number seed.
- = 0: Set the seed to a default value before each run
- = 1: Initialize the seed to a default value only before the
- first run
- = 2: Like 1, but use the seed values on the next line
- If line 11 was 2:
- line 12: INTEGER array, dimension (4)
- Four integer values for the random number seed.
- lines 12-EOF: The remaining lines occur in sets of 1 or 2 and allow
- the user to specify the matrix types. Each line contains
- a 3-character path name in columns 1-3, and the number
- of matrix types must be the first nonblank item in columns
- 4-80. If the number of matrix types is at least 1 but is
- less than the maximum number of possible types, a second
- line will be read to get the numbers of the matrix types to
- be used. For example,
- NEP 21
- requests all of the matrix types for the nonsymmetric
- eigenvalue problem, while
- NEP 4
- 9 10 11 12
- requests only matrices of type 9, 10, 11, and 12.
- The valid 3-character path names are 'NEP' or 'SHS' for the
- nonsymmetric eigenvalue routines.
- -----------------------------------------------------------------------
- SEP or DSG input file:
- line 2: NN, INTEGER
- Number of values of N.
- line 3: NVAL, INTEGER array, dimension (NN)
- The values for the matrix dimension N.
- line 4: NPARMS, INTEGER
- Number of values of the parameters NB, NBMIN, and NX.
- line 5: NBVAL, INTEGER array, dimension (NPARMS)
- The values for the blocksize NB.
- line 6: NBMIN, INTEGER array, dimension (NPARMS)
- The values for the minimum blocksize NBMIN.
- line 7: NXVAL, INTEGER array, dimension (NPARMS)
- The values for the crossover point NX.
- line 8: THRESH
- Threshold value for the test ratios. Information will be
- printed about each test for which the test ratio is greater
- than or equal to the threshold.
- line 9: TSTCHK, LOGICAL
- Flag indicating whether or not to test the LAPACK routines.
- line 10: TSTDRV, LOGICAL
- Flag indicating whether or not to test the driver routines.
- line 11: TSTERR, LOGICAL
- Flag indicating whether or not to test the error exits for
- the LAPACK routines and driver routines.
- line 12: NEWSD, INTEGER
- A code indicating how to set the random number seed.
- = 0: Set the seed to a default value before each run
- = 1: Initialize the seed to a default value only before the
- first run
- = 2: Like 1, but use the seed values on the next line
- If line 12 was 2:
- line 13: INTEGER array, dimension (4)
- Four integer values for the random number seed.
- lines 13-EOF: Lines specifying matrix types, as for NEP.
- The 3-character path names are 'SEP' or 'SST' for the
- symmetric eigenvalue routines and driver routines, and
- 'DSG' for the routines for the symmetric generalized
- eigenvalue problem.
- -----------------------------------------------------------------------
- SVD input file:
- line 2: NN, INTEGER
- Number of values of M and N.
- line 3: MVAL, INTEGER array, dimension (NN)
- The values for the matrix row dimension M.
- line 4: NVAL, INTEGER array, dimension (NN)
- The values for the matrix column dimension N.
- line 5: NPARMS, INTEGER
- Number of values of the parameter NB, NBMIN, NX, and NRHS.
- line 6: NBVAL, INTEGER array, dimension (NPARMS)
- The values for the blocksize NB.
- line 7: NBMIN, INTEGER array, dimension (NPARMS)
- The values for the minimum blocksize NBMIN.
- line 8: NXVAL, INTEGER array, dimension (NPARMS)
- The values for the crossover point NX.
- line 9: NSVAL, INTEGER array, dimension (NPARMS)
- The values for the number of right hand sides NRHS.
- line 10: THRESH
- Threshold value for the test ratios. Information will be
- printed about each test for which the test ratio is greater
- than or equal to the threshold.
- line 11: TSTCHK, LOGICAL
- Flag indicating whether or not to test the LAPACK routines.
- line 12: TSTDRV, LOGICAL
- Flag indicating whether or not to test the driver routines.
- line 13: TSTERR, LOGICAL
- Flag indicating whether or not to test the error exits for
- the LAPACK routines and driver routines.
- line 14: NEWSD, INTEGER
- A code indicating how to set the random number seed.
- = 0: Set the seed to a default value before each run
- = 1: Initialize the seed to a default value only before the
- first run
- = 2: Like 1, but use the seed values on the next line
- If line 14 was 2:
- line 15: INTEGER array, dimension (4)
- Four integer values for the random number seed.
- lines 15-EOF: Lines specifying matrix types, as for NEP.
- The 3-character path names are 'SVD' or 'SBD' for both the
- SVD routines and the SVD driver routines.
- -----------------------------------------------------------------------
- DEV and DES data files:
- line 1: 'DEV' or 'DES' in columns 1 to 3.
- line 2: NSIZES, INTEGER
- Number of sizes of matrices to use. Should be at least 0
- and at most 20. If NSIZES = 0, no testing is done
- (although the remaining 3 lines are still read).
- line 3: NN, INTEGER array, dimension(NSIZES)
- Dimensions of matrices to be tested.
- line 4: NB, NBMIN, NX, NS, NBCOL, INTEGERs
- These integer parameters determine how blocking is done
- (see ILAENV for details)
- NB : block size
- NBMIN : minimum block size
- NX : minimum dimension for blocking
- NS : number of shifts in xHSEQR
- NBCOL : minimum column dimension for blocking
- line 5: THRESH, REAL
- The test threshold against which computed residuals are
- compared. Should generally be in the range from 10. to 20.
- If it is 0., all test case data will be printed.
- line 6: TSTERR, LOGICAL
- Flag indicating whether or not to test the error exits.
- line 7: NEWSD, INTEGER
- A code indicating how to set the random number seed.
- = 0: Set the seed to a default value before each run
- = 1: Initialize the seed to a default value only before the
- first run
- = 2: Like 1, but use the seed values on the next line
- If line 7 was 2:
- line 8: INTEGER array, dimension (4)
- Four integer values for the random number seed.
- lines 9 and following: Lines specifying matrix types, as for NEP.
- The 3-character path name is 'DEV' to test SGEEV, or
- 'DES' to test SGEES.
- -----------------------------------------------------------------------
- The DVX data has two parts. The first part is identical to DEV,
- and the second part consists of test matrices with precomputed
- solutions.
- line 1: 'DVX' in columns 1-3.
- line 2: NSIZES, INTEGER
- If NSIZES = 0, no testing of randomly generated examples
- is done, but any precomputed examples are tested.
- line 3: NN, INTEGER array, dimension(NSIZES)
- line 4: NB, NBMIN, NX, NS, NBCOL, INTEGERs
- line 5: THRESH, REAL
- line 6: TSTERR, LOGICAL
- line 7: NEWSD, INTEGER
- If line 7 was 2:
- line 8: INTEGER array, dimension (4)
- lines 9 and following: The first line contains 'DVX' in columns 1-3
- followed by the number of matrix types, possibly with
- a second line to specify certain matrix types.
- If the number of matrix types = 0, no testing of randomly
- generated examples is done, but any precomputed examples
- are tested.
- remaining lines : Each matrix is stored on 1+2*N lines, where N is
- its dimension. The first line contains the dimension (a
- single integer). The next N lines contain the matrix, one
- row per line. The last N lines correspond to each
- eigenvalue. Each of these last N lines contains 4 real
- values: the real part of the eigenvalue, the imaginary
- part of the eigenvalue, the reciprocal condition number of
- the eigenvalues, and the reciprocal condition number of the
- eigenvector. The end of data is indicated by dimension N=0.
- Even if no data is to be tested, there must be at least one
- line containing N=0.
- -----------------------------------------------------------------------
- The DSX data is like DVX. The first part is identical to DEV, and the
- second part consists of test matrices with precomputed solutions.
- line 1: 'DSX' in columns 1-3.
- line 2: NSIZES, INTEGER
- If NSIZES = 0, no testing of randomly generated examples
- is done, but any precomputed examples are tested.
- line 3: NN, INTEGER array, dimension(NSIZES)
- line 4: NB, NBMIN, NX, NS, NBCOL, INTEGERs
- line 5: THRESH, REAL
- line 6: TSTERR, LOGICAL
- line 7: NEWSD, INTEGER
- If line 7 was 2:
- line 8: INTEGER array, dimension (4)
- lines 9 and following: The first line contains 'DSX' in columns 1-3
- followed by the number of matrix types, possibly with
- a second line to specify certain matrix types.
- If the number of matrix types = 0, no testing of randomly
- generated examples is done, but any precomputed examples
- are tested.
- remaining lines : Each matrix is stored on 3+N lines, where N is its
- dimension. The first line contains the dimension N and the
- dimension M of an invariant subspace. The second line
- contains M integers, identifying the eigenvalues in the
- invariant subspace (by their position in a list of
- eigenvalues ordered by increasing real part). The next N
- lines contain the matrix. The last line contains the
- reciprocal condition number for the average of the selected
- eigenvalues, and the reciprocal condition number for the
- corresponding right invariant subspace. The end of data is
- indicated by a line containing N=0 and M=0. Even if no data
- is to be tested, there must be at least one line containing
- N=0 and M=0.
- -----------------------------------------------------------------------
- DGG input file:
- line 2: NN, INTEGER
- Number of values of N.
- line 3: NVAL, INTEGER array, dimension (NN)
- The values for the matrix dimension N.
- line 4: NPARMS, INTEGER
- Number of values of the parameters NB, NBMIN, NS, MAXB, and
- NBCOL.
- line 5: NBVAL, INTEGER array, dimension (NPARMS)
- The values for the blocksize NB.
- line 6: NBMIN, INTEGER array, dimension (NPARMS)
- The values for NBMIN, the minimum row dimension for blocks.
- line 7: NSVAL, INTEGER array, dimension (NPARMS)
- The values for the number of shifts.
- line 8: MXBVAL, INTEGER array, dimension (NPARMS)
- The values for MAXB, used in determining minimum blocksize.
- line 9: NBCOL, INTEGER array, dimension (NPARMS)
- The values for NBCOL, the minimum column dimension for
- blocks.
- line 10: THRESH
- Threshold value for the test ratios. Information will be
- printed about each test for which the test ratio is greater
- than or equal to the threshold.
- line 11: TSTCHK, LOGICAL
- Flag indicating whether or not to test the LAPACK routines.
- line 12: TSTDRV, LOGICAL
- Flag indicating whether or not to test the driver routines.
- line 13: TSTERR, LOGICAL
- Flag indicating whether or not to test the error exits for
- the LAPACK routines and driver routines.
- line 14: NEWSD, INTEGER
- A code indicating how to set the random number seed.
- = 0: Set the seed to a default value before each run
- = 1: Initialize the seed to a default value only before the
- first run
- = 2: Like 1, but use the seed values on the next line
- If line 14 was 2:
- line 15: INTEGER array, dimension (4)
- Four integer values for the random number seed.
- lines 15-EOF: Lines specifying matrix types, as for NEP.
- The 3-character path name is 'DGG' for the generalized
- eigenvalue problem routines and driver routines.
- -----------------------------------------------------------------------
- DGS and DGV input files:
- line 1: 'DGS' or 'DGV' in columns 1 to 3.
- line 2: NN, INTEGER
- Number of values of N.
- line 3: NVAL, INTEGER array, dimension(NN)
- Dimensions of matrices to be tested.
- line 4: NB, NBMIN, NX, NS, NBCOL, INTEGERs
- These integer parameters determine how blocking is done
- (see ILAENV for details)
- NB : block size
- NBMIN : minimum block size
- NX : minimum dimension for blocking
- NS : number of shifts in xHGEQR
- NBCOL : minimum column dimension for blocking
- line 5: THRESH, REAL
- The test threshold against which computed residuals are
- compared. Should generally be in the range from 10. to 20.
- If it is 0., all test case data will be printed.
- line 6: TSTERR, LOGICAL
- Flag indicating whether or not to test the error exits.
- line 7: NEWSD, INTEGER
- A code indicating how to set the random number seed.
- = 0: Set the seed to a default value before each run
- = 1: Initialize the seed to a default value only before the
- first run
- = 2: Like 1, but use the seed values on the next line
- If line 17 was 2:
- line 7: INTEGER array, dimension (4)
- Four integer values for the random number seed.
- lines 7-EOF: Lines specifying matrix types, as for NEP.
- The 3-character path name is 'DGS' for the generalized
- eigenvalue problem routines and driver routines.
- -----------------------------------------------------------------------
- DXV input files:
- line 1: 'DXV' in columns 1 to 3.
- line 2: N, INTEGER
- Value of N.
- line 3: NB, NBMIN, NX, NS, NBCOL, INTEGERs
- These integer parameters determine how blocking is done
- (see ILAENV for details)
- NB : block size
- NBMIN : minimum block size
- NX : minimum dimension for blocking
- NS : number of shifts in xHGEQR
- NBCOL : minimum column dimension for blocking
- line 4: THRESH, REAL
- The test threshold against which computed residuals are
- compared. Should generally be in the range from 10. to 20.
- Information will be printed about each test for which the
- test ratio is greater than or equal to the threshold.
- line 5: TSTERR, LOGICAL
- Flag indicating whether or not to test the error exits for
- the LAPACK routines and driver routines.
- line 6: NEWSD, INTEGER
- A code indicating how to set the random number seed.
- = 0: Set the seed to a default value before each run
- = 1: Initialize the seed to a default value only before the
- first run
- = 2: Like 1, but use the seed values on the next line
- If line 6 was 2:
- line 7: INTEGER array, dimension (4)
- Four integer values for the random number seed.
- If line 2 was 0:
- line 7-EOF: Precomputed examples are tested.
- remaining lines : Each example is stored on 3+2*N lines, where N is
- its dimension. The first line contains the dimension (a
- single integer). The next N lines contain the matrix A, one
- row per line. The next N lines contain the matrix B. The
- next line contains the reciprocals of the eigenvalue
- condition numbers. The last line contains the reciprocals of
- the eigenvector condition numbers. The end of data is
- indicated by dimension N=0. Even if no data is to be tested,
- there must be at least one line containing N=0.
- -----------------------------------------------------------------------
- DGX input files:
- line 1: 'DGX' in columns 1 to 3.
- line 2: N, INTEGER
- Value of N.
- line 3: NB, NBMIN, NX, NS, NBCOL, INTEGERs
- These integer parameters determine how blocking is done
- (see ILAENV for details)
- NB : block size
- NBMIN : minimum block size
- NX : minimum dimension for blocking
- NS : number of shifts in xHGEQR
- NBCOL : minimum column dimension for blocking
- line 4: THRESH, REAL
- The test threshold against which computed residuals are
- compared. Should generally be in the range from 10. to 20.
- Information will be printed about each test for which the
- test ratio is greater than or equal to the threshold.
- line 5: TSTERR, LOGICAL
- Flag indicating whether or not to test the error exits for
- the LAPACK routines and driver routines.
- line 6: NEWSD, INTEGER
- A code indicating how to set the random number seed.
- = 0: Set the seed to a default value before each run
- = 1: Initialize the seed to a default value only before the
- first run
- = 2: Like 1, but use the seed values on the ne…
Large files files are truncated, but you can click here to view the full file