PageRenderTime 62ms CodeModel.GetById 25ms RepoModel.GetById 1ms app.codeStats 0ms

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

#
FORTRAN Legacy | 2455 lines | 978 code | 0 blank | 1477 comment | 0 complexity | e9380d89d5c4fd0b795f8e6c759bb6ba MD5 | raw file
Possible License(s): BSD-3-Clause

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

    

  1. *> \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.

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