PageRenderTime 46ms CodeModel.GetById 16ms RepoModel.GetById 1ms app.codeStats 0ms

/contrib/clapack-3.2.1-CMAKE/TESTING/EIG/dchkee.c

http://github.com/openmeeg/openmeeg
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

    

  1. /* dchkee.f -- translated by f2c (version 20061008).
    2. 

  2.    You must link the resulting object file with libf2c:
    3. 

  3. 	on Microsoft Windows system, link with libf2c.lib;
    4. 

  4. 	on Linux or Unix systems, link with .../path/to/libf2c.a -lm
    5. 

  5. 	or, if you install libf2c.a in a standard place, with -lf2c -lm
    6. 

  6. 	-- in that order, at the end of the command line, as in
    7. 

  7. 		cc *.o -lf2c -lm
    8. 

  8. 	Source for libf2c is in /netlib/f2c/libf2c.zip, e.g.,
    9. 

    10. 

  10. 		http://www.netlib.org/f2c/libf2c.zip
    11. 

  11. */
    12. 

  12. 

  13. #include "f2c.h"
    14. 

  14. #include "blaswrap.h"
    15. 

  15. 

  16. /* Common Block Declarations */
    17. 

  17. 

  18. struct {
    19. 

  19.     integer nproc, nshift, maxb;
    20. 

  20. } cenvir_;
    21. 

  21. 

  22. #define cenvir_1 cenvir_
    23. 

  23. 

  24. struct {
    25. 

  25.     integer infot, nunit;
    26. 

  26.     logical ok, lerr;
    27. 

  27. } infoc_;
    28. 

  28. 

  29. #define infoc_1 infoc_
    30. 

  30. 

  31. struct {
    32. 

  32.     char srnamt[32];
    33. 

  33. } srnamc_;
    34. 

  34. 

  35. #define srnamc_1 srnamc_
    36. 

  36. 

  37. struct {
    38. 

  38.     integer selopt, seldim;
    39. 

  39.     logical selval[20];
    40. 

  40.     doublereal selwr[20], selwi[20];
    41. 

  41. } sslct_;
    42. 

  42. 

  43. #define sslct_1 sslct_
    44. 

  44. 

  45. struct {
    46. 

  46.     integer iparms[100];
    47. 

  47. } zlaenv_;
    48. 

  48. 

  49. #define zlaenv_1 zlaenv_
    50. 

  50. 

  51. /* Table of constant values */
    52. 

  52. 

  53. static integer c__1 = 1;
    54. 

  54. static integer c__3 = 3;
    55. 

  55. static integer c__5 = 5;
    56. 

  56. static integer c__6 = 6;
    57. 

  57. static integer c__12 = 12;
    58. 

  58. static integer c__11 = 11;
    59. 

  59. static integer c__13 = 13;
    60. 

  60. static integer c__2 = 2;
    61. 

  61. static integer c__14 = 14;
    62. 

  62. static integer c__0 = 0;
    63. 

  63. static integer c__15 = 15;
    64. 

  64. static integer c__16 = 16;
    65. 

  65. static integer c__20 = 20;
    66. 

  66. static integer c__132 = 132;
    67. 

  67. static integer c__4 = 4;
    68. 

  68. static integer c__8 = 8;
    69. 

  69. static integer c__87781 = 87781;
    70. 

  70. static integer c__9 = 9;
    71. 

  71. static integer c__25 = 25;
    72. 

  72. static integer c__89760 = 89760;
    73. 

  73. static integer c__18 = 18;
    74. 

  74. static integer c__400 = 400;
    75. 

  75. static integer c__89758 = 89758;
    76. 

  76. static integer c__264 = 264;
    77. 

  77. 

  78. /* Main program */ int MAIN__(void)
    79. 

  79. {
    80. 

  80.     /* Initialized data */
    81. 

  81. 

  82.     static char intstr[10] = "0123456789";
    83. 

  83.     static integer ioldsd[4] = { 0,0,0,1 };
    84. 

  84. 

  85.     /* Format strings */
    86. 

  86.     static char fmt_9987[] = "(\002 Tests of the Nonsymmetric Eigenvalue Pro"
    87. 

  87. 	    "blem routines\002)";
    88. 

  88.     static char fmt_9986[] = "(\002 Tests of the Symmetric Eigenvalue Proble"
    89. 

  89. 	    "m routines\002)";
    90. 

  90.     static char fmt_9985[] = "(\002 Tests of the Singular Value Decompositio"
    91. 

  91. 	    "n routines\002)";
    92. 

  92.     static char fmt_9979[] = "(/\002 Tests of the Nonsymmetric Eigenvalue Pr"
    93. 

  93. 	    "oblem Driver\002,/\002    DGEEV (eigenvalues and eigevectors)"
    94. 

  94. 	    "\002)";
    95. 

  95.     static char fmt_9978[] = "(/\002 Tests of the Nonsymmetric Eigenvalue Pr"
    96. 

  96. 	    "oblem Driver\002,/\002    DGEES (Schur form)\002)";
    97. 

  97.     static char fmt_9977[] = "(/\002 Tests of the Nonsymmetric Eigenvalue Pr"
    98. 

  98. 	    "oblem Expert\002,\002 Driver\002,/\002    DGEEVX (eigenvalues, e"
    99. 

  99. 	    "igenvectors and\002,\002 condition numbers)\002)";
    100. 

  100.     static char fmt_9976[] = "(/\002 Tests of the Nonsymmetric Eigenvalue Pr"
    101. 

  101. 	    "oblem Expert\002,\002 Driver\002,/\002    DGEESX (Schur form and"
    102. 

  102. 	    " condition\002,\002 numbers)\002)";
    103. 

  103.     static char fmt_9975[] = "(/\002 Tests of the Generalized Nonsymmetric E"
    104. 

  104. 	    "igenvalue \002,\002Problem routines\002)";
    105. 

  105.     static char fmt_9964[] = "(/\002 Tests of the Generalized Nonsymmetric E"
    106. 

  106. 	    "igenvalue \002,\002Problem Driver DGGES\002)";
    107. 

  107.     static char fmt_9965[] = "(/\002 Tests of the Generalized Nonsymmetric E"
    108. 

  108. 	    "igenvalue \002,\002Problem Expert Driver DGGESX\002)";
    109. 

  109.     static char fmt_9963[] = "(/\002 Tests of the Generalized Nonsymmetric E"
    110. 

  110. 	    "igenvalue \002,\002Problem Driver DGGEV\002)";
    111. 

  111.     static char fmt_9962[] = "(/\002 Tests of the Generalized Nonsymmetric E"
    112. 

  112. 	    "igenvalue \002,\002Problem Expert Driver DGGEVX\002)";
    113. 

  113.     static char fmt_9974[] = "(\002 Tests of DSBTRD\002,/\002 (reduction of "
    114. 

  114. 	    "a symmetric band \002,\002matrix to tridiagonal form)\002)";
    115. 

  115.     static char fmt_9967[] = "(\002 Tests of DGBBRD\002,/\002 (reduction of "
    116. 

  116. 	    "a general band \002,\002matrix to real bidiagonal form)\002)";
    117. 

  117.     static char fmt_9971[] = "(/\002 Tests of the Generalized Linear Regress"
    118. 

  118. 	    "ion Model \002,\002routines\002)";
    119. 

  119.     static char fmt_9970[] = "(/\002 Tests of the Generalized QR and RQ rout"
    120. 

  120. 	    "ines\002)";
    121. 

  121.     static char fmt_9969[] = "(/\002 Tests of the Generalized Singular Valu"
    122. 

  122. 	    "e\002,\002 Decomposition routines\002)";
    123. 

  123.     static char fmt_9968[] = "(/\002 Tests of the Linear Least Squares routi"
    124. 

  124. 	    "nes\002)";
    125. 

  125.     static char fmt_9992[] = "(1x,a3,\002:  Unrecognized path name\002)";
    126. 

  126.     static char fmt_9972[] = "(/\002 LAPACK VERSION \002,i1,\002.\002,i1,"
    127. 

  127. 	    "\002.\002,i1)";
    128. 

  128.     static char fmt_9984[] = "(/\002 The following parameter values will be "
    129. 

  129. 	    "used:\002)";
    130. 

  130.     static char fmt_9989[] = "(\002 Invalid input value: \002,a,\002=\002,"
    131. 

  131. 	    "i6,\002; must be >=\002,i6)";
    132. 

  132.     static char fmt_9988[] = "(\002 Invalid input value: \002,a,\002=\002,"
    133. 

  133. 	    "i6,\002; must be <=\002,i6)";
    134. 

  134.     static char fmt_9983[] = "(4x,a,10i6,/10x,10i6)";
    135. 

  135.     static char fmt_9981[] = "(\002 Relative machine \002,a,\002 is taken to"
    136. 

  136. 	    " be\002,d16.6)";
    137. 

  137.     static char fmt_9982[] = "(/\002 Routines pass computational tests if te"
    138. 

  138. 	    "st ratio is \002,\002less than\002,f8.2,/)";
    139. 

  139.     static char fmt_9999[] = "(/\002 Execution not attempted due to input er"
    140. 

  140. 	    "rors\002)";
    141. 

  141.     static char fmt_9991[] = "(//\002 *** Invalid integer value in column"
    142. 

  142. 	    " \002,i2,\002 of input\002,\002 line:\002,/a79)";
    143. 

  143.     static char fmt_9990[] = "(//1x,a3,\002 routines were not tested\002)";
    144. 

  144.     static char fmt_9961[] = "(//1x,a3,\002:  NB =\002,i4,\002, NBMIN =\002,"
    145. 

  145. 	    "i4,\002, NX =\002,i4,\002, INMIN=\002,i4,\002, INWIN =\002,i4"
    146. 

  146. 	    ",\002, INIBL =\002,i4,\002, ISHFTS =\002,i4,\002, IACC22 =\002,i"
    147. 

  147. 	    "4)";
    148. 

  148.     static char fmt_9980[] = "(\002 *** Error code from \002,a,\002 = \002,i"
    149. 

  149. 	    "4)";
    150. 

  150.     static char fmt_9997[] = "(//1x,a3,\002:  NB =\002,i4,\002, NBMIN =\002,"
    151. 

  151. 	    "i4,\002, NX =\002,i4)";
    152. 

  152.     static char fmt_9995[] = "(//1x,a3,\002:  NB =\002,i4,\002, NBMIN =\002,"
    153. 

  153. 	    "i4,\002, NX =\002,i4,\002, NRHS =\002,i4)";
    154. 

  154.     static char fmt_9973[] = "(/1x,71(\002-\002))";
    155. 

  155.     static char fmt_9996[] = "(//1x,a3,\002:  NB =\002,i4,\002, NBMIN =\002,"
    156. 

  156. 	    "i4,\002, NS =\002,i4,\002, MAXB =\002,i4,\002, NBCOL =\002,i4)";
    157. 

  157.     static char fmt_9966[] = "(//1x,a3,\002:  NRHS =\002,i4)";
    158. 

  158.     static char fmt_9994[] = "(//\002 End of tests\002)";
    159. 

  159.     static char fmt_9993[] = "(\002 Total time used = \002,f12.2,\002 seco"
    160. 

  160. 	    "nds\002,/)";
    161. 

  161. 

  162.     /* System generated locals */
    163. 

  163.     integer i__1, i__2, i__3, i__4;
    164. 

  164.     doublereal d__1;
    165. 

  165.     cilist ci__1;
    166. 

  166. 

  167.     /* Builtin functions */
    168. 

  168.     integer s_rsfe(cilist *), do_fio(integer *, char *, ftnlen), e_rsfe(void);
    169. 

  169.     /* Subroutine */ int s_copy(char *, char *, ftnlen, ftnlen);
    170. 

  170.     integer s_cmp(char *, char *, ftnlen, ftnlen), s_wsfe(cilist *), e_wsfe(
    171. 

  171. 	    void), s_rsle(cilist *), do_lio(integer *, integer *, char *, 
    172. 

  172. 	    ftnlen), e_rsle(void), s_wsle(cilist *), e_wsle(void);
    173. 

  173.     /* Subroutine */ int s_stop(char *, ftnlen);
    174. 

  174.     integer i_len(char *, ftnlen);
    175. 

  175. 

  176.     /* Local variables */
    177. 

  177.     doublereal a[243936]	/* was [17424][14] */, b[87120]	/* was [17424]
    178. 

  178. 	    [5] */, c__[160000]	/* was [400][400] */, d__[1584]	/* was [132][
    179. 

  179. 	    12] */;
    180. 

  180.     integer i__, k;
    181. 

  181.     doublereal x[660];
    182. 

  182.     char c1[1], c3[3];
    183. 

  183.     integer i1;
    184. 

  184.     doublereal s1, s2;
    185. 

  185.     integer ic, nk, nn, vers_patch__, vers_major__, vers_minor__;
    186. 

  186.     logical dbb, dbk, dgg, dbl, dgk, dgl, dsb, des, dgs, dev, glm, dgv, nep, 
    187. 

  187. 	    lse, dgx, sep;
    188. 

  188.     doublereal eps;
    189. 

  189.     logical gqr, svd, dsx, gsv, dvx, dxv;
    190. 

  190.     char line[80];
    191. 

  191.     doublereal taua[132];
    192. 

  192.     integer info;
    193. 

  193.     char path[3];
    194. 

  194.     integer kval[20], lenp, mval[20], nval[20];
    195. 

  195.     doublereal taub[132];
    196. 

  196.     integer pval[20], itmp, nrhs;
    197. 

  197.     doublereal work[87781];
    198. 

  198.     integer iacc22[20];
    199. 

  199.     logical fatal;
    200. 

  200.     integer iseed[4], nbcol[20], inibl[20], nbval[20], nbmin[20];
    201. 

  201.     char vname[32];
    202. 

  202.     integer inmin[20], newsd, nsval[20], inwin[20], nxval[20], iwork[89760];
    203. 

  203.     extern /* Subroutine */ int dchkbb_(integer *, integer *, integer *, 
    204. 

  204. 	    integer *, integer *, integer *, logical *, integer *, integer *, 
    205. 

  205. 	    doublereal *, integer *, doublereal *, integer *, doublereal *, 
    206. 

  206. 	    integer *, doublereal *, doublereal *, doublereal *, integer *, 
    207. 

  207. 	    doublereal *, integer *, doublereal *, integer *, doublereal *, 
    208. 

  208. 	    doublereal *, integer *, doublereal *, integer *), dchkbd_(
    209. 

  209. 	    integer *, integer *, integer *, integer *, logical *, integer *, 
    210. 

  210. 	    integer *, doublereal *, doublereal *, integer *, doublereal *, 
    211. 

  211. 	    doublereal *, doublereal *, doublereal *, doublereal *, integer *, 
    212. 

  212. 	     doublereal *, doublereal *, doublereal *, integer *, doublereal *
    213. 

  213. , integer *, doublereal *, doublereal *, doublereal *, integer *, 
    214. 

  214. 	    integer *, integer *, integer *), dchkec_(doublereal *, logical *, 
    215. 

  215. 	     integer *, integer *), dchkbk_(integer *, integer *), dchkbl_(
    216. 

  216. 	    integer *, integer *);
    217. 

  217.     extern doublereal dlamch_(char *);
    218. 

  218.     extern /* Subroutine */ int dchkgg_(integer *, integer *, integer *, 
    219. 

  219. 	    logical *, integer *, doublereal *, logical *, doublereal *, 
    220. 

  220. 	    integer *, doublereal *, integer *, doublereal *, doublereal *, 
    221. 

  221. 	    doublereal *, doublereal *, doublereal *, doublereal *, 
    222. 

  222. 	    doublereal *, doublereal *, integer *, doublereal *, doublereal *, 
    223. 

  223. 	     doublereal *, doublereal *, doublereal *, doublereal *, 
    224. 

  224. 	    doublereal *, doublereal *, doublereal *, doublereal *, 
    225. 

  225. 	    doublereal *, doublereal *, integer *, logical *, doublereal *, 
    226. 

  226. 	    integer *), dchkgk_(integer *, integer *), dchkgl_(integer *, 
    227. 

  227. 	    integer *), dchksb_(integer *, integer *, integer *, integer *, 
    228. 

  228. 	    integer *, logical *, integer *, doublereal *, integer *, 
    229. 

  229. 	    doublereal *, integer *, doublereal *, doublereal *, doublereal *, 
    230. 

  230. 	     integer *, doublereal *, integer *, doublereal *, integer *);
    231. 

  231.     extern doublereal dsecnd_(void);
    232. 

  232.     extern /* Subroutine */ int dckglm_(integer *, integer *, integer *, 
    233. 

  233. 	    integer *, integer *, integer *, doublereal *, integer *, 
    234. 

  234. 	    doublereal *, doublereal *, doublereal *, doublereal *, 
    235. 

  235. 	    doublereal *, doublereal *, doublereal *, integer *, integer *, 
    236. 

  236. 	    integer *), derrbd_(char *, integer *), dchkhs_(integer *, 
    237. 

  237. 	     integer *, integer *, logical *, integer *, doublereal *, 
    238. 

  238. 	    integer *, doublereal *, integer *, doublereal *, doublereal *, 
    239. 

  239. 	    doublereal *, doublereal *, integer *, doublereal *, doublereal *, 
    240. 

  240. 	     doublereal *, doublereal *, doublereal *, doublereal *, 
    241. 

  241. 	    doublereal *, doublereal *, doublereal *, doublereal *, 
    242. 

  242. 	    doublereal *, doublereal *, doublereal *, integer *, integer *, 
    243. 

  243. 	    logical *, doublereal *, integer *), alareq_(char *, integer *, 
    244. 

  244. 	    logical *, integer *, integer *, integer *), dcklse_(
    245. 

  245. 	    integer *, integer *, integer *, integer *, integer *, integer *, 
    246. 

  246. 	    doublereal *, integer *, doublereal *, doublereal *, doublereal *, 
    247. 

  247. 	     doublereal *, doublereal *, doublereal *, doublereal *, integer *
    248. 

  248. , integer *, integer *), ddrvbd_(integer *, integer *, integer *, 
    249. 

  249. 	    integer *, logical *, integer *, doublereal *, doublereal *, 
    250. 

  250. 	    integer *, doublereal *, integer *, doublereal *, integer *, 
    251. 

  251. 	    doublereal *, doublereal *, doublereal *, doublereal *, 
    252. 

  252. 	    doublereal *, doublereal *, doublereal *, integer *, integer *, 
    253. 

  253. 	    integer *, integer *), ddrges_(integer *, integer *, integer *, 
    254. 

  254. 	    logical *, integer *, doublereal *, integer *, doublereal *, 
    255. 

  255. 	    integer *, doublereal *, doublereal *, doublereal *, doublereal *, 
    256. 

  256. 	     integer *, doublereal *, doublereal *, doublereal *, doublereal *
    257. 

  257. , doublereal *, integer *, doublereal *, logical *, integer *), 
    258. 

  258. 	    derred_(char *, integer *), derrgg_(char *, integer *), dckgqr_(integer *, integer *, integer *, integer *, 
    259. 

  259. 	    integer *, integer *, integer *, integer *, doublereal *, integer 
    260. 

  260. 	    *, doublereal *, doublereal *, doublereal *, doublereal *, 
    261. 

  261. 	    doublereal *, doublereal *, doublereal *, doublereal *, 
    262. 

  262. 	    doublereal *, doublereal *, doublereal *, doublereal *, 
    263. 

  263. 	    doublereal *, integer *, integer *, integer *), ddrgev_(integer *, 
    264. 

  264. 	     integer *, integer *, logical *, integer *, doublereal *, 
    265. 

  265. 	    integer *, doublereal *, integer *, doublereal *, doublereal *, 
    266. 

  266. 	    doublereal *, doublereal *, integer *, doublereal *, doublereal *, 
    267. 

  267. 	     integer *, doublereal *, doublereal *, doublereal *, doublereal *
    268. 

  268. , doublereal *, doublereal *, doublereal *, integer *, doublereal 
    269. 

  269. 	    *, integer *), ddrvgg_(integer *, integer *, integer *, logical *, 
    270. 

  270. 	     integer *, doublereal *, doublereal *, integer *, doublereal *, 
    271. 

  271. 	    integer *, doublereal *, doublereal *, doublereal *, doublereal *, 
    272. 

  272. 	     doublereal *, doublereal *, integer *, doublereal *, doublereal *
    273. 

  273. , doublereal *, doublereal *, doublereal *, doublereal *, 
    274. 

  274. 	    doublereal *, doublereal *, doublereal *, doublereal *, integer *, 
    275. 

  275. 	     doublereal *, integer *);
    276. 

  276.     extern logical lsamen_(integer *, char *, char *);
    277. 

  277.     extern /* Subroutine */ int dchkst_(integer *, integer *, integer *, 
    278. 

  278. 	    logical *, integer *, doublereal *, integer *, doublereal *, 
    279. 

  279. 	    integer *, doublereal *, doublereal *, doublereal *, doublereal *, 
    280. 

  280. 	     doublereal *, doublereal *, doublereal *, doublereal *, 
    281. 

  281. 	    doublereal *, doublereal *, doublereal *, doublereal *, 
    282. 

  282. 	    doublereal *, integer *, doublereal *, doublereal *, doublereal *, 
    283. 

  283. 	     doublereal *, doublereal *, integer *, integer *, integer *, 
    284. 

  284. 	    doublereal *, integer *), dckgsv_(integer *, integer *, integer *, 
    285. 

  285. 	     integer *, integer *, integer *, doublereal *, integer *, 
    286. 

  286. 	    doublereal *, doublereal *, doublereal *, doublereal *, 
    287. 

  287. 	    doublereal *, doublereal *, doublereal *, doublereal *, 
    288. 

  288. 	    doublereal *, doublereal *, integer *, doublereal *, doublereal *, 
    289. 

  289. 	     integer *, integer *, integer *), ilaver_(integer *, integer *, 
    290. 

  290. 	    integer *), ddrves_(integer *, integer *, integer *, logical *, 
    291. 

  291. 	    integer *, doublereal *, integer *, doublereal *, integer *, 
    292. 

  292. 	    doublereal *, doublereal *, doublereal *, doublereal *, 
    293. 

  293. 	    doublereal *, doublereal *, doublereal *, integer *, doublereal *, 
    294. 

  294. 	     doublereal *, integer *, integer *, logical *, integer *), 
    295. 

  295. 	    derrhs_(char *, integer *);
    296. 

  296.     integer mxbval[20];
    297. 

  297.     extern /* Subroutine */ int ddrvev_(integer *, integer *, integer *, 
    298. 

  298. 	    logical *, integer *, doublereal *, integer *, doublereal *, 
    299. 

  299. 	    integer *, doublereal *, doublereal *, doublereal *, doublereal *, 
    300. 

  300. 	     doublereal *, doublereal *, integer *, doublereal *, integer *, 
    301. 

  301. 	    doublereal *, integer *, doublereal *, doublereal *, integer *, 
    302. 

  302. 	    integer *, integer *), ddrgsx_(integer *, integer *, doublereal *, 
    303. 

  303. 	     integer *, integer *, doublereal *, integer *, doublereal *, 
    304. 

  304. 	    doublereal *, doublereal *, doublereal *, doublereal *, 
    305. 

  305. 	    doublereal *, doublereal *, doublereal *, doublereal *, integer *, 
    306. 

  306. 	     doublereal *, doublereal *, integer *, integer *, integer *, 
    307. 

  307. 	    logical *, integer *), ddrvsg_(integer *, integer *, integer *, 
    308. 

  308. 	    logical *, integer *, doublereal *, integer *, doublereal *, 
    309. 

  309. 	    integer *, doublereal *, integer *, doublereal *, doublereal *, 
    310. 

  310. 	    integer *, doublereal *, doublereal *, doublereal *, doublereal *, 
    311. 

  311. 	     doublereal *, integer *, integer *, integer *, doublereal *, 
    312. 

  312. 	    integer *);
    313. 

  313.     logical tstdif;
    314. 

  314.     doublereal thresh;
    315. 

  315.     extern /* Subroutine */ int ddrgvx_(integer *, doublereal *, integer *, 
    316. 

  316. 	    integer *, doublereal *, integer *, doublereal *, doublereal *, 
    317. 

  317. 	    doublereal *, doublereal *, doublereal *, doublereal *, 
    318. 

  318. 	    doublereal *, doublereal *, integer *, integer *, doublereal *, 
    319. 

  319. 	    doublereal *, doublereal *, doublereal *, doublereal *, 
    320. 

  320. 	    doublereal *, doublereal *, integer *, integer *, integer *, 
    321. 

  321. 	    doublereal *, logical *, integer *);
    322. 

  322.     logical tstchk;
    323. 

  323.     integer nparms, ishfts[20];
    324. 

  324.     extern /* Subroutine */ int derrst_(char *, integer *);
    325. 

  325.     logical dotype[30], logwrk[132];
    326. 

  326.     doublereal thrshn;
    327. 

  327.     extern /* Subroutine */ int ddrvst_(integer *, integer *, integer *, 
    328. 

  328. 	    logical *, integer *, doublereal *, integer *, doublereal *, 
    329. 

  329. 	    integer *, doublereal *, doublereal *, doublereal *, doublereal *, 
    330. 

  330. 	     doublereal *, doublereal *, doublereal *, doublereal *, 
    331. 

  331. 	    doublereal *, integer *, doublereal *, doublereal *, doublereal *, 
    332. 

  332. 	     doublereal *, integer *, integer *, integer *, doublereal *, 
    333. 

  333. 	    integer *), xlaenv_(integer *, integer *), ddrvsx_(integer *, 
    334. 

  334. 	    integer *, integer *, logical *, integer *, doublereal *, integer 
    335. 

  335. 	    *, integer *, doublereal *, integer *, doublereal *, doublereal *, 
    336. 

  336. 	     doublereal *, doublereal *, doublereal *, doublereal *, 
    337. 

  337. 	    doublereal *, doublereal *, doublereal *, integer *, doublereal *, 
    338. 

  338. 	     doublereal *, doublereal *, integer *, integer *, logical *, 
    339. 

  339. 	    integer *), ddrvvx_(integer *, integer *, integer *, logical *, 
    340. 

  340. 	    integer *, doublereal *, integer *, integer *, doublereal *, 
    341. 

  341. 	    integer *, doublereal *, doublereal *, doublereal *, doublereal *, 
    342. 

  342. 	     doublereal *, doublereal *, integer *, doublereal *, integer *, 
    343. 

  343. 	    doublereal *, integer *, doublereal *, doublereal *, doublereal *, 
    344. 

  344. 	     doublereal *, doublereal *, doublereal *, doublereal *, 
    345. 

  345. 	    doublereal *, doublereal *, doublereal *, integer *, integer *, 
    346. 

  346. 	    integer *);
    347. 

  347.     doublereal result[500];
    348. 

  348.     integer maxtyp;
    349. 

  349.     logical tsterr;
    350. 

  350.     integer ntypes;
    351. 

  351.     logical tstdrv;
    352. 

  352. 

  353.     /* Fortran I/O blocks */
    354. 

  354.     static cilist io___29 = { 0, 6, 0, fmt_9987, 0 };
    355. 

  355.     static cilist io___30 = { 0, 6, 0, fmt_9986, 0 };
    356. 

  356.     static cilist io___31 = { 0, 6, 0, fmt_9985, 0 };
    357. 

  357.     static cilist io___32 = { 0, 6, 0, fmt_9979, 0 };
    358. 

  358.     static cilist io___33 = { 0, 6, 0, fmt_9978, 0 };
    359. 

  359.     static cilist io___34 = { 0, 6, 0, fmt_9977, 0 };
    360. 

  360.     static cilist io___35 = { 0, 6, 0, fmt_9976, 0 };
    361. 

  361.     static cilist io___36 = { 0, 6, 0, fmt_9975, 0 };
    362. 

  362.     static cilist io___37 = { 0, 6, 0, fmt_9964, 0 };
    363. 

  363.     static cilist io___38 = { 0, 6, 0, fmt_9965, 0 };
    364. 

  364.     static cilist io___39 = { 0, 6, 0, fmt_9963, 0 };
    365. 

  365.     static cilist io___40 = { 0, 6, 0, fmt_9962, 0 };
    366. 

  366.     static cilist io___41 = { 0, 6, 0, fmt_9974, 0 };
    367. 

  367.     static cilist io___42 = { 0, 6, 0, fmt_9967, 0 };
    368. 

  368.     static cilist io___43 = { 0, 6, 0, fmt_9971, 0 };
    369. 

  369.     static cilist io___44 = { 0, 6, 0, fmt_9970, 0 };
    370. 

  370.     static cilist io___45 = { 0, 6, 0, fmt_9969, 0 };
    371. 

  371.     static cilist io___46 = { 0, 6, 0, fmt_9968, 0 };
    372. 

  372.     static cilist io___47 = { 0, 5, 0, 0, 0 };
    373. 

  373.     static cilist io___50 = { 0, 6, 0, fmt_9992, 0 };
    374. 

  374.     static cilist io___54 = { 0, 6, 0, fmt_9972, 0 };
    375. 

  375.     static cilist io___55 = { 0, 6, 0, fmt_9984, 0 };
    376. 

  376.     static cilist io___56 = { 0, 5, 0, 0, 0 };
    377. 

  377.     static cilist io___58 = { 0, 6, 0, fmt_9989, 0 };
    378. 

  378.     static cilist io___59 = { 0, 6, 0, fmt_9988, 0 };
    379. 

  379.     static cilist io___60 = { 0, 5, 0, 0, 0 };
    380. 

  380.     static cilist io___64 = { 0, 6, 0, fmt_9989, 0 };
    381. 

  381.     static cilist io___65 = { 0, 6, 0, fmt_9988, 0 };
    382. 

  382.     static cilist io___66 = { 0, 6, 0, fmt_9983, 0 };
    383. 

  383.     static cilist io___67 = { 0, 5, 0, 0, 0 };
    384. 

  384.     static cilist io___69 = { 0, 6, 0, fmt_9989, 0 };
    385. 

  385.     static cilist io___70 = { 0, 6, 0, fmt_9988, 0 };
    386. 

  386.     static cilist io___71 = { 0, 6, 0, fmt_9983, 0 };
    387. 

  387.     static cilist io___72 = { 0, 5, 0, 0, 0 };
    388. 

  388.     static cilist io___74 = { 0, 6, 0, fmt_9989, 0 };
    389. 

  389.     static cilist io___75 = { 0, 6, 0, fmt_9988, 0 };
    390. 

  390.     static cilist io___76 = { 0, 6, 0, fmt_9983, 0 };
    391. 

  391.     static cilist io___77 = { 0, 6, 0, fmt_9983, 0 };
    392. 

  392.     static cilist io___78 = { 0, 5, 0, 0, 0 };
    393. 

  393.     static cilist io___80 = { 0, 5, 0, 0, 0 };
    394. 

  394.     static cilist io___82 = { 0, 6, 0, fmt_9989, 0 };
    395. 

  395.     static cilist io___83 = { 0, 6, 0, fmt_9988, 0 };
    396. 

  396.     static cilist io___84 = { 0, 6, 0, fmt_9983, 0 };
    397. 

  397.     static cilist io___85 = { 0, 5, 0, 0, 0 };
    398. 

  398.     static cilist io___94 = { 0, 6, 0, fmt_9989, 0 };
    399. 

  399.     static cilist io___95 = { 0, 6, 0, fmt_9989, 0 };
    400. 

  400.     static cilist io___96 = { 0, 6, 0, fmt_9989, 0 };
    401. 

  401.     static cilist io___97 = { 0, 6, 0, fmt_9989, 0 };
    402. 

  402.     static cilist io___98 = { 0, 6, 0, fmt_9989, 0 };
    403. 

  403.     static cilist io___99 = { 0, 6, 0, fmt_9989, 0 };
    404. 

  404.     static cilist io___100 = { 0, 6, 0, fmt_9989, 0 };
    405. 

  405.     static cilist io___101 = { 0, 6, 0, fmt_9989, 0 };
    406. 

  406.     static cilist io___102 = { 0, 6, 0, fmt_9983, 0 };
    407. 

  407.     static cilist io___103 = { 0, 6, 0, fmt_9983, 0 };
    408. 

  408.     static cilist io___104 = { 0, 6, 0, fmt_9983, 0 };
    409. 

  409.     static cilist io___105 = { 0, 6, 0, fmt_9983, 0 };
    410. 

  410.     static cilist io___106 = { 0, 6, 0, fmt_9983, 0 };
    411. 

  411.     static cilist io___107 = { 0, 6, 0, fmt_9983, 0 };
    412. 

  412.     static cilist io___108 = { 0, 6, 0, fmt_9983, 0 };
    413. 

  413.     static cilist io___109 = { 0, 6, 0, fmt_9983, 0 };
    414. 

  414.     static cilist io___110 = { 0, 5, 0, 0, 0 };
    415. 

  415.     static cilist io___113 = { 0, 6, 0, fmt_9989, 0 };
    416. 

  416.     static cilist io___114 = { 0, 6, 0, fmt_9989, 0 };
    417. 

  417.     static cilist io___115 = { 0, 6, 0, fmt_9989, 0 };
    418. 

  418.     static cilist io___116 = { 0, 6, 0, fmt_9989, 0 };
    419. 

  419.     static cilist io___117 = { 0, 6, 0, fmt_9989, 0 };
    420. 

  420.     static cilist io___118 = { 0, 6, 0, fmt_9983, 0 };
    421. 

  421.     static cilist io___119 = { 0, 6, 0, fmt_9983, 0 };
    422. 

  422.     static cilist io___120 = { 0, 6, 0, fmt_9983, 0 };
    423. 

  423.     static cilist io___121 = { 0, 6, 0, fmt_9983, 0 };
    424. 

  424.     static cilist io___122 = { 0, 6, 0, fmt_9983, 0 };
    425. 

  425.     static cilist io___123 = { 0, 5, 0, 0, 0 };
    426. 

  426.     static cilist io___125 = { 0, 6, 0, fmt_9989, 0 };
    427. 

  427.     static cilist io___126 = { 0, 6, 0, fmt_9988, 0 };
    428. 

  428.     static cilist io___127 = { 0, 5, 0, 0, 0 };
    429. 

  429.     static cilist io___128 = { 0, 6, 0, fmt_9989, 0 };
    430. 

  430.     static cilist io___129 = { 0, 6, 0, fmt_9988, 0 };
    431. 

  431.     static cilist io___130 = { 0, 6, 0, fmt_9983, 0 };
    432. 

  432.     static cilist io___131 = { 0, 5, 0, 0, 0 };
    433. 

  433.     static cilist io___132 = { 0, 6, 0, fmt_9989, 0 };
    434. 

  434.     static cilist io___133 = { 0, 6, 0, fmt_9988, 0 };
    435. 

  435.     static cilist io___134 = { 0, 6, 0, fmt_9983, 0 };
    436. 

  436.     static cilist io___135 = { 0, 5, 0, 0, 0 };
    437. 

  437.     static cilist io___136 = { 0, 6, 0, fmt_9989, 0 };
    438. 

  438.     static cilist io___137 = { 0, 6, 0, fmt_9988, 0 };
    439. 

  439.     static cilist io___138 = { 0, 6, 0, fmt_9983, 0 };
    440. 

  440.     static cilist io___139 = { 0, 5, 0, 0, 0 };
    441. 

  441.     static cilist io___140 = { 0, 6, 0, fmt_9989, 0 };
    442. 

  442.     static cilist io___141 = { 0, 6, 0, fmt_9988, 0 };
    443. 

  443.     static cilist io___142 = { 0, 6, 0, fmt_9983, 0 };
    444. 

  444.     static cilist io___143 = { 0, 5, 0, 0, 0 };
    445. 

  445.     static cilist io___144 = { 0, 6, 0, fmt_9989, 0 };
    446. 

  446.     static cilist io___145 = { 0, 6, 0, fmt_9988, 0 };
    447. 

  447.     static cilist io___146 = { 0, 6, 0, fmt_9983, 0 };
    448. 

  448.     static cilist io___147 = { 0, 5, 0, 0, 0 };
    449. 

  449.     static cilist io___148 = { 0, 6, 0, fmt_9989, 0 };
    450. 

  450.     static cilist io___149 = { 0, 6, 0, fmt_9983, 0 };
    451. 

  451.     static cilist io___150 = { 0, 5, 0, 0, 0 };
    452. 

  452.     static cilist io___151 = { 0, 6, 0, fmt_9989, 0 };
    453. 

  453.     static cilist io___152 = { 0, 6, 0, fmt_9983, 0 };
    454. 

  454.     static cilist io___153 = { 0, 5, 0, 0, 0 };
    455. 

  455.     static cilist io___154 = { 0, 6, 0, fmt_9989, 0 };
    456. 

  456.     static cilist io___155 = { 0, 6, 0, fmt_9983, 0 };
    457. 

  457.     static cilist io___156 = { 0, 5, 0, 0, 0 };
    458. 

  458.     static cilist io___157 = { 0, 6, 0, fmt_9989, 0 };
    459. 

  459.     static cilist io___158 = { 0, 6, 0, fmt_9983, 0 };
    460. 

  460.     static cilist io___159 = { 0, 5, 0, 0, 0 };
    461. 

  461.     static cilist io___160 = { 0, 6, 0, fmt_9989, 0 };
    462. 

  462.     static cilist io___161 = { 0, 6, 0, fmt_9983, 0 };
    463. 

  463.     static cilist io___162 = { 0, 5, 0, 0, 0 };
    464. 

  464.     static cilist io___164 = { 0, 6, 0, fmt_9989, 0 };
    465. 

  465.     static cilist io___165 = { 0, 6, 0, fmt_9988, 0 };
    466. 

  466.     static cilist io___166 = { 0, 6, 0, fmt_9983, 0 };
    467. 

  467.     static cilist io___167 = { 0, 6, 0, 0, 0 };
    468. 

  468.     static cilist io___169 = { 0, 6, 0, fmt_9981, 0 };
    469. 

  469.     static cilist io___170 = { 0, 6, 0, fmt_9981, 0 };
    470. 

  470.     static cilist io___171 = { 0, 6, 0, fmt_9981, 0 };
    471. 

  471.     static cilist io___172 = { 0, 5, 0, 0, 0 };
    472. 

  472.     static cilist io___173 = { 0, 6, 0, fmt_9982, 0 };
    473. 

  473.     static cilist io___174 = { 0, 5, 0, 0, 0 };
    474. 

  474.     static cilist io___176 = { 0, 5, 0, 0, 0 };
    475. 

  475.     static cilist io___178 = { 0, 5, 0, 0, 0 };
    476. 

  476.     static cilist io___179 = { 0, 5, 0, 0, 0 };
    477. 

  477.     static cilist io___181 = { 0, 5, 0, 0, 0 };
    478. 

  478.     static cilist io___183 = { 0, 6, 0, fmt_9999, 0 };
    479. 

  479.     static cilist io___192 = { 0, 6, 0, fmt_9991, 0 };
    480. 

  480.     static cilist io___193 = { 0, 6, 0, fmt_9990, 0 };
    481. 

  481.     static cilist io___196 = { 0, 6, 0, fmt_9961, 0 };
    482. 

  482.     static cilist io___204 = { 0, 6, 0, fmt_9980, 0 };
    483. 

  483.     static cilist io___205 = { 0, 6, 0, fmt_9997, 0 };
    484. 

  484.     static cilist io___206 = { 0, 6, 0, fmt_9980, 0 };
    485. 

  485.     static cilist io___207 = { 0, 6, 0, fmt_9980, 0 };
    486. 

  486.     static cilist io___208 = { 0, 6, 0, fmt_9997, 0 };
    487. 

  487.     static cilist io___209 = { 0, 6, 0, fmt_9980, 0 };
    488. 

  488.     static cilist io___211 = { 0, 6, 0, fmt_9995, 0 };
    489. 

  489.     static cilist io___212 = { 0, 6, 0, fmt_9980, 0 };
    490. 

  490.     static cilist io___213 = { 0, 6, 0, fmt_9990, 0 };
    491. 

  491.     static cilist io___214 = { 0, 6, 0, fmt_9980, 0 };
    492. 

  492.     static cilist io___215 = { 0, 6, 0, fmt_9973, 0 };
    493. 

  493.     static cilist io___216 = { 0, 6, 0, fmt_9990, 0 };
    494. 

  494.     static cilist io___217 = { 0, 6, 0, fmt_9980, 0 };
    495. 

  495.     static cilist io___218 = { 0, 6, 0, fmt_9973, 0 };
    496. 

  496.     static cilist io___219 = { 0, 6, 0, fmt_9990, 0 };
    497. 

  497.     static cilist io___220 = { 0, 6, 0, fmt_9980, 0 };
    498. 

  498.     static cilist io___221 = { 0, 6, 0, fmt_9973, 0 };
    499. 

  499.     static cilist io___222 = { 0, 6, 0, fmt_9990, 0 };
    500. 

  500.     static cilist io___223 = { 0, 6, 0, fmt_9980, 0 };
    501. 

  501.     static cilist io___224 = { 0, 6, 0, fmt_9973, 0 };
    502. 

  502.     static cilist io___225 = { 0, 6, 0, fmt_9996, 0 };
    503. 

  503.     static cilist io___228 = { 0, 6, 0, fmt_9980, 0 };
    504. 

  504.     static cilist io___229 = { 0, 6, 0, fmt_9980, 0 };
    505. 

  505.     static cilist io___230 = { 0, 6, 0, fmt_9990, 0 };
    506. 

  506.     static cilist io___231 = { 0, 6, 0, fmt_9980, 0 };
    507. 

  507.     static cilist io___232 = { 0, 6, 0, fmt_9973, 0 };
    508. 

  508.     static cilist io___233 = { 0, 6, 0, fmt_9990, 0 };
    509. 

  509.     static cilist io___235 = { 0, 6, 0, fmt_9980, 0 };
    510. 

  510.     static cilist io___236 = { 0, 6, 0, fmt_9973, 0 };
    511. 

  511.     static cilist io___237 = { 0, 6, 0, fmt_9990, 0 };
    512. 

  512.     static cilist io___238 = { 0, 6, 0, fmt_9980, 0 };
    513. 

  513.     static cilist io___239 = { 0, 6, 0, fmt_9973, 0 };
    514. 

  514.     static cilist io___240 = { 0, 6, 0, fmt_9990, 0 };
    515. 

  515.     static cilist io___241 = { 0, 6, 0, fmt_9980, 0 };
    516. 

  516.     static cilist io___242 = { 0, 6, 0, fmt_9973, 0 };
    517. 

  517.     static cilist io___243 = { 0, 6, 0, fmt_9980, 0 };
    518. 

  518.     static cilist io___244 = { 0, 6, 0, fmt_9966, 0 };
    519. 

  519.     static cilist io___245 = { 0, 6, 0, fmt_9980, 0 };
    520. 

  520.     static cilist io___248 = { 0, 6, 0, fmt_9980, 0 };
    521. 

  521.     static cilist io___251 = { 0, 6, 0, fmt_9980, 0 };
    522. 

  522.     static cilist io___252 = { 0, 6, 0, fmt_9980, 0 };
    523. 

  523.     static cilist io___253 = { 0, 6, 0, fmt_9980, 0 };
    524. 

  524.     static cilist io___254 = { 0, 6, 0, 0, 0 };
    525. 

  525.     static cilist io___255 = { 0, 6, 0, 0, 0 };
    526. 

  526.     static cilist io___256 = { 0, 6, 0, fmt_9992, 0 };
    527. 

  527.     static cilist io___257 = { 0, 6, 0, fmt_9994, 0 };
    528. 

  528.     static cilist io___259 = { 0, 6, 0, fmt_9993, 0 };
    529. 

  529. 

  530. 

  531. 

  532. /*  -- LAPACK test routine (version 3.1.1) -- */
    533. 

  533. /*     Univ. of Tennessee, Univ. of California Berkeley and NAG Ltd.. */
    534. 

  534. /*     January 2007 */
    535. 

  535. 

  536. /*  Purpose */
    537. 

  537. /*  ======= */
    538. 

  538. 

  539. /*  DCHKEE tests the DOUBLE PRECISION LAPACK subroutines for the matrix */
    540. 

  540. /*  eigenvalue problem.  The test paths in this version are */
    541. 

  541. 

  542. /*  NEP (Nonsymmetric Eigenvalue Problem): */
    543. 

  543. /*      Test DGEHRD, DORGHR, DHSEQR, DTREVC, DHSEIN, and DORMHR */
    544. 

  544. 

  545. /*  SEP (Symmetric Eigenvalue Problem): */
    546. 

  546. /*      Test DSYTRD, DORGTR, DSTEQR, DSTERF, DSTEIN, DSTEDC, */
    547. 

  547. /*      and drivers DSYEV(X), DSBEV(X), DSPEV(X), DSTEV(X), */
    548. 

  548. /*                  DSYEVD,   DSBEVD,   DSPEVD,   DSTEVD */
    549. 

  549. 

  550. /*  SVD (Singular Value Decomposition): */
    551. 

  551. /*      Test DGEBRD, DORGBR, DBDSQR, DBDSDC */
    552. 

  552. /*      and the drivers DGESVD, DGESDD */
    553. 

  553. 

  554. /*  DEV (Nonsymmetric Eigenvalue/eigenvector Driver): */
    555. 

  555. /*      Test DGEEV */
    556. 

  556. 

  557. /*  DES (Nonsymmetric Schur form Driver): */
    558. 

  558. /*      Test DGEES */
    559. 

  559. 

  560. /*  DVX (Nonsymmetric Eigenvalue/eigenvector Expert Driver): */
    561. 

  561. /*      Test DGEEVX */
    562. 

  562. 

  563. /*  DSX (Nonsymmetric Schur form Expert Driver): */
    564. 

  564. /*      Test DGEESX */
    565. 

  565. 

  566. /*  DGG (Generalized Nonsymmetric Eigenvalue Problem): */
    567. 

  567. /*      Test DGGHRD, DGGBAL, DGGBAK, DHGEQZ, and DTGEVC */
    568. 

  568. /*      and the driver routines DGEGS and DGEGV */
    569. 

  569. 

  570. /*  DGS (Generalized Nonsymmetric Schur form Driver): */
    571. 

  571. /*      Test DGGES */
    572. 

  572. 

  573. /*  DGV (Generalized Nonsymmetric Eigenvalue/eigenvector Driver): */
    574. 

  574. /*      Test DGGEV */
    575. 

  575. 

  576. /*  DGX (Generalized Nonsymmetric Schur form Expert Driver): */
    577. 

  577. /*      Test DGGESX */
    578. 

  578. 

  579. /*  DXV (Generalized Nonsymmetric Eigenvalue/eigenvector Expert Driver): */
    580. 

  580. /*      Test DGGEVX */
    581. 

  581. 

  582. /*  DSG (Symmetric Generalized Eigenvalue Problem): */
    583. 

  583. /*      Test DSYGST, DSYGV, DSYGVD, DSYGVX, DSPGST, DSPGV, DSPGVD, */
    584. 

  584. /*      DSPGVX, DSBGST, DSBGV, DSBGVD, and DSBGVX */
    585. 

  585. 

  586. /*  DSB (Symmetric Band Eigenvalue Problem): */
    587. 

  587. /*      Test DSBTRD */
    588. 

  588. 

  589. /*  DBB (Band Singular Value Decomposition): */
    590. 

  590. /*      Test DGBBRD */
    591. 

  591. 

  592. /*  DEC (Eigencondition estimation): */
    593. 

  593. /*      Test DLALN2, DLASY2, DLAEQU, DLAEXC, DTRSYL, DTREXC, DTRSNA, */
    594. 

  594. /*      DTRSEN, and DLAQTR */
    595. 

  595. 

  596. /*  DBL (Balancing a general matrix) */
    597. 

  597. /*      Test DGEBAL */
    598. 

  598. 

  599. /*  DBK (Back transformation on a balanced matrix) */
    600. 

  600. /*      Test DGEBAK */
    601. 

  601. 

  602. /*  DGL (Balancing a matrix pair) */
    603. 

  603. /*      Test DGGBAL */
    604. 

  604. 

  605. /*  DGK (Back transformation on a matrix pair) */
    606. 

  606. /*      Test DGGBAK */
    607. 

  607. 

  608. /*  GLM (Generalized Linear Regression Model): */
    609. 

  609. /*      Tests DGGGLM */
    610. 

  610. 

  611. /*  GQR (Generalized QR and RQ factorizations): */
    612. 

  612. /*      Tests DGGQRF and DGGRQF */
    613. 

  613. 

  614. /*  GSV (Generalized Singular Value Decomposition): */
    615. 

  615. /*      Tests DGGSVD, DGGSVP, DTGSJA, DLAGS2, DLAPLL, and DLAPMT */
    616. 

  616. 

  617. /*  LSE (Constrained Linear Least Squares): */
    618. 

  618. /*      Tests DGGLSE */
    619. 

  619. 

  620. /*  Each test path has a different set of inputs, but the data sets for */
    621. 

  621. /*  the driver routines xEV, xES, xVX, and xSX can be concatenated in a */
    622. 

  622. /*  single input file.  The first line of input should contain one of the */
    623. 

  623. /*  3-character path names in columns 1-3.  The number of remaining lines */
    624. 

  624. /*  depends on what is found on the first line. */
    625. 

  625. 

  626. /*  The number of matrix types used in testing is often controllable from */
    627. 

  627. /*  the input file.  The number of matrix types for each path, and the */
    628. 

  628. /*  test routine that describes them, is as follows: */
    629. 

  629. 

  630. /*  Path name(s)  Types    Test routine */
    631. 

  631. 

  632. /*  DHS or NEP      21     DCHKHS */
    633. 

  633. /*  DST or SEP      21     DCHKST (routines) */
    634. 

  634. /*                  18     DDRVST (drivers) */
    635. 

  635. /*  DBD or SVD      16     DCHKBD (routines) */
    636. 

  636. /*                   5     DDRVBD (drivers) */
    637. 

  637. /*  DEV             21     DDRVEV */
    638. 

  638. /*  DES             21     DDRVES */
    639. 

  639. /*  DVX             21     DDRVVX */
    640. 

  640. /*  DSX             21     DDRVSX */
    641. 

  641. /*  DGG             26     DCHKGG (routines) */
    642. 

  642. /*                  26     DDRVGG (drivers) */
    643. 

  643. /*  DGS             26     DDRGES */
    644. 

  644. /*  DGX              5     DDRGSX */
    645. 

  645. /*  DGV             26     DDRGEV */
    646. 

  646. /*  DXV              2     DDRGVX */
    647. 

  647. /*  DSG             21     DDRVSG */
    648. 

  648. /*  DSB             15     DCHKSB */
    649. 

  649. /*  DBB             15     DCHKBB */
    650. 

  650. /*  DEC              -     DCHKEC */
    651. 

  651. /*  DBL              -     DCHKBL */
    652. 

  652. /*  DBK              -     DCHKBK */
    653. 

  653. /*  DGL              -     DCHKGL */
    654. 

  654. /*  DGK              -     DCHKGK */
    655. 

  655. /*  GLM              8     DCKGLM */
    656. 

  656. /*  GQR              8     DCKGQR */
    657. 

  657. /*  GSV              8     DCKGSV */
    658. 

  658. /*  LSE              8     DCKLSE */
    659. 

  659. 

  660. /* ----------------------------------------------------------------------- */
    661. 

  661. 

  662. /*  NEP input file: */
    663. 

  663. 

  664. /*  line 2:  NN, INTEGER */
    665. 

  665. /*           Number of values of N. */
    666. 

  666. 

  667. /*  line 3:  NVAL, INTEGER array, dimension (NN) */
    668. 

  668. /*           The values for the matrix dimension N. */
    669. 

  669. 

  670. /*  line 4:  NPARMS, INTEGER */
    671. 

  671. /*           Number of values of the parameters NB, NBMIN, NX, NS, and */
    672. 

  672. /*           MAXB. */
    673. 

  673. 

  674. /*  line 5:  NBVAL, INTEGER array, dimension (NPARMS) */
    675. 

  675. /*           The values for the blocksize NB. */
    676. 

  676. 

  677. /*  line 6:  NBMIN, INTEGER array, dimension (NPARMS) */
    678. 

  678. /*           The values for the minimum blocksize NBMIN. */
    679. 

  679. 

  680. /*  line 7:  NXVAL, INTEGER array, dimension (NPARMS) */
    681. 

  681. /*           The values for the crossover point NX. */
    682. 

  682. 

  683. /*  line 8:  INMIN, INTEGER array, dimension (NPARMS) */
    684. 

  684. /*           LAHQR vs TTQRE crossover point, >= 11 */
    685. 

  685. 

  686. /*  line 9:  INWIN, INTEGER array, dimension (NPARMS) */
    687. 

  687. /*           recommended deflation window size */
    688. 

  688. 

  689. /*  line 10: INIBL, INTEGER array, dimension (NPARMS) */
    690. 

  690. /*           nibble crossover point */
    691. 

  691. 

  692. /*  line 11: ISHFTS, INTEGER array, dimension (NPARMS) */
    693. 

  693. /*           number of simultaneous shifts) */
    694. 

  694. 

  695. /*  line 12: IACC22, INTEGER array, dimension (NPARMS) */
    696. 

  696. /*           select structured matrix multiply: 0, 1 or 2) */
    697. 

  697. 

  698. /*  line 13: THRESH */
    699. 

  699. /*           Threshold value for the test ratios.  Information will be */
    700. 

  700. /*           printed about each test for which the test ratio is greater */
    701. 

  701. /*           than or equal to the threshold.  To have all of the test */
    702. 

  702. /*           ratios printed, use THRESH = 0.0 . */
    703. 

  703. 

  704. /*  line 14: NEWSD, INTEGER */
    705. 

  705. /*           A code indicating how to set the random number seed. */
    706. 

  706. /*           = 0:  Set the seed to a default value before each run */
    707. 

  707. /*           = 1:  Initialize the seed to a default value only before the */
    708. 

  708. /*                 first run */
    709. 

  709. /*           = 2:  Like 1, but use the seed values on the next line */
    710. 

  710. 

  711. /*  If line 14 was 2: */
    712. 

  712. 

  713. /*  line 15: INTEGER array, dimension (4) */
    714. 

  714. /*           Four integer values for the random number seed. */
    715. 

  715. 

  716. /*  lines 15-EOF:  The remaining lines occur in sets of 1 or 2 and allow */
    717. 

  717. /*           the user to specify the matrix types.  Each line contains */
    718. 

  718. /*           a 3-character path name in columns 1-3, and the number */
    719. 

  719. /*           of matrix types must be the first nonblank item in columns */
    720. 

  720. /*           4-80.  If the number of matrix types is at least 1 but is */
    721. 

  721. /*           less than the maximum number of possible types, a second */
    722. 

  722. /*           line will be read to get the numbers of the matrix types to */
    723. 

  723. /*           be used.  For example, */
    724. 

  724. /*  NEP 21 */
    725. 

  725. /*           requests all of the matrix types for the nonsymmetric */
    726. 

  726. /*           eigenvalue problem, while */
    727. 

  727. /*  NEP  4 */
    728. 

  728. /*  9 10 11 12 */
    729. 

  729. /*           requests only matrices of type 9, 10, 11, and 12. */
    730. 

  730. 

  731. /*           The valid 3-character path names are 'NEP' or 'SHS' for the */
    732. 

  732. /*           nonsymmetric eigenvalue routines. */
    733. 

  733. 

  734. /* ----------------------------------------------------------------------- */
    735. 

  735. 

  736. /*  SEP or DSG input file: */
    737. 

  737. 

  738. /*  line 2:  NN, INTEGER */
    739. 

  739. /*           Number of values of N. */
    740. 

  740. 

  741. /*  line 3:  NVAL, INTEGER array, dimension (NN) */
    742. 

  742. /*           The values for the matrix dimension N. */
    743. 

  743. 

  744. /*  line 4:  NPARMS, INTEGER */
    745. 

  745. /*           Number of values of the parameters NB, NBMIN, and NX. */
    746. 

  746. 

  747. /*  line 5:  NBVAL, INTEGER array, dimension (NPARMS) */
    748. 

  748. /*           The values for the blocksize NB. */
    749. 

  749. 

  750. /*  line 6:  NBMIN, INTEGER array, dimension (NPARMS) */
    751. 

  751. /*           The values for the minimum blocksize NBMIN. */
    752. 

  752. 

  753. /*  line 7:  NXVAL, INTEGER array, dimension (NPARMS) */
    754. 

  754. /*           The values for the crossover point NX. */
    755. 

  755. 

  756. /*  line 8:  THRESH */
    757. 

  757. /*           Threshold value for the test ratios.  Information will be */
    758. 

  758. /*           printed about each test for which the test ratio is greater */
    759. 

  759. /*           than or equal to the threshold. */
    760. 

  760. 

  761. /*  line 9:  TSTCHK, LOGICAL */
    762. 

  762. /*           Flag indicating whether or not to test the LAPACK routines. */
    763. 

  763. 

  764. /*  line 10: TSTDRV, LOGICAL */
    765. 

  765. /*           Flag indicating whether or not to test the driver routines. */
    766. 

  766. 

  767. /*  line 11: TSTERR, LOGICAL */
    768. 

  768. /*           Flag indicating whether or not to test the error exits for */
    769. 

  769. /*           the LAPACK routines and driver routines. */
    770. 

  770. 

  771. /*  line 12: NEWSD, INTEGER */
    772. 

  772. /*           A code indicating how to set the random number seed. */
    773. 

  773. /*           = 0:  Set the seed to a default value before each run */
    774. 

  774. /*           = 1:  Initialize the seed to a default value only before the */
    775. 

  775. /*                 first run */
    776. 

  776. /*           = 2:  Like 1, but use the seed values on the next line */
    777. 

  777. 

  778. /*  If line 12 was 2: */
    779. 

  779. 

  780. /*  line 13: INTEGER array, dimension (4) */
    781. 

  781. /*           Four integer values for the random number seed. */
    782. 

  782. 

  783. /*  lines 13-EOF:  Lines specifying matrix types, as for NEP. */
    784. 

  784. /*           The 3-character path names are 'SEP' or 'SST' for the */
    785. 

  785. /*           symmetric eigenvalue routines and driver routines, and */
    786. 

  786. /*           'DSG' for the routines for the symmetric generalized */
    787. 

  787. /*           eigenvalue problem. */
    788. 

  788. 

  789. /* ----------------------------------------------------------------------- */
    790. 

  790. 

  791. /*  SVD input file: */
    792. 

  792. 

  793. /*  line 2:  NN, INTEGER */
    794. 

  794. /*           Number of values of M and N. */
    795. 

  795. 

  796. /*  line 3:  MVAL, INTEGER array, dimension (NN) */
    797. 

  797. /*           The values for the matrix row dimension M. */
    798. 

  798. 

  799. /*  line 4:  NVAL, INTEGER array, dimension (NN) */
    800. 

  800. /*           The values for the matrix column dimension N. */
    801. 

  801. 

  802. /*  line 5:  NPARMS, INTEGER */
    803. 

  803. /*           Number of values of the parameter NB, NBMIN, NX, and NRHS. */
    804. 

  804. 

  805. /*  line 6:  NBVAL, INTEGER array, dimension (NPARMS) */
    806. 

  806. /*           The values for the blocksize NB. */
    807. 

  807. 

  808. /*  line 7:  NBMIN, INTEGER array, dimension (NPARMS) */
    809. 

  809. /*           The values for the minimum blocksize NBMIN. */
    810. 

  810. 

  811. /*  line 8:  NXVAL, INTEGER array, dimension (NPARMS) */
    812. 

  812. /*           The values for the crossover point NX. */
    813. 

  813. 

  814. /*  line 9:  NSVAL, INTEGER array, dimension (NPARMS) */
    815. 

  815. /*           The values for the number of right hand sides NRHS. */
    816. 

  816. 

  817. /*  line 10: THRESH */
    818. 

  818. /*           Threshold value for the test ratios.  Information will be */
    819. 

  819. /*           printed about each test for which the test ratio is greater */
    820. 

  820. /*           than or equal to the threshold. */
    821. 

  821. 

  822. /*  line 11: TSTCHK, LOGICAL */
    823. 

  823. /*           Flag indicating whether or not to test the LAPACK routines. */
    824. 

  824. 

  825. /*  line 12: TSTDRV, LOGICAL */
    826. 

  826. /*           Flag indicating whether or not to test the driver routines. */
    827. 

  827. 

  828. /*  line 13: TSTERR, LOGICAL */
    829. 

  829. /*           Flag indicating whether or not to test the error exits for */
    830. 

  830. /*           the LAPACK routines and driver routines. */
    831. 

  831. 

  832. /*  line 14: NEWSD, INTEGER */
    833. 

  833. /*           A code indicating how to set the random number seed. */
    834. 

  834. /*           = 0:  Set the seed to a default value before each run */
    835. 

  835. /*           = 1:  Initialize the seed to a default value only before the */
    836. 

  836. /*                 first run */
    837. 

  837. /*           = 2:  Like 1, but use the seed values on the next line */
    838. 

  838. 

  839. /*  If line 14 was 2: */
    840. 

  840. 

  841. /*  line 15: INTEGER array, dimension (4) */
    842. 

  842. /*           Four integer values for the random number seed. */
    843. 

  843. 

  844. /*  lines 15-EOF:  Lines specifying matrix types, as for NEP. */
    845. 

  845. /*           The 3-character path names are 'SVD' or 'SBD' for both the */
    846. 

  846. /*           SVD routines and the SVD driver routines. */
    847. 

  847. 

  848. /* ----------------------------------------------------------------------- */
    849. 

  849. 

  850. /*  DEV and DES data files: */
    851. 

  851. 

  852. /*  line 1:  'DEV' or 'DES' in columns 1 to 3. */
    853. 

  853. 

  854. /*  line 2:  NSIZES, INTEGER */
    855. 

  855. /*           Number of sizes of matrices to use. Should be at least 0 */
    856. 

  856. /*           and at most 20. If NSIZES = 0, no testing is done */
    857. 

  857. /*           (although the remaining  3 lines are still read). */
    858. 

  858. 

  859. /*  line 3:  NN, INTEGER array, dimension(NSIZES) */
    860. 

  860. /*           Dimensions of matrices to be tested. */
    861. 

  861. 

  862. /*  line 4:  NB, NBMIN, NX, NS, NBCOL, INTEGERs */
    863. 

  863. /*           These integer parameters determine how blocking is done */
    864. 

  864. /*           (see ILAENV for details) */
    865. 

  865. /*           NB     : block size */
    866. 

  866. /*           NBMIN  : minimum block size */
    867. 

  867. /*           NX     : minimum dimension for blocking */
    868. 

  868. /*           NS     : number of shifts in xHSEQR */
    869. 

  869. /*           NBCOL  : minimum column dimension for blocking */
    870. 

  870. 

  871. /*  line 5:  THRESH, REAL */
    872. 

  872. /*           The test threshold against which computed residuals are */
    873. 

  873. /*           compared. Should generally be in the range from 10. to 20. */
    874. 

  874. /*           If it is 0., all test case data will be printed. */
    875. 

  875. 

  876. /*  line 6:  TSTERR, LOGICAL */
    877. 

  877. /*           Flag indicating whether or not to test the error exits. */
    878. 

  878. 

  879. /*  line 7:  NEWSD, INTEGER */
    880. 

  880. /*           A code indicating how to set the random number seed. */
    881. 

  881. /*           = 0:  Set the seed to a default value before each run */
    882. 

  882. /*           = 1:  Initialize the seed to a default value only before the */
    883. 

  883. /*                 first run */
    884. 

  884. /*           = 2:  Like 1, but use the seed values on the next line */
    885. 

  885. 

  886. /*  If line 7 was 2: */
    887. 

  887. 

  888. /*  line 8:  INTEGER array, dimension (4) */
    889. 

  889. /*           Four integer values for the random number seed. */
    890. 

  890. 

  891. /*  lines 9 and following:  Lines specifying matrix types, as for NEP. */
    892. 

  892. /*           The 3-character path name is 'DEV' to test SGEEV, or */
    893. 

  893. /*           'DES' to test SGEES. */
    894. 

  894. 

  895. /* ----------------------------------------------------------------------- */
    896. 

  896. 

  897. /*  The DVX data has two parts. The first part is identical to DEV, */
    898. 

  898. /*  and the second part consists of test matrices with precomputed */
    899. 

  899. /*  solutions. */
    900. 

  900. 

  901. /*  line 1:  'DVX' in columns 1-3. */
    902. 

  902. 

  903. /*  line 2:  NSIZES, INTEGER */
    904. 

  904. /*           If NSIZES = 0, no testing of randomly generated examples */
    905. 

  905. /*           is done, but any precomputed examples are tested. */
    906. 

  906. 

  907. /*  line 3:  NN, INTEGER array, dimension(NSIZES) */
    908. 

  908. 

  909. /*  line 4:  NB, NBMIN, NX, NS, NBCOL, INTEGERs */
    910. 

  910. 

  911. /*  line 5:  THRESH, REAL */
    912. 

  912. 

  913. /*  line 6:  TSTERR, LOGICAL */
    914. 

  914. 

  915. /*  line 7:  NEWSD, INTEGER */
    916. 

  916. 

  917. /*  If line 7 was 2: */
    918. 

  918. 

  919. /*  line 8:  INTEGER array, dimension (4) */
    920. 

  920. 

  921. /*  lines 9 and following: The first line contains 'DVX' in columns 1-3 */
    922. 

  922. /*           followed by the number of matrix types, possibly with */
    923. 

  923. /*           a second line to specify certain matrix types. */
    924. 

  924. /*           If the number of matrix types = 0, no testing of randomly */
    925. 

  925. /*           generated examples is done, but any precomputed examples */
    926. 

  926. /*           are tested. */
    927. 

  927. 

  928. /*  remaining lines : Each matrix is stored on 1+2*N lines, where N is */
    929. 

  929. /*           its dimension. The first line contains the dimension (a */
    930. 

  930. /*           single integer). The next N lines contain the matrix, one */
    931. 

  931. /*           row per line. The last N lines correspond to each */
    932. 

  932. /*           eigenvalue. Each of these last N lines contains 4 real */
    933. 

  933. /*           values: the real part of the eigenvalue, the imaginary */
    934. 

  934. /*           part of the eigenvalue, the reciprocal condition number of */
    935. 

  935. /*           the eigenvalues, and the reciprocal condition number of the */
    936. 

  936. /*           eigenvector.  The end of data is indicated by dimension N=0. */
    937. 

  937. /*           Even if no data is to be tested, there must be at least one */
    938. 

  938. /*           line containing N=0. */
    939. 

  939. 

  940. /* ----------------------------------------------------------------------- */
    941. 

  941. 

  942. /*  The DSX data is like DVX. The first part is identical to DEV, and the */
    943. 

  943. /*  second part consists of test matrices with precomputed solutions. */
    944. 

  944. 

  945. /*  line 1:  'DSX' in columns 1-3. */
    946. 

  946. 

  947. /*  line 2:  NSIZES, INTEGER */
    948. 

  948. /*           If NSIZES = 0, no testing of randomly generated examples */
    949. 

  949. /*           is done, but any precomputed examples are tested. */
    950. 

  950. 

  951. /*  line 3:  NN, INTEGER array, dimension(NSIZES) */
    952. 

  952. 

  953. /*  line 4:  NB, NBMIN, NX, NS, NBCOL, INTEGERs */
    954. 

  954. 

  955. /*  line 5:  THRESH, REAL */
    956. 

  956. 

  957. /*  line 6:  TSTERR, LOGICAL */
    958. 

  958. 

  959. /*  line 7:  NEWSD, INTEGER */
    960. 

  960. 

  961. /*  If line 7 was 2: */
    962. 

  962. 

  963. /*  line 8:  INTEGER array, dimension (4) */
    964. 

  964. 

  965. /*  lines 9 and following: The first line contains 'DSX' in columns 1-3 */
    966. 

  966. /*           followed by the number of matrix types, possibly with */
    967. 

  967. /*           a second line to specify certain matrix types. */
    968. 

  968. /*           If the number of matrix types = 0, no testing of randomly */
    969. 

  969. /*           generated examples is done, but any precomputed examples */
    970. 

  970. /*           are tested. */
    971. 

  971. 

  972. /*  remaining lines : Each matrix is stored on 3+N lines, where N is its */
    973. 

  973. /*           dimension. The first line contains the dimension N and the */
    974. 

  974. /*           dimension M of an invariant subspace. The second line */
    975. 

  975. /*           contains M integers, identifying the eigenvalues in the */
    976. 

  976. /*           invariant subspace (by their position in a list of */
    977. 

  977. /*           eigenvalues ordered by increasing real part). The next N */
    978. 

  978. /*           lines contain the matrix. The last line contains the */
    979. 

  979. /*           reciprocal condition number for the average of the selected */
    980. 

  980. /*           eigenvalues, and the reciprocal condition number for the */
    981. 

  981. /*           corresponding right invariant subspace. The end of data is */
    982. 

  982. /*           indicated by a line containing N=0 and M=0. Even if no data */
    983. 

  983. /*           is to be tested, there must be at least one line containing */
    984. 

  984. /*           N=0 and M=0. */
    985. 

  985. 

  986. /* ----------------------------------------------------------------------- */
    987. 

  987. 

  988. /*  DGG input file: */
    989. 

  989. 

  990. /*  line 2:  NN, INTEGER */
    991. 

  991. /*           Number of values of N. */
    992. 

  992. 

  993. /*  line 3:  NVAL, INTEGER array, dimension (NN) */
    994. 

  994. /*           The values for the matrix dimension N. */
    995. 

  995. 

  996. /*  line 4:  NPARMS, INTEGER */
    997. 

  997. /*           Number of values of the parameters NB, NBMIN, NS, MAXB, and */
    998. 

  998. /*           NBCOL. */
    999. 

  999. 

  1000. /*  line 5:  NBVAL, INTEGER array, dimension (NPARMS) */
    1001. 

  1001. /*           The values for the blocksize NB. */
    1002. 

  1002. 

  1003. /*  line 6:  NBMIN, INTEGER array, dimension (NPARMS) */
    1004. 

  1004. /*           The values for NBMIN, the minimum row dimension for blocks. */
    1005. 

  1005. 

  1006. /*  line 7:  NSVAL, INTEGER array, dimension (NPARMS) */
    1007. 

  1007. /*           The values for the number of shifts. */
    1008. 

  1008. 

  1009. /*  line 8:  MXBVAL, INTEGER array, dimension (NPARMS) */
    1010. 

  1010. /*           The values for MAXB, used in determining minimum blocksize. */
    1011. 

  1011. 

  1012. /*  line 9:  NBCOL, INTEGER array, dimension (NPARMS) */
    1013. 

  1013. /*           The values for NBCOL, the minimum column dimension for */
    1014. 

  1014. /*           blocks. */
    1015. 

  1015. 

  1016. /*  line 10: THRESH */
    1017. 

  1017. /*           Threshold value for the test ratios.  Information will be */
    1018. 

  1018. /*           printed about each test for which the test ratio is greater */
    1019. 

  1019. /*           than or equal to the threshold. */
    1020. 

  1020. 

  1021. /*  line 11: TSTCHK, LOGICAL */
    1022. 

  1022. /*           Flag indicating whether or not to test the LAPACK routines. */
    1023. 

  1023. 

  1024. /*  line 12: TSTDRV, LOGICAL */
    1025. 

  1025. /*           Flag indicating whether or not to test the driver routines. */
    1026. 

  1026. 

  1027. /*  line 13: TSTERR, LOGICAL */
    1028. 

  1028. /*           Flag indicating whether or not to test the error exits for */
    1029. 

  1029. /*           the LAPACK routines and driver routines. */
    1030. 

  1030. 

  1031. /*  line 14: NEWSD, INTEGER */
    1032. 

  1032. /*           A code indicating how to set the random number seed. */
    1033. 

  1033. /*           = 0:  Set the seed to a default value before each run */
    1034. 

  1034. /*           = 1:  Initialize the seed to a default value only before the */
    1035. 

  1035. /*                 first run */
    1036. 

  1036. /*           = 2:  Like 1, but use the seed values on the next line */
    1037. 

  1037. 

  1038. /*  If line 14 was 2: */
    1039. 

  1039. 

  1040. /*  line 15: INTEGER array, dimension (4) */
    1041. 

  1041. /*           Four integer values for the random number seed. */
    1042. 

  1042. 

  1043. /*  lines 15-EOF:  Lines specifying matrix types, as for NEP. */
    1044. 

  1044. /*           The 3-character path name is 'DGG' for the generalized */
    1045. 

  1045. /*           eigenvalue problem routines and driver routines. */
    1046. 

  1046. 

  1047. /* ----------------------------------------------------------------------- */
    1048. 

  1048. 

  1049. /*  DGS and DGV input files: */
    1050. 

  1050. 

  1051. /*  line 1:  'DGS' or 'DGV' in columns 1 to 3. */
    1052. 

  1052. 

  1053. /*  line 2:  NN, INTEGER */
    1054. 

  1054. /*           Number of values of N. */
    1055. 

  1055. 

  1056. /*  line 3:  NVAL, INTEGER array, dimension(NN) */
    1057. 

  1057. /*           Dimensions of matrices to be tested. */
    1058. 

  1058. 

  1059. /*  line 4:  NB, NBMIN, NX, NS, NBCOL, INTEGERs */
    1060. 

  1060. /*           These integer parameters determine how blocking is done */
    1061. 

  1061. /*           (see ILAENV for details) */
    1062. 

  1062. /*           NB     : block size */
    1063. 

  1063. /*           NBMIN  : minimum block size */
    1064. 

  1064. /*           NX     : minimum dimension for blocking */
    1065. 

  1065. /*           NS     : number of shifts in xHGEQR */
    1066. 

  1066. /*           NBCOL  : minimum column dimension for blocking */
    1067. 

  1067. 

  1068. /*  line 5:  THRESH, REAL */
    1069. 

  1069. /*           The test threshold against which computed residuals are */
    1070. 

  1070. /*           compared. Should generally be in the range from 10. to 20. */
    1071. 

  1071. /*           If it is 0., all test case data will be printed. */
    1072. 

  1072. 

  1073. /*  line 6:  TSTERR, LOGICAL */
    1074. 

  1074. /*           Flag indicating whether or not to test the error exits. */
    1075. 

  1075. 

  1076. /*  line 7:  NEWSD, INTEGER */
    1077. 

  1077. /*           A code indicating how to set the random number seed. */
    1078. 

  1078. /*           = 0:  Set the seed to a default value before each run */
    1079. 

  1079. /*           = 1:  Initialize the seed to a default value only before the */
    1080. 

  1080. /*                 first run */
    1081. 

  1081. /*           = 2:  Like 1, but use the seed values on the next line */
    1082. 

  1082. 

  1083. /*  If line 17 was 2: */
    1084. 

  1084. 

  1085. /*  line 7:  INTEGER array, dimension (4) */
    1086. 

  1086. /*           Four integer values for the random number seed. */
    1087. 

  1087. 

  1088. /*  lines 7-EOF:  Lines specifying matrix types, as for NEP. */
    1089. 

  1089. /*           The 3-character path name is 'DGS' for the generalized */
    1090. 

  1090. /*           eigenvalue problem routines and driver routines. */
    1091. 

  1091. 

  1092. /* ----------------------------------------------------------------------- */
    1093. 

  1093. 

  1094. /*  DXV input files: */
    1095. 

  1095. 

  1096. /*  line 1:  'DXV' in columns 1 to 3. */
    1097. 

  1097. 

  1098. /*  line 2:  N, INTEGER */
    1099. 

  1099. /*           Value of N. */
    1100. 

  1100. 

  1101. /*  line 3:  NB, NBMIN, NX, NS, NBCOL, INTEGERs */
    1102. 

  1102. /*           These integer parameters determine how blocking is done */
    1103. 

  1103. /*           (see ILAENV for details) */
    1104. 

  1104. /*           NB     : block size */
    1105. 

  1105. /*           NBMIN  : minimum block size */
    1106. 

  1106. /*           NX     : minimum dimension for blocking */
    1107. 

  1107. /*           NS     : number of shifts in xHGEQR */
    1108. 

  1108. /*           NBCOL  : minimum column dimension for blocking */
    1109. 

  1109. 

  1110. /*  line 4:  THRESH, REAL */
    1111. 

  1111. /*           The test threshold against which computed residuals are */
    1112. 

  1112. /*           compared. Should generally be in the range from 10. to 20. */
    1113. 

  1113. /*           Information will be printed about each test for which the */
    1114. 

  1114. /*           test ratio is greater than or equal to the threshold. */
    1115. 

  1115. 

  1116. /*  line 5:  TSTERR, LOGICAL */
    1117. 

  1117. /*           Flag indicating whether or not to test the error exits for */
    1118. 

  1118. /*           the LAPACK routines and driver routines. */
    1119. 

  1119. 

  1120. /*  line 6:  NEWSD, INTEGER */
    1121. 

  1121. /*           A code indicating how to set the random number seed. */
    1122. 

  1122. /*           = 0:  Set the seed to a default value before each run */
    1123. 

  1123. /*           = 1:  Initialize the seed to a default value only before the */
    1124. 

  1124. /*                 first run */
    1125. 

  1125. /*           = 2:  Like 1, but use the seed values on the next line */
    1126. 

  1126. 

  1127. /*  If line 6 was 2: */
    1128. 

  1128. 

  1129. /*  line 7: INTEGER array, dimension (4) */
    1130. 

  1130. /*           Four integer values for the random number seed. */
    1131. 

  1131. 

  1132. /*  If line 2 was 0: */
    1133. 

  1133. 

  1134. /*  line 7-EOF…
    1135.

Large files files are truncated, but you can click here to view the full file