PageRenderTime 63ms CodeModel.GetById 21ms RepoModel.GetById 0ms app.codeStats 1ms

/lapack-netlib/TESTING/EIG/dchkee.f

http://github.com/xianyi/OpenBLAS
FORTRAN Legacy | 2507 lines | 1008 code | 0 blank | 1499 comment | 0 complexity | 93ea655f0fe2f4195637cdbd8a1d43a0 MD5 | raw file
Possible License(s): BSD-3-Clause, LGPL-2.0

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

    

  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 DGGHD3, DGGBAL, DGGBAK, DHGEQZ, and DTGEVC
    48. 

  48. *>
    49. 

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

  50. *>     Test DGGES
    51. 

  51. *>
    52. 

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

  53. *>     Test DGGEV
    54. 

  54. *>
    55. 

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

  56. *>     Test DGGESX
    57. 

  57. *>
    58. 

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

  59. *>     Test DGGEVX
    60. 

  60. *>
    61. 

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

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

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

  64. *>
    65. 

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

  66. *>     Test DSBTRD
    67. 

  67. *>
    68. 

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

  69. *>     Test DGBBRD
    70. 

  70. *>
    71. 

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

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

  73. *>     DTRSEN, and DLAQTR
    74. 

  74. *>
    75. 

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

  76. *>     Test DGEBAL
    77. 

  77. *>
    78. 

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

  79. *>     Test DGEBAK
    80. 

  80. *>
    81. 

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

  82. *>     Test DGGBAL
    83. 

  83. *>
    84. 

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

  85. *>     Test DGGBAK
    86. 

  86. *>
    87. 

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

  88. *>     Tests DGGGLM
    89. 

  89. *>
    90. 

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

  91. *>     Tests DGGQRF and DGGRQF
    92. 

  92. *>
    93. 

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

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

  95. *>
    96. 

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

  97. *>     Tests DORCSD
    98. 

  98. *>
    99. 

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

  100. *>     Tests DGGLSE
    101. 

  101. *>
    102. 

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

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

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

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

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

  107. *>
    108. 

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

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

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

  111. *>
    112. 

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

  113. *>
    114. 

  114. *> DHS or NEP      21     DCHKHS
    115. 

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

  116. *>                 18     DDRVST (drivers)
    117. 

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

  118. *>                  5     DDRVBD (drivers)
    119. 

  119. *> DEV             21     DDRVEV
    120. 

  120. *> DES             21     DDRVES
    121. 

  121. *> DVX             21     DDRVVX
    122. 

  122. *> DSX             21     DDRVSX
    123. 

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

  124. *> DGS             26     DDRGES
    125. 

  125. *> DGX              5     DDRGSX
    126. 

  126. *> DGV             26     DDRGEV
    127. 

  127. *> DXV              2     DDRGVX
    128. 

  128. *> DSG             21     DDRVSG
    129. 

  129. *> DSB             15     DCHKSB
    130. 

  130. *> DBB             15     DCHKBB
    131. 

  131. *> DEC              -     DCHKEC
    132. 

  132. *> DBL              -     DCHKBL
    133. 

  133. *> DBK              -     DCHKBK
    134. 

  134. *> DGL              -     DCHKGL
    135. 

  135. *> DGK              -     DCHKGK
    136. 

  136. *> GLM              8     DCKGLM
    137. 

  137. *> GQR              8     DCKGQR
    138. 

  138. *> GSV              8     DCKGSV
    139. 

  139. *> CSD              3     DCKCSD
    140. 

  140. *> LSE              8     DCKLSE
    141. 

  141. *>
    142. 

  142. *>-----------------------------------------------------------------------
    143. 

  143. *>
    144. 

  144. *> NEP input file:
    145. 

  145. *>
    146. 

  146. *> line 2:  NN, INTEGER
    147. 

  147. *>          Number of values of N.
    148. 

  148. *>
    149. 

  149. *> line 3:  NVAL, INTEGER array, dimension (NN)
    150. 

  150. *>          The values for the matrix dimension N.
    151. 

  151. *>
    152. 

  152. *> line 4:  NPARMS, INTEGER
    153. 

  153. *>          Number of values of the parameters NB, NBMIN, NX, NS, and
    154. 

  154. *>          MAXB.
    155. 

  155. *>
    156. 

  156. *> line 5:  NBVAL, INTEGER array, dimension (NPARMS)
    157. 

  157. *>          The values for the blocksize NB.
    158. 

  158. *>
    159. 

  159. *> line 6:  NBMIN, INTEGER array, dimension (NPARMS)
    160. 

  160. *>          The values for the minimum blocksize NBMIN.
    161. 

  161. *>
    162. 

  162. *> line 7:  NXVAL, INTEGER array, dimension (NPARMS)
    163. 

  163. *>          The values for the crossover point NX.
    164. 

  164. *>
    165. 

  165. *> line 8:  INMIN, INTEGER array, dimension (NPARMS)
    166. 

  166. *>          LAHQR vs TTQRE crossover point, >= 11
    167. 

  167. *>
    168. 

  168. *> line 9:  INWIN, INTEGER array, dimension (NPARMS)
    169. 

  169. *>          recommended deflation window size
    170. 

  170. *>
    171. 

  171. *> line 10: INIBL, INTEGER array, dimension (NPARMS)
    172. 

  172. *>          nibble crossover point
    173. 

  173. *>
    174. 

  174. *> line 11: ISHFTS, INTEGER array, dimension (NPARMS)
    175. 

  175. *>          number of simultaneous shifts)
    176. 

  176. *>
    177. 

  177. *> line 12: IACC22, INTEGER array, dimension (NPARMS)
    178. 

  178. *>          select structured matrix multiply: 0, 1 or 2)
    179. 

  179. *>
    180. 

  180. *> line 13: THRESH
    181. 

  181. *>          Threshold value for the test ratios.  Information will be
    182. 

  182. *>          printed about each test for which the test ratio is greater
    183. 

  183. *>          than or equal to the threshold.  To have all of the test
    184. 

  184. *>          ratios printed, use THRESH = 0.0 .
    185. 

  185. *>
    186. 

  186. *> line 14: NEWSD, INTEGER
    187. 

  187. *>          A code indicating how to set the random number seed.
    188. 

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

  189. *>          = 1:  Initialize the seed to a default value only before the
    190. 

  190. *>                first run
    191. 

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

  192. *>
    193. 

  193. *> If line 14 was 2:
    194. 

  194. *>
    195. 

  195. *> line 15: INTEGER array, dimension (4)
    196. 

  196. *>          Four integer values for the random number seed.
    197. 

  197. *>
    198. 

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

  199. *>          the user to specify the matrix types.  Each line contains
    200. 

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

  201. *>          of matrix types must be the first nonblank item in columns
    202. 

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

  203. *>          less than the maximum number of possible types, a second
    204. 

  204. *>          line will be read to get the numbers of the matrix types to
    205. 

  205. *>          be used.  For example,
    206. 

  206. *> NEP 21
    207. 

  207. *>          requests all of the matrix types for the nonsymmetric
    208. 

  208. *>          eigenvalue problem, while
    209. 

  209. *> NEP  4
    210. 

  210. *> 9 10 11 12
    211. 

  211. *>          requests only matrices of type 9, 10, 11, and 12.
    212. 

  212. *>
    213. 

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

  214. *>          nonsymmetric eigenvalue routines.
    215. 

  215. *>
    216. 

  216. *>-----------------------------------------------------------------------
    217. 

  217. *>
    218. 

  218. *> SEP or DSG input file:
    219. 

  219. *>
    220. 

  220. *> line 2:  NN, INTEGER
    221. 

  221. *>          Number of values of N.
    222. 

  222. *>
    223. 

  223. *> line 3:  NVAL, INTEGER array, dimension (NN)
    224. 

  224. *>          The values for the matrix dimension N.
    225. 

  225. *>
    226. 

  226. *> line 4:  NPARMS, INTEGER
    227. 

  227. *>          Number of values of the parameters NB, NBMIN, and NX.
    228. 

  228. *>
    229. 

  229. *> line 5:  NBVAL, INTEGER array, dimension (NPARMS)
    230. 

  230. *>          The values for the blocksize NB.
    231. 

  231. *>
    232. 

  232. *> line 6:  NBMIN, INTEGER array, dimension (NPARMS)
    233. 

  233. *>          The values for the minimum blocksize NBMIN.
    234. 

  234. *>
    235. 

  235. *> line 7:  NXVAL, INTEGER array, dimension (NPARMS)
    236. 

  236. *>          The values for the crossover point NX.
    237. 

  237. *>
    238. 

  238. *> line 8:  THRESH
    239. 

  239. *>          Threshold value for the test ratios.  Information will be
    240. 

  240. *>          printed about each test for which the test ratio is greater
    241. 

  241. *>          than or equal to the threshold.
    242. 

  242. *>
    243. 

  243. *> line 9:  TSTCHK, LOGICAL
    244. 

  244. *>          Flag indicating whether or not to test the LAPACK routines.
    245. 

  245. *>
    246. 

  246. *> line 10: TSTDRV, LOGICAL
    247. 

  247. *>          Flag indicating whether or not to test the driver routines.
    248. 

  248. *>
    249. 

  249. *> line 11: TSTERR, LOGICAL
    250. 

  250. *>          Flag indicating whether or not to test the error exits for
    251. 

  251. *>          the LAPACK routines and driver routines.
    252. 

  252. *>
    253. 

  253. *> line 12: NEWSD, INTEGER
    254. 

  254. *>          A code indicating how to set the random number seed.
    255. 

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

  256. *>          = 1:  Initialize the seed to a default value only before the
    257. 

  257. *>                first run
    258. 

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

  259. *>
    260. 

  260. *> If line 12 was 2:
    261. 

  261. *>
    262. 

  262. *> line 13: INTEGER array, dimension (4)
    263. 

  263. *>          Four integer values for the random number seed.
    264. 

  264. *>
    265. 

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

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

  267. *>          symmetric eigenvalue routines and driver routines, and
    268. 

  268. *>          'DSG' for the routines for the symmetric generalized
    269. 

  269. *>          eigenvalue problem.
    270. 

  270. *>
    271. 

  271. *>-----------------------------------------------------------------------
    272. 

  272. *>
    273. 

  273. *> SVD input file:
    274. 

  274. *>
    275. 

  275. *> line 2:  NN, INTEGER
    276. 

  276. *>          Number of values of M and N.
    277. 

  277. *>
    278. 

  278. *> line 3:  MVAL, INTEGER array, dimension (NN)
    279. 

  279. *>          The values for the matrix row dimension M.
    280. 

  280. *>
    281. 

  281. *> line 4:  NVAL, INTEGER array, dimension (NN)
    282. 

  282. *>          The values for the matrix column dimension N.
    283. 

  283. *>
    284. 

  284. *> line 5:  NPARMS, INTEGER
    285. 

  285. *>          Number of values of the parameter NB, NBMIN, NX, and NRHS.
    286. 

  286. *>
    287. 

  287. *> line 6:  NBVAL, INTEGER array, dimension (NPARMS)
    288. 

  288. *>          The values for the blocksize NB.
    289. 

  289. *>
    290. 

  290. *> line 7:  NBMIN, INTEGER array, dimension (NPARMS)
    291. 

  291. *>          The values for the minimum blocksize NBMIN.
    292. 

  292. *>
    293. 

  293. *> line 8:  NXVAL, INTEGER array, dimension (NPARMS)
    294. 

  294. *>          The values for the crossover point NX.
    295. 

  295. *>
    296. 

  296. *> line 9:  NSVAL, INTEGER array, dimension (NPARMS)
    297. 

  297. *>          The values for the number of right hand sides NRHS.
    298. 

  298. *>
    299. 

  299. *> line 10: THRESH
    300. 

  300. *>          Threshold value for the test ratios.  Information will be
    301. 

  301. *>          printed about each test for which the test ratio is greater
    302. 

  302. *>          than or equal to the threshold.
    303. 

  303. *>
    304. 

  304. *> line 11: TSTCHK, LOGICAL
    305. 

  305. *>          Flag indicating whether or not to test the LAPACK routines.
    306. 

  306. *>
    307. 

  307. *> line 12: TSTDRV, LOGICAL
    308. 

  308. *>          Flag indicating whether or not to test the driver routines.
    309. 

  309. *>
    310. 

  310. *> line 13: TSTERR, LOGICAL
    311. 

  311. *>          Flag indicating whether or not to test the error exits for
    312. 

  312. *>          the LAPACK routines and driver routines.
    313. 

  313. *>
    314. 

  314. *> line 14: NEWSD, INTEGER
    315. 

  315. *>          A code indicating how to set the random number seed.
    316. 

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

  317. *>          = 1:  Initialize the seed to a default value only before the
    318. 

  318. *>                first run
    319. 

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

  320. *>
    321. 

  321. *> If line 14 was 2:
    322. 

  322. *>
    323. 

  323. *> line 15: INTEGER array, dimension (4)
    324. 

  324. *>          Four integer values for the random number seed.
    325. 

  325. *>
    326. 

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

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

  328. *>          SVD routines and the SVD driver routines.
    329. 

  329. *>
    330. 

  330. *>-----------------------------------------------------------------------
    331. 

  331. *>
    332. 

  332. *> DEV and DES data files:
    333. 

  333. *>
    334. 

  334. *> line 1:  'DEV' or 'DES' in columns 1 to 3.
    335. 

  335. *>
    336. 

  336. *> line 2:  NSIZES, INTEGER
    337. 

  337. *>          Number of sizes of matrices to use. Should be at least 0
    338. 

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

  339. *>          (although the remaining  3 lines are still read).
    340. 

  340. *>
    341. 

  341. *> line 3:  NN, INTEGER array, dimension(NSIZES)
    342. 

  342. *>          Dimensions of matrices to be tested.
    343. 

  343. *>
    344. 

  344. *> line 4:  NB, NBMIN, NX, NS, NBCOL, INTEGERs
    345. 

  345. *>          These integer parameters determine how blocking is done
    346. 

  346. *>          (see ILAENV for details)
    347. 

  347. *>          NB     : block size
    348. 

  348. *>          NBMIN  : minimum block size
    349. 

  349. *>          NX     : minimum dimension for blocking
    350. 

  350. *>          NS     : number of shifts in xHSEQR
    351. 

  351. *>          NBCOL  : minimum column dimension for blocking
    352. 

  352. *>
    353. 

  353. *> line 5:  THRESH, REAL
    354. 

  354. *>          The test threshold against which computed residuals are
    355. 

  355. *>          compared. Should generally be in the range from 10. to 20.
    356. 

  356. *>          If it is 0., all test case data will be printed.
    357. 

  357. *>
    358. 

  358. *> line 6:  TSTERR, LOGICAL
    359. 

  359. *>          Flag indicating whether or not to test the error exits.
    360. 

  360. *>
    361. 

  361. *> line 7:  NEWSD, INTEGER
    362. 

  362. *>          A code indicating how to set the random number seed.
    363. 

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

  364. *>          = 1:  Initialize the seed to a default value only before the
    365. 

  365. *>                first run
    366. 

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

  367. *>
    368. 

  368. *> If line 7 was 2:
    369. 

  369. *>
    370. 

  370. *> line 8:  INTEGER array, dimension (4)
    371. 

  371. *>          Four integer values for the random number seed.
    372. 

  372. *>
    373. 

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

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

  375. *>          'DES' to test SGEES.
    376. 

  376. *>
    377. 

  377. *>-----------------------------------------------------------------------
    378. 

  378. *>
    379. 

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

  380. *> and the second part consists of test matrices with precomputed
    381. 

  381. *> solutions.
    382. 

  382. *>
    383. 

  383. *> line 1:  'DVX' in columns 1-3.
    384. 

  384. *>
    385. 

  385. *> line 2:  NSIZES, INTEGER
    386. 

  386. *>          If NSIZES = 0, no testing of randomly generated examples
    387. 

  387. *>          is done, but any precomputed examples are tested.
    388. 

  388. *>
    389. 

  389. *> line 3:  NN, INTEGER array, dimension(NSIZES)
    390. 

  390. *>
    391. 

  391. *> line 4:  NB, NBMIN, NX, NS, NBCOL, INTEGERs
    392. 

  392. *>
    393. 

  393. *> line 5:  THRESH, REAL
    394. 

  394. *>
    395. 

  395. *> line 6:  TSTERR, LOGICAL
    396. 

  396. *>
    397. 

  397. *> line 7:  NEWSD, INTEGER
    398. 

  398. *>
    399. 

  399. *> If line 7 was 2:
    400. 

  400. *>
    401. 

  401. *> line 8:  INTEGER array, dimension (4)
    402. 

  402. *>
    403. 

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

  404. *>          followed by the number of matrix types, possibly with
    405. 

  405. *>          a second line to specify certain matrix types.
    406. 

  406. *>          If the number of matrix types = 0, no testing of randomly
    407. 

  407. *>          generated examples is done, but any precomputed examples
    408. 

  408. *>          are tested.
    409. 

  409. *>
    410. 

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

  411. *>          its dimension. The first line contains the dimension (a
    412. 

  412. *>          single integer). The next N lines contain the matrix, one
    413. 

  413. *>          row per line. The last N lines correspond to each
    414. 

  414. *>          eigenvalue. Each of these last N lines contains 4 real
    415. 

  415. *>          values: the real part of the eigenvalue, the imaginary
    416. 

  416. *>          part of the eigenvalue, the reciprocal condition number of
    417. 

  417. *>          the eigenvalues, and the reciprocal condition number of the
    418. 

  418. *>          eigenvector.  The end of data is indicated by dimension N=0.
    419. 

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

  420. *>          line containing N=0.
    421. 

  421. *>
    422. 

  422. *>-----------------------------------------------------------------------
    423. 

  423. *>
    424. 

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

  425. *> second part consists of test matrices with precomputed solutions.
    426. 

  426. *>
    427. 

  427. *> line 1:  'DSX' in columns 1-3.
    428. 

  428. *>
    429. 

  429. *> line 2:  NSIZES, INTEGER
    430. 

  430. *>          If NSIZES = 0, no testing of randomly generated examples
    431. 

  431. *>          is done, but any precomputed examples are tested.
    432. 

  432. *>
    433. 

  433. *> line 3:  NN, INTEGER array, dimension(NSIZES)
    434. 

  434. *>
    435. 

  435. *> line 4:  NB, NBMIN, NX, NS, NBCOL, INTEGERs
    436. 

  436. *>
    437. 

  437. *> line 5:  THRESH, REAL
    438. 

  438. *>
    439. 

  439. *> line 6:  TSTERR, LOGICAL
    440. 

  440. *>
    441. 

  441. *> line 7:  NEWSD, INTEGER
    442. 

  442. *>
    443. 

  443. *> If line 7 was 2:
    444. 

  444. *>
    445. 

  445. *> line 8:  INTEGER array, dimension (4)
    446. 

  446. *>
    447. 

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

  448. *>          followed by the number of matrix types, possibly with
    449. 

  449. *>          a second line to specify certain matrix types.
    450. 

  450. *>          If the number of matrix types = 0, no testing of randomly
    451. 

  451. *>          generated examples is done, but any precomputed examples
    452. 

  452. *>          are tested.
    453. 

  453. *>
    454. 

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

  455. *>          dimension. The first line contains the dimension N and the
    456. 

  456. *>          dimension M of an invariant subspace. The second line
    457. 

  457. *>          contains M integers, identifying the eigenvalues in the
    458. 

  458. *>          invariant subspace (by their position in a list of
    459. 

  459. *>          eigenvalues ordered by increasing real part). The next N
    460. 

  460. *>          lines contain the matrix. The last line contains the
    461. 

  461. *>          reciprocal condition number for the average of the selected
    462. 

  462. *>          eigenvalues, and the reciprocal condition number for the
    463. 

  463. *>          corresponding right invariant subspace. The end of data is
    464. 

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

  465. *>          is to be tested, there must be at least one line containing
    466. 

  466. *>          N=0 and M=0.
    467. 

  467. *>
    468. 

  468. *>-----------------------------------------------------------------------
    469. 

  469. *>
    470. 

  470. *> DGG input file:
    471. 

  471. *>
    472. 

  472. *> line 2:  NN, INTEGER
    473. 

  473. *>          Number of values of N.
    474. 

  474. *>
    475. 

  475. *> line 3:  NVAL, INTEGER array, dimension (NN)
    476. 

  476. *>          The values for the matrix dimension N.
    477. 

  477. *>
    478. 

  478. *> line 4:  NPARMS, INTEGER
    479. 

  479. *>          Number of values of the parameters NB, NBMIN, NS, MAXB, and
    480. 

  480. *>          NBCOL.
    481. 

  481. *>
    482. 

  482. *> line 5:  NBVAL, INTEGER array, dimension (NPARMS)
    483. 

  483. *>          The values for the blocksize NB.
    484. 

  484. *>
    485. 

  485. *> line 6:  NBMIN, INTEGER array, dimension (NPARMS)
    486. 

  486. *>          The values for NBMIN, the minimum row dimension for blocks.
    487. 

  487. *>
    488. 

  488. *> line 7:  NSVAL, INTEGER array, dimension (NPARMS)
    489. 

  489. *>          The values for the number of shifts.
    490. 

  490. *>
    491. 

  491. *> line 8:  MXBVAL, INTEGER array, dimension (NPARMS)
    492. 

  492. *>          The values for MAXB, used in determining minimum blocksize.
    493. 

  493. *>
    494. 

  494. *> line 9:  IACC22, INTEGER array, dimension (NPARMS)
    495. 

  495. *>          select structured matrix multiply: 1 or 2)
    496. 

  496. *>
    497. 

  497. *> line 10: NBCOL, INTEGER array, dimension (NPARMS)
    498. 

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

  499. *>          blocks.
    500. 

  500. *>
    501. 

  501. *> line 11: THRESH
    502. 

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

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

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

  505. *>
    506. 

  506. *> line 12: TSTCHK, LOGICAL
    507. 

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

  508. *>
    509. 

  509. *> line 13: TSTDRV, LOGICAL
    510. 

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

  511. *>
    512. 

  512. *> line 14: TSTERR, LOGICAL
    513. 

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

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

  515. *>
    516. 

  516. *> line 15: NEWSD, INTEGER
    517. 

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

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

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

  520. *>                first run
    521. 

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

  522. *>
    523. 

  523. *> If line 15 was 2:
    524. 

  524. *>
    525. 

  525. *> line 16: INTEGER array, dimension (4)
    526. 

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

  527. *>
    528. 

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

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

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

  531. *>
    532. 

  532. *>-----------------------------------------------------------------------
    533. 

  533. *>
    534. 

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

  535. *>
    536. 

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

  537. *>
    538. 

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

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

  540. *>
    541. 

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

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

  543. *>
    544. 

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

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

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

  547. *>          NB     : block size
    548. 

  548. *>          NBMIN  : minimum block size
    549. 

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

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

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

  552. *>
    553. 

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

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

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

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

  557. *>
    558. 

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

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

  560. *>
    561. 

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

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

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

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

  565. *>                first run
    566. 

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

  567. *>
    568. 

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

  569. *>
    570. 

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

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

  572. *>
    573. 

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

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

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

  576. *>
    577. 

  577. *>-----------------------------------------------------------------------
    578. 

  578. *>
    579. 

  579. *> DXV input files:
    580. 

  580. *>
    581. 

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

  582. *>
    583. 

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

  584. *>          Value of N.
    585. 

  585. *>
    586. 

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

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

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

  589. *>          NB     : block size
    590. 

  590. *>          NBMIN  : minimum block size
    591. 

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

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

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

  594. *>
    595. 

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

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

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

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

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

  600. *>
    601. 

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

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

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

  604. *>
    605. 

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

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

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

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

  609. *>                first run
    610. 

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

  611. *>
    612. 

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

  613. *>
    614. 

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

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

  616. *>
    617. 

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

  618. *>
    619. 

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

  620. *>
    621. 

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

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

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

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

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

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

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

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

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

  630. *>
    631. 

  631. *>-----------------------------------------------------------------------
    632. 

  632. *>
    633. 

  633. *> DGX input files:
    634. 

  634. *>
    635. 

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

  636. *>
    637. 

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

  638. *>          Value of N.
    639. 

  639. *>
    640. 

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

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

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

  643. *>          NB     : block size
    644. 

  644. *>          NBMIN  : minimum block size
    645. 

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

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

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

  648. *>
    649. 

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

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

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

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

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

  654. *>
    655. 

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

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

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

  658. *>
    659. 

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

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

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

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

  663. *>                first run
    664. 

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

  665. *>
    666. 

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

  667. *>
    668. 

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

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

  670. *>
    671. 

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

  672. *>
    673. 

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

  674. *>
    675. 

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

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

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

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

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

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

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

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

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

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

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

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

  687. *>
    688. 

  688. *>-----------------------------------------------------------------------
    689. 

  689. *>
    690. 

  690. *> DSB input file:
    691. 

  691. *>
    692. 

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

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

  694. *>
    695. 

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

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

  697. *>
    698. 

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

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

  700. *>
    701. 

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

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

  703. *>
    704. 

  704. *> line 6:  THRESH
    705. 

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

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

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

  708. *>
    709. 

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

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

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

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

  713. *>                first run
    714. 

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

  715. *>
    716. 

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

  717. *>
    718. 

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

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

  720. *>
    721. 

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

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

  723. *>
    724. 

  724. *>-----------------------------------------------------------------------
    725. 

  725. *>
    726. 

  726. *> DBB input file:
    727. 

  727. *>
    728. 

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

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

  730. *>
    731. 

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

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

  733. *>
    734. 

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

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

  736. *>
    737. 

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

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

  739. *>
    740. 

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

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

  742. *>
    743. 

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

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

  745. *>
    746. 

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

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

  748. *>
    749. 

  749. *> line 8:  THRESH
    750. 

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

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

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

  753. *>
    754. 

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

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

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

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

  758. *>                first run
    759. 

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

  760. *>
    761. 

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

  762. *>
    763. 

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

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

  765. *>
    766. 

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

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

  768. *>
    769. 

  769. *>-----------------------------------------------------------------------
    770. 

  770. *>
    771. 

  771. *> DEC input file:
    772. 

  772. *>
    773. 

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

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

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

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

  777. *>
    778. 

  778. *> lines  3-EOF:
    779. 

  779. *>
    780. 

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

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

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

  783. *>
    784. 

  784. *>-----------------------------------------------------------------------
    785. 

  785. *>
    786. 

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

  787. *>
    788. 

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

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

  790. *>
    791. 

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

  792. *>
    793. 

  793. *>-----------------------------------------------------------------------
    794. 

  794. *>
    795. 

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

  796. *>
    797. 

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

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

  799. *>
    800. 

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

  801. *>
    802. 

  802. *>-----------------------------------------------------------------------
    803. 

  803. *>
    804. 

  804. *> GLM data file:
    805. 

  805. *>
    806. 

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

  807. *>
    808. 

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

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

  810. *>
    811. 

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

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

  813. *>
    814. 

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

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

  816. *>
    817. 

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

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

  819. *>
    820. 

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

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

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

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

  824. *>
    825. 

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

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

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

  828. *>
    829. 

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

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

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

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

  833. *>                first run
    834. 

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

  835. *>
    836. 

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

  837. *>
    838. 

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

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

  840. *>
    841. 

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

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

  843. *>          linear regression model routines.
    844. 

  844. *>
    845. 

  845. *>-----------------------------------------------------------------------
    846. 

  846. *>
    847. 

  847. *> GQR data file:
    848. 

  848. *>
    849. 

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

  850. *>
    851. 

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

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

  853. *>
    854. 

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

  855. *>          Values of M.
    856. 

  856. *>
    857. 

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

  858. *>          Values of P.
    859. 

  859. *>
    860. 

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

  861. *>          Values of N.
    862. 

  862. *>
    863. 

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

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

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

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

  867. *>
    868. 

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

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

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

  871. *>
    872. 

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

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

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

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

  876. *>                first run
    877. 

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

  878. *>
    879. 

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

  880. *>
    881. 

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

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

  883. *>
    884. 

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

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

  886. *>          QR and RQ routines.
    887. 

  887. *>
    888. 

  888. *>-----------------------------------------------------------------------
    889. 

  889. *>
    890. 

  890. *> GSV data file:
    891. 

  891. *>
    892. 

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

  893. *>
    894. 

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

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

  896. *>
    897. 

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

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

  899. *>
    900. 

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

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

  902. *>
    903. 

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

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

  905. *>
    906. 

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

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

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

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

  910. *>
    911. 

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

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

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

  914. *>
    915. 

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

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

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

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

  919. *>                first run
    920. 

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

  921. *>
    922. 

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

  923. *>
    924. 

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

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

  926. *>
    927. 

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

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

  929. *>          SVD routines.
    930. 

  930. *>
    931. 

  931. *>-----------------------------------------------------------------------
    932. 

  932. *>
    933. 

  933. *> CSD data file:
    934. 

  934. *>
    935. 

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

  936. *>
    937. 

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

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

  939. *>
    940. 

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

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

  942. *>
    943. 

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

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

  945. *>
    946. 

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

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

  948. *>
    949. 

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

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

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

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

  953. *>
    954. 

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

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

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

  957. *>
    958. 

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

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

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

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

  962. *>                first run
    963. 

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

  964. *>
    965. 

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

  966. *>
    967. 

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

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

  969. *>
    970. 

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

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

  972. *>
    973. 

  973. *>-----------------------------------------------------------------------
    974. 

  974. *>
    975. 

  975. *> LSE data file:
    976. 

  976. *>
    977. 

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

  978. *>
    979. 

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

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

  981. *>
    982. 

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

  983. *>          Values of M.
    984. 

  984. *>
    985. 

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

  986. *>          Values of P.
    987. 

  987. *>
    988. 

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

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

  990. *>
    991. 

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

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

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

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

  995. *>
    996. 

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

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

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

  999. *>
    1000. 

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

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

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

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

  1004. *>                first run
    1005. 

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

  1006. *>
    1007. 

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

  1008. *>
    1009. 

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

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

  1011. *>
    1012. 

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

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

  1014. *>          SVD routines.
    1015. 

  1015. *>
    1016. 

  1016. *>-----------------------------------------------------------------------
    1017. 

  1017. *>
    1018. 

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

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

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

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

  1022. *> \endverbatim
    1023. 

  1023. *
    1024. 

  1024. *  Arguments:
    1025. 

  1025. *  ==========
    1026. 

  1026. *
    1027. 

  1027. *
    1028. 

  1028. *  Authors:
    1029. 

  1029. *  ========
    1030. 

  1030. *
    1031. 

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

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

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

  1034. *> \author NAG Ltd.
    1035. 

  1035. *
    1036. 

  1036. *> \date June 2016
    1037. 

  1037. *
    1038. 

  1038. *> \ingroup double_eig
    1039. 

  1039. *
    1040. 

  1040. *  =====================================================================
    1041. 

  1041.       PROGRAM DCHKEE
    1042. 

  1042. *
    1043. 

  1043. *  -- LAPACK test routine (version 3.7.0) --
    1044. 

  1044. *  -- LAPACK is a software package provided by Univ. of Tennessee,    --
    1045. 

  1045. *  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
    1046. 

  1046. *     June 2016
    1047. 

  1047. *
    1048. 

  1048. *  =====================================================================
    1049. 

  1049. *
    1050. 

  1050. *     .. Parameters ..
    1051. 

  1051.       INTEGER            NMAX
    1052. 

  1052.       PARAMETER          ( NMAX = 132 )
    1053. 

  1053.       INTEGER            NCMAX
    1054. 

  1054.       PARAMETER          ( NCMAX = 20 )
    1055. 

  1055.       INTEGER            NEED
    1056. 

  1056.       PARAMETER          ( NEED = 14 )
    1057. 

  1057.       INTEGER            LWORK
    1058. 

  1058.       PARAMETER          ( LWORK = NMAX*( 5*NMAX+5 )+1 )
    1059. 

  1059.       INTEGER            LIWORK
    1060. 

  1060.       PARAMETER          ( LIWORK = NMAX*( 5*NMAX+20 ) )
    1061. 

  1061.       INTEGER            MAXIN
    1062. 

  1062.       PARAMETER          ( MAXIN = 20 )
    1063. 

  1063.       INTEGER            MAXT
    1064. 

  1064.       PARAMETER          ( MAXT = 30 )
    1065. 

  1065.       INTEGER            NIN, NOUT
    1066. 

  1066.       PARAMETER          ( NIN = 5, NOUT = 6 )
    1067. 

  1067. *     ..
    1068. 

  1068. *     .. Local Scalars ..
    1069. 

  1069.       LOGICAL            CSD, DBB, DGG, DSB, FATAL, GLM, GQR, GSV, LSE,
    1070. 

  1070.      $                   NEP, DBK, DBL, SEP, DES, DEV, DGK, DGL, DGS,
    1071. 

  1071.      $                   DGV, DGX, DSX, SVD, DVX, DXV, TSTCHK, TSTDIF,
    1072. 

  1072.      $                   TSTDRV, TSTERR
    1073. 

  1073.       CHARACTER          C1
    1074. 

  1074.       CHARACTER*3        C3, PATH
    1075. 

  1075.       CHARACTER*32       VNAME
    1076. 

  1076.       CHARACTER*10       INTSTR
    1077. 

  1077.       CHARACTER*80       LINE
    1078. 

  1078.       INTEGER            I, I1, IC, INFO, ITMP, K, LENP, MAXTYP, NEWSD,
    1079. 

  1079.      $                   NK, NN, NPARMS, NRHS, NTYPES,
    1080. 

  1080.      $                   VERS_MAJOR, VERS_MINOR, VERS_PATCH
    1081. 

  1081.       DOUBLE PRECISION   EPS, S1, S2, THRESH, THRSHN
    1082. 

  1082. *     ..
    1083. 

  1083. *     .. Local Arrays ..
    1084. 

  1084.       LOGICAL            DOTYPE( MAXT ), LOGWRK( NMAX )
    1085. 

  1085.       INTEGER            IOLDSD( 4 ), ISEED( 4 ), IWORK( LIWORK ),
    1086. 

  1086.      $                   KVAL( MAXIN ), MVAL( MAXIN ), MXBVAL( MAXIN ),
    1087. 

  1087.      $                   NBCOL( MAXIN ), NBMIN( MAXIN ), NBVAL( MAXIN ),
    1088. 

  1088.      $                   NSVAL( MAXIN ), NVAL( MAXIN ), NXVAL( MAXIN ),
    1089. 

  1089.      $                   PVAL( MAXIN )
    1090. 

  1090.       INTEGER            INMIN( MAXIN ), INWIN( MAXIN ), INIBL( MAXIN ),
    1091. 

  1091.      $                   ISHFTS( MAXIN ), IACC22( MAXIN )
    1092. 

  1092.       DOUBLE PRECISION   A( NMAX*NMAX, NEED ), B( NMAX*NMAX, 5 ),
    1093. 

  1093.      $                   C( NCMAX*NCMAX, NCMAX*NCMAX ), D( NMAX, 12 ),
    1094. 

  1094.      $                   RESULT( 500 ), TAUA( NMAX ), TAUB( NMAX ),
    1095. 

  1095.      $                   WORK( LWORK ), X( 5*NMAX )
    1096. 

  1096. *     ..
    1097. 

  1097. *     .. External Functions ..
    1098. 

  1098.       LOGICAL            LSAMEN
    1099. 

  1099.       DOUBLE PRECISION   DLAMCH, DSECND
    1100. 

  1100.       EXTERNAL           LSAMEN, DLAMCH, DSECND
    1101. 

  1101. *     ..
    1102. 

  1102. *     .. External Subroutines ..
    1103. 

  1103.       EXTERNAL           ALAREQ, DCHKBB, DCHKBD, DCHKBK, DCHKBL, DCHKEC,
    1104. 

  1104.      $                   DCHKGG, DCHKGK, DCHKGL, DCHKHS, DCHKSB, DCHKST,
    1105. 

  1105.      $                   DCKCSD, DCKGLM, DCKGQR, DCKGSV, DCKLSE, DDRGES,
    1106. 

  1106.      $                   DDRGEV, DDRGSX, DDRGVX, DDRVBD, DDRVES, DDRVEV,
    1107. 

  1107.      $                   DDRVSG, DDRVST, DDRVSX, DDRVVX, DERRBD,
    1108. 

  1108.      $                   DERRED, DERRGG, DERRHS, DERRST, ILAVER, XLAENV,
    1109. 

  1109.      $                   DDRGES3, DDRGEV3, 
    1110. 

  1110.      $                   DCHKST2STG, DDRVST2STG, DCHKSB2STG, DDRVSG2STG
    1111. 

  1111. *     ..
    1112. 

  1112. *     .. Intrinsic Functions ..
    1113. 

  1113.       INTRINSIC          LEN, MIN
    1114. 

  1114. *     ..
    1115. 

  1115. *     .. Scalars in Common ..
    1116. 

  1116.       LOGICAL            LERR, OK
    1117. 

  1117.       CHARACTER*32       SRNAMT
    1118. 

  1118.       INTEGER            INFOT, MAXB, NPROC, NSHIFT, NUNIT, SELDIM,
    1119. 

  1119.      $                   SELOPT
    1120. 

  1120. *     ..
    1121. 

  1121. *     .. Arrays in Common ..
    1122. 

  1122.       LOGICAL            SELVAL( 20 )
    1123. 

  1123.       INTEGER            IPARMS( 100 )
    1124. 

  1124.       DOUBLE PRECISION   SELWI( 20 ), SELWR( 20 )
    1125. 

  1125. *     ..
    1126. 

  1126. *     .. Common blocks ..
    1127. 

  1127.       COMMON             / CENVIR / NPROC, NSHIFT, MAXB
    1128. 

  1128.       COMMON             / INFOC / INFOT, NUNIT, OK, LERR
    1129. 

  1129.       COMMON             / SRNAMC / SRNAMT
    1130. 

  1130.       COMMON             / SSLCT / SELOPT, SELDIM, SELVAL, SELWR, SELWI
    1131. 

  1131.       COMMON             / CLAENV / IPARMS
    1132. 

  1132. *     ..
    1133. 

  1133. *     .. Data statements ..
    1134. 

  1134.       DATA               INTSTR / '0123456789' /
    1135. 

  1135.       DATA               IOLDSD / 0, 0, 0, 1 /
    1136. 

  1136. *     ..
    1137. 

  1137. *     .. Executable Statements ..
    1138. 

  1138. *
    1139. 

  1139.       A = 0.0
    1140. 

  1140.       B = 0.0
    1141. 

  1141.       C = 0.0
    1142. 

  1142.       D = 0.0
    1143. 

  1143.       S1 = DSECND( )
    1144. 

  1144.       FATAL = .FALSE.
    1145. 

  1145.       NUNIT = NOUT
    1146. 

  1146. *
    1147. 

  1147. *     Return to here to read multiple sets of data
    1148. 

  1148. *
    1149. 

  1149.    10 CONTINUE
    1150. 

  1150. *
    1151. 

  1151. *     Read the first line and set the 3-character test path
    1152. 

  1152. *
    1153. 

  1153.       READ( NIN, FMT = '(A80)', END = 380 )LINE
    1154. 

  1154.       PATH = LINE( 1: 3 )
    1155. 

  1155.       NEP = LSAMEN( 3, PATH, 'NEP' ) .OR. LSAMEN( 3, PATH, 'DHS' )
    1156. 

  1156.       SEP = LSAMEN( 3, PATH, 'SEP' ) .OR. LSAMEN( 3, PATH, 'DST' ) .OR.
    1157. 

  1157.      $      LSAMEN( 3, PATH, 'DSG' ) .OR. LSAMEN( 3, PATH, 'SE2' )
    1158. 

  1158.       SVD = LSAMEN( 3, PATH, 'SVD' ) .OR. LSAMEN( 3, PATH, 'DBD' )
    1159. 

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

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

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

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

  1163.       DGG = LSAMEN( 3, PATH, 'DGG' )
    1164. 

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

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

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

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

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

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

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

  1171.       GQR = LSAMEN( 3, PATH, 'GQR' ) .OR. LSAMEN( 3, PATH, 'GRQ' )
    1172. 

  1172.       GSV = LSAMEN( 3, PATH, 'GSV' )
    1173. 

  1173.       CSD = LSAMEN( 3, PATH, 'CSD' )
    1174. 

  1174.       LSE = LSAMEN( 3, PATH, 'LSE' )
    1175. 

  1175.       DBL = LSAMEN( 3, PATH, 'DBL' )
    1176. 

  1176.       DBK = LSAMEN( 3, PATH, 'DBK' )
    1177. 

  1177.       DGL = LSAMEN( 3, PATH, 'DGL' )
    1178. 

  1178.       DGK = LSAMEN( 3, PATH, 'DGK' )
    1179. 

  1179. *
    1180. 

  1180. *     Report values of parameters.
    1181. 

  1181. *
    1182. 

  1182.       IF( PATH.EQ.'   ' ) THEN
    1183. 

  1183.          GO TO 10
    1184. 

  1184.       ELSE IF( NEP ) THEN
    1185. 

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

  1186.       ELSE IF( SEP ) THEN
    1187. 

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

  1188.       ELSE IF( SVD ) THEN
    1189. 

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

  1190.       ELSE IF( DEV ) THEN
    1191. 

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

  1192.       ELSE IF( DES ) THEN
    1193. 

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

  1194.       ELSE IF( DVX ) THEN
    1195. 

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

  1196.       ELSE IF( DSX ) THEN
    1197. 

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

  1198.       ELSE IF( DGG ) THEN
    1199. 

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

  1200.       ELSE IF( DGS ) THEN
    1201. 

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

  1202.       ELSE IF( DGX ) THEN
    1203. 

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

  1204.       ELSE IF( DGV ) THEN
    1205. 

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

  1206.       ELSE IF( DXV ) THEN
    1207. 

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

  1208.       ELSE IF( DSB ) THEN
    1209. 

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

  1210.       ELSE IF( DBB ) THEN
    1211. 

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

  1212.       ELSE IF( GLM ) THEN
    1213. 

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

  1214.       ELSE IF( GQR ) THEN
    1215. 

  1215.          WRITE( NOUT, FMT = 9970 )
    1216. 

  1216.       ELSE IF( GSV ) THEN
    1217. 

  1217.          WRITE( NOUT, FMT = 9969 )
    1218. 

  1218.       ELSE IF( CSD ) THEN
    1219. 

  1219.          WRITE( NOUT, FMT = 9960 )
    1220. 

  1220.       ELSE IF( LSE ) THEN
    1221. 

  1221.          WRITE( NOUT, FMT = 9968 )
    1222. 

  1222.       ELSE IF( DBL ) THEN
    1223. 

  1223. *
    1224. 

  1224. *        DGEBAL:  Balancing
    1225. 

  1225. *
    1226. 

  1226.          CALL DCHKBL( NIN, NOUT )
    1227. 

  1227.          GO TO 10
    1228. 

  1228.       ELSE IF( DBK ) THEN
    1229. 

  1229. *
    1230. 

  1230. *        DGEBAK:  Back transformation
    1231. 

  1231. *
    1232. 

  1232.          CALL DCHKBK( NIN, NOUT )
    1233. 

  1233.          GO TO 10
    1234. 

  1234.       ELSE IF( DGL ) THEN
    1235. 

  1235. *
    1236. 

  1236. *        DGGBAL:  Balancing
    1237. 

  1237. *
    1238. 

  1238.          CALL DCHKGL( NIN, NOUT )
    1239. 

  1239.          GO TO 10
    1240. 

  1240.       ELSE IF( DGK ) THEN
    1241. 

  1241. *
    1242. 

  1242. *        DGGBAK:  Back transformation
    1243. 

  1243. *
    1244. 

  1244.          CALL DCHKGK( NIN, NOUT )
    1245. 

  1245.          GO TO 10
    1246. 

  1246.       ELSE IF( LSAMEN( 3, PATH, 'DEC' ) ) THEN
    1247. 

  1247. *
    1248. 

  1248. *        DEC:  Eigencondition estimation
    1249. 

  1249. *
    1250. 

  1250.          READ( NIN, FMT = * )THRESH
    1251. 

  1251.          CALL XLAENV( 1, 1 )
    1252. 

  1252.          CALL XLAENV( 12, 11 )
    1253. 

  1253.          CALL XLAENV( 13, 2 )
    1254. 

  1254.          CALL XLAENV( 14, 0 )
    1255. 

  1255.          CALL XLAENV( 15, 2 )
    1256. 

  1256.          CALL XLAENV( 16, 2 )
    1257. 

  1257.          TSTERR = .TRUE.
    1258. 

  1258.          CALL DCHKEC( THRESH, TSTERR, NIN, NOUT )
    1259. 

  1259.          GO TO 10
    1260. 

  1260.       ELSE
    1261. 

  1261.          WRITE( NOUT, FMT = 9992 )PATH
    1262. 

  1262.          GO TO 10
    1263. 

  1263.       END IF
    1264. 

  1264.       CALL ILAVER( VERS_MAJOR, VERS_MINOR, VERS_PATCH )
    1265. 

  1265.       WRITE( NOUT, FMT = 9972 ) VERS_MAJOR, VERS_MINOR, VERS_PATCH
    1266. 

  1266.       WRITE( NOUT, FMT = 9984 )
    1267. 

  1267. *
    1268. 

  1268. *     Read the number of values of M, P, and N.
    1269. 

  1269. *
    1270. 

  1270.       READ( NIN, FMT = * )NN
    1271. 

  1271.       IF( NN.LT.0 ) THEN
    1272. 

  1272.          WRITE( NOUT, FMT = 9989 )'   NN ', NN, 1
    1273. 

  1273.          NN = 0
    1274. 

  1274.          FATAL = .TRUE.
    1275. 

  1275.       ELSE IF( NN.GT.MAXIN ) THEN
    1276. 

  1276.          WRITE( NOUT, FMT = 9988 )'   NN ', NN, MAXIN
    1277. 

  1277.          NN = 0
    1278. 

  1278.          FATAL = .TRUE.
    1279. 

  1279.       END IF
    1280. 

  1280. *
    1281. 

  1281. *     Read the values of M
    1282. 

  1282. *
    1283. 

  1283.       IF( .NOT.( DGX .OR. DXV ) ) THEN
    1284. 

  1284.          READ( NIN, FMT = * )( MVAL( I ), I = 1, NN )
    1285. 

  1285.          IF( SVD ) THEN
    1286. 

  1286.             VNAME = '    M '
    1287. 

  1287.          ELSE
    1288. 

  1288.             VNAME = '    N '
    1289. 

  1289.          END IF
    1290. 

  1290.          DO 20 I = 1, NN
    1291. 

  1291.             IF( MVAL( I ).LT.0 ) THEN
    1292. 

  1292.                WRITE( NOUT, FMT = 9989 )VNAME, MVAL( I ), 0
    1293. 

  1293.                FATAL = .TRUE.
    1294. 

  1294.             ELSE IF( MVAL( I ).GT.NMAX ) THEN
    1295. 

  1295.                WRITE( NOUT, FMT = 9988 )VNAME, MVAL( I ), NMAX
    1296. 

  1296.                FATAL = .TRUE.
    1297. 

  1297.             END IF
    1298. 

  1298.    20    CONTINUE
    1299. 

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

  1300.       END IF
    1301. 

  1301. *
    1302. 

  1302. *     Read the values of P
    1303. 

  1303. *
    1304. 

  1304.       IF( GLM .OR. GQR .OR. GSV .OR. CSD .OR. LSE ) THEN
    1305. 

  1305.          READ( NIN, FMT = * )( PVAL( I ), I = 1, NN )
    1306. 

  1306.          DO 30 I = 1, NN
    1307. 

  1307.             IF( PVAL( I ).LT.0 ) THEN
    1308. 

  1308.                WRITE( NOUT, FMT = 9989 )' P  ', PVAL( I ), 0
    1309. 

  1309.                FATAL = .TRUE.
    1310. 

  1310.             ELSE IF( PVAL( I ).GT.NMAX ) THEN
    1311. 

  1311.                WRITE( NOUT, FMT = 9988 )' P  ', PVAL( I ), NMAX
    1312. 

  1312.                FATAL = .TRUE.
    1313. 

  1313.             END IF
    1314. 

  1314.    30    CONTINUE
    1315. 

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

  1316.       END IF
    1317. 

  1317. *
    1318. 

  1318. *     Read the values of N
    1319. 

  1319. *
    1320. 

  1320.       IF( SVD .OR. DBB .OR. GLM .OR. GQR .OR. GSV .OR. CSD .OR.
    1321. 

  1321.      $    LSE ) THEN
    1322. 

  1322.          READ( NIN, FMT = * )( NVAL( I ), I = 1, NN )
    1323. 

  1323.          DO 40 I = 1, NN
    1324. 

  1324.             IF( NVAL( I ).LT.0 ) THEN
    1325. 

  1325.                WRITE( NOUT, FMT = 9989 )'    N ', NVAL( I ), 0
    1326. 

  1326.                FATAL = .TRUE.
    1327. 

  1327.             ELSE IF( NVAL( I ).GT.NMAX ) THEN
    1328. 

  1328.                WRITE( NOUT, FMT = 9988 )'    N ', NVAL( I ), NMAX
    1329. 

  1329.                FATAL = .TRUE.
    1330. 

  1330.             END IF
    1331. 

  1331.    40    CONTINUE
    1332. 

  1332.       ELSE
    1333. 

  1333.          DO 50 I = 1, NN
    1334. 

  1334.             NVAL( I ) = MVAL( I )
    1335. 

  1335.    50    CONTINUE
    1336. 

  1336.       END IF
    1337. 

  1337.       IF( .NOT.( DGX .OR. DXV ) ) THEN
    1338. 

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

  1339.       ELSE
    1340. 

  1340.          WRITE( NOUT, FMT = 9983 )'N:    ', NN
    1341. 

  1341.       END IF
    1342. 

  1342. *
    1343. 

  1343. *     Read the number of values of K, followed by the values of K
    1344. 

  1344. *
    1345. 

  1345.       IF( DSB .OR. DBB ) THEN
    1346. 

  1346.          READ( NIN, FMT = * )NK
    1347. 

  1347.          READ( NIN, FMT = * )( KVAL( I ), I = 1, NK )
    1348. 

  1348.          DO 60 I = 1, NK
    1349. 

  1349.             IF( KVAL( I ).LT.0 ) THEN
    1350. 

  1350.                WRITE( NOUT, FMT = 9989 )'    K ', KVAL( I ), 0
    1351. 

  1351.                FATAL = .TRUE.
    1352. 

  1352.             ELSE IF( KVAL( I ).GT.NMAX ) THEN
    1353. 

  1353.                WRITE( NOUT, FMT = 9988 )'    K ', KVAL( I ), NMAX
    1354. 

  1354.                FATAL = .TRUE.
    1355. 

  1355.             END IF
    1356. 

  1356.    60    CONTINUE
    1357. 

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

  1358.       END IF
    1359. 

  1359. *
    1360. 

  1360.       IF( DEV .OR. DES .OR. DVX .OR. DSX ) THEN
    1361. 

  1361. *
    1362. 

  1362. *        For the nonsymmetric QR driver routines, only one set of
    1363. 

  1363. *        parameters is allowed.
    1364. 

  1364. *
    1365. 

  1365.          READ( NIN, FMT = * )NBVAL( 1 ), NBMIN( 1 ), NXVAL( 1 ),
    1366. 

  1366.      $      INMIN( 1 ), INWIN( 1 ), INIBL(1), ISHFTS(1), IACC22(1)
    1367. 

  1367.          IF( NBVAL( 1 ).LT.1 ) THEN
    1368. 

  1368.             WRITE( NOUT, FMT = 9989 )'   NB ', NBVAL( 1 ), 1
    1369. 

  1369.             FATAL = .TRUE.
    1370. 

  1370.          ELSE IF( NBMIN( 1 ).LT.1 ) THEN
    1371. 

  1371.             WRITE( NOUT, FMT = 9989 )'NBMIN ', NBMIN( 1 ), 1
    1372. 

  1372.             FATAL = .TRUE.
    1373. 

  1373.          ELSE IF( NXVAL( 1 ).LT.1 ) THEN
    1374. 

  1374.             WRITE( NOUT, FMT = 9989 )'   NX ', NXVAL( 1 ), 1
    1375. 

  1375.             FATAL = .TRUE.
    1376. 

  1376.          ELSE IF( INMIN( 1 ).LT.1 ) THEN
    1377. 

  1377.             WRITE( NOUT, FMT = 9989 )'   INMIN ', INMIN( 1 ), 1
    1378. 

  1378.             FATAL = .TRUE.
    1379. 

  1379.          ELSE IF( INWIN( 1 ).LT.1 ) THEN
    1380. 

  1380.             WRITE( NOUT, FMT = 9989 )'   INWIN ', INWIN( 1 ), 1
    1381. 

  1381.             FATAL = .TRUE.
    1382. 

  1382.          ELSE IF( INIBL( 1 ).LT.1 ) THEN
    1383. 

  1383.             WRITE( NOUT, FMT = 9989 )'   INIBL ', INIBL( 1 ), 1
    1384. 

  1384.             FATAL = .TRUE.
    1385. 

  1385.          ELSE IF( ISHFTS( 1 ).LT.1 ) THEN
    1386. 

  1386.             WRITE( NOUT, FMT = 9989 )'   ISHFTS ', ISHFTS( 1 ), 1
    1387. 

  1387.             FATAL = .TRUE.
    1388. 

  1388.          ELSE IF( IACC22( 1 ).LT.0 ) THEN
    1389. 

  1389.             WRITE( NOUT, FMT = 9989 )'   IACC22 ', IACC22( 1 ), 0
    1390. 

  1390.             FATAL = .TRUE.
    1391. 

  1391.          END IF
    1392. 

  1392.          CALL XLAENV( 1, NBVAL( 1 ) )
    1393. 

  1393.          CALL XLAENV( 2, NBMIN( 1 ) )
    1394. 

  1394.          CALL XLAENV( 3, NXVAL( 1 ) )
    1395. 

  1395.          CALL XLAENV(12, MAX( 11, INMIN( 1 ) ) )
    1396. 

  1396.          CALL XLAENV(13, INWIN( 1 ) )
    1397. 

  1397.          CALL XLAENV(14, INIBL( 1 ) )
    1398. 

  1398.          CALL XLAENV(15, ISHFTS( 1 ) )
    1399. 

  1399.          CALL XLAENV(16, IACC22( 1 ) )
    1400. 

  1400.          WRITE( NOUT, FMT = 9983 )'NB:   ', NBVAL( 1 )
    1401. 

  1401.          WRITE( NOUT, FMT = 9983 )'NBMIN:', NBMIN( 1 )
    1402. 

  1402.          WRITE( NOUT, FMT = 9983 )'NX:   ', NXVAL( 1 )
    1403. 

  1403.          WRITE( NOUT, FMT = 9983 )'INMIN:   ', INMIN( 1 )
    1404. 

  1404.          WRITE( NOUT, FMT = 9983 )'INWIN: ', INWIN( 1 )
    1405. 

  1405.          WRITE( NOUT, FMT = 9983 )'INIBL: ', INIBL( 1 )
    1406. 

  1406.          WRITE( NOUT, FMT = 9983 )'ISHFTS: ', ISHFTS( 1 )
    1407. 

  1407.          WRITE( NOUT, FMT = 9983 )'IACC22: ', IACC22( 1 )
    1408. 

  1408. *
    1409. 

  1409.       ELSEIF( DGS .OR. DGX .OR. DGV .OR.  DXV ) THEN
    1410. 

  1410. *
    1411. 

  1411. *      …
    1412.

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