/contrib/clapack-3.2.1-CMAKE/TESTING/EIG/dchkee.c
C | 1560 lines | 522 code | 301 blank | 737 comment | 10 complexity | a84eb8bdc9d3367b0e007ad98cd92073 MD5 | raw file
Possible License(s): BSD-3-Clause
Large files files are truncated, but you can click here to view the full file
- /* dchkee.f -- translated by f2c (version 20061008).
- You must link the resulting object file with libf2c:
- on Microsoft Windows system, link with libf2c.lib;
- on Linux or Unix systems, link with .../path/to/libf2c.a -lm
- or, if you install libf2c.a in a standard place, with -lf2c -lm
- -- in that order, at the end of the command line, as in
- cc *.o -lf2c -lm
- Source for libf2c is in /netlib/f2c/libf2c.zip, e.g.,
- http://www.netlib.org/f2c/libf2c.zip
- */
- #include "f2c.h"
- #include "blaswrap.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[32];
- } 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__12 = 12;
- static integer c__11 = 11;
- static integer c__13 = 13;
- static integer c__2 = 2;
- static integer c__14 = 14;
- static integer c__0 = 0;
- static integer c__15 = 15;
- static integer c__16 = 16;
- static integer c__20 = 20;
- static integer c__132 = 132;
- 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 */ int 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 \002,i1,\002.\002,i1,"
- "\002.\002,i1)";
- static char fmt_9984[] = "(/\002 The following parameter values will be "
- "used:\002)";
- static char fmt_9989[] = "(\002 Invalid input value: \002,a,\002=\002,"
- "i6,\002; must be >=\002,i6)";
- static char fmt_9988[] = "(\002 Invalid input value: \002,a,\002=\002,"
- "i6,\002; must be <=\002,i6)";
- static char fmt_9983[] = "(4x,a,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_9961[] = "(//1x,a3,\002: NB =\002,i4,\002, NBMIN =\002,"
- "i4,\002, NX =\002,i4,\002, INMIN=\002,i4,\002, INWIN =\002,i4"
- ",\002, INIBL =\002,i4,\002, ISHFTS =\002,i4,\002, IACC22 =\002,i"
- "4)";
- static char fmt_9980[] = "(\002 *** Error code from \002,a,\002 = \002,i"
- "4)";
- 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, i__2, i__3, i__4;
- doublereal d__1;
- cilist ci__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 */
- doublereal a[243936] /* was [17424][14] */, b[87120] /* was [17424]
- [5] */, c__[160000] /* was [400][400] */, d__[1584] /* was [132][
- 12] */;
- integer i__, k;
- doublereal x[660];
- char c1[1], c3[3];
- integer i1;
- doublereal s1, s2;
- integer ic, nk, nn, vers_patch__, vers_major__, vers_minor__;
- logical dbb, dbk, dgg, dbl, dgk, dgl, dsb, des, dgs, dev, glm, dgv, nep,
- lse, dgx, sep;
- doublereal eps;
- logical gqr, svd, dsx, gsv, dvx, dxv;
- char line[80];
- doublereal taua[132];
- integer info;
- char path[3];
- integer kval[20], lenp, mval[20], nval[20];
- doublereal taub[132];
- integer pval[20], itmp, nrhs;
- doublereal work[87781];
- integer iacc22[20];
- logical fatal;
- integer iseed[4], nbcol[20], inibl[20], nbval[20], nbmin[20];
- char vname[32];
- integer inmin[20], newsd, nsval[20], inwin[20], nxval[20], iwork[89760];
- 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 *);
- 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 *), 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 *), ilaver_(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 *);
- 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 *);
- logical tstdif;
- 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 *);
- logical tstchk;
- integer nparms, ishfts[20];
- extern /* Subroutine */ int derrst_(char *, integer *);
- logical dotype[30], logwrk[132];
- 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 *), xlaenv_(integer *, 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 *);
- doublereal result[500];
- integer maxtyp;
- logical tsterr;
- integer ntypes;
- logical tstdrv;
- /* Fortran I/O blocks */
- static cilist io___29 = { 0, 6, 0, fmt_9987, 0 };
- static cilist io___30 = { 0, 6, 0, fmt_9986, 0 };
- static cilist io___31 = { 0, 6, 0, fmt_9985, 0 };
- static cilist io___32 = { 0, 6, 0, fmt_9979, 0 };
- static cilist io___33 = { 0, 6, 0, fmt_9978, 0 };
- static cilist io___34 = { 0, 6, 0, fmt_9977, 0 };
- static cilist io___35 = { 0, 6, 0, fmt_9976, 0 };
- static cilist io___36 = { 0, 6, 0, fmt_9975, 0 };
- static cilist io___37 = { 0, 6, 0, fmt_9964, 0 };
- static cilist io___38 = { 0, 6, 0, fmt_9965, 0 };
- static cilist io___39 = { 0, 6, 0, fmt_9963, 0 };
- static cilist io___40 = { 0, 6, 0, fmt_9962, 0 };
- static cilist io___41 = { 0, 6, 0, fmt_9974, 0 };
- static cilist io___42 = { 0, 6, 0, fmt_9967, 0 };
- static cilist io___43 = { 0, 6, 0, fmt_9971, 0 };
- static cilist io___44 = { 0, 6, 0, fmt_9970, 0 };
- static cilist io___45 = { 0, 6, 0, fmt_9969, 0 };
- static cilist io___46 = { 0, 6, 0, fmt_9968, 0 };
- static cilist io___47 = { 0, 5, 0, 0, 0 };
- static cilist io___50 = { 0, 6, 0, fmt_9992, 0 };
- static cilist io___54 = { 0, 6, 0, fmt_9972, 0 };
- static cilist io___55 = { 0, 6, 0, fmt_9984, 0 };
- static cilist io___56 = { 0, 5, 0, 0, 0 };
- static cilist io___58 = { 0, 6, 0, fmt_9989, 0 };
- static cilist io___59 = { 0, 6, 0, fmt_9988, 0 };
- static cilist io___60 = { 0, 5, 0, 0, 0 };
- static cilist io___64 = { 0, 6, 0, fmt_9989, 0 };
- static cilist io___65 = { 0, 6, 0, fmt_9988, 0 };
- static cilist io___66 = { 0, 6, 0, fmt_9983, 0 };
- static cilist io___67 = { 0, 5, 0, 0, 0 };
- static cilist io___69 = { 0, 6, 0, fmt_9989, 0 };
- static cilist io___70 = { 0, 6, 0, fmt_9988, 0 };
- static cilist io___71 = { 0, 6, 0, fmt_9983, 0 };
- static cilist io___72 = { 0, 5, 0, 0, 0 };
- static cilist io___74 = { 0, 6, 0, fmt_9989, 0 };
- static cilist io___75 = { 0, 6, 0, fmt_9988, 0 };
- static cilist io___76 = { 0, 6, 0, fmt_9983, 0 };
- static cilist io___77 = { 0, 6, 0, fmt_9983, 0 };
- static cilist io___78 = { 0, 5, 0, 0, 0 };
- static cilist io___80 = { 0, 5, 0, 0, 0 };
- static cilist io___82 = { 0, 6, 0, fmt_9989, 0 };
- static cilist io___83 = { 0, 6, 0, fmt_9988, 0 };
- static cilist io___84 = { 0, 6, 0, fmt_9983, 0 };
- static cilist io___85 = { 0, 5, 0, 0, 0 };
- static cilist io___94 = { 0, 6, 0, fmt_9989, 0 };
- static cilist io___95 = { 0, 6, 0, fmt_9989, 0 };
- static cilist io___96 = { 0, 6, 0, fmt_9989, 0 };
- static cilist io___97 = { 0, 6, 0, fmt_9989, 0 };
- static cilist io___98 = { 0, 6, 0, fmt_9989, 0 };
- static cilist io___99 = { 0, 6, 0, fmt_9989, 0 };
- static cilist io___100 = { 0, 6, 0, fmt_9989, 0 };
- static cilist io___101 = { 0, 6, 0, fmt_9989, 0 };
- static cilist io___102 = { 0, 6, 0, fmt_9983, 0 };
- static cilist io___103 = { 0, 6, 0, fmt_9983, 0 };
- static cilist io___104 = { 0, 6, 0, fmt_9983, 0 };
- static cilist io___105 = { 0, 6, 0, fmt_9983, 0 };
- static cilist io___106 = { 0, 6, 0, fmt_9983, 0 };
- static cilist io___107 = { 0, 6, 0, fmt_9983, 0 };
- static cilist io___108 = { 0, 6, 0, fmt_9983, 0 };
- static cilist io___109 = { 0, 6, 0, fmt_9983, 0 };
- static cilist io___110 = { 0, 5, 0, 0, 0 };
- static cilist io___113 = { 0, 6, 0, fmt_9989, 0 };
- static cilist io___114 = { 0, 6, 0, fmt_9989, 0 };
- static cilist io___115 = { 0, 6, 0, fmt_9989, 0 };
- static cilist io___116 = { 0, 6, 0, fmt_9989, 0 };
- static cilist io___117 = { 0, 6, 0, fmt_9989, 0 };
- static cilist io___118 = { 0, 6, 0, fmt_9983, 0 };
- static cilist io___119 = { 0, 6, 0, fmt_9983, 0 };
- static cilist io___120 = { 0, 6, 0, fmt_9983, 0 };
- static cilist io___121 = { 0, 6, 0, fmt_9983, 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, 5, 0, 0, 0 };
- static cilist io___128 = { 0, 6, 0, fmt_9989, 0 };
- static cilist io___129 = { 0, 6, 0, fmt_9988, 0 };
- static cilist io___130 = { 0, 6, 0, fmt_9983, 0 };
- static cilist io___131 = { 0, 5, 0, 0, 0 };
- static cilist io___132 = { 0, 6, 0, fmt_9989, 0 };
- static cilist io___133 = { 0, 6, 0, fmt_9988, 0 };
- static cilist io___134 = { 0, 6, 0, fmt_9983, 0 };
- static cilist io___135 = { 0, 5, 0, 0, 0 };
- static cilist io___136 = { 0, 6, 0, fmt_9989, 0 };
- static cilist io___137 = { 0, 6, 0, fmt_9988, 0 };
- static cilist io___138 = { 0, 6, 0, fmt_9983, 0 };
- static cilist io___139 = { 0, 5, 0, 0, 0 };
- static cilist io___140 = { 0, 6, 0, fmt_9989, 0 };
- static cilist io___141 = { 0, 6, 0, fmt_9988, 0 };
- static cilist io___142 = { 0, 6, 0, fmt_9983, 0 };
- static cilist io___143 = { 0, 5, 0, 0, 0 };
- static cilist io___144 = { 0, 6, 0, fmt_9989, 0 };
- static cilist io___145 = { 0, 6, 0, fmt_9988, 0 };
- static cilist io___146 = { 0, 6, 0, fmt_9983, 0 };
- static cilist io___147 = { 0, 5, 0, 0, 0 };
- static cilist io___148 = { 0, 6, 0, fmt_9989, 0 };
- static cilist io___149 = { 0, 6, 0, fmt_9983, 0 };
- static cilist io___150 = { 0, 5, 0, 0, 0 };
- static cilist io___151 = { 0, 6, 0, fmt_9989, 0 };
- static cilist io___152 = { 0, 6, 0, fmt_9983, 0 };
- static cilist io___153 = { 0, 5, 0, 0, 0 };
- static cilist io___154 = { 0, 6, 0, fmt_9989, 0 };
- static cilist io___155 = { 0, 6, 0, fmt_9983, 0 };
- static cilist io___156 = { 0, 5, 0, 0, 0 };
- static cilist io___157 = { 0, 6, 0, fmt_9989, 0 };
- static cilist io___158 = { 0, 6, 0, fmt_9983, 0 };
- static cilist io___159 = { 0, 5, 0, 0, 0 };
- static cilist io___160 = { 0, 6, 0, fmt_9989, 0 };
- static cilist io___161 = { 0, 6, 0, fmt_9983, 0 };
- static cilist io___162 = { 0, 5, 0, 0, 0 };
- static cilist io___164 = { 0, 6, 0, fmt_9989, 0 };
- static cilist io___165 = { 0, 6, 0, fmt_9988, 0 };
- static cilist io___166 = { 0, 6, 0, fmt_9983, 0 };
- static cilist io___167 = { 0, 6, 0, 0, 0 };
- static cilist io___169 = { 0, 6, 0, fmt_9981, 0 };
- static cilist io___170 = { 0, 6, 0, fmt_9981, 0 };
- static cilist io___171 = { 0, 6, 0, fmt_9981, 0 };
- static cilist io___172 = { 0, 5, 0, 0, 0 };
- static cilist io___173 = { 0, 6, 0, fmt_9982, 0 };
- static cilist io___174 = { 0, 5, 0, 0, 0 };
- static cilist io___176 = { 0, 5, 0, 0, 0 };
- static cilist io___178 = { 0, 5, 0, 0, 0 };
- static cilist io___179 = { 0, 5, 0, 0, 0 };
- static cilist io___181 = { 0, 5, 0, 0, 0 };
- static cilist io___183 = { 0, 6, 0, fmt_9999, 0 };
- static cilist io___192 = { 0, 6, 0, fmt_9991, 0 };
- static cilist io___193 = { 0, 6, 0, fmt_9990, 0 };
- static cilist io___196 = { 0, 6, 0, fmt_9961, 0 };
- static cilist io___204 = { 0, 6, 0, fmt_9980, 0 };
- static cilist io___205 = { 0, 6, 0, fmt_9997, 0 };
- static cilist io___206 = { 0, 6, 0, fmt_9980, 0 };
- static cilist io___207 = { 0, 6, 0, fmt_9980, 0 };
- static cilist io___208 = { 0, 6, 0, fmt_9997, 0 };
- static cilist io___209 = { 0, 6, 0, fmt_9980, 0 };
- static cilist io___211 = { 0, 6, 0, fmt_9995, 0 };
- static cilist io___212 = { 0, 6, 0, fmt_9980, 0 };
- static cilist io___213 = { 0, 6, 0, fmt_9990, 0 };
- static cilist io___214 = { 0, 6, 0, fmt_9980, 0 };
- static cilist io___215 = { 0, 6, 0, fmt_9973, 0 };
- static cilist io___216 = { 0, 6, 0, fmt_9990, 0 };
- static cilist io___217 = { 0, 6, 0, fmt_9980, 0 };
- static cilist io___218 = { 0, 6, 0, fmt_9973, 0 };
- static cilist io___219 = { 0, 6, 0, fmt_9990, 0 };
- static cilist io___220 = { 0, 6, 0, fmt_9980, 0 };
- static cilist io___221 = { 0, 6, 0, fmt_9973, 0 };
- static cilist io___222 = { 0, 6, 0, fmt_9990, 0 };
- static cilist io___223 = { 0, 6, 0, fmt_9980, 0 };
- static cilist io___224 = { 0, 6, 0, fmt_9973, 0 };
- static cilist io___225 = { 0, 6, 0, fmt_9996, 0 };
- static cilist io___228 = { 0, 6, 0, fmt_9980, 0 };
- static cilist io___229 = { 0, 6, 0, fmt_9980, 0 };
- static cilist io___230 = { 0, 6, 0, fmt_9990, 0 };
- static cilist io___231 = { 0, 6, 0, fmt_9980, 0 };
- static cilist io___232 = { 0, 6, 0, fmt_9973, 0 };
- static cilist io___233 = { 0, 6, 0, fmt_9990, 0 };
- static cilist io___235 = { 0, 6, 0, fmt_9980, 0 };
- static cilist io___236 = { 0, 6, 0, fmt_9973, 0 };
- static cilist io___237 = { 0, 6, 0, fmt_9990, 0 };
- static cilist io___238 = { 0, 6, 0, fmt_9980, 0 };
- static cilist io___239 = { 0, 6, 0, fmt_9973, 0 };
- static cilist io___240 = { 0, 6, 0, fmt_9990, 0 };
- static cilist io___241 = { 0, 6, 0, fmt_9980, 0 };
- static cilist io___242 = { 0, 6, 0, fmt_9973, 0 };
- static cilist io___243 = { 0, 6, 0, fmt_9980, 0 };
- static cilist io___244 = { 0, 6, 0, fmt_9966, 0 };
- static cilist io___245 = { 0, 6, 0, fmt_9980, 0 };
- static cilist io___248 = { 0, 6, 0, fmt_9980, 0 };
- static cilist io___251 = { 0, 6, 0, fmt_9980, 0 };
- static cilist io___252 = { 0, 6, 0, fmt_9980, 0 };
- static cilist io___253 = { 0, 6, 0, fmt_9980, 0 };
- static cilist io___254 = { 0, 6, 0, 0, 0 };
- static cilist io___255 = { 0, 6, 0, 0, 0 };
- static cilist io___256 = { 0, 6, 0, fmt_9992, 0 };
- static cilist io___257 = { 0, 6, 0, fmt_9994, 0 };
- static cilist io___259 = { 0, 6, 0, fmt_9993, 0 };
- /* -- LAPACK test routine (version 3.1.1) -- */
- /* Univ. of Tennessee, Univ. of California Berkeley and NAG Ltd.. */
- /* January 2007 */
- /* 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: INMIN, INTEGER array, dimension (NPARMS) */
- /* LAHQR vs TTQRE crossover point, >= 11 */
- /* line 9: INWIN, INTEGER array, dimension (NPARMS) */
- /* recommended deflation window size */
- /* line 10: INIBL, INTEGER array, dimension (NPARMS) */
- /* nibble crossover point */
- /* line 11: ISHFTS, INTEGER array, dimension (NPARMS) */
- /* number of simultaneous shifts) */
- /* line 12: IACC22, INTEGER array, dimension (NPARMS) */
- /* select structured matrix multiply: 0, 1 or 2) */
- /* line 13: 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 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: 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…
Large files files are truncated, but you can click here to view the full file