PageRenderTime 67ms CodeModel.GetById 22ms RepoModel.GetById 1ms app.codeStats 0ms

/structural/supporting_matl_misc_and_old/support/LAPACK/lapack-3.4.0/TESTING/EIG/dchkee.f

https://bitbucket.org/slawton/windturbinemdo-stevebitb
FORTRAN Legacy | 2455 lines | 978 code | 0 blank | 1477 comment | 0 complexity | e9380d89d5c4fd0b795f8e6c759bb6ba MD5 | raw file
Possible License(s): GPL-2.0, BSD-3-Clause 
    

  1. *> \brief \b DCHKEE
    2. 

  2. *
    3. 

  3. *  =========== DOCUMENTATION ===========
    4. 

  4. *
    5. 

  5. * Online html documentation available at 
    6. 

  6. *            http://www.netlib.org/lapack/explore-html/ 
    7. 

  7. *
    8. 

  8. *  Definition:
    9. 

  9. *  ===========
    10. 

  10. *
    11. 

  11. *       PROGRAM DCHKEE
    12. 

  12. * 
    13. 

  13. *
    14. 

  14. *> \par Purpose:
    15. 

  15. *  =============
    16. 

  16. *>
    17. 

  17. *> \verbatim
    18. 

  18. *>
    19. 

  19. *> DCHKEE tests the DOUBLE PRECISION LAPACK subroutines for the matrix
    20. 

  20. *> eigenvalue problem.  The test paths in this version are
    21. 

  21. *>
    22. 

  22. *> NEP (Nonsymmetric Eigenvalue Problem):
    23. 

  23. *>     Test DGEHRD, DORGHR, DHSEQR, DTREVC, DHSEIN, and DORMHR
    24. 

  24. *>
    25. 

  25. *> SEP (Symmetric Eigenvalue Problem):
    26. 

  26. *>     Test DSYTRD, DORGTR, DSTEQR, DSTERF, DSTEIN, DSTEDC,
    27. 

  27. *>     and drivers DSYEV(X), DSBEV(X), DSPEV(X), DSTEV(X),
    28. 

  28. *>                 DSYEVD,   DSBEVD,   DSPEVD,   DSTEVD
    29. 

  29. *>
    30. 

  30. *> SVD (Singular Value Decomposition):
    31. 

  31. *>     Test DGEBRD, DORGBR, DBDSQR, DBDSDC
    32. 

  32. *>     and the drivers DGESVD, DGESDD
    33. 

  33. *>
    34. 

  34. *> DEV (Nonsymmetric Eigenvalue/eigenvector Driver):
    35. 

  35. *>     Test DGEEV
    36. 

  36. *>
    37. 

  37. *> DES (Nonsymmetric Schur form Driver):
    38. 

  38. *>     Test DGEES
    39. 

  39. *>
    40. 

  40. *> DVX (Nonsymmetric Eigenvalue/eigenvector Expert Driver):
    41. 

  41. *>     Test DGEEVX
    42. 

  42. *>
    43. 

  43. *> DSX (Nonsymmetric Schur form Expert Driver):
    44. 

  44. *>     Test DGEESX
    45. 

  45. *>
    46. 

  46. *> DGG (Generalized Nonsymmetric Eigenvalue Problem):
    47. 

  47. *>     Test DGGHRD, DGGBAL, DGGBAK, DHGEQZ, and DTGEVC
    48. 

  48. *>     and the driver routines DGEGS and DGEGV
    49. 

  49. *>
    50. 

  50. *> DGS (Generalized Nonsymmetric Schur form Driver):
    51. 

  51. *>     Test DGGES
    52. 

  52. *>
    53. 

  53. *> DGV (Generalized Nonsymmetric Eigenvalue/eigenvector Driver):
    54. 

  54. *>     Test DGGEV
    55. 

  55. *>
    56. 

  56. *> DGX (Generalized Nonsymmetric Schur form Expert Driver):
    57. 

  57. *>     Test DGGESX
    58. 

  58. *>
    59. 

  59. *> DXV (Generalized Nonsymmetric Eigenvalue/eigenvector Expert Driver):
    60. 

  60. *>     Test DGGEVX
    61. 

  61. *>
    62. 

  62. *> DSG (Symmetric Generalized Eigenvalue Problem):
    63. 

  63. *>     Test DSYGST, DSYGV, DSYGVD, DSYGVX, DSPGST, DSPGV, DSPGVD,
    64. 

  64. *>     DSPGVX, DSBGST, DSBGV, DSBGVD, and DSBGVX
    65. 

  65. *>
    66. 

  66. *> DSB (Symmetric Band Eigenvalue Problem):
    67. 

  67. *>     Test DSBTRD
    68. 

  68. *>
    69. 

  69. *> DBB (Band Singular Value Decomposition):
    70. 

  70. *>     Test DGBBRD
    71. 

  71. *>
    72. 

  72. *> DEC (Eigencondition estimation):
    73. 

  73. *>     Test DLALN2, DLASY2, DLAEQU, DLAEXC, DTRSYL, DTREXC, DTRSNA,
    74. 

  74. *>     DTRSEN, and DLAQTR
    75. 

  75. *>
    76. 

  76. *> DBL (Balancing a general matrix)
    77. 

  77. *>     Test DGEBAL
    78. 

  78. *>
    79. 

  79. *> DBK (Back transformation on a balanced matrix)
    80. 

  80. *>     Test DGEBAK
    81. 

  81. *>
    82. 

  82. *> DGL (Balancing a matrix pair)
    83. 

  83. *>     Test DGGBAL
    84. 

  84. *>
    85. 

  85. *> DGK (Back transformation on a matrix pair)
    86. 

  86. *>     Test DGGBAK
    87. 

  87. *>
    88. 

  88. *> GLM (Generalized Linear Regression Model):
    89. 

  89. *>     Tests DGGGLM
    90. 

  90. *>
    91. 

  91. *> GQR (Generalized QR and RQ factorizations):
    92. 

  92. *>     Tests DGGQRF and DGGRQF
    93. 

  93. *>
    94. 

  94. *> GSV (Generalized Singular Value Decomposition):
    95. 

  95. *>     Tests DGGSVD, DGGSVP, DTGSJA, DLAGS2, DLAPLL, and DLAPMT
    96. 

  96. *>
    97. 

  97. *> CSD (CS decomposition):
    98. 

  98. *>     Tests DORCSD
    99. 

  99. *>
    100. 

  100. *> LSE (Constrained Linear Least Squares):
    101. 

  101. *>     Tests DGGLSE
    102. 

  102. *>
    103. 

  103. *> Each test path has a different set of inputs, but the data sets for
    104. 

  104. *> the driver routines xEV, xES, xVX, and xSX can be concatenated in a
    105. 

  105. *> single input file.  The first line of input should contain one of the
    106. 

  106. *> 3-character path names in columns 1-3.  The number of remaining lines
    107. 

  107. *> depends on what is found on the first line.
    108. 

  108. *>
    109. 

  109. *> The number of matrix types used in testing is often controllable from
    110. 

  110. *> the input file.  The number of matrix types for each path, and the
    111. 

  111. *> test routine that describes them, is as follows:
    112. 

  112. *>
    113. 

  113. *> Path name(s)  Types    Test routine
    114. 

  114. *>
    115. 

  115. *> DHS or NEP      21     DCHKHS
    116. 

  116. *> DST or SEP      21     DCHKST (routines)
    117. 

  117. *>                 18     DDRVST (drivers)
    118. 

  118. *> DBD or SVD      16     DCHKBD (routines)
    119. 

  119. *>                  5     DDRVBD (drivers)
    120. 

  120. *> DEV             21     DDRVEV
    121. 

  121. *> DES             21     DDRVES
    122. 

  122. *> DVX             21     DDRVVX
    123. 

  123. *> DSX             21     DDRVSX
    124. 

  124. *> DGG             26     DCHKGG (routines)
    125. 

  125. *>                 26     DDRVGG (drivers)
    126. 

  126. *> DGS             26     DDRGES
    127. 

  127. *> DGX              5     DDRGSX
    128. 

  128. *> DGV             26     DDRGEV
    129. 

  129. *> DXV              2     DDRGVX
    130. 

  130. *> DSG             21     DDRVSG
    131. 

  131. *> DSB             15     DCHKSB
    132. 

  132. *> DBB             15     DCHKBB
    133. 

  133. *> DEC              -     DCHKEC
    134. 

  134. *> DBL              -     DCHKBL
    135. 

  135. *> DBK              -     DCHKBK
    136. 

  136. *> DGL              -     DCHKGL
    137. 

  137. *> DGK              -     DCHKGK
    138. 

  138. *> GLM              8     DCKGLM
    139. 

  139. *> GQR              8     DCKGQR
    140. 

  140. *> GSV              8     DCKGSV
    141. 

  141. *> CSD              3     DCKCSD
    142. 

  142. *> LSE              8     DCKLSE
    143. 

  143. *>
    144. 

  144. *>-----------------------------------------------------------------------
    145. 

  145. *>
    146. 

  146. *> NEP input file:
    147. 

  147. *>
    148. 

  148. *> line 2:  NN, INTEGER
    149. 

  149. *>          Number of values of N.
    150. 

  150. *>
    151. 

  151. *> line 3:  NVAL, INTEGER array, dimension (NN)
    152. 

  152. *>          The values for the matrix dimension N.
    153. 

  153. *>
    154. 

  154. *> line 4:  NPARMS, INTEGER
    155. 

  155. *>          Number of values of the parameters NB, NBMIN, NX, NS, and
    156. 

  156. *>          MAXB.
    157. 

  157. *>
    158. 

  158. *> line 5:  NBVAL, INTEGER array, dimension (NPARMS)
    159. 

  159. *>          The values for the blocksize NB.
    160. 

  160. *>
    161. 

  161. *> line 6:  NBMIN, INTEGER array, dimension (NPARMS)
    162. 

  162. *>          The values for the minimum blocksize NBMIN.
    163. 

  163. *>
    164. 

  164. *> line 7:  NXVAL, INTEGER array, dimension (NPARMS)
    165. 

  165. *>          The values for the crossover point NX.
    166. 

  166. *>
    167. 

  167. *> line 8:  INMIN, INTEGER array, dimension (NPARMS)
    168. 

  168. *>          LAHQR vs TTQRE crossover point, >= 11
    169. 

  169. *>
    170. 

  170. *> line 9:  INWIN, INTEGER array, dimension (NPARMS)
    171. 

  171. *>          recommended deflation window size
    172. 

  172. *>
    173. 

  173. *> line 10: INIBL, INTEGER array, dimension (NPARMS)
    174. 

  174. *>          nibble crossover point
    175. 

  175. *>
    176. 

  176. *> line 11: ISHFTS, INTEGER array, dimension (NPARMS)
    177. 

  177. *>          number of simultaneous shifts)
    178. 

  178. *>
    179. 

  179. *> line 12: IACC22, INTEGER array, dimension (NPARMS)
    180. 

  180. *>          select structured matrix multiply: 0, 1 or 2)
    181. 

  181. *>
    182. 

  182. *> line 13: THRESH
    183. 

  183. *>          Threshold value for the test ratios.  Information will be
    184. 

  184. *>          printed about each test for which the test ratio is greater
    185. 

  185. *>          than or equal to the threshold.  To have all of the test
    186. 

  186. *>          ratios printed, use THRESH = 0.0 .
    187. 

  187. *>
    188. 

  188. *> line 14: NEWSD, INTEGER
    189. 

  189. *>          A code indicating how to set the random number seed.
    190. 

  190. *>          = 0:  Set the seed to a default value before each run
    191. 

  191. *>          = 1:  Initialize the seed to a default value only before the
    192. 

  192. *>                first run
    193. 

  193. *>          = 2:  Like 1, but use the seed values on the next line
    194. 

  194. *>
    195. 

  195. *> If line 14 was 2:
    196. 

  196. *>
    197. 

  197. *> line 15: INTEGER array, dimension (4)
    198. 

  198. *>          Four integer values for the random number seed.
    199. 

  199. *>
    200. 

  200. *> lines 15-EOF:  The remaining lines occur in sets of 1 or 2 and allow
    201. 

  201. *>          the user to specify the matrix types.  Each line contains
    202. 

  202. *>          a 3-character path name in columns 1-3, and the number
    203. 

  203. *>          of matrix types must be the first nonblank item in columns
    204. 

  204. *>          4-80.  If the number of matrix types is at least 1 but is
    205. 

  205. *>          less than the maximum number of possible types, a second
    206. 

  206. *>          line will be read to get the numbers of the matrix types to
    207. 

  207. *>          be used.  For example,
    208. 

  208. *> NEP 21
    209. 

  209. *>          requests all of the matrix types for the nonsymmetric
    210. 

  210. *>          eigenvalue problem, while
    211. 

  211. *> NEP  4
    212. 

  212. *> 9 10 11 12
    213. 

  213. *>          requests only matrices of type 9, 10, 11, and 12.
    214. 

  214. *>
    215. 

  215. *>          The valid 3-character path names are 'NEP' or 'SHS' for the
    216. 

  216. *>          nonsymmetric eigenvalue routines.
    217. 

  217. *>
    218. 

  218. *>-----------------------------------------------------------------------
    219. 

  219. *>
    220. 

  220. *> SEP or DSG input file:
    221. 

  221. *>
    222. 

  222. *> line 2:  NN, INTEGER
    223. 

  223. *>          Number of values of N.
    224. 

  224. *>
    225. 

  225. *> line 3:  NVAL, INTEGER array, dimension (NN)
    226. 

  226. *>          The values for the matrix dimension N.
    227. 

  227. *>
    228. 

  228. *> line 4:  NPARMS, INTEGER
    229. 

  229. *>          Number of values of the parameters NB, NBMIN, and NX.
    230. 

  230. *>
    231. 

  231. *> line 5:  NBVAL, INTEGER array, dimension (NPARMS)
    232. 

  232. *>          The values for the blocksize NB.
    233. 

  233. *>
    234. 

  234. *> line 6:  NBMIN, INTEGER array, dimension (NPARMS)
    235. 

  235. *>          The values for the minimum blocksize NBMIN.
    236. 

  236. *>
    237. 

  237. *> line 7:  NXVAL, INTEGER array, dimension (NPARMS)
    238. 

  238. *>          The values for the crossover point NX.
    239. 

  239. *>
    240. 

  240. *> line 8:  THRESH
    241. 

  241. *>          Threshold value for the test ratios.  Information will be
    242. 

  242. *>          printed about each test for which the test ratio is greater
    243. 

  243. *>          than or equal to the threshold.
    244. 

  244. *>
    245. 

  245. *> line 9:  TSTCHK, LOGICAL
    246. 

  246. *>          Flag indicating whether or not to test the LAPACK routines.
    247. 

  247. *>
    248. 

  248. *> line 10: TSTDRV, LOGICAL
    249. 

  249. *>          Flag indicating whether or not to test the driver routines.
    250. 

  250. *>
    251. 

  251. *> line 11: TSTERR, LOGICAL
    252. 

  252. *>          Flag indicating whether or not to test the error exits for
    253. 

  253. *>          the LAPACK routines and driver routines.
    254. 

  254. *>
    255. 

  255. *> line 12: NEWSD, INTEGER
    256. 

  256. *>          A code indicating how to set the random number seed.
    257. 

  257. *>          = 0:  Set the seed to a default value before each run
    258. 

  258. *>          = 1:  Initialize the seed to a default value only before the
    259. 

  259. *>                first run
    260. 

  260. *>          = 2:  Like 1, but use the seed values on the next line
    261. 

  261. *>
    262. 

  262. *> If line 12 was 2:
    263. 

  263. *>
    264. 

  264. *> line 13: INTEGER array, dimension (4)
    265. 

  265. *>          Four integer values for the random number seed.
    266. 

  266. *>
    267. 

  267. *> lines 13-EOF:  Lines specifying matrix types, as for NEP.
    268. 

  268. *>          The 3-character path names are 'SEP' or 'SST' for the
    269. 

  269. *>          symmetric eigenvalue routines and driver routines, and
    270. 

  270. *>          'DSG' for the routines for the symmetric generalized
    271. 

  271. *>          eigenvalue problem.
    272. 

  272. *>
    273. 

  273. *>-----------------------------------------------------------------------
    274. 

  274. *>
    275. 

  275. *> SVD input file:
    276. 

  276. *>
    277. 

  277. *> line 2:  NN, INTEGER
    278. 

  278. *>          Number of values of M and N.
    279. 

  279. *>
    280. 

  280. *> line 3:  MVAL, INTEGER array, dimension (NN)
    281. 

  281. *>          The values for the matrix row dimension M.
    282. 

  282. *>
    283. 

  283. *> line 4:  NVAL, INTEGER array, dimension (NN)
    284. 

  284. *>          The values for the matrix column dimension N.
    285. 

  285. *>
    286. 

  286. *> line 5:  NPARMS, INTEGER
    287. 

  287. *>          Number of values of the parameter NB, NBMIN, NX, and NRHS.
    288. 

  288. *>
    289. 

  289. *> line 6:  NBVAL, INTEGER array, dimension (NPARMS)
    290. 

  290. *>          The values for the blocksize NB.
    291. 

  291. *>
    292. 

  292. *> line 7:  NBMIN, INTEGER array, dimension (NPARMS)
    293. 

  293. *>          The values for the minimum blocksize NBMIN.
    294. 

  294. *>
    295. 

  295. *> line 8:  NXVAL, INTEGER array, dimension (NPARMS)
    296. 

  296. *>          The values for the crossover point NX.
    297. 

  297. *>
    298. 

  298. *> line 9:  NSVAL, INTEGER array, dimension (NPARMS)
    299. 

  299. *>          The values for the number of right hand sides NRHS.
    300. 

  300. *>
    301. 

  301. *> line 10: THRESH
    302. 

  302. *>          Threshold value for the test ratios.  Information will be
    303. 

  303. *>          printed about each test for which the test ratio is greater
    304. 

  304. *>          than or equal to the threshold.
    305. 

  305. *>
    306. 

  306. *> line 11: TSTCHK, LOGICAL
    307. 

  307. *>          Flag indicating whether or not to test the LAPACK routines.
    308. 

  308. *>
    309. 

  309. *> line 12: TSTDRV, LOGICAL
    310. 

  310. *>          Flag indicating whether or not to test the driver routines.
    311. 

  311. *>
    312. 

  312. *> line 13: TSTERR, LOGICAL
    313. 

  313. *>          Flag indicating whether or not to test the error exits for
    314. 

  314. *>          the LAPACK routines and driver routines.
    315. 

  315. *>
    316. 

  316. *> line 14: NEWSD, INTEGER
    317. 

  317. *>          A code indicating how to set the random number seed.
    318. 

  318. *>          = 0:  Set the seed to a default value before each run
    319. 

  319. *>          = 1:  Initialize the seed to a default value only before the
    320. 

  320. *>                first run
    321. 

  321. *>          = 2:  Like 1, but use the seed values on the next line
    322. 

  322. *>
    323. 

  323. *> If line 14 was 2:
    324. 

  324. *>
    325. 

  325. *> line 15: INTEGER array, dimension (4)
    326. 

  326. *>          Four integer values for the random number seed.
    327. 

  327. *>
    328. 

  328. *> lines 15-EOF:  Lines specifying matrix types, as for NEP.
    329. 

  329. *>          The 3-character path names are 'SVD' or 'SBD' for both the
    330. 

  330. *>          SVD routines and the SVD driver routines.
    331. 

  331. *>
    332. 

  332. *>-----------------------------------------------------------------------
    333. 

  333. *>
    334. 

  334. *> DEV and DES data files:
    335. 

  335. *>
    336. 

  336. *> line 1:  'DEV' or 'DES' in columns 1 to 3.
    337. 

  337. *>
    338. 

  338. *> line 2:  NSIZES, INTEGER
    339. 

  339. *>          Number of sizes of matrices to use. Should be at least 0
    340. 

  340. *>          and at most 20. If NSIZES = 0, no testing is done
    341. 

  341. *>          (although the remaining  3 lines are still read).
    342. 

  342. *>
    343. 

  343. *> line 3:  NN, INTEGER array, dimension(NSIZES)
    344. 

  344. *>          Dimensions of matrices to be tested.
    345. 

  345. *>
    346. 

  346. *> line 4:  NB, NBMIN, NX, NS, NBCOL, INTEGERs
    347. 

  347. *>          These integer parameters determine how blocking is done
    348. 

  348. *>          (see ILAENV for details)
    349. 

  349. *>          NB     : block size
    350. 

  350. *>          NBMIN  : minimum block size
    351. 

  351. *>          NX     : minimum dimension for blocking
    352. 

  352. *>          NS     : number of shifts in xHSEQR
    353. 

  353. *>          NBCOL  : minimum column dimension for blocking
    354. 

  354. *>
    355. 

  355. *> line 5:  THRESH, REAL
    356. 

  356. *>          The test threshold against which computed residuals are
    357. 

  357. *>          compared. Should generally be in the range from 10. to 20.
    358. 

  358. *>          If it is 0., all test case data will be printed.
    359. 

  359. *>
    360. 

  360. *> line 6:  TSTERR, LOGICAL
    361. 

  361. *>          Flag indicating whether or not to test the error exits.
    362. 

  362. *>
    363. 

  363. *> line 7:  NEWSD, INTEGER
    364. 

  364. *>          A code indicating how to set the random number seed.
    365. 

  365. *>          = 0:  Set the seed to a default value before each run
    366. 

  366. *>          = 1:  Initialize the seed to a default value only before the
    367. 

  367. *>                first run
    368. 

  368. *>          = 2:  Like 1, but use the seed values on the next line
    369. 

  369. *>
    370. 

  370. *> If line 7 was 2:
    371. 

  371. *>
    372. 

  372. *> line 8:  INTEGER array, dimension (4)
    373. 

  373. *>          Four integer values for the random number seed.
    374. 

  374. *>
    375. 

  375. *> lines 9 and following:  Lines specifying matrix types, as for NEP.
    376. 

  376. *>          The 3-character path name is 'DEV' to test SGEEV, or
    377. 

  377. *>          'DES' to test SGEES.
    378. 

  378. *>
    379. 

  379. *>-----------------------------------------------------------------------
    380. 

  380. *>
    381. 

  381. *> The DVX data has two parts. The first part is identical to DEV,
    382. 

  382. *> and the second part consists of test matrices with precomputed
    383. 

  383. *> solutions.
    384. 

  384. *>
    385. 

  385. *> line 1:  'DVX' in columns 1-3.
    386. 

  386. *>
    387. 

  387. *> line 2:  NSIZES, INTEGER
    388. 

  388. *>          If NSIZES = 0, no testing of randomly generated examples
    389. 

  389. *>          is done, but any precomputed examples are tested.
    390. 

  390. *>
    391. 

  391. *> line 3:  NN, INTEGER array, dimension(NSIZES)
    392. 

  392. *>
    393. 

  393. *> line 4:  NB, NBMIN, NX, NS, NBCOL, INTEGERs
    394. 

  394. *>
    395. 

  395. *> line 5:  THRESH, REAL
    396. 

  396. *>
    397. 

  397. *> line 6:  TSTERR, LOGICAL
    398. 

  398. *>
    399. 

  399. *> line 7:  NEWSD, INTEGER
    400. 

  400. *>
    401. 

  401. *> If line 7 was 2:
    402. 

  402. *>
    403. 

  403. *> line 8:  INTEGER array, dimension (4)
    404. 

  404. *>
    405. 

  405. *> lines 9 and following: The first line contains 'DVX' in columns 1-3
    406. 

  406. *>          followed by the number of matrix types, possibly with
    407. 

  407. *>          a second line to specify certain matrix types.
    408. 

  408. *>          If the number of matrix types = 0, no testing of randomly
    409. 

  409. *>          generated examples is done, but any precomputed examples
    410. 

  410. *>          are tested.
    411. 

  411. *>
    412. 

  412. *> remaining lines : Each matrix is stored on 1+2*N lines, where N is
    413. 

  413. *>          its dimension. The first line contains the dimension (a
    414. 

  414. *>          single integer). The next N lines contain the matrix, one
    415. 

  415. *>          row per line. The last N lines correspond to each
    416. 

  416. *>          eigenvalue. Each of these last N lines contains 4 real
    417. 

  417. *>          values: the real part of the eigenvalue, the imaginary
    418. 

  418. *>          part of the eigenvalue, the reciprocal condition number of
    419. 

  419. *>          the eigenvalues, and the reciprocal condition number of the
    420. 

  420. *>          eigenvector.  The end of data is indicated by dimension N=0.
    421. 

  421. *>          Even if no data is to be tested, there must be at least one
    422. 

  422. *>          line containing N=0.
    423. 

  423. *>
    424. 

  424. *>-----------------------------------------------------------------------
    425. 

  425. *>
    426. 

  426. *> The DSX data is like DVX. The first part is identical to DEV, and the
    427. 

  427. *> second part consists of test matrices with precomputed solutions.
    428. 

  428. *>
    429. 

  429. *> line 1:  'DSX' in columns 1-3.
    430. 

  430. *>
    431. 

  431. *> line 2:  NSIZES, INTEGER
    432. 

  432. *>          If NSIZES = 0, no testing of randomly generated examples
    433. 

  433. *>          is done, but any precomputed examples are tested.
    434. 

  434. *>
    435. 

  435. *> line 3:  NN, INTEGER array, dimension(NSIZES)
    436. 

  436. *>
    437. 

  437. *> line 4:  NB, NBMIN, NX, NS, NBCOL, INTEGERs
    438. 

  438. *>
    439. 

  439. *> line 5:  THRESH, REAL
    440. 

  440. *>
    441. 

  441. *> line 6:  TSTERR, LOGICAL
    442. 

  442. *>
    443. 

  443. *> line 7:  NEWSD, INTEGER
    444. 

  444. *>
    445. 

  445. *> If line 7 was 2:
    446. 

  446. *>
    447. 

  447. *> line 8:  INTEGER array, dimension (4)
    448. 

  448. *>
    449. 

  449. *> lines 9 and following: The first line contains 'DSX' in columns 1-3
    450. 

  450. *>          followed by the number of matrix types, possibly with
    451. 

  451. *>          a second line to specify certain matrix types.
    452. 

  452. *>          If the number of matrix types = 0, no testing of randomly
    453. 

  453. *>          generated examples is done, but any precomputed examples
    454. 

  454. *>          are tested.
    455. 

  455. *>
    456. 

  456. *> remaining lines : Each matrix is stored on 3+N lines, where N is its
    457. 

  457. *>          dimension. The first line contains the dimension N and the
    458. 

  458. *>          dimension M of an invariant subspace. The second line
    459. 

  459. *>          contains M integers, identifying the eigenvalues in the
    460. 

  460. *>          invariant subspace (by their position in a list of
    461. 

  461. *>          eigenvalues ordered by increasing real part). The next N
    462. 

  462. *>          lines contain the matrix. The last line contains the
    463. 

  463. *>          reciprocal condition number for the average of the selected
    464. 

  464. *>          eigenvalues, and the reciprocal condition number for the
    465. 

  465. *>          corresponding right invariant subspace. The end of data is
    466. 

  466. *>          indicated by a line containing N=0 and M=0. Even if no data
    467. 

  467. *>          is to be tested, there must be at least one line containing
    468. 

  468. *>          N=0 and M=0.
    469. 

  469. *>
    470. 

  470. *>-----------------------------------------------------------------------
    471. 

  471. *>
    472. 

  472. *> DGG input file:
    473. 

  473. *>
    474. 

  474. *> line 2:  NN, INTEGER
    475. 

  475. *>          Number of values of N.
    476. 

  476. *>
    477. 

  477. *> line 3:  NVAL, INTEGER array, dimension (NN)
    478. 

  478. *>          The values for the matrix dimension N.
    479. 

  479. *>
    480. 

  480. *> line 4:  NPARMS, INTEGER
    481. 

  481. *>          Number of values of the parameters NB, NBMIN, NS, MAXB, and
    482. 

  482. *>          NBCOL.
    483. 

  483. *>
    484. 

  484. *> line 5:  NBVAL, INTEGER array, dimension (NPARMS)
    485. 

  485. *>          The values for the blocksize NB.
    486. 

  486. *>
    487. 

  487. *> line 6:  NBMIN, INTEGER array, dimension (NPARMS)
    488. 

  488. *>          The values for NBMIN, the minimum row dimension for blocks.
    489. 

  489. *>
    490. 

  490. *> line 7:  NSVAL, INTEGER array, dimension (NPARMS)
    491. 

  491. *>          The values for the number of shifts.
    492. 

  492. *>
    493. 

  493. *> line 8:  MXBVAL, INTEGER array, dimension (NPARMS)
    494. 

  494. *>          The values for MAXB, used in determining minimum blocksize.
    495. 

  495. *>
    496. 

  496. *> line 9:  NBCOL, INTEGER array, dimension (NPARMS)
    497. 

  497. *>          The values for NBCOL, the minimum column dimension for
    498. 

  498. *>          blocks.
    499. 

  499. *>
    500. 

  500. *> line 10: THRESH
    501. 

  501. *>          Threshold value for the test ratios.  Information will be
    502. 

  502. *>          printed about each test for which the test ratio is greater
    503. 

  503. *>          than or equal to the threshold.
    504. 

  504. *>
    505. 

  505. *> line 11: TSTCHK, LOGICAL
    506. 

  506. *>          Flag indicating whether or not to test the LAPACK routines.
    507. 

  507. *>
    508. 

  508. *> line 12: TSTDRV, LOGICAL
    509. 

  509. *>          Flag indicating whether or not to test the driver routines.
    510. 

  510. *>
    511. 

  511. *> line 13: TSTERR, LOGICAL
    512. 

  512. *>          Flag indicating whether or not to test the error exits for
    513. 

  513. *>          the LAPACK routines and driver routines.
    514. 

  514. *>
    515. 

  515. *> line 14: NEWSD, INTEGER
    516. 

  516. *>          A code indicating how to set the random number seed.
    517. 

  517. *>          = 0:  Set the seed to a default value before each run
    518. 

  518. *>          = 1:  Initialize the seed to a default value only before the
    519. 

  519. *>                first run
    520. 

  520. *>          = 2:  Like 1, but use the seed values on the next line
    521. 

  521. *>
    522. 

  522. *> If line 14 was 2:
    523. 

  523. *>
    524. 

  524. *> line 15: INTEGER array, dimension (4)
    525. 

  525. *>          Four integer values for the random number seed.
    526. 

  526. *>
    527. 

  527. *> lines 15-EOF:  Lines specifying matrix types, as for NEP.
    528. 

  528. *>          The 3-character path name is 'DGG' for the generalized
    529. 

  529. *>          eigenvalue problem routines and driver routines.
    530. 

  530. *>
    531. 

  531. *>-----------------------------------------------------------------------
    532. 

  532. *>
    533. 

  533. *> DGS and DGV input files:
    534. 

  534. *>
    535. 

  535. *> line 1:  'DGS' or 'DGV' in columns 1 to 3.
    536. 

  536. *>
    537. 

  537. *> line 2:  NN, INTEGER
    538. 

  538. *>          Number of values of N.
    539. 

  539. *>
    540. 

  540. *> line 3:  NVAL, INTEGER array, dimension(NN)
    541. 

  541. *>          Dimensions of matrices to be tested.
    542. 

  542. *>
    543. 

  543. *> line 4:  NB, NBMIN, NX, NS, NBCOL, INTEGERs
    544. 

  544. *>          These integer parameters determine how blocking is done
    545. 

  545. *>          (see ILAENV for details)
    546. 

  546. *>          NB     : block size
    547. 

  547. *>          NBMIN  : minimum block size
    548. 

  548. *>          NX     : minimum dimension for blocking
    549. 

  549. *>          NS     : number of shifts in xHGEQR
    550. 

  550. *>          NBCOL  : minimum column dimension for blocking
    551. 

  551. *>
    552. 

  552. *> line 5:  THRESH, REAL
    553. 

  553. *>          The test threshold against which computed residuals are
    554. 

  554. *>          compared. Should generally be in the range from 10. to 20.
    555. 

  555. *>          If it is 0., all test case data will be printed.
    556. 

  556. *>
    557. 

  557. *> line 6:  TSTERR, LOGICAL
    558. 

  558. *>          Flag indicating whether or not to test the error exits.
    559. 

  559. *>
    560. 

  560. *> line 7:  NEWSD, INTEGER
    561. 

  561. *>          A code indicating how to set the random number seed.
    562. 

  562. *>          = 0:  Set the seed to a default value before each run
    563. 

  563. *>          = 1:  Initialize the seed to a default value only before the
    564. 

  564. *>                first run
    565. 

  565. *>          = 2:  Like 1, but use the seed values on the next line
    566. 

  566. *>
    567. 

  567. *> If line 17 was 2:
    568. 

  568. *>
    569. 

  569. *> line 7:  INTEGER array, dimension (4)
    570. 

  570. *>          Four integer values for the random number seed.
    571. 

  571. *>
    572. 

  572. *> lines 7-EOF:  Lines specifying matrix types, as for NEP.
    573. 

  573. *>          The 3-character path name is 'DGS' for the generalized
    574. 

  574. *>          eigenvalue problem routines and driver routines.
    575. 

  575. *>
    576. 

  576. *>-----------------------------------------------------------------------
    577. 

  577. *>
    578. 

  578. *> DXV input files:
    579. 

  579. *>
    580. 

  580. *> line 1:  'DXV' in columns 1 to 3.
    581. 

  581. *>
    582. 

  582. *> line 2:  N, INTEGER
    583. 

  583. *>          Value of N.
    584. 

  584. *>
    585. 

  585. *> line 3:  NB, NBMIN, NX, NS, NBCOL, INTEGERs
    586. 

  586. *>          These integer parameters determine how blocking is done
    587. 

  587. *>          (see ILAENV for details)
    588. 

  588. *>          NB     : block size
    589. 

  589. *>          NBMIN  : minimum block size
    590. 

  590. *>          NX     : minimum dimension for blocking
    591. 

  591. *>          NS     : number of shifts in xHGEQR
    592. 

  592. *>          NBCOL  : minimum column dimension for blocking
    593. 

  593. *>
    594. 

  594. *> line 4:  THRESH, REAL
    595. 

  595. *>          The test threshold against which computed residuals are
    596. 

  596. *>          compared. Should generally be in the range from 10. to 20.
    597. 

  597. *>          Information will be printed about each test for which the
    598. 

  598. *>          test ratio is greater than or equal to the threshold.
    599. 

  599. *>
    600. 

  600. *> line 5:  TSTERR, LOGICAL
    601. 

  601. *>          Flag indicating whether or not to test the error exits for
    602. 

  602. *>          the LAPACK routines and driver routines.
    603. 

  603. *>
    604. 

  604. *> line 6:  NEWSD, INTEGER
    605. 

  605. *>          A code indicating how to set the random number seed.
    606. 

  606. *>          = 0:  Set the seed to a default value before each run
    607. 

  607. *>          = 1:  Initialize the seed to a default value only before the
    608. 

  608. *>                first run
    609. 

  609. *>          = 2:  Like 1, but use the seed values on the next line
    610. 

  610. *>
    611. 

  611. *> If line 6 was 2:
    612. 

  612. *>
    613. 

  613. *> line 7: INTEGER array, dimension (4)
    614. 

  614. *>          Four integer values for the random number seed.
    615. 

  615. *>
    616. 

  616. *> If line 2 was 0:
    617. 

  617. *>
    618. 

  618. *> line 7-EOF: Precomputed examples are tested.
    619. 

  619. *>
    620. 

  620. *> remaining lines : Each example is stored on 3+2*N lines, where N is
    621. 

  621. *>          its dimension. The first line contains the dimension (a
    622. 

  622. *>          single integer). The next N lines contain the matrix A, one
    623. 

  623. *>          row per line. The next N lines contain the matrix B.  The
    624. 

  624. *>          next line contains the reciprocals of the eigenvalue
    625. 

  625. *>          condition numbers.  The last line contains the reciprocals of
    626. 

  626. *>          the eigenvector condition numbers.  The end of data is
    627. 

  627. *>          indicated by dimension N=0.  Even if no data is to be tested,
    628. 

  628. *>          there must be at least one line containing N=0.
    629. 

  629. *>
    630. 

  630. *>-----------------------------------------------------------------------
    631. 

  631. *>
    632. 

  632. *> DGX input files:
    633. 

  633. *>
    634. 

  634. *> line 1:  'DGX' in columns 1 to 3.
    635. 

  635. *>
    636. 

  636. *> line 2:  N, INTEGER
    637. 

  637. *>          Value of N.
    638. 

  638. *>
    639. 

  639. *> line 3:  NB, NBMIN, NX, NS, NBCOL, INTEGERs
    640. 

  640. *>          These integer parameters determine how blocking is done
    641. 

  641. *>          (see ILAENV for details)
    642. 

  642. *>          NB     : block size
    643. 

  643. *>          NBMIN  : minimum block size
    644. 

  644. *>          NX     : minimum dimension for blocking
    645. 

  645. *>          NS     : number of shifts in xHGEQR
    646. 

  646. *>          NBCOL  : minimum column dimension for blocking
    647. 

  647. *>
    648. 

  648. *> line 4:  THRESH, REAL
    649. 

  649. *>          The test threshold against which computed residuals are
    650. 

  650. *>          compared. Should generally be in the range from 10. to 20.
    651. 

  651. *>          Information will be printed about each test for which the
    652. 

  652. *>          test ratio is greater than or equal to the threshold.
    653. 

  653. *>
    654. 

  654. *> line 5:  TSTERR, LOGICAL
    655. 

  655. *>          Flag indicating whether or not to test the error exits for
    656. 

  656. *>          the LAPACK routines and driver routines.
    657. 

  657. *>
    658. 

  658. *> line 6:  NEWSD, INTEGER
    659. 

  659. *>          A code indicating how to set the random number seed.
    660. 

  660. *>          = 0:  Set the seed to a default value before each run
    661. 

  661. *>          = 1:  Initialize the seed to a default value only before the
    662. 

  662. *>                first run
    663. 

  663. *>          = 2:  Like 1, but use the seed values on the next line
    664. 

  664. *>
    665. 

  665. *> If line 6 was 2:
    666. 

  666. *>
    667. 

  667. *> line 7: INTEGER array, dimension (4)
    668. 

  668. *>          Four integer values for the random number seed.
    669. 

  669. *>
    670. 

  670. *> If line 2 was 0:
    671. 

  671. *>
    672. 

  672. *> line 7-EOF: Precomputed examples are tested.
    673. 

  673. *>
    674. 

  674. *> remaining lines : Each example is stored on 3+2*N lines, where N is
    675. 

  675. *>          its dimension. The first line contains the dimension (a
    676. 

  676. *>          single integer).  The next line contains an integer k such
    677. 

  677. *>          that only the last k eigenvalues will be selected and appear
    678. 

  678. *>          in the leading diagonal blocks of $A$ and $B$. The next N
    679. 

  679. *>          lines contain the matrix A, one row per line.  The next N
    680. 

  680. *>          lines contain the matrix B.  The last line contains the
    681. 

  681. *>          reciprocal of the eigenvalue cluster condition number and the
    682. 

  682. *>          reciprocal of the deflating subspace (associated with the
    683. 

  683. *>          selected eigencluster) condition number.  The end of data is
    684. 

  684. *>          indicated by dimension N=0.  Even if no data is to be tested,
    685. 

  685. *>          there must be at least one line containing N=0.
    686. 

  686. *>
    687. 

  687. *>-----------------------------------------------------------------------
    688. 

  688. *>
    689. 

  689. *> DSB input file:
    690. 

  690. *>
    691. 

  691. *> line 2:  NN, INTEGER
    692. 

  692. *>          Number of values of N.
    693. 

  693. *>
    694. 

  694. *> line 3:  NVAL, INTEGER array, dimension (NN)
    695. 

  695. *>          The values for the matrix dimension N.
    696. 

  696. *>
    697. 

  697. *> line 4:  NK, INTEGER
    698. 

  698. *>          Number of values of K.
    699. 

  699. *>
    700. 

  700. *> line 5:  KVAL, INTEGER array, dimension (NK)
    701. 

  701. *>          The values for the matrix dimension K.
    702. 

  702. *>
    703. 

  703. *> line 6:  THRESH
    704. 

  704. *>          Threshold value for the test ratios.  Information will be
    705. 

  705. *>          printed about each test for which the test ratio is greater
    706. 

  706. *>          than or equal to the threshold.
    707. 

  707. *>
    708. 

  708. *> line 7:  NEWSD, INTEGER
    709. 

  709. *>          A code indicating how to set the random number seed.
    710. 

  710. *>          = 0:  Set the seed to a default value before each run
    711. 

  711. *>          = 1:  Initialize the seed to a default value only before the
    712. 

  712. *>                first run
    713. 

  713. *>          = 2:  Like 1, but use the seed values on the next line
    714. 

  714. *>
    715. 

  715. *> If line 7 was 2:
    716. 

  716. *>
    717. 

  717. *> line 8:  INTEGER array, dimension (4)
    718. 

  718. *>          Four integer values for the random number seed.
    719. 

  719. *>
    720. 

  720. *> lines 8-EOF:  Lines specifying matrix types, as for NEP.
    721. 

  721. *>          The 3-character path name is 'DSB'.
    722. 

  722. *>
    723. 

  723. *>-----------------------------------------------------------------------
    724. 

  724. *>
    725. 

  725. *> DBB input file:
    726. 

  726. *>
    727. 

  727. *> line 2:  NN, INTEGER
    728. 

  728. *>          Number of values of M and N.
    729. 

  729. *>
    730. 

  730. *> line 3:  MVAL, INTEGER array, dimension (NN)
    731. 

  731. *>          The values for the matrix row dimension M.
    732. 

  732. *>
    733. 

  733. *> line 4:  NVAL, INTEGER array, dimension (NN)
    734. 

  734. *>          The values for the matrix column dimension N.
    735. 

  735. *>
    736. 

  736. *> line 4:  NK, INTEGER
    737. 

  737. *>          Number of values of K.
    738. 

  738. *>
    739. 

  739. *> line 5:  KVAL, INTEGER array, dimension (NK)
    740. 

  740. *>          The values for the matrix bandwidth K.
    741. 

  741. *>
    742. 

  742. *> line 6:  NPARMS, INTEGER
    743. 

  743. *>          Number of values of the parameter NRHS
    744. 

  744. *>
    745. 

  745. *> line 7:  NSVAL, INTEGER array, dimension (NPARMS)
    746. 

  746. *>          The values for the number of right hand sides NRHS.
    747. 

  747. *>
    748. 

  748. *> line 8:  THRESH
    749. 

  749. *>          Threshold value for the test ratios.  Information will be
    750. 

  750. *>          printed about each test for which the test ratio is greater
    751. 

  751. *>          than or equal to the threshold.
    752. 

  752. *>
    753. 

  753. *> line 9:  NEWSD, INTEGER
    754. 

  754. *>          A code indicating how to set the random number seed.
    755. 

  755. *>          = 0:  Set the seed to a default value before each run
    756. 

  756. *>          = 1:  Initialize the seed to a default value only before the
    757. 

  757. *>                first run
    758. 

  758. *>          = 2:  Like 1, but use the seed values on the next line
    759. 

  759. *>
    760. 

  760. *> If line 9 was 2:
    761. 

  761. *>
    762. 

  762. *> line 10: INTEGER array, dimension (4)
    763. 

  763. *>          Four integer values for the random number seed.
    764. 

  764. *>
    765. 

  765. *> lines 10-EOF:  Lines specifying matrix types, as for SVD.
    766. 

  766. *>          The 3-character path name is 'DBB'.
    767. 

  767. *>
    768. 

  768. *>-----------------------------------------------------------------------
    769. 

  769. *>
    770. 

  770. *> DEC input file:
    771. 

  771. *>
    772. 

  772. *> line  2: THRESH, REAL
    773. 

  773. *>          Threshold value for the test ratios.  Information will be
    774. 

  774. *>          printed about each test for which the test ratio is greater
    775. 

  775. *>          than or equal to the threshold.
    776. 

  776. *>
    777. 

  777. *> lines  3-EOF:
    778. 

  778. *>
    779. 

  779. *> Input for testing the eigencondition routines consists of a set of
    780. 

  780. *> specially constructed test cases and their solutions.  The data
    781. 

  781. *> format is not intended to be modified by the user.
    782. 

  782. *>
    783. 

  783. *>-----------------------------------------------------------------------
    784. 

  784. *>
    785. 

  785. *> DBL and DBK input files:
    786. 

  786. *>
    787. 

  787. *> line 1:  'DBL' in columns 1-3 to test SGEBAL, or 'DBK' in
    788. 

  788. *>          columns 1-3 to test SGEBAK.
    789. 

  789. *>
    790. 

  790. *> The remaining lines consist of specially constructed test cases.
    791. 

  791. *>
    792. 

  792. *>-----------------------------------------------------------------------
    793. 

  793. *>
    794. 

  794. *> DGL and DGK input files:
    795. 

  795. *>
    796. 

  796. *> line 1:  'DGL' in columns 1-3 to test DGGBAL, or 'DGK' in
    797. 

  797. *>          columns 1-3 to test DGGBAK.
    798. 

  798. *>
    799. 

  799. *> The remaining lines consist of specially constructed test cases.
    800. 

  800. *>
    801. 

  801. *>-----------------------------------------------------------------------
    802. 

  802. *>
    803. 

  803. *> GLM data file:
    804. 

  804. *>
    805. 

  805. *> line 1:  'GLM' in columns 1 to 3.
    806. 

  806. *>
    807. 

  807. *> line 2:  NN, INTEGER
    808. 

  808. *>          Number of values of M, P, and N.
    809. 

  809. *>
    810. 

  810. *> line 3:  MVAL, INTEGER array, dimension(NN)
    811. 

  811. *>          Values of M (row dimension).
    812. 

  812. *>
    813. 

  813. *> line 4:  PVAL, INTEGER array, dimension(NN)
    814. 

  814. *>          Values of P (row dimension).
    815. 

  815. *>
    816. 

  816. *> line 5:  NVAL, INTEGER array, dimension(NN)
    817. 

  817. *>          Values of N (column dimension), note M <= N <= M+P.
    818. 

  818. *>
    819. 

  819. *> line 6:  THRESH, REAL
    820. 

  820. *>          Threshold value for the test ratios.  Information will be
    821. 

  821. *>          printed about each test for which the test ratio is greater
    822. 

  822. *>          than or equal to the threshold.
    823. 

  823. *>
    824. 

  824. *> line 7:  TSTERR, LOGICAL
    825. 

  825. *>          Flag indicating whether or not to test the error exits for
    826. 

  826. *>          the LAPACK routines and driver routines.
    827. 

  827. *>
    828. 

  828. *> line 8:  NEWSD, INTEGER
    829. 

  829. *>          A code indicating how to set the random number seed.
    830. 

  830. *>          = 0:  Set the seed to a default value before each run
    831. 

  831. *>          = 1:  Initialize the seed to a default value only before the
    832. 

  832. *>                first run
    833. 

  833. *>          = 2:  Like 1, but use the seed values on the next line
    834. 

  834. *>
    835. 

  835. *> If line 8 was 2:
    836. 

  836. *>
    837. 

  837. *> line 9:  INTEGER array, dimension (4)
    838. 

  838. *>          Four integer values for the random number seed.
    839. 

  839. *>
    840. 

  840. *> lines 9-EOF:  Lines specifying matrix types, as for NEP.
    841. 

  841. *>          The 3-character path name is 'GLM' for the generalized
    842. 

  842. *>          linear regression model routines.
    843. 

  843. *>
    844. 

  844. *>-----------------------------------------------------------------------
    845. 

  845. *>
    846. 

  846. *> GQR data file:
    847. 

  847. *>
    848. 

  848. *> line 1:  'GQR' in columns 1 to 3.
    849. 

  849. *>
    850. 

  850. *> line 2:  NN, INTEGER
    851. 

  851. *>          Number of values of M, P, and N.
    852. 

  852. *>
    853. 

  853. *> line 3:  MVAL, INTEGER array, dimension(NN)
    854. 

  854. *>          Values of M.
    855. 

  855. *>
    856. 

  856. *> line 4:  PVAL, INTEGER array, dimension(NN)
    857. 

  857. *>          Values of P.
    858. 

  858. *>
    859. 

  859. *> line 5:  NVAL, INTEGER array, dimension(NN)
    860. 

  860. *>          Values of N.
    861. 

  861. *>
    862. 

  862. *> line 6:  THRESH, REAL
    863. 

  863. *>          Threshold value for the test ratios.  Information will be
    864. 

  864. *>          printed about each test for which the test ratio is greater
    865. 

  865. *>          than or equal to the threshold.
    866. 

  866. *>
    867. 

  867. *> line 7:  TSTERR, LOGICAL
    868. 

  868. *>          Flag indicating whether or not to test the error exits for
    869. 

  869. *>          the LAPACK routines and driver routines.
    870. 

  870. *>
    871. 

  871. *> line 8:  NEWSD, INTEGER
    872. 

  872. *>          A code indicating how to set the random number seed.
    873. 

  873. *>          = 0:  Set the seed to a default value before each run
    874. 

  874. *>          = 1:  Initialize the seed to a default value only before the
    875. 

  875. *>                first run
    876. 

  876. *>          = 2:  Like 1, but use the seed values on the next line
    877. 

  877. *>
    878. 

  878. *> If line 8 was 2:
    879. 

  879. *>
    880. 

  880. *> line 9:  INTEGER array, dimension (4)
    881. 

  881. *>          Four integer values for the random number seed.
    882. 

  882. *>
    883. 

  883. *> lines 9-EOF:  Lines specifying matrix types, as for NEP.
    884. 

  884. *>          The 3-character path name is 'GQR' for the generalized
    885. 

  885. *>          QR and RQ routines.
    886. 

  886. *>
    887. 

  887. *>-----------------------------------------------------------------------
    888. 

  888. *>
    889. 

  889. *> GSV data file:
    890. 

  890. *>
    891. 

  891. *> line 1:  'GSV' in columns 1 to 3.
    892. 

  892. *>
    893. 

  893. *> line 2:  NN, INTEGER
    894. 

  894. *>          Number of values of M, P, and N.
    895. 

  895. *>
    896. 

  896. *> line 3:  MVAL, INTEGER array, dimension(NN)
    897. 

  897. *>          Values of M (row dimension).
    898. 

  898. *>
    899. 

  899. *> line 4:  PVAL, INTEGER array, dimension(NN)
    900. 

  900. *>          Values of P (row dimension).
    901. 

  901. *>
    902. 

  902. *> line 5:  NVAL, INTEGER array, dimension(NN)
    903. 

  903. *>          Values of N (column dimension).
    904. 

  904. *>
    905. 

  905. *> line 6:  THRESH, REAL
    906. 

  906. *>          Threshold value for the test ratios.  Information will be
    907. 

  907. *>          printed about each test for which the test ratio is greater
    908. 

  908. *>          than or equal to the threshold.
    909. 

  909. *>
    910. 

  910. *> line 7:  TSTERR, LOGICAL
    911. 

  911. *>          Flag indicating whether or not to test the error exits for
    912. 

  912. *>          the LAPACK routines and driver routines.
    913. 

  913. *>
    914. 

  914. *> line 8:  NEWSD, INTEGER
    915. 

  915. *>          A code indicating how to set the random number seed.
    916. 

  916. *>          = 0:  Set the seed to a default value before each run
    917. 

  917. *>          = 1:  Initialize the seed to a default value only before the
    918. 

  918. *>                first run
    919. 

  919. *>          = 2:  Like 1, but use the seed values on the next line
    920. 

  920. *>
    921. 

  921. *> If line 8 was 2:
    922. 

  922. *>
    923. 

  923. *> line 9:  INTEGER array, dimension (4)
    924. 

  924. *>          Four integer values for the random number seed.
    925. 

  925. *>
    926. 

  926. *> lines 9-EOF:  Lines specifying matrix types, as for NEP.
    927. 

  927. *>          The 3-character path name is 'GSV' for the generalized
    928. 

  928. *>          SVD routines.
    929. 

  929. *>
    930. 

  930. *>-----------------------------------------------------------------------
    931. 

  931. *>
    932. 

  932. *> CSD data file:
    933. 

  933. *>
    934. 

  934. *> line 1:  'CSD' in columns 1 to 3.
    935. 

  935. *>
    936. 

  936. *> line 2:  NM, INTEGER
    937. 

  937. *>          Number of values of M, P, and N.
    938. 

  938. *>
    939. 

  939. *> line 3:  MVAL, INTEGER array, dimension(NM)
    940. 

  940. *>          Values of M (row and column dimension of orthogonal matrix).
    941. 

  941. *>
    942. 

  942. *> line 4:  PVAL, INTEGER array, dimension(NM)
    943. 

  943. *>          Values of P (row dimension of top-left block).
    944. 

  944. *>
    945. 

  945. *> line 5:  NVAL, INTEGER array, dimension(NM)
    946. 

  946. *>          Values of N (column dimension of top-left block).
    947. 

  947. *>
    948. 

  948. *> line 6:  THRESH, REAL
    949. 

  949. *>          Threshold value for the test ratios.  Information will be
    950. 

  950. *>          printed about each test for which the test ratio is greater
    951. 

  951. *>          than or equal to the threshold.
    952. 

  952. *>
    953. 

  953. *> line 7:  TSTERR, LOGICAL
    954. 

  954. *>          Flag indicating whether or not to test the error exits for
    955. 

  955. *>          the LAPACK routines and driver routines.
    956. 

  956. *>
    957. 

  957. *> line 8:  NEWSD, INTEGER
    958. 

  958. *>          A code indicating how to set the random number seed.
    959. 

  959. *>          = 0:  Set the seed to a default value before each run
    960. 

  960. *>          = 1:  Initialize the seed to a default value only before the
    961. 

  961. *>                first run
    962. 

  962. *>          = 2:  Like 1, but use the seed values on the next line
    963. 

  963. *>
    964. 

  964. *> If line 8 was 2:
    965. 

  965. *>
    966. 

  966. *> line 9:  INTEGER array, dimension (4)
    967. 

  967. *>          Four integer values for the random number seed.
    968. 

  968. *>
    969. 

  969. *> lines 9-EOF:  Lines specifying matrix types, as for NEP.
    970. 

  970. *>          The 3-character path name is 'CSD' for the CSD routine.
    971. 

  971. *>
    972. 

  972. *>-----------------------------------------------------------------------
    973. 

  973. *>
    974. 

  974. *> LSE data file:
    975. 

  975. *>
    976. 

  976. *> line 1:  'LSE' in columns 1 to 3.
    977. 

  977. *>
    978. 

  978. *> line 2:  NN, INTEGER
    979. 

  979. *>          Number of values of M, P, and N.
    980. 

  980. *>
    981. 

  981. *> line 3:  MVAL, INTEGER array, dimension(NN)
    982. 

  982. *>          Values of M.
    983. 

  983. *>
    984. 

  984. *> line 4:  PVAL, INTEGER array, dimension(NN)
    985. 

  985. *>          Values of P.
    986. 

  986. *>
    987. 

  987. *> line 5:  NVAL, INTEGER array, dimension(NN)
    988. 

  988. *>          Values of N, note P <= N <= P+M.
    989. 

  989. *>
    990. 

  990. *> line 6:  THRESH, REAL
    991. 

  991. *>          Threshold value for the test ratios.  Information will be
    992. 

  992. *>          printed about each test for which the test ratio is greater
    993. 

  993. *>          than or equal to the threshold.
    994. 

  994. *>
    995. 

  995. *> line 7:  TSTERR, LOGICAL
    996. 

  996. *>          Flag indicating whether or not to test the error exits for
    997. 

  997. *>          the LAPACK routines and driver routines.
    998. 

  998. *>
    999. 

  999. *> line 8:  NEWSD, INTEGER
    1000. 

  1000. *>          A code indicating how to set the random number seed.
    1001. 

  1001. *>          = 0:  Set the seed to a default value before each run
    1002. 

  1002. *>          = 1:  Initialize the seed to a default value only before the
    1003. 

  1003. *>                first run
    1004. 

  1004. *>          = 2:  Like 1, but use the seed values on the next line
    1005. 

  1005. *>
    1006. 

  1006. *> If line 8 was 2:
    1007. 

  1007. *>
    1008. 

  1008. *> line 9:  INTEGER array, dimension (4)
    1009. 

  1009. *>          Four integer values for the random number seed.
    1010. 

  1010. *>
    1011. 

  1011. *> lines 9-EOF:  Lines specifying matrix types, as for NEP.
    1012. 

  1012. *>          The 3-character path name is 'GSV' for the generalized
    1013. 

  1013. *>          SVD routines.
    1014. 

  1014. *>
    1015. 

  1015. *>-----------------------------------------------------------------------
    1016. 

  1016. *>
    1017. 

  1017. *> NMAX is currently set to 132 and must be at least 12 for some of the
    1018. 

  1018. *> precomputed examples, and LWORK = NMAX*(5*NMAX+5)+1 in the parameter
    1019. 

  1019. *> statements below.  For SVD, we assume NRHS may be as big as N.  The
    1020. 

  1020. *> parameter NEED is set to 14 to allow for 14 N-by-N matrices for DGG.
    1021. 

  1021. *> \endverbatim
    1022. 

  1022. *
    1023. 

  1023. *  Arguments:
    1024. 

  1024. *  ==========
    1025. 

  1025. *
    1026. 

  1026. *
    1027. 

  1027. *  Authors:
    1028. 

  1028. *  ========
    1029. 

  1029. *
    1030. 

  1030. *> \author Univ. of Tennessee 
    1031. 

  1031. *> \author Univ. of California Berkeley 
    1032. 

  1032. *> \author Univ. of Colorado Denver 
    1033. 

  1033. *> \author NAG Ltd. 
    1034. 

  1034. *
    1035. 

  1035. *> \date November 2011
    1036. 

  1036. *
    1037. 

  1037. *> \ingroup double_eig
    1038. 

  1038. *
    1039. 

  1039. *  =====================================================================      PROGRAM DCHKEE
    1040. 

  1040. *
    1041. 

  1041. *  -- LAPACK test routine (version 3.4.0) --
    1042. 

  1042. *  -- LAPACK is a software package provided by Univ. of Tennessee,    --
    1043. 

  1043. *  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
    1044. 

  1044. *     November 2011
    1045. 

  1045. *
    1046. 

  1046. *  =====================================================================
    1047. 

  1047. *
    1048. 

  1048. *     .. Parameters ..
    1049. 

  1049.       INTEGER            NMAX
    1050. 

  1050.       PARAMETER          ( NMAX = 132 )
    1051. 

  1051.       INTEGER            NCMAX
    1052. 

  1052.       PARAMETER          ( NCMAX = 20 )
    1053. 

  1053.       INTEGER            NEED
    1054. 

  1054.       PARAMETER          ( NEED = 14 )
    1055. 

  1055.       INTEGER            LWORK
    1056. 

  1056.       PARAMETER          ( LWORK = NMAX*( 5*NMAX+5 )+1 )
    1057. 

  1057.       INTEGER            LIWORK
    1058. 

  1058.       PARAMETER          ( LIWORK = NMAX*( 5*NMAX+20 ) )
    1059. 

  1059.       INTEGER            MAXIN
    1060. 

  1060.       PARAMETER          ( MAXIN = 20 )
    1061. 

  1061.       INTEGER            MAXT
    1062. 

  1062.       PARAMETER          ( MAXT = 30 )
    1063. 

  1063.       INTEGER            NIN, NOUT
    1064. 

  1064.       PARAMETER          ( NIN = 5, NOUT = 6 )
    1065. 

  1065. *     ..
    1066. 

  1066. *     .. Local Scalars ..
    1067. 

  1067.       LOGICAL            CSD, DBB, DGG, DSB, FATAL, GLM, GQR, GSV, LSE,
    1068. 

  1068.      $                   NEP, DBK, DBL, SEP, DES, DEV, DGK, DGL, DGS,
    1069. 

  1069.      $                   DGV, DGX, DSX, SVD, DVX, DXV, TSTCHK, TSTDIF,
    1070. 

  1070.      $                   TSTDRV, TSTERR
    1071. 

  1071.       CHARACTER          C1
    1072. 

  1072.       CHARACTER*3        C3, PATH
    1073. 

  1073.       CHARACTER*32       VNAME
    1074. 

  1074.       CHARACTER*10       INTSTR
    1075. 

  1075.       CHARACTER*80       LINE
    1076. 

  1076.       INTEGER            I, I1, IC, INFO, ITMP, K, LENP, MAXTYP, NEWSD,
    1077. 

  1077.      $                   NK, NN, NPARMS, NRHS, NTYPES,
    1078. 

  1078.      $                   VERS_MAJOR, VERS_MINOR, VERS_PATCH 
    1079. 

  1079.       DOUBLE PRECISION   EPS, S1, S2, THRESH, THRSHN
    1080. 

  1080. *     ..
    1081. 

  1081. *     .. Local Arrays ..
    1082. 

  1082.       LOGICAL            DOTYPE( MAXT ), LOGWRK( NMAX )
    1083. 

  1083.       INTEGER            IOLDSD( 4 ), ISEED( 4 ), IWORK( LIWORK ),
    1084. 

  1084.      $                   KVAL( MAXIN ), MVAL( MAXIN ), MXBVAL( MAXIN ),
    1085. 

  1085.      $                   NBCOL( MAXIN ), NBMIN( MAXIN ), NBVAL( MAXIN ),
    1086. 

  1086.      $                   NSVAL( MAXIN ), NVAL( MAXIN ), NXVAL( MAXIN ),
    1087. 

  1087.      $                   PVAL( MAXIN )
    1088. 

  1088.       INTEGER            INMIN( MAXIN ), INWIN( MAXIN ), INIBL( MAXIN ),
    1089. 

  1089.      $                   ISHFTS( MAXIN ), IACC22( MAXIN )
    1090. 

  1090.       DOUBLE PRECISION   A( NMAX*NMAX, NEED ), B( NMAX*NMAX, 5 ),
    1091. 

  1091.      $                   C( NCMAX*NCMAX, NCMAX*NCMAX ), D( NMAX, 12 ),
    1092. 

  1092.      $                   RESULT( 500 ), TAUA( NMAX ), TAUB( NMAX ),
    1093. 

  1093.      $                   WORK( LWORK ), X( 5*NMAX )
    1094. 

  1094. *     ..
    1095. 

  1095. *     .. External Functions ..
    1096. 

  1096.       LOGICAL            LSAMEN
    1097. 

  1097.       DOUBLE PRECISION   DLAMCH, DSECND
    1098. 

  1098.       EXTERNAL           LSAMEN, DLAMCH, DSECND
    1099. 

  1099. *     ..
    1100. 

  1100. *     .. External Subroutines ..
    1101. 

  1101.       EXTERNAL           ALAREQ, DCHKBB, DCHKBD, DCHKBK, DCHKBL, DCHKEC,
    1102. 

  1102.      $                   DCHKGG, DCHKGK, DCHKGL, DCHKHS, DCHKSB, DCHKST,
    1103. 

  1103.      $                   DCKCSD, DCKGLM, DCKGQR, DCKGSV, DCKLSE, DDRGES,
    1104. 

  1104.      $                   DDRGEV, DDRGSX, DDRGVX, DDRVBD, DDRVES, DDRVEV,
    1105. 

  1105.      $                   DDRVGG, DDRVSG, DDRVST, DDRVSX, DDRVVX, DERRBD,
    1106. 

  1106.      $                   DERRED, DERRGG, DERRHS, DERRST, ILAVER, XLAENV
    1107. 

  1107. *     ..
    1108. 

  1108. *     .. Intrinsic Functions ..
    1109. 

  1109.       INTRINSIC          LEN, MIN
    1110. 

  1110. *     ..
    1111. 

  1111. *     .. Scalars in Common ..
    1112. 

  1112.       LOGICAL            LERR, OK
    1113. 

  1113.       CHARACTER*32       SRNAMT
    1114. 

  1114.       INTEGER            INFOT, MAXB, NPROC, NSHIFT, NUNIT, SELDIM,
    1115. 

  1115.      $                   SELOPT
    1116. 

  1116. *     ..
    1117. 

  1117. *     .. Arrays in Common ..
    1118. 

  1118.       LOGICAL            SELVAL( 20 )
    1119. 

  1119.       INTEGER            IPARMS( 100 )
    1120. 

  1120.       DOUBLE PRECISION   SELWI( 20 ), SELWR( 20 )
    1121. 

  1121. *     ..
    1122. 

  1122. *     .. Common blocks ..
    1123. 

  1123.       COMMON             / CENVIR / NPROC, NSHIFT, MAXB
    1124. 

  1124.       COMMON             / INFOC / INFOT, NUNIT, OK, LERR
    1125. 

  1125.       COMMON             / SRNAMC / SRNAMT
    1126. 

  1126.       COMMON             / SSLCT / SELOPT, SELDIM, SELVAL, SELWR, SELWI
    1127. 

  1127.       COMMON             / CLAENV / IPARMS
    1128. 

  1128. *     ..
    1129. 

  1129. *     .. Data statements ..
    1130. 

  1130.       DATA               INTSTR / '0123456789' /
    1131. 

  1131.       DATA               IOLDSD / 0, 0, 0, 1 /
    1132. 

  1132. *     ..
    1133. 

  1133. *     .. Executable Statements ..
    1134. 

  1134. *
    1135. 

  1135.       S1 = DSECND( )
    1136. 

  1136.       FATAL = .FALSE.
    1137. 

  1137.       NUNIT = NOUT
    1138. 

  1138. *
    1139. 

  1139. *     Return to here to read multiple sets of data
    1140. 

  1140. *
    1141. 

  1141.    10 CONTINUE
    1142. 

  1142. *
    1143. 

  1143. *     Read the first line and set the 3-character test path
    1144. 

  1144. *
    1145. 

  1145.       READ( NIN, FMT = '(A80)', END = 380 )LINE
    1146. 

  1146.       PATH = LINE( 1: 3 )
    1147. 

  1147.       NEP = LSAMEN( 3, PATH, 'NEP' ) .OR. LSAMEN( 3, PATH, 'DHS' )
    1148. 

  1148.       SEP = LSAMEN( 3, PATH, 'SEP' ) .OR. LSAMEN( 3, PATH, 'DST' ) .OR.
    1149. 

  1149.      $      LSAMEN( 3, PATH, 'DSG' )
    1150. 

  1150.       SVD = LSAMEN( 3, PATH, 'SVD' ) .OR. LSAMEN( 3, PATH, 'DBD' )
    1151. 

  1151.       DEV = LSAMEN( 3, PATH, 'DEV' )
    1152. 

  1152.       DES = LSAMEN( 3, PATH, 'DES' )
    1153. 

  1153.       DVX = LSAMEN( 3, PATH, 'DVX' )
    1154. 

  1154.       DSX = LSAMEN( 3, PATH, 'DSX' )
    1155. 

  1155.       DGG = LSAMEN( 3, PATH, 'DGG' )
    1156. 

  1156.       DGS = LSAMEN( 3, PATH, 'DGS' )
    1157. 

  1157.       DGX = LSAMEN( 3, PATH, 'DGX' )
    1158. 

  1158.       DGV = LSAMEN( 3, PATH, 'DGV' )
    1159. 

  1159.       DXV = LSAMEN( 3, PATH, 'DXV' )
    1160. 

  1160.       DSB = LSAMEN( 3, PATH, 'DSB' )
    1161. 

  1161.       DBB = LSAMEN( 3, PATH, 'DBB' )
    1162. 

  1162.       GLM = LSAMEN( 3, PATH, 'GLM' )
    1163. 

  1163.       GQR = LSAMEN( 3, PATH, 'GQR' ) .OR. LSAMEN( 3, PATH, 'GRQ' )
    1164. 

  1164.       GSV = LSAMEN( 3, PATH, 'GSV' )
    1165. 

  1165.       CSD = LSAMEN( 3, PATH, 'CSD' )
    1166. 

  1166.       LSE = LSAMEN( 3, PATH, 'LSE' )
    1167. 

  1167.       DBL = LSAMEN( 3, PATH, 'DBL' )
    1168. 

  1168.       DBK = LSAMEN( 3, PATH, 'DBK' )
    1169. 

  1169.       DGL = LSAMEN( 3, PATH, 'DGL' )
    1170. 

  1170.       DGK = LSAMEN( 3, PATH, 'DGK' )
    1171. 

  1171. *
    1172. 

  1172. *     Report values of parameters.
    1173. 

  1173. *
    1174. 

  1174.       IF( PATH.EQ.'   ' ) THEN
    1175. 

  1175.          GO TO 10
    1176. 

  1176.       ELSE IF( NEP ) THEN
    1177. 

  1177.          WRITE( NOUT, FMT = 9987 )
    1178. 

  1178.       ELSE IF( SEP ) THEN
    1179. 

  1179.          WRITE( NOUT, FMT = 9986 )
    1180. 

  1180.       ELSE IF( SVD ) THEN
    1181. 

  1181.          WRITE( NOUT, FMT = 9985 )
    1182. 

  1182.       ELSE IF( DEV ) THEN
    1183. 

  1183.          WRITE( NOUT, FMT = 9979 )
    1184. 

  1184.       ELSE IF( DES ) THEN
    1185. 

  1185.          WRITE( NOUT, FMT = 9978 )
    1186. 

  1186.       ELSE IF( DVX ) THEN
    1187. 

  1187.          WRITE( NOUT, FMT = 9977 )
    1188. 

  1188.       ELSE IF( DSX ) THEN
    1189. 

  1189.          WRITE( NOUT, FMT = 9976 )
    1190. 

  1190.       ELSE IF( DGG ) THEN
    1191. 

  1191.          WRITE( NOUT, FMT = 9975 )
    1192. 

  1192.       ELSE IF( DGS ) THEN
    1193. 

  1193.          WRITE( NOUT, FMT = 9964 )
    1194. 

  1194.       ELSE IF( DGX ) THEN
    1195. 

  1195.          WRITE( NOUT, FMT = 9965 )
    1196. 

  1196.       ELSE IF( DGV ) THEN
    1197. 

  1197.          WRITE( NOUT, FMT = 9963 )
    1198. 

  1198.       ELSE IF( DXV ) THEN
    1199. 

  1199.          WRITE( NOUT, FMT = 9962 )
    1200. 

  1200.       ELSE IF( DSB ) THEN
    1201. 

  1201.          WRITE( NOUT, FMT = 9974 )
    1202. 

  1202.       ELSE IF( DBB ) THEN
    1203. 

  1203.          WRITE( NOUT, FMT = 9967 )
    1204. 

  1204.       ELSE IF( GLM ) THEN
    1205. 

  1205.          WRITE( NOUT, FMT = 9971 )
    1206. 

  1206.       ELSE IF( GQR ) THEN
    1207. 

  1207.          WRITE( NOUT, FMT = 9970 )
    1208. 

  1208.       ELSE IF( GSV ) THEN
    1209. 

  1209.          WRITE( NOUT, FMT = 9969 )
    1210. 

  1210.       ELSE IF( CSD ) THEN
    1211. 

  1211.          WRITE( NOUT, FMT = 9960 )
    1212. 

  1212.       ELSE IF( LSE ) THEN
    1213. 

  1213.          WRITE( NOUT, FMT = 9968 )
    1214. 

  1214.       ELSE IF( DBL ) THEN
    1215. 

  1215. *
    1216. 

  1216. *        DGEBAL:  Balancing
    1217. 

  1217. *
    1218. 

  1218.          CALL DCHKBL( NIN, NOUT )
    1219. 

  1219.          GO TO 10
    1220. 

  1220.       ELSE IF( DBK ) THEN
    1221. 

  1221. *
    1222. 

  1222. *        DGEBAK:  Back transformation
    1223. 

  1223. *
    1224. 

  1224.          CALL DCHKBK( NIN, NOUT )
    1225. 

  1225.          GO TO 10
    1226. 

  1226.       ELSE IF( DGL ) THEN
    1227. 

  1227. *
    1228. 

  1228. *        DGGBAL:  Balancing
    1229. 

  1229. *
    1230. 

  1230.          CALL DCHKGL( NIN, NOUT )
    1231. 

  1231.          GO TO 10
    1232. 

  1232.       ELSE IF( DGK ) THEN
    1233. 

  1233. *
    1234. 

  1234. *        DGGBAK:  Back transformation
    1235. 

  1235. *
    1236. 

  1236.          CALL DCHKGK( NIN, NOUT )
    1237. 

  1237.          GO TO 10
    1238. 

  1238.       ELSE IF( LSAMEN( 3, PATH, 'DEC' ) ) THEN
    1239. 

  1239. *
    1240. 

  1240. *        DEC:  Eigencondition estimation
    1241. 

  1241. *
    1242. 

  1242.          READ( NIN, FMT = * )THRESH
    1243. 

  1243.          CALL XLAENV( 1, 1 )
    1244. 

  1244.          CALL XLAENV( 12, 11 )
    1245. 

  1245.          CALL XLAENV( 13, 2 )
    1246. 

  1246.          CALL XLAENV( 14, 0 )
    1247. 

  1247.          CALL XLAENV( 15, 2 )
    1248. 

  1248.          CALL XLAENV( 16, 2 )
    1249. 

  1249.          TSTERR = .TRUE.
    1250. 

  1250.          CALL DCHKEC( THRESH, TSTERR, NIN, NOUT )
    1251. 

  1251.          GO TO 10
    1252. 

  1252.       ELSE
    1253. 

  1253.          WRITE( NOUT, FMT = 9992 )PATH
    1254. 

  1254.          GO TO 10
    1255. 

  1255.       END IF
    1256. 

  1256.       CALL ILAVER( VERS_MAJOR, VERS_MINOR, VERS_PATCH )
    1257. 

  1257.       WRITE( NOUT, FMT = 9972 ) VERS_MAJOR, VERS_MINOR, VERS_PATCH
    1258. 

  1258.       WRITE( NOUT, FMT = 9984 )
    1259. 

  1259. *
    1260. 

  1260. *     Read the number of values of M, P, and N.
    1261. 

  1261. *
    1262. 

  1262.       READ( NIN, FMT = * )NN
    1263. 

  1263.       IF( NN.LT.0 ) THEN
    1264. 

  1264.          WRITE( NOUT, FMT = 9989 )'   NN ', NN, 1
    1265. 

  1265.          NN = 0
    1266. 

  1266.          FATAL = .TRUE.
    1267. 

  1267.       ELSE IF( NN.GT.MAXIN ) THEN
    1268. 

  1268.          WRITE( NOUT, FMT = 9988 )'   NN ', NN, MAXIN
    1269. 

  1269.          NN = 0
    1270. 

  1270.          FATAL = .TRUE.
    1271. 

  1271.       END IF
    1272. 

  1272. *
    1273. 

  1273. *     Read the values of M
    1274. 

  1274. *
    1275. 

  1275.       IF( .NOT.( DGX .OR. DXV ) ) THEN
    1276. 

  1276.          READ( NIN, FMT = * )( MVAL( I ), I = 1, NN )
    1277. 

  1277.          IF( SVD ) THEN
    1278. 

  1278.             VNAME = '    M '
    1279. 

  1279.          ELSE
    1280. 

  1280.             VNAME = '    N '
    1281. 

  1281.          END IF
    1282. 

  1282.          DO 20 I = 1, NN
    1283. 

  1283.             IF( MVAL( I ).LT.0 ) THEN
    1284. 

  1284.                WRITE( NOUT, FMT = 9989 )VNAME, MVAL( I ), 0
    1285. 

  1285.                FATAL = .TRUE.
    1286. 

  1286.             ELSE IF( MVAL( I ).GT.NMAX ) THEN
    1287. 

  1287.                WRITE( NOUT, FMT = 9988 )VNAME, MVAL( I ), NMAX
    1288. 

  1288.                FATAL = .TRUE.
    1289. 

  1289.             END IF
    1290. 

  1290.    20    CONTINUE
    1291. 

  1291.          WRITE( NOUT, FMT = 9983 )'M:    ', ( MVAL( I ), I = 1, NN )
    1292. 

  1292.       END IF
    1293. 

  1293. *
    1294. 

  1294. *     Read the values of P
    1295. 

  1295. *
    1296. 

  1296.       IF( GLM .OR. GQR .OR. GSV .OR. CSD .OR. LSE ) THEN
    1297. 

  1297.          READ( NIN, FMT = * )( PVAL( I ), I = 1, NN )
    1298. 

  1298.          DO 30 I = 1, NN
    1299. 

  1299.             IF( PVAL( I ).LT.0 ) THEN
    1300. 

  1300.                WRITE( NOUT, FMT = 9989 )' P  ', PVAL( I ), 0
    1301. 

  1301.                FATAL = .TRUE.
    1302. 

  1302.             ELSE IF( PVAL( I ).GT.NMAX ) THEN
    1303. 

  1303.                WRITE( NOUT, FMT = 9988 )' P  ', PVAL( I ), NMAX
    1304. 

  1304.                FATAL = .TRUE.
    1305. 

  1305.             END IF
    1306. 

  1306.    30    CONTINUE
    1307. 

  1307.          WRITE( NOUT, FMT = 9983 )'P:    ', ( PVAL( I ), I = 1, NN )
    1308. 

  1308.       END IF
    1309. 

  1309. *
    1310. 

  1310. *     Read the values of N
    1311. 

  1311. *
    1312. 

  1312.       IF( SVD .OR. DBB .OR. GLM .OR. GQR .OR. GSV .OR. CSD .OR.
    1313. 

  1313.      $    LSE ) THEN
    1314. 

  1314.          READ( NIN, FMT = * )( NVAL( I ), I = 1, NN )
    1315. 

  1315.          DO 40 I = 1, NN
    1316. 

  1316.             IF( NVAL( I ).LT.0 ) THEN
    1317. 

  1317.                WRITE( NOUT, FMT = 9989 )'    N ', NVAL( I ), 0
    1318. 

  1318.                FATAL = .TRUE.
    1319. 

  1319.             ELSE IF( NVAL( I ).GT.NMAX ) THEN
    1320. 

  1320.                WRITE( NOUT, FMT = 9988 )'    N ', NVAL( I ), NMAX
    1321. 

  1321.                FATAL = .TRUE.
    1322. 

  1322.             END IF
    1323. 

  1323.    40    CONTINUE
    1324. 

  1324.       ELSE
    1325. 

  1325.          DO 50 I = 1, NN
    1326. 

  1326.             NVAL( I ) = MVAL( I )
    1327. 

  1327.    50    CONTINUE
    1328. 

  1328.       END IF
    1329. 

  1329.       IF( .NOT.( DGX .OR. DXV ) ) THEN
    1330. 

  1330.          WRITE( NOUT, FMT = 9983 )'N:    ', ( NVAL( I ), I = 1, NN )
    1331. 

  1331.       ELSE
    1332. 

  1332.          WRITE( NOUT, FMT = 9983 )'N:    ', NN
    1333. 

  1333.       END IF
    1334. 

  1334. *
    1335. 

  1335. *     Read the number of values of K, followed by the values of K
    1336. 

  1336. *
    1337. 

  1337.       IF( DSB .OR. DBB ) THEN
    1338. 

  1338.          READ( NIN, FMT = * )NK
    1339. 

  1339.          READ( NIN, FMT = * )( KVAL( I ), I = 1, NK )
    1340. 

  1340.          DO 60 I = 1, NK
    1341. 

  1341.             IF( KVAL( I ).LT.0 ) THEN
    1342. 

  1342.                WRITE( NOUT, FMT = 9989 )'    K ', KVAL( I ), 0
    1343. 

  1343.                FATAL = .TRUE.
    1344. 

  1344.             ELSE IF( KVAL( I ).GT.NMAX ) THEN
    1345. 

  1345.                WRITE( NOUT, FMT = 9988 )'    K ', KVAL( I ), NMAX
    1346. 

  1346.                FATAL = .TRUE.
    1347. 

  1347.             END IF
    1348. 

  1348.    60    CONTINUE
    1349. 

  1349.          WRITE( NOUT, FMT = 9983 )'K:    ', ( KVAL( I ), I = 1, NK )
    1350. 

  1350.       END IF
    1351. 

  1351. *
    1352. 

  1352.       IF( DEV .OR. DES .OR. DVX .OR. DSX ) THEN
    1353. 

  1353. *
    1354. 

  1354. *        For the nonsymmetric QR driver routines, only one set of
    1355. 

  1355. *        parameters is allowed.
    1356. 

  1356. *
    1357. 

  1357.          READ( NIN, FMT = * )NBVAL( 1 ), NBMIN( 1 ), NXVAL( 1 ),
    1358. 

  1358.      $      INMIN( 1 ), INWIN( 1 ), INIBL(1), ISHFTS(1), IACC22(1)
    1359. 

  1359.          IF( NBVAL( 1 ).LT.1 ) THEN
    1360. 

  1360.             WRITE( NOUT, FMT = 9989 )'   NB ', NBVAL( 1 ), 1
    1361. 

  1361.             FATAL = .TRUE.
    1362. 

  1362.          ELSE IF( NBMIN( 1 ).LT.1 ) THEN
    1363. 

  1363.             WRITE( NOUT, FMT = 9989 )'NBMIN ', NBMIN( 1 ), 1
    1364. 

  1364.             FATAL = .TRUE.
    1365. 

  1365.          ELSE IF( NXVAL( 1 ).LT.1 ) THEN
    1366. 

  1366.             WRITE( NOUT, FMT = 9989 )'   NX ', NXVAL( 1 ), 1
    1367. 

  1367.             FATAL = .TRUE.
    1368. 

  1368.          ELSE IF( INMIN( 1 ).LT.1 ) THEN
    1369. 

  1369.             WRITE( NOUT, FMT = 9989 )'   INMIN ', INMIN( 1 ), 1
    1370. 

  1370.             FATAL = .TRUE.
    1371. 

  1371.          ELSE IF( INWIN( 1 ).LT.1 ) THEN
    1372. 

  1372.             WRITE( NOUT, FMT = 9989 )'   INWIN ', INWIN( 1 ), 1
    1373. 

  1373.             FATAL = .TRUE.
    1374. 

  1374.          ELSE IF( INIBL( 1 ).LT.1 ) THEN
    1375. 

  1375.             WRITE( NOUT, FMT = 9989 )'   INIBL ', INIBL( 1 ), 1
    1376. 

  1376.             FATAL = .TRUE.
    1377. 

  1377.          ELSE IF( ISHFTS( 1 ).LT.1 ) THEN
    1378. 

  1378.             WRITE( NOUT, FMT = 9989 )'   ISHFTS ', ISHFTS( 1 ), 1
    1379. 

  1379.             FATAL = .TRUE.
    1380. 

  1380.          ELSE IF( IACC22( 1 ).LT.0 ) THEN
    1381. 

  1381.             WRITE( NOUT, FMT = 9989 )'   IACC22 ', IACC22( 1 ), 0
    1382. 

  1382.             FATAL = .TRUE.
    1383. 

  1383.          END IF
    1384. 

  1384.          CALL XLAENV( 1, NBVAL( 1 ) )
    1385. 

  1385.          CALL XLAENV( 2, NBMIN( 1 ) )
    1386. 

  1386.          CALL XLAENV( 3, NXVAL( 1 ) )
    1387. 

  1387.          CALL XLAENV(12, MAX( 11, INMIN( 1 ) ) )
    1388. 

  1388.          CALL XLAENV(13, INWIN( 1 ) )
    1389. 

  1389.          CALL XLAENV(14, INIBL( 1 ) )
    1390. 

  1390.          CALL XLAENV(15, ISHFTS( 1 ) )
    1391. 

  1391.          CALL XLAENV(16, IACC22( 1 ) )
    1392. 

  1392.          WRITE( NOUT, FMT = 9983 )'NB:   ', NBVAL( 1 )
    1393. 

  1393.          WRITE( NOUT, FMT = 9983 )'NBMIN:', NBMIN( 1 )
    1394. 

  1394.          WRITE( NOUT, FMT = 9983 )'NX:   ', NXVAL( 1 )
    1395. 

  1395.          WRITE( NOUT, FMT = 9983 )'INMIN:   ', INMIN( 1 )
    1396. 

  1396.          WRITE( NOUT, FMT = 9983 )'INWIN: ', INWIN( 1 )
    1397. 

  1397.          WRITE( NOUT, FMT = 9983 )'INIBL: ', INIBL( 1 )
    1398. 

  1398.          WRITE( NOUT, FMT = 9983 )'ISHFTS: ', ISHFTS( 1 )
    1399. 

  1399.          WRITE( NOUT, FMT = 9983 )'IACC22: ', IACC22( 1 )
    1400. 

  1400. *
    1401. 

  1401.       ELSEIF( DGS .OR. DGX .OR. DGV .OR.  DXV ) THEN
    1402. 

  1402. *
    1403. 

  1403. *        For the nonsymmetric generalized driver routines, only one set
    1404. 

  1404. *        of parameters is allowed.
    1405. 

  1405. *
    1406. 

  1406.          READ( NIN, FMT = * )NBVAL( 1 ), NBMIN( 1 ), NXVAL( 1 ),
    1407. 

  1407.      $      NSVAL( 1 ), MXBVAL( 1 )
    1408. 

  1408.          IF( NBVAL( 1 ).LT.1 ) THEN
    1409. 

  1409.             WRITE( NOUT, FMT = 9989 )'   NB ', NBVAL( 1 ), 1
    1410. 

  1410.             FATAL = .TRUE.
    1411. 

  1411.          ELSE IF( NBMIN( 1 ).LT.1 ) THEN
    1412. 

  1412.             WRITE( NOUT, FMT = 9989 )'NBMIN ', NBMIN( 1 ), 1
    1413. 

  1413.             FATAL = .TRUE.
    1414. 

  1414.          ELSE IF( NXVAL( 1 ).LT.1 ) THEN
    1415. 

  1415.             WRITE( NOUT, FMT = 9989 )'   NX ', NXVAL( 1 ), 1
    1416. 

  1416.             FATAL = .TRUE.
    1417. 

  1417.          ELSE IF( NSVAL( 1 ).LT.2 ) THEN
    1418. 

  1418.             WRITE( NOUT, FMT = 9989 )'   NS ', NSVAL( 1 ), 2
    1419. 

  1419.             FATAL = .TRUE.
    1420. 

  1420.          ELSE IF( MXBVAL( 1 ).LT.1 ) THEN
    1421. 

  1421.             WRITE( NOUT, FMT = 9989 )' MAXB ', MXBVAL( 1 ), 1
    1422. 

  1422.             FATAL = .TRUE.
    1423. 

  1423.          END IF
    1424. 

  1424.          CALL XLAENV( 1, NBVAL( 1 ) )
    1425. 

  1425.          CALL XLAENV( 2, NBMIN( 1 ) )
    1426. 

  1426.          CALL XLAENV( 3, NXVAL( 1 ) )
    1427. 

  1427.          CALL XLAENV( 4, NSVAL( 1 ) )
    1428. 

  1428.          CALL XLAENV( 8, MXBVAL( 1 ) )
    1429. 

  1429.          WRITE( NOUT, FMT = 9983 )'NB:   ', NBVAL( 1 )
    1430. 

  1430.          WRITE( NOUT, FMT = 9983 )'NBMIN:', NBMIN( 1 )
    1431. 

  1431.          WRITE( NOUT, FMT = 9983 )'NX:   ', NXVAL( 1 )
    1432. 

  1432.          WRITE( NOUT, FMT = 9983 )'NS:   ', NSVAL( 1 )
    1433. 

  1433.          WRITE( NOUT, FMT = 9983 )'MAXB: ', MXBVAL( 1 )
    1434. 

  1434. *
    1435. 

  1435.       ELSE IF( .NOT.DSB .AND. .NOT.GLM .AND. .NOT.GQR .AND. .NOT.
    1436. 

  1436.      $         GSV .AND. .NOT.CSD .AND. .NOT.LSE ) THEN
    1437. 

  1437. *
    1438. 

  1438. *        For the other paths, the number of parameters can be varied
    1439. 

  1439. *        from the input file.  Read the number of parameter values.
    1440. 

  1440. *
    1441. 

  1441.          READ( NIN, FMT = * )NPARMS
    1442. 

  1442.          IF( NPARMS.LT.1 ) THEN
    1443. 

  1443.             WRITE( NOUT, FMT = 9989 )'NPARMS', NPARMS, 1
    1444. 

  1444.             NPARMS = 0
    1445. 

  1445.             FATAL = .TRUE.
    1446. 

  1446.          ELSE IF( NPARMS.GT.MAXIN ) THEN
    1447. 

  1447.             WRITE( NOUT, FMT = 9988 )'NPARMS', NPARMS, MAXIN
    1448. 

  1448.             NPARMS = 0
    1449. 

  1449.             FATAL = .TRUE.
    1450. 

  1450.          END IF
    1451. 

  1451. *
    1452. 

  1452. *        Read the values of NB
    1453. 

  1453. *
    1454. 

  1454.          IF( .NOT.DBB ) THEN
    1455. 

  1455.             READ( NIN, FMT = * )( NBVAL( I ), I = 1, NPARMS )
    1456. 

  1456.             DO 70 I = 1, NPARMS
    1457. 

  1457.                IF( NBVAL( I ).LT.0 ) THEN
    1458. 

  1458.                   WRITE( NOUT, FMT = 9989 )'   NB ', NBVAL( I ), 0
    1459. 

  1459.                   FATAL = .TRUE.
    1460. 

  1460.                ELSE IF( NBVAL( I ).GT.NMAX ) THEN
    1461. 

  1461.                   WRITE( NOUT, FMT = 9988 )'   NB ', NBVAL( I ), NMAX
    1462. 

  1462.                   FATAL = .TRUE.
    1463. 

  1463.                END IF
    1464. 

  1464.    70       CONTINUE
    1465. 

  1465.             WRITE( NOUT, FMT = 9983 )'NB:   ',
    1466. 

  1466.      $         ( NBVAL( I ), I = 1, NPARMS )
    1467. 

  1467.          END IF
    1468. 

  1468. *
    1469. 

  1469. *        Read the values of NBMIN
    1470. 

  1470. *
    1471. 

  1471.          IF( NEP .OR. SEP .OR. SVD .OR. DGG ) THEN
    1472. 

  1472.             READ( NIN, FMT = * )( NBMIN( I ), I = 1, NPARMS )
    1473. 

  1473.             DO 80 I = 1, NPARMS
    1474. 

  1474.                IF( NBMIN( I ).LT.0 ) THEN
    1475. 

  1475.                   WRITE( NOUT, FMT = 9989 )'NBMIN ', NBMIN( I ), 0
    1476. 

  1476.                   FATAL = .TRUE.
    1477. 

  1477.                ELSE IF( NBMIN( I ).GT.NMAX ) THEN
    1478. 

  1478.                   WRITE( NOUT, FMT = 9988 )'NBMIN ', NBMIN( I ), NMAX
    1479. 

  1479.                   FATAL = .TRUE.
    1480. 

  1480.                END IF
    1481. 

  1481.    80       CONTINUE
    1482. 

  1482.             WRITE( NOUT, FMT = 9983 )'NBMIN:',
    1483. 

  1483.      $         ( NBMIN( I ), I = 1, NPARMS )
    1484. 

  1484.          ELSE
    1485. 

  1485.             DO 90 I = 1, NPARMS
    1486. 

  1486.                NBMIN( I ) = 1
    1487. 

  1487.    90       CONTINUE
    1488. 

  1488.          END IF
    1489. 

  1489. *
    1490. 

  1490. *        Read the values of NX
    1491. 

  1491. *
    1492. 

  1492.          IF( NEP .OR. SEP .OR. SVD ) THEN
    1493. 

  1493.             READ( NIN, FMT = * )( NXVAL( I ), I = 1, NPARMS )
    1494. 

  1494.             DO 100 I = 1, NPARMS
    1495. 

  1495.                IF( NXVAL( I ).LT.0 ) THEN
    1496. 

  1496.                   WRITE( NOUT, FMT = 9989 )'   NX ', NXVAL( I ), 0
    1497. 

  1497.                   FATAL = .TRUE.
    1498. 

  1498.                ELSE IF( NXVAL( I ).GT.NMAX ) THEN
    1499. 

  1499.                   WRITE( NOUT, FMT = 9988 )'   NX ', NXVAL( I ), NMAX
    1500. 

  1500.                   FATAL = .TRUE.
    1501. 

  1501.                END IF
    1502. 

  1502.   100       CONTINUE
    1503. 

  1503.             WRITE( NOUT, FMT = 9983 )'NX:   ',
    1504. 

  1504.      $         ( NXVAL( I ), I = 1, NPARMS )
    1505. 

  1505.          ELSE
    1506. 

  1506.             DO 110 I = 1, NPARMS
    1507. 

  1507.                NXVAL( I ) = 1
    1508. 

  1508.   110       CONTINUE
    1509. 

  1509.          END IF
    1510. 

  1510. *
    1511. 

  1511. *        Read the values of NSHIFT (if DGG) or NRHS (if SVD
    1512. 

  1512. *        or DBB).
    1513. 

  1513. *
    1514. 

  1514.          IF( SVD .OR. DBB .OR. DGG ) THEN
    1515. 

  1515.             READ( NIN, FMT = * )( NSVAL( I ), I = 1, NPARMS )
    1516. 

  1516.             DO 120 I = 1, NPARMS
    1517. 

  1517.                IF( NSVAL( I ).LT.0 ) THEN
    1518. 

  1518.                   WRITE( NOUT, FMT = 9989 )'   NS ', NSVAL( I ), 0
    1519. 

  1519.                   FATAL = .TRUE.
    1520. 

  1520.                ELSE IF( NSVAL( I ).GT.NMAX ) THEN
    1521. 

  1521.                   WRITE( NOUT, FMT = 9988 )'   NS ', NSVAL( I ), NMAX
    1522. 

  1522.                   FATAL = .TRUE.
    1523. 

  1523.                END IF
    1524. 

  1524.   120       CONTINUE
    1525. 

  1525.             WRITE( NOUT, FMT = 9983 )'NS:   ',
    1526. 

  1526.      $         ( NSVAL( I ), I = 1, NPARMS )
    1527. 

  1527.          ELSE
    1528. 

  1528.             DO 130 I = 1, NPARMS
    1529. 

  1529.                NSVAL( I ) = 1
    1530. 

  1530.   130       CONTINUE
    1531. 

  1531.          END IF
    1532. 

  1532. *
    1533. 

  1533. *        Read the values for MAXB.
    1534. 

  1534. *
    1535. 

  1535.          IF( DGG ) THEN
    1536. 

  1536.             READ( NIN, FMT = * )( MXBVAL( I ), I = 1, NPARMS )
    1537. 

  1537.             DO 140 I = 1, NPARMS
    1538. 

  1538.                IF( MXBVAL( I ).LT.0 ) THEN
    1539. 

  1539.                   WRITE( NOUT, FMT = 9989 )' MAXB ', MXBVAL( I ), 0
    1540. 

  1540.                   FATAL = .TRUE.
    1541. 

  1541.                ELSE IF( MXBVAL( I ).GT.NMAX ) THEN
    1542. 

  1542.                   WRITE( NOUT, FMT = 9988 )' MAXB ', MXBVAL( I ), NMAX
    1543. 

  1543.                   FATAL = .TRUE.
    1544. 

  1544.                END IF
    1545. 

  1545.   140       CONTINUE
    1546. 

  1546.             WRITE( NOUT, FMT = 9983 )'MAXB: ',
    1547. 

  1547.      $         ( MXBVAL( I ), I = 1, NPARMS )
    1548. 

  1548.          ELSE
    1549. 

  1549.             DO 150 I = 1, NPARMS
    1550. 

  1550.                MXBVAL( I ) = 1
    1551. 

  1551.   150       CONTINUE
    1552. 

  1552.          END IF
    1553. 

  1553. *
    1554. 

  1554. *        Read the values for INMIN.
    1555. 

  1555. *
    1556. 

  1556.          IF( NEP ) THEN
    1557. 

  1557.             READ( NIN, FMT = * )( INMIN( I ), I = 1, NPARMS )
    1558. 

  1558.             DO 540 I = 1, NPARMS
    1559. 

  1559.                IF( INMIN( I ).LT.0 ) THEN
    1560. 

  1560.                   WRITE( NOUT, FMT = 9989 )' INMIN ', INMIN( I ), 0
    1561. 

  1561.                   FATAL = .TRUE.
    1562. 

  1562.                END IF
    1563. 

  1563.   540       CONTINUE
    1564. 

  1564.             WRITE( NOUT, FMT = 9983 )'INMIN: ',
    1565. 

  1565.      $         ( INMIN( I ), I = 1, NPARMS )
    1566. 

  1566.          ELSE
    1567. 

  1567.             DO 550 I = 1, NPARMS
    1568. 

  1568.                INMIN( I ) = 1
    1569. 

  1569.   550       CONTINUE
    1570. 

  1570.          END IF
    1571. 

  1571. *
    1572. 

  1572. *        Read the values for INWIN.
    1573. 

  1573. *
    1574. 

  1574.          IF( NEP ) THEN
    1575. 

  1575.             READ( NIN, FMT = * )( INWIN( I ), I = 1, NPARMS )
    1576. 

  1576.             DO 560 I = 1, NPARMS
    1577. 

  1577.                IF( INWIN( I ).LT.0 ) THEN
    1578. 

  1578.                   WRITE( NOUT, FMT = 9989 )' INWIN ', INWIN( I ), 0
    1579. 

  1579.                   FATAL = .TRUE.
    1580. 

  1580.                END IF
    1581. 

  1581.   560       CONTINUE
    1582. 

  1582.             WRITE( NOUT, FMT = 9983 )'INWIN: ',
    1583. 

  1583.      $         ( INWIN( I ), I = 1, NPARMS )
    1584. 

  1584.          ELSE
    1585. 

  1585.             DO 570 I = 1, NPARMS
    1586. 

  1586.                INWIN( I ) = 1
    1587. 

  1587.   570       CONTINUE
    1588. 

  1588.          END IF
    1589. 

  1589. *
    1590. 

  1590. *        Read the values for INIBL.
    1591. 

  1591. *
    1592. 

  1592.          IF( NEP ) THEN
    1593. 

  1593.             READ( NIN, FMT = * )( INIBL( I ), I = 1, NPARMS )
    1594. 

  1594.             DO 580 I = 1, NPARMS
    1595. 

  1595.                IF( INIBL( I ).LT.0 ) THEN
    1596. 

  1596.                   WRITE( NOUT, FMT = 9989 )' INIBL ', INIBL( I ), 0
    1597. 

  1597.                   FATAL = .TRUE.
    1598. 

  1598.                END IF
    1599. 

  1599.   580       CONTINUE
    1600. 

  1600.             WRITE( NOUT, FMT = 9983 )'INIBL: ',
    1601. 

  1601.      $         ( INIBL( I ), I = 1, NPARMS )
    1602. 

  1602.          ELSE
    1603. 

  1603.             DO 590 I = 1, NPARMS
    1604. 

  1604.                INIBL( I ) = 1
    1605. 

  1605.   590       CONTINUE
    1606. 

  1606.          END IF
    1607. 

  1607. *
    1608. 

  1608. *        Read the values for ISHFTS.
    1609. 

  1609. *
    1610. 

  1610.          IF( NEP ) THEN
    1611. 

  1611.             READ( NIN, FMT = * )( ISHFTS( I ), I = 1, NPARMS )
    1612. 

  1612.             DO 600 I = 1, NPARMS
    1613. 

  1613.                IF( ISHFTS( I ).LT.0 ) THEN
    1614. 

  1614.                   WRITE( NOUT, FMT = 9989 )' ISHFTS ', ISHFTS( I ), 0
    1615. 

  1615.                   FATAL = .TRUE.
    1616. 

  1616.                END IF
    1617. 

  1617.   600       CONTINUE
    1618. 

  1618.             WRITE( NOUT, FMT = 9983 )'ISHFTS: ',
    1619. 

  1619.      $         ( ISHFTS( I ), I = 1, NPARMS )
    1620. 

  1620.          ELSE
    1621. 

  1621.             DO 610 I = 1, NPARMS
    1622. 

  1622.                ISHFTS( I ) = 1
    1623. 

  1623.   610       CONTINUE
    1624. 

  1624.          END IF
    1625. 

  1625. *
    1626. 

  1626. *        Read the values for IACC22.
    1627. 

  1627. *
    1628. 

  1628.          IF( NEP ) THEN
    1629. 

  1629.             READ( NIN, FMT = * )( IACC22( I ), I = 1, NPARMS )
    1630. 

  1630.             DO 620 I = 1, NPARMS
    1631. 

  1631.                IF( IACC22( I ).LT.0 ) THEN
    1632. 

  1632.                   WRITE( NOUT, FMT = 9989 )' IACC22 ', IACC22( I ), 0
    1633. 

  1633.                   FATAL = .TRUE.
    1634. 

  1634.                END IF
    1635. 

  1635.   620       CONTINUE
    1636. 

  1636.             WRITE( NOUT, FMT = 9983 )'IACC22: ',
    1637. 

  1637.      $         ( IACC22( I ), I = 1, NPARMS )
    1638. 

  1638.          ELSE
    1639. 

  1639.             DO 630 I = 1, NPARMS
    1640. 

  1640.                IACC22( I ) = 1
    1641. 

  1641.   630       CONTINUE
    1642. 

  1642.          END IF
    1643. 

  1643. *
    1644. 

  1644. *        Read the values for NBCOL.
    1645. 

  1645. *
    1646. 

  1646.          IF( DGG ) THEN
    1647. 

  1647.             READ( NIN, FMT = * )( NBCOL( I ), I = 1, NPARMS )
    1648. 

  1648.             DO 160 I = 1, NPARMS
    1649. 

  1649.                IF( NBCOL( I ).LT.0 ) THEN
    1650. 

  1650.                   WRITE( NOUT, FMT = 9989 )'NBCOL ', NBCOL( I ), 0
    1651. 

  1651.                   FATAL = .TRUE.
    1652. 

  1652.                ELSE IF( NBCOL( I ).GT.NMAX ) THEN
    1653. 

  1653.                   WRITE( NOUT, FMT = 9988 )'NBCOL ', NBCOL( I ), NMAX
    1654. 

  1654.                   FATAL = .TRUE.
    1655. 

  1655.                END IF
    1656. 

  1656.   160       CONTINUE
    1657. 

  1657.             WRITE( NOUT, FMT = 9983 )'NBCOL:',
    1658. 

  1658.      $         ( NBCOL( I ), I = 1, NPARMS )
    1659. 

  1659.          ELSE
    1660. 

  1660.             DO 170 I = 1, NPARMS
    1661. 

  1661.                NBCOL( I ) = 1
    1662. 

  1662.   170       CONTINUE
    1663. 

  1663.          END IF
    1664. 

  1664.       END IF
    1665. 

  1665. *
    1666. 

  1666. *     Calculate and print the machine dependent constants.
    1667. 

  1667. *
    1668. 

  1668.       WRITE( NOUT, FMT = * )
    1669. 

  1669.       EPS = DLAMCH( 'Underflow threshold' )
    1670. 

  1670.       WRITE( NOUT, FMT = 9981 )'underflow', EPS
    1671. 

  1671.       EPS = DLAMCH( 'Overflow threshold' )
    1672. 

  1672.       WRITE( NOUT, FMT = 9981 )'overflow ', EPS
    1673. 

  1673.       EPS = DLAMCH( 'Epsilon' )
    1674. 

  1674.       WRITE( NOUT, FMT = 9981 )'precision', EPS
    1675. 

  1675. *
    1676. 

  1676. *     Read the threshold value for the test ratios.
    1677. 

  1677. *
    1678. 

  1678.       READ( NIN, FMT = * )THRESH
    1679. 

  1679.       WRITE( NOUT, FMT = 9982 )THRESH
    1680. 

  1680.       IF( SEP .OR. SVD .OR. DGG ) THEN
    1681. 

  1681. *
    1682. 

  1682. *        Read the flag that indicates whether to test LAPACK routines.
    1683. 

  1683. *
    1684. 

  1684.          READ( NIN, FMT = * )TSTCHK
    1685. 

  1685. *
    1686. 

  1686. *        Read the flag that indicates whether to test driver routines.
    1687. 

  1687. *
    1688. 

  1688.          READ( NIN, FMT = * )TSTDRV
    1689. 

  1689.       END IF
    1690. 

  1690. *
    1691. 

  1691. *     Read the flag that indicates whether to test the error exits.
    1692. 

  1692. *
    1693. 

  1693.       READ( NIN, FMT = * )TSTERR
    1694. 

  1694. *
    1695. 

  1695. *     Read the code describing how to set the random number seed.
    1696. 

  1696. *
    1697. 

  1697.       READ( NIN, FMT = * )NEWSD
    1698. 

  1698. *
    1699. 

  1699. *     If NEWSD = 2, read another line with 4 integers for the seed.
    1700. 

  1700. *
    1701. 

  1701.       IF( NEWSD.EQ.2 )
    1702. 

  1702.      $   READ( NIN, FMT = * )( IOLDSD( I ), I = 1, 4 )
    1703. 

  1703. *
    1704. 

  1704.       DO 180 I = 1, 4
    1705. 

  1705.          ISEED( I ) = IOLDSD( I )
    1706. 

  1706.   180 CONTINUE
    1707. 

  1707. *
    1708. 

  1708.       IF( FATAL ) THEN
    1709. 

  1709.          WRITE( NOUT, FMT = 9999 )
    1710. 

  1710.          STOP
    1711. 

  1711.       END IF
    1712. 

  1712. *
    1713. 

  1713. *     Read the input lines indicating the test path and its parameters.
    1714. 

  1714. *     The first three characters indicate the test path, and the number
    1715. 

  1715. *     of test matrix types must be the first nonblank item in columns
    1716. 

  1716. *     4-80.
    1717. 

  1717. *
    1718. 

  1718.   190 CONTINUE
    1719. 

  1719. *
    1720. 

  1720.       IF( .NOT.( DGX .OR. DXV ) ) THEN
    1721. 

  1721. *
    1722. 

  1722.   200    CONTINUE
    1723. 

  1723.          READ( NIN, FMT = '(A80)', END = 380 )LINE
    1724. 

  1724.          C3 = LINE( 1: 3 )
    1725. 

  1725.          LENP = LEN( LINE )
    1726. 

  1726.          I = 3
    1727. 

  1727.          ITMP = 0
    1728. 

  1728.          I1 = 0
    1729. 

  1729.   210    CONTINUE
    1730. 

  1730.          I = I + 1
    1731. 

  1731.          IF( I.GT.LENP ) THEN
    1732. 

  1732.             IF( I1.GT.0 ) THEN
    1733. 

  1733.                GO TO 240
    1734. 

  1734.             ELSE
    1735. 

  1735.                NTYPES = MAXT
    1736. 

  1736.                GO TO 240
    1737. 

  1737.             END IF
    1738. 

  1738.          END IF
    1739. 

  1739.          IF( LINE( I: I ).NE.' ' .AND. LINE( I: I ).NE.',' ) THEN
    1740. 

  1740.             I1 = I
    1741. 

  1741.             C1 = LINE( I1: I1 )
    1742. 

  1742. *
    1743. 

  1743. *        Check that a valid integer was read
    1744. 

  1744. *
    1745. 

  1745.             DO 220 K = 1, 10
    1746. 

  1746.                IF( C1.EQ.INTSTR( K: K ) ) THEN
    1747. 

  1747.                   IC = K - 1
    1748. 

  1748.                   GO TO 230
    1749. 

  1749.                END IF
    1750. 

  1750.   220       CONTINUE
    1751. 

  1751.             WRITE( NOUT, FMT = 9991 )I, LINE
    1752. 

  1752.             GO TO 200
    1753. 

  1753.   230       CONTINUE
    1754. 

  1754.             ITMP = 10*ITMP + IC
    1755. 

  1755.             GO TO 210
    1756. 

  1756.          ELSE IF( I1.GT.0 ) THEN
    1757. 

  1757.             GO TO 240
    1758. 

  1758.          ELSE
    1759. 

  1759.             GO TO 210
    1760. 

  1760.          END IF
    1761. 

  1761.   240    CONTINUE
    1762. 

  1762.          NTYPES = ITMP
    1763. 

  1763. *
    1764. 

  1764. *     Skip the tests if NTYPES is <= 0.
    1765. 

  1765. *
    1766. 

  1766.          IF( .NOT.( DEV .OR. DES .OR. DVX .OR. DSX .OR. DGV .OR.
    1767. 

  1767.      $       DGS ) .AND. NTYPES.LE.0 ) THEN
    1768. 

  1768.             WRITE( NOUT, FMT = 9990 )C3
    1769. 

  1769.             GO TO 200
    1770. 

  1770.          END IF
    1771. 

  1771. *
    1772. 

  1772.       ELSE
    1773. 

  1773.          IF( DXV )
    1774. 

  1774.      $      C3 = 'DXV'
    1775. 

  1775.          IF( DGX )
    1776. 

  1776.      $      C3 = 'DGX'
    1777. 

  1777.       END IF
    1778. 

  1778. *
    1779. 

  1779. *     Reset the random number seed.
    1780. 

  1780. *
    1781. 

  1781.       IF( NEWSD.EQ.0 ) THEN
    1782. 

  1782.          DO 250 K = 1, 4
    1783. 

  1783.             ISEED( K ) = IOLDSD( K )
    1784. 

  1784.   250    CONTINUE
    1785. 

  1785.       END IF
    1786. 

  1786. *
    1787. 

  1787.       IF( LSAMEN( 3, C3, 'DHS' ) .OR. LSAMEN( 3, C3, 'NEP' ) ) THEN
    1788. 

  1788. *
    1789. 

  1789. *        -------------------------------------
    1790. 

  1790. *        NEP:  Nonsymmetric Eigenvalue Problem
    1791. 

  1791. *        -------------------------------------
    1792. 

  1792. *        Vary the parameters
    1793. 

  1793. *           NB    = block size
    1794. 

  1794. *           NBMIN = minimum block size
    1795. 

  1795. *           NX    = crossover point
    1796. 

  1796. *           NS    = number of shifts
    1797. 

  1797. *           MAXB  = minimum submatrix size
    1798. 

  1798. *
    1799. 

  1799.          MAXTYP = 21
    1800. 

  1800.          NTYPES = MIN( MAXTYP, NTYPES )
    1801. 

  1801.          CALL ALAREQ( C3, NTYPES, DOTYPE, MAXTYP, NIN, NOUT )
    1802. 

  1802.          CALL XLAENV( 1, 1 )
    1803. 

  1803.          IF( TSTERR )
    1804. 

  1804.      $      CALL DERRHS( 'DHSEQR', NOUT )
    1805. 

  1805.          DO 270 I = 1, NPARMS
    1806. 

  1806.             CALL XLAENV( 1, NBVAL( I ) )
    1807. 

  1807.             CALL XLAENV( 2, NBMIN( I ) )
    1808. 

  1808.             CALL XLAENV( 3, NXVAL( I ) )
    1809. 

  1809.             CALL XLAENV(12, MAX( 11, INMIN( I ) ) )
    1810. 

  1810.             CALL XLAENV(13, INWIN( I ) )
    1811. 

  1811.             CALL XLAENV(14, INIBL( I ) )
    1812. 

  1812.             CALL XLAENV(15, ISHFTS( I ) )
    1813. 

  1813.             CALL XLAENV(16, IACC22( I ) )
    1814. 

  1814. *
    1815. 

  1815.             IF( NEWSD.EQ.0 ) THEN
    1816. 

  1816.                DO 260 K = 1, 4
    1817. 

  1817.                   ISEED( K ) = IOLDSD( K )
    1818. 

  1818.   260          CONTINUE
    1819. 

  1819.             END IF
    1820. 

  1820.             WRITE( NOUT, FMT = 9961 )C3, NBVAL( I ), NBMIN( I ),
    1821. 

  1821.      $         NXVAL( I ), MAX( 11, INMIN(I)),
    1822. 

  1822.      $         INWIN( I ), INIBL( I ), ISHFTS( I ), IACC22( I )
    1823. 

  1823.             CALL DCHKHS( NN, NVAL, MAXTYP, DOTYPE, ISEED, THRESH, NOUT,
    1824. 

  1824.      $                   A( 1, 1 ), NMAX, A( 1, 2 ), A( 1, 3 ),
    1825. 

  1825.      $                   A( 1, 4 ), A( 1, 5 ), NMAX, A( 1, 6 ),
    1826. 

  1826.      $                   A( 1, 7 ), D( 1, 1 ), D( 1, 2 ), D( 1, 3 ),
    1827. 

  1827.      $                   D( 1, 4 ), A( 1, 8 ), A( 1, 9 ), A( 1, 10 ),
    1828. 

  1828.      $                   A( 1, 11 ), A( 1, 12 ), D( 1, 5 ), WORK, LWORK,
    1829. 

  1829.      $                   IWORK, LOGWRK, RESULT, INFO )
    1830. 

  1830.             IF( INFO.NE.0 )
    1831. 

  1831.      $         WRITE( NOUT, FMT = 9980 )'DCHKHS', INFO
    1832. 

  1832.   270    CONTINUE
    1833. 

  1833. *
    1834. 

  1834.       ELSE IF( LSAMEN( 3, C3, 'DST' ) .OR. LSAMEN( 3, C3, 'SEP' ) ) THEN
    1835. 

  1835. *
    1836. 

  1836. *        ----------------------------------
    1837. 

  1837. *        SEP:  Symmetric Eigenvalue Problem
    1838. 

  1838. *        ----------------------------------
    1839. 

  1839. *        Vary the parameters
    1840. 

  1840. *           NB    = block size
    1841. 

  1841. *           NBMIN = minimum block size
    1842. 

  1842. *           NX    = crossover point
    1843. 

  1843. *
    1844. 

  1844.          MAXTYP = 21
    1845. 

  1845.          NTYPES = MIN( MAXTYP, NTYPES )
    1846. 

  1846.          CALL ALAREQ( C3, NTYPES, DOTYPE, MAXTYP, NIN, NOUT )
    1847. 

  1847.          CALL XLAENV( 1, 1 )
    1848. 

  1848.          CALL XLAENV( 9, 25 )
    1849. 

  1849.          IF( TSTERR )
    1850. 

  1850.      $      CALL DERRST( 'DST', NOUT )
    1851. 

  1851.          DO 290 I = 1, NPARMS
    1852. 

  1852.             CALL XLAENV( 1, NBVAL( I ) )
    1853. 

  1853.             CALL XLAENV( 2, NBMIN( I ) )
    1854. 

  1854.             CALL XLAENV( 3, NXVAL( I ) )
    1855. 

  1855. *
    1856. 

  1856.             IF( NEWSD.EQ.0 ) THEN
    1857. 

  1857.                DO 280 K = 1, 4
    1858. 

  1858.                   ISEED( K ) = IOLDSD( K )
    1859. 

  1859.   280          CONTINUE
    1860. 

  1860.             END IF
    1861. 

  1861.             WRITE( NOUT, FMT = 9997 )C3, NBVAL( I ), NBMIN( I ),
    1862. 

  1862.      $         NXVAL( I )
    1863. 

  1863.             IF( TSTCHK ) THEN
    1864. 

  1864.                CALL DCHKST( NN, NVAL, MAXTYP, DOTYPE, ISEED, THRESH,
    1865. 

  1865.      $                      NOUT, A( 1, 1 ), NMAX, A( 1, 2 ), D( 1, 1 ),
    1866. 

  1866.      $                      D( 1, 2 ), D( 1, 3 ), D( 1, 4 ), D( 1, 5 ),
    1867. 

  1867.      $                      D( 1, 6 ), D( 1, 7 ), D( 1, 8 ), D( 1, 9 ),
    1868. 

  1868.      $                      D( 1, 10 ), D( 1, 11 ), A( 1, 3 ), NMAX,
    1869. 

  1869.      $                      A( 1, 4 ), A( 1, 5 ), D( 1, 12 ), A( 1, 6 ),
    1870. 

  1870.      $                      WORK, LWORK, IWORK, LIWORK, RESULT, INFO )
    1871. 

  1871.                IF( INFO.NE.0 )
    1872. 

  1872.      $            WRITE( NOUT, FMT = 9980 )'DCHKST', INFO
    1873. 

  1873.             END IF
    1874. 

  1874.             IF( TSTDRV ) THEN
    1875. 

  1875.                CALL DDRVST( NN, NVAL, 18, DOTYPE, ISEED, THRESH, NOUT,
    1876. 

  1876.      $                      A( 1, 1 ), NMAX, D( 1, 3 ), D( 1, 4 ),
    1877. 

  1877.      $                      D( 1, 5 ), D( 1, 6 ), D( 1, 8 ), D( 1, 9 ),
    1878. 

  1878.      $                      D( 1, 10 ), D( 1, 11 ), A( 1, 2 ), NMAX,
    1879. 

  1879.      $                      A( 1, 3 ), D( 1, 12 ), A( 1, 4 ), WORK,
    1880. 

  1880.      $                      LWORK, IWORK, LIWORK, RESULT, INFO )
    1881. 

  1881.                IF( INFO.NE.0 )
    1882. 

  1882.      $            WRITE( NOUT, FMT = 9980 )'DDRVST', INFO
    1883. 

  1883.             END IF
    1884. 

  1884.   290    CONTINUE
    1885. 

  1885. *
    1886. 

  1886.       ELSE IF( LSAMEN( 3, C3, 'DSG' ) ) THEN
    1887. 

  1887. *
    1888. 

  1888. *        ----------------------------------------------
    1889. 

  1889. *        DSG:  Symmetric Generalized Eigenvalue Problem
    1890. 

  1890. *        ----------------------------------------------
    1891. 

  1891. *        Vary the parameters
    1892. 

  1892. *           NB    = block size
    1893. 

  1893. *           NBMIN = minimum block size
    1894. 

  1894. *           NX    = crossover point
    1895. 

  1895. *
    1896. 

  1896.          MAXTYP = 21
    1897. 

  1897.          NTYPES = MIN( MAXTYP, NTYPES )
    1898. 

  1898.          CALL ALAREQ( C3, NTYPES, DOTYPE, MAXTYP, NIN, NOUT )
    1899. 

  1899.          CALL XLAENV( 9, 25 )
    1900. 

  1900.          DO 310 I = 1, NPARMS
    1901. 

  1901.             CALL XLAENV( 1, NBVAL( I ) )
    1902. 

  1902.             CALL XLAENV( 2, NBMIN( I ) )
    1903. 

  1903.             CALL XLAENV( 3, NXVAL( I ) )
    1904. 

  1904. *
    1905. 

  1905.             IF( NEWSD.EQ.0 ) THEN
    1906. 

  1906.                DO 300 K = 1, 4
    1907. 

  1907.                   ISEED( K ) = IOLDSD( K )
    1908. 

  1908.   300          CONTINUE
    1909. 

  1909.             END IF
    1910. 

  1910.             WRITE( NOUT, FMT = 9997 )C3, NBVAL( I ), NBMIN( I ),
    1911. 

  1911.      $         NXVAL( I )
    1912. 

  1912.             IF( TSTCHK ) THEN
    1913. 

  1913.                CALL DDRVSG( NN, NVAL, MAXTYP, DOTYPE, ISEED, THRESH,
    1914. 

  1914.      $                      NOUT, A( 1, 1 ), NMAX, A( 1, 2 ), NMAX,
    1915. 

  1915.      $                      D( 1, 3 ), A( 1, 3 ), NMAX, A( 1, 4 ),
    1916. 

  1916.      $                      A( 1, 5 ), A( 1, 6 ), A( 1, 7 ), WORK,
    1917. 

  1917.      $                      LWORK, IWORK, LIWORK, RESULT, INFO )
    1918. 

  1918.                IF( INFO.NE.0 )
    1919. 

  1919.      $            WRITE( NOUT, FMT = 9980 )'DDRVSG', INFO
    1920. 

  1920.             END IF
    1921. 

  1921.   310    CONTINUE
    1922. 

  1922. *
    1923. 

  1923.       ELSE IF( LSAMEN( 3, C3, 'DBD' ) .OR. LSAMEN( 3, C3, 'SVD' ) ) THEN
    1924. 

  1924. *
    1925. 

  1925. *        ----------------------------------
    1926. 

  1926. *        SVD:  Singular Value Decomposition
    1927. 

  1927. *        ----------------------------------
    1928. 

  1928. *        Vary the parameters
    1929. 

  1929. *           NB    = block size
    1930. 

  1930. *           NBMIN = minimum block size
    1931. 

  1931. *           NX    = crossover point
    1932. 

  1932. *           NRHS  = number of right hand sides
    1933. 

  1933. *
    1934. 

  1934.          MAXTYP = 16
    1935. 

  1935.          NTYPES = MIN( MAXTYP, NTYPES )
    1936. 

  1936.          CALL ALAREQ( C3, NTYPES, DOTYPE, MAXTYP, NIN, NOUT )
    1937. 

  1937.          CALL XLAENV( 1, 1 )
    1938. 

  1938.          CALL XLAENV( 9, 25 )
    1939. 

  1939. *
    1940. 

  1940. *        Test the error exits
    1941. 

  1941. *
    1942. 

  1942.          IF( TSTERR .AND. TSTCHK )
    1943. 

  1943.      $      CALL DERRBD( 'DBD', NOUT )
    1944. 

  1944.          IF( TSTERR .AND. TSTDRV )
    1945. 

  1945.      $      CALL DERRED( 'DBD', NOUT )
    1946. 

  1946. *
    1947. 

  1947.          DO 330 I = 1, NPARMS
    1948. 

  1948.             NRHS = NSVAL( I )
    1949. 

  1949.             CALL XLAENV( 1, NBVAL( I ) )
    1950. 

  1950.             CALL XLAENV( 2, NBMIN( I ) )
    1951. 

  1951.             CALL XLAENV( 3, NXVAL( I ) )
    1952. 

  1952.             IF( NEWSD.EQ.0 ) THEN
    1953. 

  1953.                DO 320 K = 1, 4
    1954. 

  1954.                   ISEED( K ) = IOLDSD( K )
    1955. 

  1955.   320          CONTINUE
    1956. 

  1956.             END IF
    1957. 

  1957.             WRITE( NOUT, FMT = 9995 )C3, NBVAL( I ), NBMIN( I ),
    1958. 

  1958.      $         NXVAL( I ), NRHS
    1959. 

  1959.             IF( TSTCHK ) THEN
    1960. 

  1960.                CALL DCHKBD( NN, MVAL, NVAL, MAXTYP, DOTYPE, NRHS, ISEED,
    1961. 

  1961.      $                      THRESH, A( 1, 1 ), NMAX, D( 1, 1 ),
    1962. 

  1962.      $                      D( 1, 2 ), D( 1, 3 ), D( 1, 4 ), A( 1, 2 ),
    1963. 

  1963.      $                      NMAX, A( 1, 3 ), A( 1, 4 ), A( 1, 5 ), NMAX,
    1964. 

  1964.      $                      A( 1, 6 ), NMAX, A( 1, 7 ), A( 1, 8 ), WORK,
    1965. 

  1965.      $                      LWORK, IWORK, NOUT, INFO )
    1966. 

  1966.                IF( INFO.NE.0 )
    1967. 

  1967.      $            WRITE( NOUT, FMT = 9980 )'DCHKBD', INFO
    1968. 

  1968.             END IF
    1969. 

  1969.             IF( TSTDRV )
    1970. 

  1970.      $         CALL DDRVBD( NN, MVAL, NVAL, MAXTYP, DOTYPE, ISEED,
    1971. 

  1971.      $                      THRESH, A( 1, 1 ), NMAX, A( 1, 2 ), NMAX,
    1972. 

  1972.      $                      A( 1, 3 ), NMAX, A( 1, 4 ), A( 1, 5 ),
    1973. 

  1973.      $                      A( 1, 6 ), D( 1, 1 ), D( 1, 2 ), D( 1, 3 ),
    1974. 

  1974.      $                      WORK, LWORK, IWORK, NOUT, INFO )
    1975. 

  1975.   330    CONTINUE
    1976. 

  1976. *
    1977. 

  1977.       ELSE IF( LSAMEN( 3, C3, 'DEV' ) ) THEN
    1978. 

  1978. *
    1979. 

  1979. *        --------------------------------------------
    1980. 

  1980. *        DEV:  Nonsymmetric Eigenvalue Problem Driver
    1981. 

  1981. *              DGEEV (eigenvalues and eigenvectors)
    1982. 

  1982. *        --------------------------------------------
    1983. 

  1983. *
    1984. 

  1984.          MAXTYP = 21
    1985. 

  1985.          NTYPES = MIN( MAXTYP, NTYPES )
    1986. 

  1986.          IF( NTYPES.LE.0 ) THEN
    1987. 

  1987.             WRITE( NOUT, FMT = 9990 )C3
    1988. 

  1988.          ELSE
    1989. 

  1989.             IF( TSTERR )
    1990. 

  1990.      $         CALL DERRED( C3, NOUT )
    1991. 

  1991.             CALL ALAREQ( C3, NTYPES, DOTYPE, MAXTYP, NIN, NOUT )
    1992. 

  1992.             CALL DDRVEV( NN, NVAL, NTYPES, DOTYPE, ISEED, THRESH, NOUT,
    1993. 

  1993.      $                   A( 1, 1 ), NMAX, A( 1, 2 ), D( 1, 1 ),
    1994. 

  1994.      $                   D( 1, 2 ), D( 1, 3 ), D( 1, 4 ), A( 1, 3 ),
    1995. 

  1995.      $                   NMAX, A( 1, 4 ), NMAX, A( 1, 5 ), NMAX, RESULT,
    1996. 

  1996.      $                   WORK, LWORK, IWORK, INFO )
    1997. 

  1997.             IF( INFO.NE.0 )
    1998. 

  1998.      $         WRITE( NOUT, FMT = 9980 )'DGEEV', INFO
    1999. 

  1999.          END IF
    2000. 

  2000.          WRITE( NOUT, FMT = 9973 )
    2001. 

  2001.          GO TO 10
    2002. 

  2002. *
    2003. 

  2003.       ELSE IF( LSAMEN( 3, C3, 'DES' ) ) THEN
    2004. 

  2004. *
    2005. 

  2005. *        --------------------------------------------
    2006. 

  2006. *        DES:  Nonsymmetric Eigenvalue Problem Driver
    2007. 

  2007. *              DGEES (Schur form)
    2008. 

  2008. *        --------------------------------------------
    2009. 

  2009. *
    2010. 

  2010.          MAXTYP = 21
    2011. 

  2011.          NTYPES = MIN( MAXTYP, NTYPES )
    2012. 

  2012.          IF( NTYPES.LE.0 ) THEN
    2013. 

  2013.             WRITE( NOUT, FMT = 9990 )C3
    2014. 

  2014.          ELSE
    2015. 

  2015.             IF( TSTERR )
    2016. 

  2016.      $         CALL DERRED( C3, NOUT )
    2017. 

  2017.             CALL ALAREQ( C3, NTYPES, DOTYPE, MAXTYP, NIN, NOUT )
    2018. 

  2018.             CALL DDRVES( NN, NVAL, NTYPES, DOTYPE, ISEED, THRESH, NOUT,
    2019. 

  2019.      $                   A( 1, 1 ), NMAX, A( 1, 2 ), A( 1, 3 ),
    2020. 

  2020.      $                   D( 1, 1 ), D( 1, 2 ), D( 1, 3 ), D( 1, 4 ),
    2021. 

  2021.      $                   A( 1, 4 ), NMAX, RESULT, WORK, LWORK, IWORK,
    2022. 

  2022.      $                   LOGWRK, INFO )
    2023. 

  2023.             IF( INFO.NE.0 )
    2024. 

  2024.      $         WRITE( NOUT, FMT = 9980 )'DGEES', INFO
    2025. 

  2025.          END IF
    2026. 

  2026.          WRITE( NOUT, FMT = 9973 )
    2027. 

  2027.          GO TO 10
    2028. 

  2028. *
    2029. 

  2029.       ELSE IF( LSAMEN( 3, C3, 'DVX' ) ) THEN
    2030. 

  2030. *
    2031. 

  2031. *        --------------------------------------------------------------
    2032. 

  2032. *        DVX:  Nonsymmetric Eigenvalue Problem Expert Driver
    2033. 

  2033. *              DGEEVX (eigenvalues, eigenvectors and condition numbers)
    2034. 

  2034. *        --------------------------------------------------------------
    2035. 

  2035. *
    2036. 

  2036.          MAXTYP = 21
    2037. 

  2037.          NTYPES = MIN( MAXTYP, NTYPES )
    2038. 

  2038.          IF( NTYPES.LT.0 ) THEN
    2039. 

  2039.             WRITE( NOUT, FMT = 9990 )C3
    2040. 

  2040.          ELSE
    2041. 

  2041.             IF( TSTERR )
    2042. 

  2042.      $         CALL DERRED( C3, NOUT )
    2043. 

  2043.             CALL ALAREQ( C3, NTYPES, DOTYPE, MAXTYP, NIN, NOUT )
    2044. 

  2044.             CALL DDRVVX( NN, NVAL, NTYPES, DOTYPE, ISEED, THRESH, NIN,
    2045. 

  2045.      $                   NOUT, A( 1, 1 ), NMAX, A( 1, 2 ), D( 1, 1 ),
    2046. 

  2046.      $                   D( 1, 2 ), D( 1, 3 ), D( 1, 4 ), A( 1, 3 ),
    2047. 

  2047.      $                   NMAX, A( 1, 4 ), NMAX, A( 1, 5 ), NMAX,
    2048. 

  2048.      $                   D( 1, 5 ), D( 1, 6 ), D( 1, 7 ), D( 1, 8 ),
    2049. 

  2049.      $                   D( 1, 9 ), D( 1, 10 ), D( 1, 11 ), D( 1, 12 ),
    2050. 

  2050.      $                   RESULT, WORK, LWORK, IWORK, INFO )
    2051. 

  2051.             IF( INFO.NE.0 )
    2052. 

  2052.      $         WRITE( NOUT, FMT = 9980 )'DGEEVX', INFO
    2053. 

  2053.          END IF
    2054. 

  2054.          WRITE( NOUT, FMT = 9973 )
    2055. 

  2055.          GO TO 10
    2056. 

  2056. *
    2057. 

  2057.       ELSE IF( LSAMEN( 3, C3, 'DSX' ) ) THEN
    2058. 

  2058. *
    2059. 

  2059. *        ---------------------------------------------------
    2060. 

  2060. *        DSX:  Nonsymmetric Eigenvalue Problem Expert Driver
    2061. 

  2061. *              DGEESX (Schur form and condition numbers)
    2062. 

  2062. *        ---------------------------------------------------
    2063. 

  2063. *
    2064. 

  2064.          MAXTYP = 21
    2065. 

  2065.          NTYPES = MIN( MAXTYP, NTYPES )
    2066. 

  2066.          IF( NTYPES.LT.0 ) THEN
    2067. 

  2067.             WRITE( NOUT, FMT = 9990 )C3
    2068. 

  2068.          ELSE
    2069. 

  2069.             IF( TSTERR )
    2070. 

  2070.      $         CALL DERRED( C3, NOUT )
    2071. 

  2071.             CALL ALAREQ( C3, NTYPES, DOTYPE, MAXTYP, NIN, NOUT )
    2072. 

  2072.             CALL DDRVSX( NN, NVAL, NTYPES, DOTYPE, ISEED, THRESH, NIN,
    2073. 

  2073.      $                   NOUT, A( 1, 1 ), NMAX, A( 1, 2 ), A( 1, 3 ),
    2074. 

  2074.      $                   D( 1, 1 ), D( 1, 2 ), D( 1, 3 ), D( 1, 4 ),
    2075. 

  2075.      $                   D( 1, 5 ), D( 1, 6 ), A( 1, 4 ), NMAX,
    2076. 

  2076.      $                   A( 1, 5 ), RESULT, WORK, LWORK, IWORK, LOGWRK,
    2077. 

  2077.      $                   INFO )
    2078. 

  2078.             IF( INFO.NE.0 )
    2079. 

  2079.      $         WRITE( NOUT, FMT = 9980 )'DGEESX', INFO
    2080. 

  2080.          END IF
    2081. 

  2081.          WRITE( NOUT, FMT = 9973 )
    2082. 

  2082.          GO TO 10
    2083. 

  2083. *
    2084. 

  2084.       ELSE IF( LSAMEN( 3, C3, 'DGG' ) ) THEN
    2085. 

  2085. *
    2086. 

  2086. *        -------------------------------------------------
    2087. 

  2087. *        DGG:  Generalized Nonsymmetric Eigenvalue Problem
    2088. 

  2088. *        -------------------------------------------------
    2089. 

  2089. *        Vary the parameters
    2090. 

  2090. *           NB    = block size
    2091. 

  2091. *           NBMIN = minimum block size
    2092. 

  2092. *           NS    = number of shifts
    2093. 

  2093. *           MAXB  = minimum submatrix size
    2094. 

  2094. *           NBCOL = minimum column dimension for blocks
    2095. 

  2095. *
    2096. 

  2096.          MAXTYP = 26
    2097. 

  2097.          NTYPES = MIN( MAXTYP, NTYPES )
    2098. 

  2098.          CALL ALAREQ( C3, NTYPES, DOTYPE, MAXTYP, NIN, NOUT )
    2099. 

  2099.          IF( TSTCHK .AND. TSTERR )
    2100. 

  2100.      $      CALL DERRGG( C3, NOUT )
    2101. 

  2101.          DO 350 I = 1, NPARMS
    2102. 

  2102.             CALL XLAENV( 1, NBVAL( I ) )
    2103. 

  2103.             CALL XLAENV( 2, NBMIN( I ) )
    2104. 

  2104.             CALL XLAENV( 4, NSVAL( I ) )
    2105. 

  2105.             CALL XLAENV( 8, MXBVAL( I ) )
    2106. 

  2106.             CALL XLAENV( 5, NBCOL( I ) )
    2107. 

  2107. *
    2108. 

  2108.             IF( NEWSD.EQ.0 ) THEN
    2109. 

  2109.                DO 340 K = 1, 4
    2110. 

  2110.                   ISEED( K ) = IOLDSD( K )
    2111. 

  2111.   340          CONTINUE
    2112. 

  2112.             END IF
    2113. 

  2113.             WRITE( NOUT, FMT = 9996 )C3, NBVAL( I ), NBMIN( I ),
    2114. 

  2114.      $         NSVAL( I ), MXBVAL( I ), NBCOL( I )
    2115. 

  2115.             TSTDIF = .FALSE.
    2116. 

  2116.             THRSHN = 10.D0
    2117. 

  2117.             IF( TSTCHK ) THEN
    2118. 

  2118.                CALL DCHKGG( NN, NVAL, MAXTYP, DOTYPE, ISEED, THRESH,
    2119. 

  2119.      $                      TSTDIF, THRSHN, NOUT, A( 1, 1 ), NMAX,
    2120. 

  2120.      $                      A( 1, 2 ), A( 1, 3 ), A( 1, 4 ), A( 1, 5 ),
    2121. 

  2121.      $                      A( 1, 6 ), A( 1, 7 ), A( 1, 8 ), A( 1, 9 ),
    2122. 

  2122.      $                      NMAX, A( 1, 10 ), A( 1, 11 ), A( 1, 12 ),
    2123. 

  2123.      $                      D( 1, 1 ), D( 1, 2 ), D( 1, 3 ), D( 1, 4 ),
    2124. 

  2124.      $                      D( 1, 5 ), D( 1, 6 ), A( 1, 13 ),
    2125. 

  2125.      $                      A( 1, 14 ), WORK, LWORK, LOGWRK, RESULT,
    2126. 

  2126.      $                      INFO )
    2127. 

  2127.                IF( INFO.NE.0 )
    2128. 

  2128.      $            WRITE( NOUT, FMT = 9980 )'DCHKGG', INFO
    2129. 

  2129.             END IF
    2130. 

  2130.             CALL XLAENV( 1, 1 )
    2131. 

  2131.             IF( TSTDRV ) THEN
    2132. 

  2132.                CALL DDRVGG( NN, NVAL, MAXTYP, DOTYPE, ISEED, THRESH,
    2133. 

  2133.      $                      THRSHN, NOUT, A( 1, 1 ), NMAX, A( 1, 2 ),
    2134. 

  2134.      $                      A( 1, 3 ), A( 1, 4 ), A( 1, 5 ), A( 1, 6 ),
    2135. 

  2135.      $                      A( 1, 7 ), NMAX, A( 1, 8 ), D( 1, 1 ),
    2136. 

  2136.      $                      D( 1, 2 ), D( 1, 3 ), D( 1, 4 ), D( 1, 5 ),
    2137. 

  2137.      $                      D( 1, 6 ), A( 1, 13 ), A( 1, 14 ), WORK,
    2138. 

  2138.      $                      LWORK, RESULT, INFO )
    2139. 

  2139.                IF( INFO.NE.0 )
    2140. 

  2140.      $            WRITE( NOUT, FMT = 9980 )'DDRVGG', INFO
    2141. 

  2141.             END IF
    2142. 

  2142.   350    CONTINUE
    2143. 

  2143. *
    2144. 

  2144.       ELSE IF( LSAMEN( 3, C3, 'DGS' ) ) THEN
    2145. 

  2145. *
    2146. 

  2146. *        -------------------------------------------------
    2147. 

  2147. *        DGS:  Generalized Nonsymmetric Eigenvalue Problem
    2148. 

  2148. *              DGGES (Schur form)
    2149. 

  2149. *        -------------------------------------------------
    2150. 

  2150. *
    2151. 

  2151.          MAXTYP = 26
    2152. 

  2152.          NTYPES = MIN( MAXTYP, NTYPES )
    2153. 

  2153.          IF( NTYPES.LE.0 ) THEN
    2154. 

  2154.             WRITE( NOUT, FMT = 9990 )C3
    2155. 

  2155.          ELSE
    2156. 

  2156.             IF( TSTERR )
    2157. 

  2157.      $         CALL DERRGG( C3, NOUT )
    2158. 

  2158.             CALL ALAREQ( C3, NTYPES, DOTYPE, MAXTYP, NIN, NOUT )
    2159. 

  2159.             CALL DDRGES( NN, NVAL, MAXTYP, DOTYPE, ISEED, THRESH, NOUT,
    2160. 

  2160.      $                   A( 1, 1 ), NMAX, A( 1, 2 ), A( 1, 3 ),
    2161. 

  2161.      $                   A( 1, 4 ), A( 1, 7 ), NMAX, A( 1, 8 ),
    2162. 

  2162.      $                   D( 1, 1 ), D( 1, 2 ), D( 1, 3 ), WORK, LWORK,
    2163. 

  2163.      $                   RESULT, LOGWRK, INFO )
    2164. 

  2164. *
    2165. 

  2165.             IF( INFO.NE.0 )
    2166. 

  2166.      $         WRITE( NOUT, FMT = 9980 )'DDRGES', INFO
    2167. 

  2167.          END IF
    2168. 

  2168.          WRITE( NOUT, FMT = 9973 )
    2169. 

  2169.          GO TO 10
    2170. 

  2170. *
    2171. 

  2171.       ELSE IF( DGX ) THEN
    2172. 

  2172. *
    2173. 

  2173. *        -------------------------------------------------
    2174. 

  2174. *        DGX:  Generalized Nonsymmetric Eigenvalue Problem
    2175. 

  2175. *              DGGESX (Schur form and condition numbers)
    2176. 

  2176. *        -------------------------------------------------
    2177. 

  2177. *
    2178. 

  2178.          MAXTYP = 5
    2179. 

  2179.          NTYPES = MAXTYP
    2180. 

  2180.          IF( NN.LT.0 ) THEN
    2181. 

  2181.             WRITE( NOUT, FMT = 9990 )C3
    2182. 

  2182.          ELSE
    2183. 

  2183.             IF( TSTERR )
    2184. 

  2184.      $         CALL DERRGG( C3, NOUT )
    2185. 

  2185.             CALL ALAREQ( C3, NTYPES, DOTYPE, MAXTYP, NIN, NOUT )
    2186. 

  2186.             CALL XLAENV( 5, 2 )
    2187. 

  2187.             CALL DDRGSX( NN, NCMAX, THRESH, NIN, NOUT, A( 1, 1 ), NMAX,
    2188. 

  2188.      $                   A( 1, 2 ), A( 1, 3 ), A( 1, 4 ), A( 1, 5 ),
    2189. 

  2189.      $                   A( 1, 6 ), D( 1, 1 ), D( 1, 2 ), D( 1, 3 ),
    2190. 

  2190.      $                   C( 1, 1 ), NCMAX*NCMAX, A( 1, 12 ), WORK,
    2191. 

  2191.      $                   LWORK, IWORK, LIWORK, LOGWRK, INFO )
    2192. 

  2192.             IF( INFO.NE.0 )
    2193. 

  2193.      $         WRITE( NOUT, FMT = 9980 )'DDRGSX', INFO
    2194. 

  2194.          END IF
    2195. 

  2195.          WRITE( NOUT, FMT = 9973 )
    2196. 

  2196.          GO TO 10
    2197. 

  2197. *
    2198. 

  2198.       ELSE IF( LSAMEN( 3, C3, 'DGV' ) ) THEN
    2199. 

  2199. *
    2200. 

  2200. *        -------------------------------------------------
    2201. 

  2201. *        DGV:  Generalized Nonsymmetric Eigenvalue Problem
    2202. 

  2202. *              DGGEV (Eigenvalue/vector form)
    2203. 

  2203. *        -------------------------------------------------
    2204. 

  2204. *
    2205. 

  2205.          MAXTYP = 26
    2206. 

  2206.          NTYPES = MIN( MAXTYP, NTYPES )
    2207. 

  2207.          IF( NTYPES.LE.0 ) THEN
    2208. 

  2208.             WRITE( NOUT, FMT = 9990 )C3
    2209. 

  2209.          ELSE
    2210. 

  2210.             IF( TSTERR )
    2211. 

  2211.      $         CALL DERRGG( C3, NOUT )
    2212. 

  2212.             CALL ALAREQ( C3, NTYPES, DOTYPE, MAXTYP, NIN, NOUT )
    2213. 

  2213.             CALL DDRGEV( NN, NVAL, MAXTYP, DOTYPE, ISEED, THRESH, NOUT,
    2214. 

  2214.      $                   A( 1, 1 ), NMAX, A( 1, 2 ), A( 1, 3 ),
    2215. 

  2215.      $                   A( 1, 4 ), A( 1, 7 ), NMAX, A( 1, 8 ),
    2216. 

  2216.      $                   A( 1, 9 ), NMAX, D( 1, 1 ), D( 1, 2 ),
    2217. 

  2217.      $                   D( 1, 3 ), D( 1, 4 ), D( 1, 5 ), D( 1, 6 ),
    2218. 

  2218.      $                   WORK, LWORK, RESULT, INFO )
    2219. 

  2219.             IF( INFO.NE.0 )
    2220. 

  2220.      $         WRITE( NOUT, FMT = 9980 )'DDRGEV', INFO
    2221. 

  2221.          END IF
    2222. 

  2222.          WRITE( NOUT, FMT = 9973 )
    2223. 

  2223.          GO TO 10
    2224. 

  2224. *
    2225. 

  2225.       ELSE IF( DXV ) THEN
    2226. 

  2226. *
    2227. 

  2227. *        -------------------------------------------------
    2228. 

  2228. *        DXV:  Generalized Nonsymmetric Eigenvalue Problem
    2229. 

  2229. *              DGGEVX (eigenvalue/vector with condition numbers)
    2230. 

  2230. *        -------------------------------------------------
    2231. 

  2231. *
    2232. 

  2232.          MAXTYP = 2
    2233. 

  2233.          NTYPES = MAXTYP
    2234. 

  2234.          IF( NN.LT.0 ) THEN
    2235. 

  2235.             WRITE( NOUT, FMT = 9990 )C3
    2236. 

  2236.          ELSE
    2237. 

  2237.             IF( TSTERR )
    2238. 

  2238.      $         CALL DERRGG( C3, NOUT )
    2239. 

  2239.             CALL ALAREQ( C3, NTYPES, DOTYPE, MAXTYP, NIN, NOUT )
    2240. 

  2240.             CALL DDRGVX( NN, THRESH, NIN, NOUT, A( 1, 1 ), NMAX,
    2241. 

  2241.      $                   A( 1, 2 ), A( 1, 3 ), A( 1, 4 ), D( 1, 1 ),
    2242. 

  2242.      $                   D( 1, 2 ), D( 1, 3 ), A( 1, 5 ), A( 1, 6 ),
    2243. 

  2243.      $                   IWORK( 1 ), IWORK( 2 ), D( 1, 4 ), D( 1, 5 ),
    2244. 

  2244.      $                   D( 1, 6 ), D( 1, 7 ), D( 1, 8 ), D( 1, 9 ),
    2245. 

  2245.      $                   WORK, LWORK, IWORK( 3 ), LIWORK-2, RESULT,
    2246. 

  2246.      $                   LOGWRK, INFO )
    2247. 

  2247. *
    2248. 

  2248.             IF( INFO.NE.0 )
    2249. 

  2249.      $         WRITE( NOUT, FMT = 9980 )'DDRGVX', INFO
    2250. 

  2250.          END IF
    2251. 

  2251.          WRITE( NOUT, FMT = 9973 )
    2252. 

  2252.          GO TO 10
    2253. 

  2253. *
    2254. 

  2254.       ELSE IF( LSAMEN( 3, C3, 'DSB' ) ) THEN
    2255. 

  2255. *
    2256. 

  2256. *        ------------------------------
    2257. 

  2257. *        DSB:  Symmetric Band Reduction
    2258. 

  2258. *        ------------------------------
    2259. 

  2259. *
    2260. 

  2260.          MAXTYP = 15
    2261. 

  2261.          NTYPES = MIN( MAXTYP, NTYPES )
    2262. 

  2262.          CALL ALAREQ( C3, NTYPES, DOTYPE, MAXTYP, NIN, NOUT )
    2263. 

  2263.          IF( TSTERR )
    2264. 

  2264.      $      CALL DERRST( 'DSB', NOUT )
    2265. 

  2265.          CALL DCHKSB( NN, NVAL, NK, KVAL, MAXTYP, DOTYPE, ISEED, THRESH,
    2266. 

  2266.      $                NOUT, A( 1, 1 ), NMAX, D( 1, 1 ), D( 1, 2 ),
    2267. 

  2267.      $                A( 1, 2 ), NMAX, WORK, LWORK, RESULT, INFO )
    2268. 

  2268.          IF( INFO.NE.0 )
    2269. 

  2269.      $      WRITE( NOUT, FMT = 9980 )'DCHKSB', INFO
    2270. 

  2270. *
    2271. 

  2271.       ELSE IF( LSAMEN( 3, C3, 'DBB' ) ) THEN
    2272. 

  2272. *
    2273. 

  2273. *        ------------------------------
    2274. 

  2274. *        DBB:  General Band Reduction
    2275. 

  2275. *        ------------------------------
    2276. 

  2276. *
    2277. 

  2277.          MAXTYP = 15
    2278. 

  2278.          NTYPES = MIN( MAXTYP, NTYPES )
    2279. 

  2279.          CALL ALAREQ( C3, NTYPES, DOTYPE, MAXTYP, NIN, NOUT )
    2280. 

  2280.          DO 370 I = 1, NPARMS
    2281. 

  2281.             NRHS = NSVAL( I )
    2282. 

  2282. *
    2283. 

  2283.             IF( NEWSD.EQ.0 ) THEN
    2284. 

  2284.                DO 360 K = 1, 4
    2285. 

  2285.                   ISEED( K ) = IOLDSD( K )
    2286. 

  2286.   360          CONTINUE
    2287. 

  2287.             END IF
    2288. 

  2288.             WRITE( NOUT, FMT = 9966 )C3, NRHS
    2289. 

  2289.             CALL DCHKBB( NN, MVAL, NVAL, NK, KVAL, MAXTYP, DOTYPE, NRHS,
    2290. 

  2290.      $                   ISEED, THRESH, NOUT, A( 1, 1 ), NMAX,
    2291. 

  2291.      $                   A( 1, 2 ), 2*NMAX, D( 1, 1 ), D( 1, 2 ),
    2292. 

  2292.      $                   A( 1, 4 ), NMAX, A( 1, 5 ), NMAX, A( 1, 6 ),
    2293. 

  2293.      $                   NMAX, A( 1, 7 ), WORK, LWORK, RESULT, INFO )
    2294. 

  2294.             IF( INFO.NE.0 )
    2295. 

  2295.      $         WRITE( NOUT, FMT = 9980 )'DCHKBB', INFO
    2296. 

  2296.   370    CONTINUE
    2297. 

  2297. *
    2298. 

  2298.       ELSE IF( LSAMEN( 3, C3, 'GLM' ) ) THEN
    2299. 

  2299. *
    2300. 

  2300. *        -----------------------------------------
    2301. 

  2301. *        GLM:  Generalized Linear Regression Model
    2302. 

  2302. *        -----------------------------------------
    2303. 

  2303. *
    2304. 

  2304.          CALL XLAENV( 1, 1 )
    2305. 

  2305.          IF( TSTERR )
    2306. 

  2306.      $      CALL DERRGG( 'GLM', NOUT )
    2307. 

  2307.          CALL DCKGLM( NN, MVAL, PVAL, NVAL, NTYPES, ISEED, THRESH, NMAX,
    2308. 

  2308.      $                A( 1, 1 ), A( 1, 2 ), B( 1, 1 ), B( 1, 2 ), X,
    2309. 

  2309.      $                WORK, D( 1, 1 ), NIN, NOUT, INFO )
    2310. 

  2310.          IF( INFO.NE.0 )
    2311. 

  2311.      $      WRITE( NOUT, FMT = 9980 )'DCKGLM', INFO
    2312. 

  2312. *
    2313. 

  2313.       ELSE IF( LSAMEN( 3, C3, 'GQR' ) ) THEN
    2314. 

  2314. *
    2315. 

  2315. *        ------------------------------------------
    2316. 

  2316. *        GQR:  Generalized QR and RQ factorizations
    2317. 

  2317. *        ------------------------------------------
    2318. 

  2318. *
    2319. 

  2319.          CALL XLAENV( 1, 1 )
    2320. 

  2320.          IF( TSTERR )
    2321. 

  2321.      $      CALL DERRGG( 'GQR', NOUT )
    2322. 

  2322.          CALL DCKGQR( NN, MVAL, NN, PVAL, NN, NVAL, NTYPES, ISEED,
    2323. 

  2323.      $                THRESH, NMAX, A( 1, 1 ), A( 1, 2 ), A( 1, 3 ),
    2324. 

  2324.      $                A( 1, 4 ), TAUA, B( 1, 1 ), B( 1, 2 ), B( 1, 3 ),
    2325. 

  2325.      $                B( 1, 4 ), B( 1, 5 ), TAUB, WORK, D( 1, 1 ), NIN,
    2326. 

  2326.      $                NOUT, INFO )
    2327. 

  2327.          IF( INFO.NE.0 )
    2328. 

  2328.      $      WRITE( NOUT, FMT = 9980 )'DCKGQR', INFO
    2329. 

  2329. *
    2330. 

  2330.       ELSE IF( LSAMEN( 3, C3, 'GSV' ) ) THEN
    2331. 

  2331. *
    2332. 

  2332. *        ----------------------------------------------
    2333. 

  2333. *        GSV:  Generalized Singular Value Decomposition
    2334. 

  2334. *        ----------------------------------------------
    2335. 

  2335. *
    2336. 

  2336.          IF( TSTERR )
    2337. 

  2337.      $      CALL DERRGG( 'GSV', NOUT )
    2338. 

  2338.          CALL DCKGSV( NN, MVAL, PVAL, NVAL, NTYPES, ISEED, THRESH, NMAX,
    2339. 

  2339.      $                A( 1, 1 ), A( 1, 2 ), B( 1, 1 ), B( 1, 2 ),
    2340. 

  2340.      $                A( 1, 3 ), B( 1, 3 ), A( 1, 4 ), TAUA, TAUB,
    2341. 

  2341.      $                B( 1, 4 ), IWORK, WORK, D( 1, 1 ), NIN, NOUT,
    2342. 

  2342.      $                INFO )
    2343. 

  2343.          IF( INFO.NE.0 )
    2344. 

  2344.      $      WRITE( NOUT, FMT = 9980 )'DCKGSV', INFO
    2345. 

  2345. *
    2346. 

  2346.       ELSE IF( LSAMEN( 3, C3, 'CSD' ) ) THEN
    2347. 

  2347. *
    2348. 

  2348. *        ----------------------------------------------
    2349. 

  2349. *        CSD:  CS Decomposition
    2350. 

  2350. *        ----------------------------------------------
    2351. 

  2351. *
    2352. 

  2352.          CALL XLAENV(1,1)
    2353. 

  2353.          IF( TSTERR )
    2354. 

  2354.      $      CALL DERRGG( 'CSD', NOUT )
    2355. 

  2355.          CALL DCKCSD( NN, MVAL, PVAL, NVAL, NTYPES, ISEED, THRESH, NMAX,
    2356. 

  2356.      $                A( 1, 1 ), A( 1, 2 ), A( 1, 3 ), A( 1, 4 ),
    2357. 

  2357.      $                A( 1, 5 ), A( 1, 6 ), A( 1, 7 ), IWORK, WORK,
    2358. 

  2358.      $                D( 1, 1 ), NIN, NOUT, INFO )
    2359. 

  2359.          IF( INFO.NE.0 )
    2360. 

  2360.      $      WRITE( NOUT, FMT = 9980 )'DCKCSD', INFO
    2361. 

  2361. *
    2362. 

  2362.       ELSE IF( LSAMEN( 3, C3, 'LSE' ) ) THEN
    2363. 

  2363. *
    2364. 

  2364. *        --------------------------------------
    2365. 

  2365. *        LSE:  Constrained Linear Least Squares
    2366. 

  2366. *        --------------------------------------
    2367. 

  2367. *
    2368. 

  2368.          CALL XLAENV( 1, 1 )
    2369. 

  2369.          IF( TSTERR )
    2370. 

  2370.      $      CALL DERRGG( 'LSE', NOUT )
    2371. 

  2371.          CALL DCKLSE( NN, MVAL, PVAL, NVAL, NTYPES, ISEED, THRESH, NMAX,
    2372. 

  2372.      $                A( 1, 1 ), A( 1, 2 ), B( 1, 1 ), B( 1, 2 ), X,
    2373. 

  2373.      $                WORK, D( 1, 1 ), NIN, NOUT, INFO )
    2374. 

  2374.          IF( INFO.NE.0 )
    2375. 

  2375.      $      WRITE( NOUT, FMT = 9980 )'DCKLSE', INFO
    2376. 

  2376. *
    2377. 

  2377.       ELSE
    2378. 

  2378.          WRITE( NOUT, FMT = * )
    2379. 

  2379.          WRITE( NOUT, FMT = * )
    2380. 

  2380.          WRITE( NOUT, FMT = 9992 )C3
    2381. 

  2381.       END IF
    2382. 

  2382.       IF( .NOT.( DGX .OR. DXV ) )
    2383. 

  2383.      $   GO TO 190
    2384. 

  2384.   380 CONTINUE
    2385. 

  2385.       WRITE( NOUT, FMT = 9994 )
    2386. 

  2386.       S2 = DSECND( )
    2387. 

  2387.       WRITE( NOUT, FMT = 9993 )S2 - S1
    2388. 

  2388. *
    2389. 

  2389.  9999 FORMAT( / ' Execution not attempted due to input errors' )
    2390. 

  2390.  9997 FORMAT( / / 1X, A3, ':  NB =', I4, ', NBMIN =', I4, ', NX =', I4 )
    2391. 

  2391.  9996 FORMAT( / / 1X, A3, ':  NB =', I4, ', NBMIN =', I4, ', NS =', I4,
    2392. 

  2392.      $      ', MAXB =', I4, ', NBCOL =', I4 )
    2393. 

  2393.  9995 FORMAT( / / 1X, A3, ':  NB =', I4, ', NBMIN =', I4, ', NX =', I4,
    2394. 

  2394.      $      ', NRHS =', I4 )
    2395. 

  2395.  9994 FORMAT( / / ' End of tests' )
    2396. 

  2396.  9993 FORMAT( ' Total time used = ', F12.2, ' seconds', / )
    2397. 

  2397.  9992 FORMAT( 1X, A3, ':  Unrecognized path name' )
    2398. 

  2398.  9991 FORMAT( / / ' *** Invalid integer value in column ', I2,
    2399. 

  2399.      $      ' of input', ' line:', / A79 )
    2400. 

  2400.  9990 FORMAT( / / 1X, A3, ' routines were not tested' )
    2401. 

  2401.  9989 FORMAT( ' Invalid input value: ', A, '=', I6, '; must be >=',
    2402. 

  2402.      $      I6 )
    2403. 

  2403.  9988 FORMAT( ' Invalid input value: ', A, '=', I6, '; must be <=',
    2404. 

  2404.      $      I6 )
    2405. 

  2405.  9987 FORMAT( ' Tests of the Nonsymmetric Eigenvalue Problem routines' )
    2406. 

  2406.  9986 FORMAT( ' Tests of the Symmetric Eigenvalue Problem routines' )
    2407. 

  2407.  9985 FORMAT( ' Tests of the Singular Value Decomposition routines' )
    2408. 

  2408.  9984 FORMAT( / ' The following parameter values will be used:' )
    2409. 

  2409.  9983 FORMAT( 4X, A, 10I6, / 10X, 10I6 )
    2410. 

  2410.  9982 FORMAT( / ' Routines pass computational tests if test ratio is ',
    2411. 

  2411.      $      'less than', F8.2, / )
    2412. 

  2412.  9981 FORMAT( ' Relative machine ', A, ' is taken to be', D16.6 )
    2413. 

  2413.  9980 FORMAT( ' *** Error code from ', A, ' = ', I4 )
    2414. 

  2414.  9979 FORMAT( / ' Tests of the Nonsymmetric Eigenvalue Problem Driver',
    2415. 

  2415.      $      / '    DGEEV (eigenvalues and eigevectors)' )
    2416. 

  2416.  9978 FORMAT( / ' Tests of the Nonsymmetric Eigenvalue Problem Driver',
    2417. 

  2417.      $      / '    DGEES (Schur form)' )
    2418. 

  2418.  9977 FORMAT( / ' Tests of the Nonsymmetric Eigenvalue Problem Expert',
    2419. 

  2419.      $      ' Driver', / '    DGEEVX (eigenvalues, eigenvectors and',
    2420. 

  2420.      $      ' condition numbers)' )
    2421. 

  2421.  9976 FORMAT( / ' Tests of the Nonsymmetric Eigenvalue Problem Expert',
    2422. 

  2422.      $      ' Driver', / '    DGEESX (Schur form and condition',
    2423. 

  2423.      $      ' numbers)' )
    2424. 

  2424.  9975 FORMAT( / ' Tests of the Generalized Nonsymmetric Eigenvalue ',
    2425. 

  2425.      $      'Problem routines' )
    2426. 

  2426.  9974 FORMAT( ' Tests of DSBTRD', / ' (reduction of a symmetric band ',
    2427. 

  2427.      $      'matrix to tridiagonal form)' )
    2428. 

  2428.  9973 FORMAT( / 1X, 71( '-' ) )
    2429. 

  2429.  9972 FORMAT( / ' LAPACK VERSION ', I1, '.', I1, '.', I1 )
    2430. 

  2430.  9971 FORMAT( / ' Tests of the Generalized Linear Regression Model ',
    2431. 

  2431.      $      'routines' )
    2432. 

  2432.  9970 FORMAT( / ' Tests of the Generalized QR and RQ routines' )
    2433. 

  2433.  9969 FORMAT( / ' Tests of the Generalized Singular Value',
    2434. 

  2434.      $      ' Decomposition routines' )
    2435. 

  2435.  9968 FORMAT( / ' Tests of the Linear Least Squares routines' )
    2436. 

  2436.  9967 FORMAT( ' Tests of DGBBRD', / ' (reduction of a general band ',
    2437. 

  2437.      $      'matrix to real bidiagonal form)' )
    2438. 

  2438.  9966 FORMAT( / / 1X, A3, ':  NRHS =', I4 )
    2439. 

  2439.  9965 FORMAT( / ' Tests of the Generalized Nonsymmetric Eigenvalue ',
    2440. 

  2440.      $      'Problem Expert Driver DGGESX' )
    2441. 

  2441.  9964 FORMAT( / ' Tests of the Generalized Nonsymmetric Eigenvalue ',
    2442. 

  2442.      $      'Problem Driver DGGES' )
    2443. 

  2443.  9963 FORMAT( / ' Tests of the Generalized Nonsymmetric Eigenvalue ',
    2444. 

  2444.      $      'Problem Driver DGGEV' )
    2445. 

  2445.  9962 FORMAT( / ' Tests of the Generalized Nonsymmetric Eigenvalue ',
    2446. 

  2446.      $      'Problem Expert Driver DGGEVX' )
    2447. 

  2447.  9961 FORMAT( / / 1X, A3, ':  NB =', I4, ', NBMIN =', I4, ', NX =', I4,
    2448. 

  2448.      $      ', INMIN=', I4, 
    2449. 

  2449.      $      ', INWIN =', I4, ', INIBL =', I4, ', ISHFTS =', I4,
    2450. 

  2450.      $      ', IACC22 =', I4)
    2451. 

  2451.  9960 FORMAT( / ' Tests of the CS Decomposition routines' )
    2452. 

  2452. *
    2453. 

  2453. *     End of DCHKEE
    2454. 

  2454. *
    2455. 

  2455.       END
    2456. 

  2456.