 |
LAPACK 3.3.1
Linear Algebra PACKage
|
|
LAPACK 3.3.1
Linear Algebra PACKage
|
00001 PROGRAM ZCHKEE 00002 * 00003 * -- LAPACK test routine (version 3.3.0) -- 00004 * Univ. of Tennessee, Univ. of California Berkeley and NAG Ltd.. 00005 * February 2007 00006 * 00007 * Purpose 00008 * ======= 00009 * 00010 * ZCHKEE tests the COMPLEX*16 LAPACK subroutines for the matrix 00011 * eigenvalue problem. The test paths in this version are 00012 * 00013 * NEP (Nonsymmetric Eigenvalue Problem): 00014 * Test ZGEHRD, ZUNGHR, ZHSEQR, ZTREVC, ZHSEIN, and ZUNMHR 00015 * 00016 * SEP (Hermitian Eigenvalue Problem): 00017 * Test ZHETRD, ZUNGTR, ZSTEQR, ZSTERF, ZSTEIN, ZSTEDC, 00018 * and drivers ZHEEV(X), ZHBEV(X), ZHPEV(X), 00019 * ZHEEVD, ZHBEVD, ZHPEVD 00020 * 00021 * SVD (Singular Value Decomposition): 00022 * Test ZGEBRD, ZUNGBR, and ZBDSQR 00023 * and the drivers ZGESVD, ZGESDD 00024 * 00025 * ZEV (Nonsymmetric Eigenvalue/eigenvector Driver): 00026 * Test ZGEEV 00027 * 00028 * ZES (Nonsymmetric Schur form Driver): 00029 * Test ZGEES 00030 * 00031 * ZVX (Nonsymmetric Eigenvalue/eigenvector Expert Driver): 00032 * Test ZGEEVX 00033 * 00034 * ZSX (Nonsymmetric Schur form Expert Driver): 00035 * Test ZGEESX 00036 * 00037 * ZGG (Generalized Nonsymmetric Eigenvalue Problem): 00038 * Test ZGGHRD, ZGGBAL, ZGGBAK, ZHGEQZ, and ZTGEVC 00039 * and the driver routines ZGEGS and ZGEGV 00040 * 00041 * ZGS (Generalized Nonsymmetric Schur form Driver): 00042 * Test ZGGES 00043 * 00044 * ZGV (Generalized Nonsymmetric Eigenvalue/eigenvector Driver): 00045 * Test ZGGEV 00046 * 00047 * ZGX (Generalized Nonsymmetric Schur form Expert Driver): 00048 * Test ZGGESX 00049 * 00050 * ZXV (Generalized Nonsymmetric Eigenvalue/eigenvector Expert Driver): 00051 * Test ZGGEVX 00052 * 00053 * ZSG (Hermitian Generalized Eigenvalue Problem): 00054 * Test ZHEGST, ZHEGV, ZHEGVD, ZHEGVX, ZHPGST, ZHPGV, ZHPGVD, 00055 * ZHPGVX, ZHBGST, ZHBGV, ZHBGVD, and ZHBGVX 00056 * 00057 * ZHB (Hermitian Band Eigenvalue Problem): 00058 * Test ZHBTRD 00059 * 00060 * ZBB (Band Singular Value Decomposition): 00061 * Test ZGBBRD 00062 * 00063 * ZEC (Eigencondition estimation): 00064 * Test ZTRSYL, ZTREXC, ZTRSNA, and ZTRSEN 00065 * 00066 * ZBL (Balancing a general matrix) 00067 * Test ZGEBAL 00068 * 00069 * ZBK (Back transformation on a balanced matrix) 00070 * Test ZGEBAK 00071 * 00072 * ZGL (Balancing a matrix pair) 00073 * Test ZGGBAL 00074 * 00075 * ZGK (Back transformation on a matrix pair) 00076 * Test ZGGBAK 00077 * 00078 * GLM (Generalized Linear Regression Model): 00079 * Tests ZGGGLM 00080 * 00081 * GQR (Generalized QR and RQ factorizations): 00082 * Tests ZGGQRF and ZGGRQF 00083 * 00084 * GSV (Generalized Singular Value Decomposition): 00085 * Tests ZGGSVD, ZGGSVP, ZTGSJA, ZLAGS2, ZLAPLL, and ZLAPMT 00086 * 00087 * CSD (CS decomposition): 00088 * Tests ZUNCSD 00089 * 00090 * LSE (Constrained Linear Least Squares): 00091 * Tests ZGGLSE 00092 * 00093 * Each test path has a different set of inputs, but the data sets for 00094 * the driver routines xEV, xES, xVX, and xSX can be concatenated in a 00095 * single input file. The first line of input should contain one of the 00096 * 3-character path names in columns 1-3. The number of remaining lines 00097 * depends on what is found on the first line. 00098 * 00099 * The number of matrix types used in testing is often controllable from 00100 * the input file. The number of matrix types for each path, and the 00101 * test routine that describes them, is as follows: 00102 * 00103 * Path name(s) Types Test routine 00104 * 00105 * ZHS or NEP 21 ZCHKHS 00106 * ZST or SEP 21 ZCHKST (routines) 00107 * 18 ZDRVST (drivers) 00108 * ZBD or SVD 16 ZCHKBD (routines) 00109 * 5 ZDRVBD (drivers) 00110 * ZEV 21 ZDRVEV 00111 * ZES 21 ZDRVES 00112 * ZVX 21 ZDRVVX 00113 * ZSX 21 ZDRVSX 00114 * ZGG 26 ZCHKGG (routines) 00115 * 26 ZDRVGG (drivers) 00116 * ZGS 26 ZDRGES 00117 * ZGX 5 ZDRGSX 00118 * ZGV 26 ZDRGEV 00119 * ZXV 2 ZDRGVX 00120 * ZSG 21 ZDRVSG 00121 * ZHB 15 ZCHKHB 00122 * ZBB 15 ZCHKBB 00123 * ZEC - ZCHKEC 00124 * ZBL - ZCHKBL 00125 * ZBK - ZCHKBK 00126 * ZGL - ZCHKGL 00127 * ZGK - ZCHKGK 00128 * GLM 8 ZCKGLM 00129 * GQR 8 ZCKGQR 00130 * GSV 8 ZCKGSV 00131 * CSD 3 ZCKCSD 00132 * LSE 8 ZCKLSE 00133 * 00134 *----------------------------------------------------------------------- 00135 * 00136 * NEP input file: 00137 * 00138 * line 2: NN, INTEGER 00139 * Number of values of N. 00140 * 00141 * line 3: NVAL, INTEGER array, dimension (NN) 00142 * The values for the matrix dimension N. 00143 * 00144 * line 4: NPARMS, INTEGER 00145 * Number of values of the parameters NB, NBMIN, NX, NS, and 00146 * MAXB. 00147 * 00148 * line 5: NBVAL, INTEGER array, dimension (NPARMS) 00149 * The values for the blocksize NB. 00150 * 00151 * line 6: NBMIN, INTEGER array, dimension (NPARMS) 00152 * The values for the minimum blocksize NBMIN. 00153 * 00154 * line 7: NXVAL, INTEGER array, dimension (NPARMS) 00155 * The values for the crossover point NX. 00156 * 00157 * line 8: INMIN, INTEGER array, dimension (NPARMS) 00158 * LAHQR vs TTQRE crossover point, >= 11 00159 * 00160 * line 9: INWIN, INTEGER array, dimension (NPARMS) 00161 * recommended deflation window size 00162 * 00163 * line 10: INIBL, INTEGER array, dimension (NPARMS) 00164 * nibble crossover point 00165 * 00166 * line 11: ISHFTS, INTEGER array, dimension (NPARMS) 00167 * number of simultaneous shifts) 00168 * 00169 * line 12: IACC22, INTEGER array, dimension (NPARMS) 00170 * select structured matrix multiply: 0, 1 or 2) 00171 * 00172 * line 13: THRESH 00173 * Threshold value for the test ratios. Information will be 00174 * printed about each test for which the test ratio is greater 00175 * than or equal to the threshold. To have all of the test 00176 * ratios printed, use THRESH = 0.0 . 00177 * 00178 * line 14: NEWSD, INTEGER 00179 * A code indicating how to set the random number seed. 00180 * = 0: Set the seed to a default value before each run 00181 * = 1: Initialize the seed to a default value only before the 00182 * first run 00183 * = 2: Like 1, but use the seed values on the next line 00184 * 00185 * If line 14 was 2: 00186 * 00187 * line 15: INTEGER array, dimension (4) 00188 * Four integer values for the random number seed. 00189 * 00190 * lines 15-EOF: The remaining lines occur in sets of 1 or 2 and allow 00191 * the user to specify the matrix types. Each line contains 00192 * a 3-character path name in columns 1-3, and the number 00193 * of matrix types must be the first nonblank item in columns 00194 * 4-80. If the number of matrix types is at least 1 but is 00195 * less than the maximum number of possible types, a second 00196 * line will be read to get the numbers of the matrix types to 00197 * be used. For example, 00198 * NEP 21 00199 * requests all of the matrix types for the nonsymmetric 00200 * eigenvalue problem, while 00201 * NEP 4 00202 * 9 10 11 12 00203 * requests only matrices of type 9, 10, 11, and 12. 00204 * 00205 * The valid 3-character path names are 'NEP' or 'ZHS' for the 00206 * nonsymmetric eigenvalue routines. 00207 * 00208 *----------------------------------------------------------------------- 00209 * 00210 * SEP or ZSG input file: 00211 * 00212 * line 2: NN, INTEGER 00213 * Number of values of N. 00214 * 00215 * line 3: NVAL, INTEGER array, dimension (NN) 00216 * The values for the matrix dimension N. 00217 * 00218 * line 4: NPARMS, INTEGER 00219 * Number of values of the parameters NB, NBMIN, and NX. 00220 * 00221 * line 5: NBVAL, INTEGER array, dimension (NPARMS) 00222 * The values for the blocksize NB. 00223 * 00224 * line 6: NBMIN, INTEGER array, dimension (NPARMS) 00225 * The values for the minimum blocksize NBMIN. 00226 * 00227 * line 7: NXVAL, INTEGER array, dimension (NPARMS) 00228 * The values for the crossover point NX. 00229 * 00230 * line 8: THRESH 00231 * Threshold value for the test ratios. Information will be 00232 * printed about each test for which the test ratio is greater 00233 * than or equal to the threshold. 00234 * 00235 * line 9: TSTCHK, LOGICAL 00236 * Flag indicating whether or not to test the LAPACK routines. 00237 * 00238 * line 10: TSTDRV, LOGICAL 00239 * Flag indicating whether or not to test the driver routines. 00240 * 00241 * line 11: TSTERR, LOGICAL 00242 * Flag indicating whether or not to test the error exits for 00243 * the LAPACK routines and driver routines. 00244 * 00245 * line 12: NEWSD, INTEGER 00246 * A code indicating how to set the random number seed. 00247 * = 0: Set the seed to a default value before each run 00248 * = 1: Initialize the seed to a default value only before the 00249 * first run 00250 * = 2: Like 1, but use the seed values on the next line 00251 * 00252 * If line 12 was 2: 00253 * 00254 * line 13: INTEGER array, dimension (4) 00255 * Four integer values for the random number seed. 00256 * 00257 * lines 13-EOF: Lines specifying matrix types, as for NEP. 00258 * The valid 3-character path names are 'SEP' or 'ZST' for the 00259 * Hermitian eigenvalue routines and driver routines, and 00260 * 'ZSG' for the routines for the Hermitian generalized 00261 * eigenvalue problem. 00262 * 00263 *----------------------------------------------------------------------- 00264 * 00265 * SVD input file: 00266 * 00267 * line 2: NN, INTEGER 00268 * Number of values of M and N. 00269 * 00270 * line 3: MVAL, INTEGER array, dimension (NN) 00271 * The values for the matrix row dimension M. 00272 * 00273 * line 4: NVAL, INTEGER array, dimension (NN) 00274 * The values for the matrix column dimension N. 00275 * 00276 * line 5: NPARMS, INTEGER 00277 * Number of values of the parameter NB, NBMIN, NX, and NRHS. 00278 * 00279 * line 6: NBVAL, INTEGER array, dimension (NPARMS) 00280 * The values for the blocksize NB. 00281 * 00282 * line 7: NBMIN, INTEGER array, dimension (NPARMS) 00283 * The values for the minimum blocksize NBMIN. 00284 * 00285 * line 8: NXVAL, INTEGER array, dimension (NPARMS) 00286 * The values for the crossover point NX. 00287 * 00288 * line 9: NSVAL, INTEGER array, dimension (NPARMS) 00289 * The values for the number of right hand sides NRHS. 00290 * 00291 * line 10: THRESH 00292 * Threshold value for the test ratios. Information will be 00293 * printed about each test for which the test ratio is greater 00294 * than or equal to the threshold. 00295 * 00296 * line 11: TSTCHK, LOGICAL 00297 * Flag indicating whether or not to test the LAPACK routines. 00298 * 00299 * line 12: TSTDRV, LOGICAL 00300 * Flag indicating whether or not to test the driver routines. 00301 * 00302 * line 13: TSTERR, LOGICAL 00303 * Flag indicating whether or not to test the error exits for 00304 * the LAPACK routines and driver routines. 00305 * 00306 * line 14: NEWSD, INTEGER 00307 * A code indicating how to set the random number seed. 00308 * = 0: Set the seed to a default value before each run 00309 * = 1: Initialize the seed to a default value only before the 00310 * first run 00311 * = 2: Like 1, but use the seed values on the next line 00312 * 00313 * If line 14 was 2: 00314 * 00315 * line 15: INTEGER array, dimension (4) 00316 * Four integer values for the random number seed. 00317 * 00318 * lines 15-EOF: Lines specifying matrix types, as for NEP. 00319 * The 3-character path names are 'SVD' or 'ZBD' for both the 00320 * SVD routines and the SVD driver routines. 00321 * 00322 *----------------------------------------------------------------------- 00323 * 00324 * ZEV and ZES data files: 00325 * 00326 * line 1: 'ZEV' or 'ZES' in columns 1 to 3. 00327 * 00328 * line 2: NSIZES, INTEGER 00329 * Number of sizes of matrices to use. Should be at least 0 00330 * and at most 20. If NSIZES = 0, no testing is done 00331 * (although the remaining 3 lines are still read). 00332 * 00333 * line 3: NN, INTEGER array, dimension(NSIZES) 00334 * Dimensions of matrices to be tested. 00335 * 00336 * line 4: NB, NBMIN, NX, NS, NBCOL, INTEGERs 00337 * These integer parameters determine how blocking is done 00338 * (see ILAENV for details) 00339 * NB : block size 00340 * NBMIN : minimum block size 00341 * NX : minimum dimension for blocking 00342 * NS : number of shifts in xHSEQR 00343 * NBCOL : minimum column dimension for blocking 00344 * 00345 * line 5: THRESH, REAL 00346 * The test threshold against which computed residuals are 00347 * compared. Should generally be in the range from 10. to 20. 00348 * If it is 0., all test case data will be printed. 00349 * 00350 * line 6: NEWSD, INTEGER 00351 * A code indicating how to set the random number seed. 00352 * = 0: Set the seed to a default value before each run 00353 * = 1: Initialize the seed to a default value only before the 00354 * first run 00355 * = 2: Like 1, but use the seed values on the next line 00356 * 00357 * If line 6 was 2: 00358 * 00359 * line 7: INTEGER array, dimension (4) 00360 * Four integer values for the random number seed. 00361 * 00362 * lines 8 and following: Lines specifying matrix types, as for NEP. 00363 * The 3-character path name is 'ZEV' to test CGEEV, or 00364 * 'ZES' to test CGEES. 00365 * 00366 *----------------------------------------------------------------------- 00367 * 00368 * The ZVX data has two parts. The first part is identical to ZEV, 00369 * and the second part consists of test matrices with precomputed 00370 * solutions. 00371 * 00372 * line 1: 'ZVX' in columns 1-3. 00373 * 00374 * line 2: NSIZES, INTEGER 00375 * If NSIZES = 0, no testing of randomly generated examples 00376 * is done, but any precomputed examples are tested. 00377 * 00378 * line 3: NN, INTEGER array, dimension(NSIZES) 00379 * 00380 * line 4: NB, NBMIN, NX, NS, NBCOL, INTEGERs 00381 * 00382 * line 5: THRESH, REAL 00383 * 00384 * line 6: NEWSD, INTEGER 00385 * 00386 * If line 6 was 2: 00387 * 00388 * line 7: INTEGER array, dimension (4) 00389 * 00390 * lines 8 and following: The first line contains 'ZVX' in columns 1-3 00391 * followed by the number of matrix types, possibly with 00392 * a second line to specify certain matrix types. 00393 * If the number of matrix types = 0, no testing of randomly 00394 * generated examples is done, but any precomputed examples 00395 * are tested. 00396 * 00397 * remaining lines : Each matrix is stored on 1+N+N**2 lines, where N is 00398 * its dimension. The first line contains the dimension N and 00399 * ISRT (two integers). ISRT indicates whether the last N lines 00400 * are sorted by increasing real part of the eigenvalue 00401 * (ISRT=0) or by increasing imaginary part (ISRT=1). The next 00402 * N**2 lines contain the matrix rowwise, one entry per line. 00403 * The last N lines correspond to each eigenvalue. Each of 00404 * these last N lines contains 4 real values: the real part of 00405 * the eigenvalues, the imaginary part of the eigenvalue, the 00406 * reciprocal condition number of the eigenvalues, and the 00407 * reciprocal condition number of the vector eigenvector. The 00408 * end of data is indicated by dimension N=0. Even if no data 00409 * is to be tested, there must be at least one line containing 00410 * N=0. 00411 * 00412 *----------------------------------------------------------------------- 00413 * 00414 * The ZSX data is like ZVX. The first part is identical to ZEV, and the 00415 * second part consists of test matrices with precomputed solutions. 00416 * 00417 * line 1: 'ZSX' in columns 1-3. 00418 * 00419 * line 2: NSIZES, INTEGER 00420 * If NSIZES = 0, no testing of randomly generated examples 00421 * is done, but any precomputed examples are tested. 00422 * 00423 * line 3: NN, INTEGER array, dimension(NSIZES) 00424 * 00425 * line 4: NB, NBMIN, NX, NS, NBCOL, INTEGERs 00426 * 00427 * line 5: THRESH, REAL 00428 * 00429 * line 6: NEWSD, INTEGER 00430 * 00431 * If line 6 was 2: 00432 * 00433 * line 7: INTEGER array, dimension (4) 00434 * 00435 * lines 8 and following: The first line contains 'ZSX' in columns 1-3 00436 * followed by the number of matrix types, possibly with 00437 * a second line to specify certain matrix types. 00438 * If the number of matrix types = 0, no testing of randomly 00439 * generated examples is done, but any precomputed examples 00440 * are tested. 00441 * 00442 * remaining lines : Each matrix is stored on 3+N**2 lines, where N is 00443 * its dimension. The first line contains the dimension N, the 00444 * dimension M of an invariant subspace, and ISRT. The second 00445 * line contains M integers, identifying the eigenvalues in the 00446 * invariant subspace (by their position in a list of 00447 * eigenvalues ordered by increasing real part (if ISRT=0) or 00448 * by increasing imaginary part (if ISRT=1)). The next N**2 00449 * lines contain the matrix rowwise. The last line contains the 00450 * reciprocal condition number for the average of the selected 00451 * eigenvalues, and the reciprocal condition number for the 00452 * corresponding right invariant subspace. The end of data in 00453 * indicated by a line containing N=0, M=0, and ISRT = 0. Even 00454 * if no data is to be tested, there must be at least one line 00455 * containing N=0, M=0 and ISRT=0. 00456 * 00457 *----------------------------------------------------------------------- 00458 * 00459 * ZGG input file: 00460 * 00461 * line 2: NN, INTEGER 00462 * Number of values of N. 00463 * 00464 * line 3: NVAL, INTEGER array, dimension (NN) 00465 * The values for the matrix dimension N. 00466 * 00467 * line 4: NPARMS, INTEGER 00468 * Number of values of the parameters NB, NBMIN, NBCOL, NS, and 00469 * MAXB. 00470 * 00471 * line 5: NBVAL, INTEGER array, dimension (NPARMS) 00472 * The values for the blocksize NB. 00473 * 00474 * line 6: NBMIN, INTEGER array, dimension (NPARMS) 00475 * The values for NBMIN, the minimum row dimension for blocks. 00476 * 00477 * line 7: NSVAL, INTEGER array, dimension (NPARMS) 00478 * The values for the number of shifts. 00479 * 00480 * line 8: MXBVAL, INTEGER array, dimension (NPARMS) 00481 * The values for MAXB, used in determining minimum blocksize. 00482 * 00483 * line 9: NBCOL, INTEGER array, dimension (NPARMS) 00484 * The values for NBCOL, the minimum column dimension for 00485 * blocks. 00486 * 00487 * line 10: THRESH 00488 * Threshold value for the test ratios. Information will be 00489 * printed about each test for which the test ratio is greater 00490 * than or equal to the threshold. 00491 * 00492 * line 11: TSTCHK, LOGICAL 00493 * Flag indicating whether or not to test the LAPACK routines. 00494 * 00495 * line 12: TSTDRV, LOGICAL 00496 * Flag indicating whether or not to test the driver routines. 00497 * 00498 * line 13: TSTERR, LOGICAL 00499 * Flag indicating whether or not to test the error exits for 00500 * the LAPACK routines and driver routines. 00501 * 00502 * line 14: NEWSD, INTEGER 00503 * A code indicating how to set the random number seed. 00504 * = 0: Set the seed to a default value before each run 00505 * = 1: Initialize the seed to a default value only before the 00506 * first run 00507 * = 2: Like 1, but use the seed values on the next line 00508 * 00509 * If line 14 was 2: 00510 * 00511 * line 15: INTEGER array, dimension (4) 00512 * Four integer values for the random number seed. 00513 * 00514 * lines 16-EOF: Lines specifying matrix types, as for NEP. 00515 * The 3-character path name is 'ZGG' for the generalized 00516 * eigenvalue problem routines and driver routines. 00517 * 00518 *----------------------------------------------------------------------- 00519 * 00520 * ZGS and ZGV input files: 00521 * 00522 * line 1: 'ZGS' or 'ZGV' in columns 1 to 3. 00523 * 00524 * line 2: NN, INTEGER 00525 * Number of values of N. 00526 * 00527 * line 3: NVAL, INTEGER array, dimension(NN) 00528 * Dimensions of matrices to be tested. 00529 * 00530 * line 4: NB, NBMIN, NX, NS, NBCOL, INTEGERs 00531 * These integer parameters determine how blocking is done 00532 * (see ILAENV for details) 00533 * NB : block size 00534 * NBMIN : minimum block size 00535 * NX : minimum dimension for blocking 00536 * NS : number of shifts in xHGEQR 00537 * NBCOL : minimum column dimension for blocking 00538 * 00539 * line 5: THRESH, REAL 00540 * The test threshold against which computed residuals are 00541 * compared. Should generally be in the range from 10. to 20. 00542 * If it is 0., all test case data will be printed. 00543 * 00544 * line 6: TSTERR, LOGICAL 00545 * Flag indicating whether or not to test the error exits. 00546 * 00547 * line 7: NEWSD, INTEGER 00548 * A code indicating how to set the random number seed. 00549 * = 0: Set the seed to a default value before each run 00550 * = 1: Initialize the seed to a default value only before the 00551 * first run 00552 * = 2: Like 1, but use the seed values on the next line 00553 * 00554 * If line 17 was 2: 00555 * 00556 * line 7: INTEGER array, dimension (4) 00557 * Four integer values for the random number seed. 00558 * 00559 * lines 7-EOF: Lines specifying matrix types, as for NEP. 00560 * The 3-character path name is 'ZGS' for the generalized 00561 * eigenvalue problem routines and driver routines. 00562 * 00563 *----------------------------------------------------------------------- 00564 * 00565 * ZGX input file: 00566 * line 1: 'ZGX' in columns 1 to 3. 00567 * 00568 * line 2: N, INTEGER 00569 * Value of N. 00570 * 00571 * line 3: NB, NBMIN, NX, NS, NBCOL, INTEGERs 00572 * These integer parameters determine how blocking is done 00573 * (see ILAENV for details) 00574 * NB : block size 00575 * NBMIN : minimum block size 00576 * NX : minimum dimension for blocking 00577 * NS : number of shifts in xHGEQR 00578 * NBCOL : minimum column dimension for blocking 00579 * 00580 * line 4: THRESH, REAL 00581 * The test threshold against which computed residuals are 00582 * compared. Should generally be in the range from 10. to 20. 00583 * Information will be printed about each test for which the 00584 * test ratio is greater than or equal to the threshold. 00585 * 00586 * line 5: TSTERR, LOGICAL 00587 * Flag indicating whether or not to test the error exits for 00588 * the LAPACK routines and driver routines. 00589 * 00590 * line 6: NEWSD, INTEGER 00591 * A code indicating how to set the random number seed. 00592 * = 0: Set the seed to a default value before each run 00593 * = 1: Initialize the seed to a default value only before the 00594 * first run 00595 * = 2: Like 1, but use the seed values on the next line 00596 * 00597 * If line 6 was 2: 00598 * 00599 * line 7: INTEGER array, dimension (4) 00600 * Four integer values for the random number seed. 00601 * 00602 * If line 2 was 0: 00603 * 00604 * line 7-EOF: Precomputed examples are tested. 00605 * 00606 * remaining lines : Each example is stored on 3+2*N*N lines, where N is 00607 * its dimension. The first line contains the dimension (a 00608 * single integer). The next line contains an integer k such 00609 * that only the last k eigenvalues will be selected and appear 00610 * in the leading diagonal blocks of $A$ and $B$. The next N*N 00611 * lines contain the matrix A, one element per line. The next N*N 00612 * lines contain the matrix B. The last line contains the 00613 * reciprocal of the eigenvalue cluster condition number and the 00614 * reciprocal of the deflating subspace (associated with the 00615 * selected eigencluster) condition number. The end of data is 00616 * indicated by dimension N=0. Even if no data is to be tested, 00617 * there must be at least one line containing N=0. 00618 * 00619 *----------------------------------------------------------------------- 00620 * 00621 * ZXV input files: 00622 * line 1: 'ZXV' in columns 1 to 3. 00623 * 00624 * line 2: N, INTEGER 00625 * Value of N. 00626 * 00627 * line 3: NB, NBMIN, NX, NS, NBCOL, INTEGERs 00628 * These integer parameters determine how blocking is done 00629 * (see ILAENV for details) 00630 * NB : block size 00631 * NBMIN : minimum block size 00632 * NX : minimum dimension for blocking 00633 * NS : number of shifts in xHGEQR 00634 * NBCOL : minimum column dimension for blocking 00635 * 00636 * line 4: THRESH, REAL 00637 * The test threshold against which computed residuals are 00638 * compared. Should generally be in the range from 10. to 20. 00639 * Information will be printed about each test for which the 00640 * test ratio is greater than or equal to the threshold. 00641 * 00642 * line 5: TSTERR, LOGICAL 00643 * Flag indicating whether or not to test the error exits for 00644 * the LAPACK routines and driver routines. 00645 * 00646 * line 6: NEWSD, INTEGER 00647 * A code indicating how to set the random number seed. 00648 * = 0: Set the seed to a default value before each run 00649 * = 1: Initialize the seed to a default value only before the 00650 * first run 00651 * = 2: Like 1, but use the seed values on the next line 00652 * 00653 * If line 6 was 2: 00654 * 00655 * line 7: INTEGER array, dimension (4) 00656 * Four integer values for the random number seed. 00657 * 00658 * If line 2 was 0: 00659 * 00660 * line 7-EOF: Precomputed examples are tested. 00661 * 00662 * remaining lines : Each example is stored on 3+2*N*N lines, where N is 00663 * its dimension. The first line contains the dimension (a 00664 * single integer). The next N*N lines contain the matrix A, one 00665 * element per line. The next N*N lines contain the matrix B. 00666 * The next line contains the reciprocals of the eigenvalue 00667 * condition numbers. The last line contains the reciprocals of 00668 * the eigenvector condition numbers. The end of data is 00669 * indicated by dimension N=0. Even if no data is to be tested, 00670 * there must be at least one line containing N=0. 00671 * 00672 *----------------------------------------------------------------------- 00673 * 00674 * ZHB input file: 00675 * 00676 * line 2: NN, INTEGER 00677 * Number of values of N. 00678 * 00679 * line 3: NVAL, INTEGER array, dimension (NN) 00680 * The values for the matrix dimension N. 00681 * 00682 * line 4: NK, INTEGER 00683 * Number of values of K. 00684 * 00685 * line 5: KVAL, INTEGER array, dimension (NK) 00686 * The values for the matrix dimension K. 00687 * 00688 * line 6: THRESH 00689 * Threshold value for the test ratios. Information will be 00690 * printed about each test for which the test ratio is greater 00691 * than or equal to the threshold. 00692 * 00693 * line 7: NEWSD, INTEGER 00694 * A code indicating how to set the random number seed. 00695 * = 0: Set the seed to a default value before each run 00696 * = 1: Initialize the seed to a default value only before the 00697 * first run 00698 * = 2: Like 1, but use the seed values on the next line 00699 * 00700 * If line 7 was 2: 00701 * 00702 * line 8: INTEGER array, dimension (4) 00703 * Four integer values for the random number seed. 00704 * 00705 * lines 8-EOF: Lines specifying matrix types, as for NEP. 00706 * The 3-character path name is 'ZHB'. 00707 * 00708 *----------------------------------------------------------------------- 00709 * 00710 * ZBB input file: 00711 * 00712 * line 2: NN, INTEGER 00713 * Number of values of M and N. 00714 * 00715 * line 3: MVAL, INTEGER array, dimension (NN) 00716 * The values for the matrix row dimension M. 00717 * 00718 * line 4: NVAL, INTEGER array, dimension (NN) 00719 * The values for the matrix column dimension N. 00720 * 00721 * line 4: NK, INTEGER 00722 * Number of values of K. 00723 * 00724 * line 5: KVAL, INTEGER array, dimension (NK) 00725 * The values for the matrix bandwidth K. 00726 * 00727 * line 6: NPARMS, INTEGER 00728 * Number of values of the parameter NRHS 00729 * 00730 * line 7: NSVAL, INTEGER array, dimension (NPARMS) 00731 * The values for the number of right hand sides NRHS. 00732 * 00733 * line 8: THRESH 00734 * Threshold value for the test ratios. Information will be 00735 * printed about each test for which the test ratio is greater 00736 * than or equal to the threshold. 00737 * 00738 * line 9: NEWSD, INTEGER 00739 * A code indicating how to set the random number seed. 00740 * = 0: Set the seed to a default value before each run 00741 * = 1: Initialize the seed to a default value only before the 00742 * first run 00743 * = 2: Like 1, but use the seed values on the next line 00744 * 00745 * If line 9 was 2: 00746 * 00747 * line 10: INTEGER array, dimension (4) 00748 * Four integer values for the random number seed. 00749 * 00750 * lines 10-EOF: Lines specifying matrix types, as for SVD. 00751 * The 3-character path name is 'ZBB'. 00752 * 00753 *----------------------------------------------------------------------- 00754 * 00755 * ZEC input file: 00756 * 00757 * line 2: THRESH, REAL 00758 * Threshold value for the test ratios. Information will be 00759 * printed about each test for which the test ratio is greater 00760 * than or equal to the threshold. 00761 * 00762 * lines 3-EOF: 00763 * 00764 * Input for testing the eigencondition routines consists of a set of 00765 * specially constructed test cases and their solutions. The data 00766 * format is not intended to be modified by the user. 00767 * 00768 *----------------------------------------------------------------------- 00769 * 00770 * ZBL and ZBK input files: 00771 * 00772 * line 1: 'ZBL' in columns 1-3 to test CGEBAL, or 'ZBK' in 00773 * columns 1-3 to test CGEBAK. 00774 * 00775 * The remaining lines consist of specially constructed test cases. 00776 * 00777 *----------------------------------------------------------------------- 00778 * 00779 * ZGL and ZGK input files: 00780 * 00781 * line 1: 'ZGL' in columns 1-3 to test ZGGBAL, or 'ZGK' in 00782 * columns 1-3 to test ZGGBAK. 00783 * 00784 * The remaining lines consist of specially constructed test cases. 00785 * 00786 *----------------------------------------------------------------------- 00787 * 00788 * GLM data file: 00789 * 00790 * line 1: 'GLM' in columns 1 to 3. 00791 * 00792 * line 2: NN, INTEGER 00793 * Number of values of M, P, and N. 00794 * 00795 * line 3: MVAL, INTEGER array, dimension(NN) 00796 * Values of M (row dimension). 00797 * 00798 * line 4: PVAL, INTEGER array, dimension(NN) 00799 * Values of P (row dimension). 00800 * 00801 * line 5: NVAL, INTEGER array, dimension(NN) 00802 * Values of N (column dimension), note M <= N <= M+P. 00803 * 00804 * line 6: THRESH, REAL 00805 * Threshold value for the test ratios. Information will be 00806 * printed about each test for which the test ratio is greater 00807 * than or equal to the threshold. 00808 * 00809 * line 7: TSTERR, LOGICAL 00810 * Flag indicating whether or not to test the error exits for 00811 * the LAPACK routines and driver routines. 00812 * 00813 * line 8: NEWSD, INTEGER 00814 * A code indicating how to set the random number seed. 00815 * = 0: Set the seed to a default value before each run 00816 * = 1: Initialize the seed to a default value only before the 00817 * first run 00818 * = 2: Like 1, but use the seed values on the next line 00819 * 00820 * If line 8 was 2: 00821 * 00822 * line 9: INTEGER array, dimension (4) 00823 * Four integer values for the random number seed. 00824 * 00825 * lines 9-EOF: Lines specifying matrix types, as for NEP. 00826 * The 3-character path name is 'GLM' for the generalized 00827 * linear regression model routines. 00828 * 00829 *----------------------------------------------------------------------- 00830 * 00831 * GQR data file: 00832 * 00833 * line 1: 'GQR' in columns 1 to 3. 00834 * 00835 * line 2: NN, INTEGER 00836 * Number of values of M, P, and N. 00837 * 00838 * line 3: MVAL, INTEGER array, dimension(NN) 00839 * Values of M. 00840 * 00841 * line 4: PVAL, INTEGER array, dimension(NN) 00842 * Values of P. 00843 * 00844 * line 5: NVAL, INTEGER array, dimension(NN) 00845 * Values of N. 00846 * 00847 * line 6: THRESH, REAL 00848 * Threshold value for the test ratios. Information will be 00849 * printed about each test for which the test ratio is greater 00850 * than or equal to the threshold. 00851 * 00852 * line 7: TSTERR, LOGICAL 00853 * Flag indicating whether or not to test the error exits for 00854 * the LAPACK routines and driver routines. 00855 * 00856 * line 8: NEWSD, INTEGER 00857 * A code indicating how to set the random number seed. 00858 * = 0: Set the seed to a default value before each run 00859 * = 1: Initialize the seed to a default value only before the 00860 * first run 00861 * = 2: Like 1, but use the seed values on the next line 00862 * 00863 * If line 8 was 2: 00864 * 00865 * line 9: INTEGER array, dimension (4) 00866 * Four integer values for the random number seed. 00867 * 00868 * lines 9-EOF: Lines specifying matrix types, as for NEP. 00869 * The 3-character path name is 'GQR' for the generalized 00870 * QR and RQ routines. 00871 * 00872 *----------------------------------------------------------------------- 00873 * 00874 * GSV data file: 00875 * 00876 * line 1: 'GSV' in columns 1 to 3. 00877 * 00878 * line 2: NN, INTEGER 00879 * Number of values of M, P, and N. 00880 * 00881 * line 3: MVAL, INTEGER array, dimension(NN) 00882 * Values of M (row dimension). 00883 * 00884 * line 4: PVAL, INTEGER array, dimension(NN) 00885 * Values of P (row dimension). 00886 * 00887 * line 5: NVAL, INTEGER array, dimension(NN) 00888 * Values of N (column dimension). 00889 * 00890 * line 6: THRESH, REAL 00891 * Threshold value for the test ratios. Information will be 00892 * printed about each test for which the test ratio is greater 00893 * than or equal to the threshold. 00894 * 00895 * line 7: TSTERR, LOGICAL 00896 * Flag indicating whether or not to test the error exits for 00897 * the LAPACK routines and driver routines. 00898 * 00899 * line 8: NEWSD, INTEGER 00900 * A code indicating how to set the random number seed. 00901 * = 0: Set the seed to a default value before each run 00902 * = 1: Initialize the seed to a default value only before the 00903 * first run 00904 * = 2: Like 1, but use the seed values on the next line 00905 * 00906 * If line 8 was 2: 00907 * 00908 * line 9: INTEGER array, dimension (4) 00909 * Four integer values for the random number seed. 00910 * 00911 * lines 9-EOF: Lines specifying matrix types, as for NEP. 00912 * The 3-character path name is 'GSV' for the generalized 00913 * SVD routines. 00914 * 00915 *----------------------------------------------------------------------- 00916 * 00917 * CSD data file: 00918 * 00919 * line 1: 'CSD' in columns 1 to 3. 00920 * 00921 * line 2: NM, INTEGER 00922 * Number of values of M, P, and N. 00923 * 00924 * line 3: MVAL, INTEGER array, dimension(NM) 00925 * Values of M (row and column dimension of orthogonal matrix). 00926 * 00927 * line 4: PVAL, INTEGER array, dimension(NM) 00928 * Values of P (row dimension of top-left block). 00929 * 00930 * line 5: NVAL, INTEGER array, dimension(NM) 00931 * Values of N (column dimension of top-left block). 00932 * 00933 * line 6: THRESH, REAL 00934 * Threshold value for the test ratios. Information will be 00935 * printed about each test for which the test ratio is greater 00936 * than or equal to the threshold. 00937 * 00938 * line 7: TSTERR, LOGICAL 00939 * Flag indicating whether or not to test the error exits for 00940 * the LAPACK routines and driver routines. 00941 * 00942 * line 8: NEWSD, INTEGER 00943 * A code indicating how to set the random number seed. 00944 * = 0: Set the seed to a default value before each run 00945 * = 1: Initialize the seed to a default value only before the 00946 * first run 00947 * = 2: Like 1, but use the seed values on the next line 00948 * 00949 * If line 8 was 2: 00950 * 00951 * line 9: INTEGER array, dimension (4) 00952 * Four integer values for the random number seed. 00953 * 00954 * lines 9-EOF: Lines specifying matrix types, as for NEP. 00955 * The 3-character path name is 'CSD' for the CSD routine. 00956 * 00957 *----------------------------------------------------------------------- 00958 * 00959 * LSE data file: 00960 * 00961 * line 1: 'LSE' in columns 1 to 3. 00962 * 00963 * line 2: NN, INTEGER 00964 * Number of values of M, P, and N. 00965 * 00966 * line 3: MVAL, INTEGER array, dimension(NN) 00967 * Values of M. 00968 * 00969 * line 4: PVAL, INTEGER array, dimension(NN) 00970 * Values of P. 00971 * 00972 * line 5: NVAL, INTEGER array, dimension(NN) 00973 * Values of N, note P <= N <= P+M. 00974 * 00975 * line 6: THRESH, REAL 00976 * Threshold value for the test ratios. Information will be 00977 * printed about each test for which the test ratio is greater 00978 * than or equal to the threshold. 00979 * 00980 * line 7: TSTERR, LOGICAL 00981 * Flag indicating whether or not to test the error exits for 00982 * the LAPACK routines and driver routines. 00983 * 00984 * line 8: NEWSD, INTEGER 00985 * A code indicating how to set the random number seed. 00986 * = 0: Set the seed to a default value before each run 00987 * = 1: Initialize the seed to a default value only before the 00988 * first run 00989 * = 2: Like 1, but use the seed values on the next line 00990 * 00991 * If line 8 was 2: 00992 * 00993 * line 9: INTEGER array, dimension (4) 00994 * Four integer values for the random number seed. 00995 * 00996 * lines 9-EOF: Lines specifying matrix types, as for NEP. 00997 * The 3-character path name is 'GSV' for the generalized 00998 * SVD routines. 00999 * 01000 *----------------------------------------------------------------------- 01001 * 01002 * NMAX is currently set to 132 and must be at least 12 for some of the 01003 * precomputed examples, and LWORK = NMAX*(5*NMAX+20) in the parameter 01004 * statements below. For SVD, we assume NRHS may be as big as N. The 01005 * parameter NEED is set to 14 to allow for 14 N-by-N matrices for ZGG. 01006 * 01007 * ===================================================================== 01008 * 01009 * .. Parameters .. 01010 INTEGER NMAX 01011 PARAMETER ( NMAX = 132 ) 01012 INTEGER NCMAX 01013 PARAMETER ( NCMAX = 20 ) 01014 INTEGER NEED 01015 PARAMETER ( NEED = 14 ) 01016 INTEGER LWORK 01017 PARAMETER ( LWORK = NMAX*( 5*NMAX+20 ) ) 01018 INTEGER LIWORK 01019 PARAMETER ( LIWORK = NMAX*( NMAX+20 ) ) 01020 INTEGER MAXIN 01021 PARAMETER ( MAXIN = 20 ) 01022 INTEGER MAXT 01023 PARAMETER ( MAXT = 30 ) 01024 INTEGER NIN, NOUT 01025 PARAMETER ( NIN = 5, NOUT = 6 ) 01026 * .. 01027 * .. Local Scalars .. 01028 LOGICAL ZBK, ZBL, ZES, ZEV, ZGK, ZGL, ZGS, ZGV, ZGX, 01029 $ ZSX, ZVX, ZXV, CSD, FATAL, GLM, GQR, GSV, LSE, 01030 $ NEP, SEP, SVD, TSTCHK, TSTDIF, TSTDRV, TSTERR, 01031 $ ZBB, ZGG, ZHB 01032 CHARACTER C1 01033 CHARACTER*3 C3, PATH 01034 CHARACTER*32 VNAME 01035 CHARACTER*10 INTSTR 01036 CHARACTER*80 LINE 01037 INTEGER I, I1, IC, INFO, ITMP, K, LENP, MAXTYP, NEWSD, 01038 $ NK, NN, NPARMS, NRHS, NTYPES, 01039 $ VERS_MAJOR, VERS_MINOR, VERS_PATCH 01040 DOUBLE PRECISION EPS, S1, S2, THRESH, THRSHN 01041 * .. 01042 * .. Local Arrays .. 01043 LOGICAL DOTYPE( MAXT ), LOGWRK( NMAX ) 01044 INTEGER IOLDSD( 4 ), ISEED( 4 ), IWORK( LIWORK ), 01045 $ KVAL( MAXIN ), MVAL( MAXIN ), MXBVAL( MAXIN ), 01046 $ NBCOL( MAXIN ), NBMIN( MAXIN ), NBVAL( MAXIN ), 01047 $ NSVAL( MAXIN ), NVAL( MAXIN ), NXVAL( MAXIN ), 01048 $ PVAL( MAXIN ) 01049 INTEGER INMIN( MAXIN ), INWIN( MAXIN ), INIBL( MAXIN ), 01050 $ ISHFTS( MAXIN ), IACC22( MAXIN ) 01051 DOUBLE PRECISION ALPHA( NMAX ), BETA( NMAX ), DR( NMAX, 12 ), 01052 $ RESULT( 500 ), RWORK( LWORK ), S( NMAX*NMAX ) 01053 COMPLEX*16 A( NMAX*NMAX, NEED ), B( NMAX*NMAX, 5 ), 01054 $ C( NCMAX*NCMAX, NCMAX*NCMAX ), DC( NMAX, 6 ), 01055 $ TAUA( NMAX ), TAUB( NMAX ), WORK( LWORK ), 01056 $ X( 5*NMAX ) 01057 * .. 01058 * .. External Functions .. 01059 LOGICAL LSAMEN 01060 DOUBLE PRECISION DLAMCH, DSECND 01061 EXTERNAL LSAMEN, DLAMCH, DSECND 01062 * .. 01063 * .. External Subroutines .. 01064 EXTERNAL ALAREQ, XLAENV, ZCHKBB, ZCHKBD, ZCHKBK, ZCHKBL, 01065 $ ZCHKEC, ZCHKGG, ZCHKGK, ZCHKGL, ZCHKHB, ZCHKHS, 01066 $ ZCHKST, ZCKCSD, ZCKGLM, ZCKGQR, ZCKGSV, ZCKLSE, 01067 $ ZDRGES, ZDRGEV, ZDRGSX, ZDRGVX, ZDRVBD, ZDRVES, 01068 $ ZDRVEV, ZDRVGG, ZDRVSG, ZDRVST, ZDRVSX, ZDRVVX, 01069 $ ZERRBD, ZERRED, ZERRGG, ZERRHS, ZERRST, ILAVER 01070 * .. 01071 * .. Intrinsic Functions .. 01072 INTRINSIC LEN, MIN 01073 * .. 01074 * .. Scalars in Common .. 01075 LOGICAL LERR, OK 01076 CHARACTER*32 SRNAMT 01077 INTEGER INFOT, MAXB, NPROC, NSHIFT, NUNIT, SELDIM, 01078 $ SELOPT 01079 * .. 01080 * .. Arrays in Common .. 01081 LOGICAL SELVAL( 20 ) 01082 INTEGER IPARMS( 100 ) 01083 DOUBLE PRECISION SELWI( 20 ), SELWR( 20 ) 01084 * .. 01085 * .. Common blocks .. 01086 COMMON / CENVIR / NPROC, NSHIFT, MAXB 01087 COMMON / INFOC / INFOT, NUNIT, OK, LERR 01088 COMMON / SRNAMC / SRNAMT 01089 COMMON / SSLCT / SELOPT, SELDIM, SELVAL, SELWR, SELWI 01090 COMMON / CLAENV / IPARMS 01091 * .. 01092 * .. Data statements .. 01093 DATA INTSTR / '0123456789' / 01094 DATA IOLDSD / 0, 0, 0, 1 / 01095 * .. 01096 * .. Executable Statements .. 01097 * 01098 S1 = DSECND( ) 01099 FATAL = .FALSE. 01100 NUNIT = NOUT 01101 * 01102 * Return to here to read multiple sets of data 01103 * 01104 10 CONTINUE 01105 * 01106 * Read the first line and set the 3-character test path 01107 * 01108 READ( NIN, FMT = '(A80)', END = 380 )LINE 01109 PATH = LINE( 1: 3 ) 01110 NEP = LSAMEN( 3, PATH, 'NEP' ) .OR. LSAMEN( 3, PATH, 'ZHS' ) 01111 SEP = LSAMEN( 3, PATH, 'SEP' ) .OR. LSAMEN( 3, PATH, 'ZST' ) .OR. 01112 $ LSAMEN( 3, PATH, 'ZSG' ) 01113 SVD = LSAMEN( 3, PATH, 'SVD' ) .OR. LSAMEN( 3, PATH, 'ZBD' ) 01114 ZEV = LSAMEN( 3, PATH, 'ZEV' ) 01115 ZES = LSAMEN( 3, PATH, 'ZES' ) 01116 ZVX = LSAMEN( 3, PATH, 'ZVX' ) 01117 ZSX = LSAMEN( 3, PATH, 'ZSX' ) 01118 ZGG = LSAMEN( 3, PATH, 'ZGG' ) 01119 ZGS = LSAMEN( 3, PATH, 'ZGS' ) 01120 ZGX = LSAMEN( 3, PATH, 'ZGX' ) 01121 ZGV = LSAMEN( 3, PATH, 'ZGV' ) 01122 ZXV = LSAMEN( 3, PATH, 'ZXV' ) 01123 ZHB = LSAMEN( 3, PATH, 'ZHB' ) 01124 ZBB = LSAMEN( 3, PATH, 'ZBB' ) 01125 GLM = LSAMEN( 3, PATH, 'GLM' ) 01126 GQR = LSAMEN( 3, PATH, 'GQR' ) .OR. LSAMEN( 3, PATH, 'GRQ' ) 01127 GSV = LSAMEN( 3, PATH, 'GSV' ) 01128 CSD = LSAMEN( 3, PATH, 'CSD' ) 01129 LSE = LSAMEN( 3, PATH, 'LSE' ) 01130 ZBL = LSAMEN( 3, PATH, 'ZBL' ) 01131 ZBK = LSAMEN( 3, PATH, 'ZBK' ) 01132 ZGL = LSAMEN( 3, PATH, 'ZGL' ) 01133 ZGK = LSAMEN( 3, PATH, 'ZGK' ) 01134 * 01135 * Report values of parameters. 01136 * 01137 IF( PATH.EQ.' ' ) THEN 01138 GO TO 10 01139 ELSE IF( NEP ) THEN 01140 WRITE( NOUT, FMT = 9987 ) 01141 ELSE IF( SEP ) THEN 01142 WRITE( NOUT, FMT = 9986 ) 01143 ELSE IF( SVD ) THEN 01144 WRITE( NOUT, FMT = 9985 ) 01145 ELSE IF( ZEV ) THEN 01146 WRITE( NOUT, FMT = 9979 ) 01147 ELSE IF( ZES ) THEN 01148 WRITE( NOUT, FMT = 9978 ) 01149 ELSE IF( ZVX ) THEN 01150 WRITE( NOUT, FMT = 9977 ) 01151 ELSE IF( ZSX ) THEN 01152 WRITE( NOUT, FMT = 9976 ) 01153 ELSE IF( ZGG ) THEN 01154 WRITE( NOUT, FMT = 9975 ) 01155 ELSE IF( ZGS ) THEN 01156 WRITE( NOUT, FMT = 9964 ) 01157 ELSE IF( ZGX ) THEN 01158 WRITE( NOUT, FMT = 9965 ) 01159 ELSE IF( ZGV ) THEN 01160 WRITE( NOUT, FMT = 9963 ) 01161 ELSE IF( ZXV ) THEN 01162 WRITE( NOUT, FMT = 9962 ) 01163 ELSE IF( ZHB ) THEN 01164 WRITE( NOUT, FMT = 9974 ) 01165 ELSE IF( ZBB ) THEN 01166 WRITE( NOUT, FMT = 9967 ) 01167 ELSE IF( GLM ) THEN 01168 WRITE( NOUT, FMT = 9971 ) 01169 ELSE IF( GQR ) THEN 01170 WRITE( NOUT, FMT = 9970 ) 01171 ELSE IF( GSV ) THEN 01172 WRITE( NOUT, FMT = 9969 ) 01173 ELSE IF( CSD ) THEN 01174 WRITE( NOUT, FMT = 9960 ) 01175 ELSE IF( LSE ) THEN 01176 WRITE( NOUT, FMT = 9968 ) 01177 ELSE IF( ZBL ) THEN 01178 * 01179 * ZGEBAL: Balancing 01180 * 01181 CALL ZCHKBL( NIN, NOUT ) 01182 GO TO 380 01183 ELSE IF( ZBK ) THEN 01184 * 01185 * ZGEBAK: Back transformation 01186 * 01187 CALL ZCHKBK( NIN, NOUT ) 01188 GO TO 380 01189 ELSE IF( ZGL ) THEN 01190 * 01191 * ZGGBAL: Balancing 01192 * 01193 CALL ZCHKGL( NIN, NOUT ) 01194 GO TO 380 01195 ELSE IF( ZGK ) THEN 01196 * 01197 * ZGGBAK: Back transformation 01198 * 01199 CALL ZCHKGK( NIN, NOUT ) 01200 GO TO 380 01201 ELSE IF( LSAMEN( 3, PATH, 'ZEC' ) ) THEN 01202 * 01203 * ZEC: Eigencondition estimation 01204 * 01205 READ( NIN, FMT = * )THRESH 01206 CALL XLAENV( 1, 1 ) 01207 TSTERR = .TRUE. 01208 CALL ZCHKEC( THRESH, TSTERR, NIN, NOUT ) 01209 GO TO 380 01210 ELSE 01211 WRITE( NOUT, FMT = 9992 )PATH 01212 GO TO 380 01213 END IF 01214 CALL ILAVER( VERS_MAJOR, VERS_MINOR, VERS_PATCH ) 01215 WRITE( NOUT, FMT = 9972 ) VERS_MAJOR, VERS_MINOR, VERS_PATCH 01216 WRITE( NOUT, FMT = 9984 ) 01217 * 01218 * Read the number of values of M, P, and N. 01219 * 01220 READ( NIN, FMT = * )NN 01221 IF( NN.LT.0 ) THEN 01222 WRITE( NOUT, FMT = 9989 )' NN ', NN, 1 01223 NN = 0 01224 FATAL = .TRUE. 01225 ELSE IF( NN.GT.MAXIN ) THEN 01226 WRITE( NOUT, FMT = 9988 )' NN ', NN, MAXIN 01227 NN = 0 01228 FATAL = .TRUE. 01229 END IF 01230 * 01231 * Read the values of M 01232 * 01233 IF( .NOT.( ZGX .OR. ZXV ) ) THEN 01234 READ( NIN, FMT = * )( MVAL( I ), I = 1, NN ) 01235 IF( SVD ) THEN 01236 VNAME = ' M ' 01237 ELSE 01238 VNAME = ' N ' 01239 END IF 01240 DO 20 I = 1, NN 01241 IF( MVAL( I ).LT.0 ) THEN 01242 WRITE( NOUT, FMT = 9989 )VNAME, MVAL( I ), 0 01243 FATAL = .TRUE. 01244 ELSE IF( MVAL( I ).GT.NMAX ) THEN 01245 WRITE( NOUT, FMT = 9988 )VNAME, MVAL( I ), NMAX 01246 FATAL = .TRUE. 01247 END IF 01248 20 CONTINUE 01249 WRITE( NOUT, FMT = 9983 )'M: ', ( MVAL( I ), I = 1, NN ) 01250 END IF 01251 * 01252 * Read the values of P 01253 * 01254 IF( GLM .OR. GQR .OR. GSV .OR. CSD .OR. LSE ) THEN 01255 READ( NIN, FMT = * )( PVAL( I ), I = 1, NN ) 01256 DO 30 I = 1, NN 01257 IF( PVAL( I ).LT.0 ) THEN 01258 WRITE( NOUT, FMT = 9989 )' P ', PVAL( I ), 0 01259 FATAL = .TRUE. 01260 ELSE IF( PVAL( I ).GT.NMAX ) THEN 01261 WRITE( NOUT, FMT = 9988 )' P ', PVAL( I ), NMAX 01262 FATAL = .TRUE. 01263 END IF 01264 30 CONTINUE 01265 WRITE( NOUT, FMT = 9983 )'P: ', ( PVAL( I ), I = 1, NN ) 01266 END IF 01267 * 01268 * Read the values of N 01269 * 01270 IF( SVD .OR. ZBB .OR. GLM .OR. GQR .OR. GSV .OR. CSD .OR. 01271 $ LSE ) THEN 01272 READ( NIN, FMT = * )( NVAL( I ), I = 1, NN ) 01273 DO 40 I = 1, NN 01274 IF( NVAL( I ).LT.0 ) THEN 01275 WRITE( NOUT, FMT = 9989 )' N ', NVAL( I ), 0 01276 FATAL = .TRUE. 01277 ELSE IF( NVAL( I ).GT.NMAX ) THEN 01278 WRITE( NOUT, FMT = 9988 )' N ', NVAL( I ), NMAX 01279 FATAL = .TRUE. 01280 END IF 01281 40 CONTINUE 01282 ELSE 01283 DO 50 I = 1, NN 01284 NVAL( I ) = MVAL( I ) 01285 50 CONTINUE 01286 END IF 01287 IF( .NOT.( ZGX .OR. ZXV ) ) THEN 01288 WRITE( NOUT, FMT = 9983 )'N: ', ( NVAL( I ), I = 1, NN ) 01289 ELSE 01290 WRITE( NOUT, FMT = 9983 )'N: ', NN 01291 END IF 01292 * 01293 * Read the number of values of K, followed by the values of K 01294 * 01295 IF( ZHB .OR. ZBB ) THEN 01296 READ( NIN, FMT = * )NK 01297 READ( NIN, FMT = * )( KVAL( I ), I = 1, NK ) 01298 DO 60 I = 1, NK 01299 IF( KVAL( I ).LT.0 ) THEN 01300 WRITE( NOUT, FMT = 9989 )' K ', KVAL( I ), 0 01301 FATAL = .TRUE. 01302 ELSE IF( KVAL( I ).GT.NMAX ) THEN 01303 WRITE( NOUT, FMT = 9988 )' K ', KVAL( I ), NMAX 01304 FATAL = .TRUE. 01305 END IF 01306 60 CONTINUE 01307 WRITE( NOUT, FMT = 9983 )'K: ', ( KVAL( I ), I = 1, NK ) 01308 END IF 01309 * 01310 IF( ZEV .OR. ZES .OR. ZVX .OR. ZSX ) THEN 01311 * 01312 * For the nonsymmetric QR driver routines, only one set of 01313 * parameters is allowed. 01314 * 01315 READ( NIN, FMT = * )NBVAL( 1 ), NBMIN( 1 ), NXVAL( 1 ), 01316 $ INMIN( 1 ), INWIN( 1 ), INIBL(1), ISHFTS(1), IACC22(1) 01317 IF( NBVAL( 1 ).LT.1 ) THEN 01318 WRITE( NOUT, FMT = 9989 )' NB ', NBVAL( 1 ), 1 01319 FATAL = .TRUE. 01320 ELSE IF( NBMIN( 1 ).LT.1 ) THEN 01321 WRITE( NOUT, FMT = 9989 )'NBMIN ', NBMIN( 1 ), 1 01322 FATAL = .TRUE. 01323 ELSE IF( NXVAL( 1 ).LT.1 ) THEN 01324 WRITE( NOUT, FMT = 9989 )' NX ', NXVAL( 1 ), 1 01325 FATAL = .TRUE. 01326 ELSE IF( INMIN( 1 ).LT.1 ) THEN 01327 WRITE( NOUT, FMT = 9989 )' INMIN ', INMIN( 1 ), 1 01328 FATAL = .TRUE. 01329 ELSE IF( INWIN( 1 ).LT.1 ) THEN 01330 WRITE( NOUT, FMT = 9989 )' INWIN ', INWIN( 1 ), 1 01331 FATAL = .TRUE. 01332 ELSE IF( INIBL( 1 ).LT.1 ) THEN 01333 WRITE( NOUT, FMT = 9989 )' INIBL ', INIBL( 1 ), 1 01334 FATAL = .TRUE. 01335 ELSE IF( ISHFTS( 1 ).LT.1 ) THEN 01336 WRITE( NOUT, FMT = 9989 )' ISHFTS ', ISHFTS( 1 ), 1 01337 FATAL = .TRUE. 01338 ELSE IF( IACC22( 1 ).LT.0 ) THEN 01339 WRITE( NOUT, FMT = 9989 )' IACC22 ', IACC22( 1 ), 0 01340 FATAL = .TRUE. 01341 END IF 01342 CALL XLAENV( 1, NBVAL( 1 ) ) 01343 CALL XLAENV( 2, NBMIN( 1 ) ) 01344 CALL XLAENV( 3, NXVAL( 1 ) ) 01345 CALL XLAENV(12, MAX( 11, INMIN( 1 ) ) ) 01346 CALL XLAENV(13, INWIN( 1 ) ) 01347 CALL XLAENV(14, INIBL( 1 ) ) 01348 CALL XLAENV(15, ISHFTS( 1 ) ) 01349 CALL XLAENV(16, IACC22( 1 ) ) 01350 WRITE( NOUT, FMT = 9983 )'NB: ', NBVAL( 1 ) 01351 WRITE( NOUT, FMT = 9983 )'NBMIN:', NBMIN( 1 ) 01352 WRITE( NOUT, FMT = 9983 )'NX: ', NXVAL( 1 ) 01353 WRITE( NOUT, FMT = 9983 )'INMIN: ', INMIN( 1 ) 01354 WRITE( NOUT, FMT = 9983 )'INWIN: ', INWIN( 1 ) 01355 WRITE( NOUT, FMT = 9983 )'INIBL: ', INIBL( 1 ) 01356 WRITE( NOUT, FMT = 9983 )'ISHFTS: ', ISHFTS( 1 ) 01357 WRITE( NOUT, FMT = 9983 )'IACC22: ', IACC22( 1 ) 01358 * 01359 ELSE IF( ZGS .OR. ZGX .OR. ZGV .OR. ZXV ) THEN 01360 * 01361 * For the nonsymmetric generalized driver routines, only one set of 01362 * parameters is allowed. 01363 * 01364 READ( NIN, FMT = * )NBVAL( 1 ), NBMIN( 1 ), NXVAL( 1 ), 01365 $ NSVAL( 1 ), MXBVAL( 1 ) 01366 IF( NBVAL( 1 ).LT.1 ) THEN 01367 WRITE( NOUT, FMT = 9989 )' NB ', NBVAL( 1 ), 1 01368 FATAL = .TRUE. 01369 ELSE IF( NBMIN( 1 ).LT.1 ) THEN 01370 WRITE( NOUT, FMT = 9989 )'NBMIN ', NBMIN( 1 ), 1 01371 FATAL = .TRUE. 01372 ELSE IF( NXVAL( 1 ).LT.1 ) THEN 01373 WRITE( NOUT, FMT = 9989 )' NX ', NXVAL( 1 ), 1 01374 FATAL = .TRUE. 01375 ELSE IF( NSVAL( 1 ).LT.2 ) THEN 01376 WRITE( NOUT, FMT = 9989 )' NS ', NSVAL( 1 ), 2 01377 FATAL = .TRUE. 01378 ELSE IF( MXBVAL( 1 ).LT.1 ) THEN 01379 WRITE( NOUT, FMT = 9989 )' MAXB ', MXBVAL( 1 ), 1 01380 FATAL = .TRUE. 01381 END IF 01382 CALL XLAENV( 1, NBVAL( 1 ) ) 01383 CALL XLAENV( 2, NBMIN( 1 ) ) 01384 CALL XLAENV( 3, NXVAL( 1 ) ) 01385 CALL XLAENV( 4, NSVAL( 1 ) ) 01386 CALL XLAENV( 8, MXBVAL( 1 ) ) 01387 WRITE( NOUT, FMT = 9983 )'NB: ', NBVAL( 1 ) 01388 WRITE( NOUT, FMT = 9983 )'NBMIN:', NBMIN( 1 ) 01389 WRITE( NOUT, FMT = 9983 )'NX: ', NXVAL( 1 ) 01390 WRITE( NOUT, FMT = 9983 )'NS: ', NSVAL( 1 ) 01391 WRITE( NOUT, FMT = 9983 )'MAXB: ', MXBVAL( 1 ) 01392 ELSE IF( .NOT.ZHB .AND. .NOT.GLM .AND. .NOT.GQR .AND. .NOT. 01393 $ GSV .AND. .NOT.CSD .AND. .NOT.LSE ) THEN 01394 * 01395 * For the other paths, the number of parameters can be varied 01396 * from the input file. Read the number of parameter values. 01397 * 01398 READ( NIN, FMT = * )NPARMS 01399 IF( NPARMS.LT.1 ) THEN 01400 WRITE( NOUT, FMT = 9989 )'NPARMS', NPARMS, 1 01401 NPARMS = 0 01402 FATAL = .TRUE. 01403 ELSE IF( NPARMS.GT.MAXIN ) THEN 01404 WRITE( NOUT, FMT = 9988 )'NPARMS', NPARMS, MAXIN 01405 NPARMS = 0 01406 FATAL = .TRUE. 01407 END IF 01408 * 01409 * Read the values of NB 01410 * 01411 IF( .NOT.ZBB ) THEN 01412 READ( NIN, FMT = * )( NBVAL( I ), I = 1, NPARMS ) 01413 DO 70 I = 1, NPARMS 01414 IF( NBVAL( I ).LT.0 ) THEN 01415 WRITE( NOUT, FMT = 9989 )' NB ', NBVAL( I ), 0 01416 FATAL = .TRUE. 01417 ELSE IF( NBVAL( I ).GT.NMAX ) THEN 01418 WRITE( NOUT, FMT = 9988 )' NB ', NBVAL( I ), NMAX 01419 FATAL = .TRUE. 01420 END IF 01421 70 CONTINUE 01422 WRITE( NOUT, FMT = 9983 )'NB: ', 01423 $ ( NBVAL( I ), I = 1, NPARMS ) 01424 END IF 01425 * 01426 * Read the values of NBMIN 01427 * 01428 IF( NEP .OR. SEP .OR. SVD .OR. ZGG ) THEN 01429 READ( NIN, FMT = * )( NBMIN( I ), I = 1, NPARMS ) 01430 DO 80 I = 1, NPARMS 01431 IF( NBMIN( I ).LT.0 ) THEN 01432 WRITE( NOUT, FMT = 9989 )'NBMIN ', NBMIN( I ), 0 01433 FATAL = .TRUE. 01434 ELSE IF( NBMIN( I ).GT.NMAX ) THEN 01435 WRITE( NOUT, FMT = 9988 )'NBMIN ', NBMIN( I ), NMAX 01436 FATAL = .TRUE. 01437 END IF 01438 80 CONTINUE 01439 WRITE( NOUT, FMT = 9983 )'NBMIN:', 01440 $ ( NBMIN( I ), I = 1, NPARMS ) 01441 ELSE 01442 DO 90 I = 1, NPARMS 01443 NBMIN( I ) = 1 01444 90 CONTINUE 01445 END IF 01446 * 01447 * Read the values of NX 01448 * 01449 IF( NEP .OR. SEP .OR. SVD ) THEN 01450 READ( NIN, FMT = * )( NXVAL( I ), I = 1, NPARMS ) 01451 DO 100 I = 1, NPARMS 01452 IF( NXVAL( I ).LT.0 ) THEN 01453 WRITE( NOUT, FMT = 9989 )' NX ', NXVAL( I ), 0 01454 FATAL = .TRUE. 01455 ELSE IF( NXVAL( I ).GT.NMAX ) THEN 01456 WRITE( NOUT, FMT = 9988 )' NX ', NXVAL( I ), NMAX 01457 FATAL = .TRUE. 01458 END IF 01459 100 CONTINUE 01460 WRITE( NOUT, FMT = 9983 )'NX: ', 01461 $ ( NXVAL( I ), I = 1, NPARMS ) 01462 ELSE 01463 DO 110 I = 1, NPARMS 01464 NXVAL( I ) = 1 01465 110 CONTINUE 01466 END IF 01467 * 01468 * Read the values of NSHIFT (if ZGG) or NRHS (if SVD 01469 * or ZBB). 01470 * 01471 IF( SVD .OR. ZBB .OR. ZGG ) THEN 01472 READ( NIN, FMT = * )( NSVAL( I ), I = 1, NPARMS ) 01473 DO 120 I = 1, NPARMS 01474 IF( NSVAL( I ).LT.0 ) THEN 01475 WRITE( NOUT, FMT = 9989 )' NS ', NSVAL( I ), 0 01476 FATAL = .TRUE. 01477 ELSE IF( NSVAL( I ).GT.NMAX ) THEN 01478 WRITE( NOUT, FMT = 9988 )' NS ', NSVAL( I ), NMAX 01479 FATAL = .TRUE. 01480 END IF 01481 120 CONTINUE 01482 WRITE( NOUT, FMT = 9983 )'NS: ', 01483 $ ( NSVAL( I ), I = 1, NPARMS ) 01484 ELSE 01485 DO 130 I = 1, NPARMS 01486 NSVAL( I ) = 1 01487 130 CONTINUE 01488 END IF 01489 * 01490 * Read the values for MAXB. 01491 * 01492 IF( ZGG ) THEN 01493 READ( NIN, FMT = * )( MXBVAL( I ), I = 1, NPARMS ) 01494 DO 140 I = 1, NPARMS 01495 IF( MXBVAL( I ).LT.0 ) THEN 01496 WRITE( NOUT, FMT = 9989 )' MAXB ', MXBVAL( I ), 0 01497 FATAL = .TRUE. 01498 ELSE IF( MXBVAL( I ).GT.NMAX ) THEN 01499 WRITE( NOUT, FMT = 9988 )' MAXB ', MXBVAL( I ), NMAX 01500 FATAL = .TRUE. 01501 END IF 01502 140 CONTINUE 01503 WRITE( NOUT, FMT = 9983 )'MAXB: ', 01504 $ ( MXBVAL( I ), I = 1, NPARMS ) 01505 ELSE 01506 DO 150 I = 1, NPARMS 01507 MXBVAL( I ) = 1 01508 150 CONTINUE 01509 END IF 01510 * 01511 * Read the values for INMIN. 01512 * 01513 IF( NEP ) THEN 01514 READ( NIN, FMT = * )( INMIN( I ), I = 1, NPARMS ) 01515 DO 540 I = 1, NPARMS 01516 IF( INMIN( I ).LT.0 ) THEN 01517 WRITE( NOUT, FMT = 9989 )' INMIN ', INMIN( I ), 0 01518 FATAL = .TRUE. 01519 END IF 01520 540 CONTINUE 01521 WRITE( NOUT, FMT = 9983 )'INMIN: ', 01522 $ ( INMIN( I ), I = 1, NPARMS ) 01523 ELSE 01524 DO 550 I = 1, NPARMS 01525 INMIN( I ) = 1 01526 550 CONTINUE 01527 END IF 01528 * 01529 * Read the values for INWIN. 01530 * 01531 IF( NEP ) THEN 01532 READ( NIN, FMT = * )( INWIN( I ), I = 1, NPARMS ) 01533 DO 560 I = 1, NPARMS 01534 IF( INWIN( I ).LT.0 ) THEN 01535 WRITE( NOUT, FMT = 9989 )' INWIN ', INWIN( I ), 0 01536 FATAL = .TRUE. 01537 END IF 01538 560 CONTINUE 01539 WRITE( NOUT, FMT = 9983 )'INWIN: ', 01540 $ ( INWIN( I ), I = 1, NPARMS ) 01541 ELSE 01542 DO 570 I = 1, NPARMS 01543 INWIN( I ) = 1 01544 570 CONTINUE 01545 END IF 01546 * 01547 * Read the values for INIBL. 01548 * 01549 IF( NEP ) THEN 01550 READ( NIN, FMT = * )( INIBL( I ), I = 1, NPARMS ) 01551 DO 580 I = 1, NPARMS 01552 IF( INIBL( I ).LT.0 ) THEN 01553 WRITE( NOUT, FMT = 9989 )' INIBL ', INIBL( I ), 0 01554 FATAL = .TRUE. 01555 END IF 01556 580 CONTINUE 01557 WRITE( NOUT, FMT = 9983 )'INIBL: ', 01558 $ ( INIBL( I ), I = 1, NPARMS ) 01559 ELSE 01560 DO 590 I = 1, NPARMS 01561 INIBL( I ) = 1 01562 590 CONTINUE 01563 END IF 01564 * 01565 * Read the values for ISHFTS. 01566 * 01567 IF( NEP ) THEN 01568 READ( NIN, FMT = * )( ISHFTS( I ), I = 1, NPARMS ) 01569 DO 600 I = 1, NPARMS 01570 IF( ISHFTS( I ).LT.0 ) THEN 01571 WRITE( NOUT, FMT = 9989 )' ISHFTS ', ISHFTS( I ), 0 01572 FATAL = .TRUE. 01573 END IF 01574 600 CONTINUE 01575 WRITE( NOUT, FMT = 9983 )'ISHFTS: ', 01576 $ ( ISHFTS( I ), I = 1, NPARMS ) 01577 ELSE 01578 DO 610 I = 1, NPARMS 01579 ISHFTS( I ) = 1 01580 610 CONTINUE 01581 END IF 01582 * 01583 * Read the values for IACC22. 01584 * 01585 IF( NEP ) THEN 01586 READ( NIN, FMT = * )( IACC22( I ), I = 1, NPARMS ) 01587 DO 620 I = 1, NPARMS 01588 IF( IACC22( I ).LT.0 ) THEN 01589 WRITE( NOUT, FMT = 9989 )' IACC22 ', IACC22( I ), 0 01590 FATAL = .TRUE. 01591 END IF 01592 620 CONTINUE 01593 WRITE( NOUT, FMT = 9983 )'IACC22: ', 01594 $ ( IACC22( I ), I = 1, NPARMS ) 01595 ELSE 01596 DO 630 I = 1, NPARMS 01597 IACC22( I ) = 1 01598 630 CONTINUE 01599 END IF 01600 * 01601 * Read the values for NBCOL. 01602 * 01603 IF( ZGG ) THEN 01604 READ( NIN, FMT = * )( NBCOL( I ), I = 1, NPARMS ) 01605 DO 160 I = 1, NPARMS 01606 IF( NBCOL( I ).LT.0 ) THEN 01607 WRITE( NOUT, FMT = 9989 )'NBCOL ', NBCOL( I ), 0 01608 FATAL = .TRUE. 01609 ELSE IF( NBCOL( I ).GT.NMAX ) THEN 01610 WRITE( NOUT, FMT = 9988 )'NBCOL ', NBCOL( I ), NMAX 01611 FATAL = .TRUE. 01612 END IF 01613 160 CONTINUE 01614 WRITE( NOUT, FMT = 9983 )'NBCOL:', 01615 $ ( NBCOL( I ), I = 1, NPARMS ) 01616 ELSE 01617 DO 170 I = 1, NPARMS 01618 NBCOL( I ) = 1 01619 170 CONTINUE 01620 END IF 01621 END IF 01622 * 01623 * Calculate and print the machine dependent constants. 01624 * 01625 WRITE( NOUT, FMT = * ) 01626 EPS = DLAMCH( 'Underflow threshold' ) 01627 WRITE( NOUT, FMT = 9981 )'underflow', EPS 01628 EPS = DLAMCH( 'Overflow threshold' ) 01629 WRITE( NOUT, FMT = 9981 )'overflow ', EPS 01630 EPS = DLAMCH( 'Epsilon' ) 01631 WRITE( NOUT, FMT = 9981 )'precision', EPS 01632 * 01633 * Read the threshold value for the test ratios. 01634 * 01635 READ( NIN, FMT = * )THRESH 01636 WRITE( NOUT, FMT = 9982 )THRESH 01637 IF( SEP .OR. SVD .OR. ZGG ) THEN 01638 * 01639 * Read the flag that indicates whether to test LAPACK routines. 01640 * 01641 READ( NIN, FMT = * )TSTCHK 01642 * 01643 * Read the flag that indicates whether to test driver routines. 01644 * 01645 READ( NIN, FMT = * )TSTDRV 01646 END IF 01647 * 01648 * Read the flag that indicates whether to test the error exits. 01649 * 01650 READ( NIN, FMT = * )TSTERR 01651 * 01652 * Read the code describing how to set the random number seed. 01653 * 01654 READ( NIN, FMT = * )NEWSD 01655 * 01656 * If NEWSD = 2, read another line with 4 integers for the seed. 01657 * 01658 IF( NEWSD.EQ.2 ) 01659 $ READ( NIN, FMT = * )( IOLDSD( I ), I = 1, 4 ) 01660 * 01661 DO 180 I = 1, 4 01662 ISEED( I ) = IOLDSD( I ) 01663 180 CONTINUE 01664 * 01665 IF( FATAL ) THEN 01666 WRITE( NOUT, FMT = 9999 ) 01667 STOP 01668 END IF 01669 * 01670 * Read the input lines indicating the test path and its parameters. 01671 * The first three characters indicate the test path, and the number 01672 * of test matrix types must be the first nonblank item in columns 01673 * 4-80. 01674 * 01675 190 CONTINUE 01676 * 01677 IF( .NOT.( ZGX .OR. ZXV ) ) THEN 01678 * 01679 200 CONTINUE 01680 READ( NIN, FMT = '(A80)', END = 380 )LINE 01681 C3 = LINE( 1: 3 ) 01682 LENP = LEN( LINE ) 01683 I = 3 01684 ITMP = 0 01685 I1 = 0 01686 210 CONTINUE 01687 I = I + 1 01688 IF( I.GT.LENP ) THEN 01689 IF( I1.GT.0 ) THEN 01690 GO TO 240 01691 ELSE 01692 NTYPES = MAXT 01693 GO TO 240 01694 END IF 01695 END IF 01696 IF( LINE( I: I ).NE.' ' .AND. LINE( I: I ).NE.',' ) THEN 01697 I1 = I 01698 C1 = LINE( I1: I1 ) 01699 * 01700 * Check that a valid integer was read 01701 * 01702 DO 220 K = 1, 10 01703 IF( C1.EQ.INTSTR( K: K ) ) THEN 01704 IC = K - 1 01705 GO TO 230 01706 END IF 01707 220 CONTINUE 01708 WRITE( NOUT, FMT = 9991 )I, LINE 01709 GO TO 200 01710 230 CONTINUE 01711 ITMP = 10*ITMP + IC 01712 GO TO 210 01713 ELSE IF( I1.GT.0 ) THEN 01714 GO TO 240 01715 ELSE 01716 GO TO 210 01717 END IF 01718 240 CONTINUE 01719 NTYPES = ITMP 01720 * 01721 * Skip the tests if NTYPES is <= 0. 01722 * 01723 IF( .NOT.( ZEV .OR. ZES .OR. ZVX .OR. ZSX .OR. ZGV .OR. 01724 $ ZGS ) .AND. NTYPES.LE.0 ) THEN 01725 WRITE( NOUT, FMT = 9990 )C3 01726 GO TO 200 01727 END IF 01728 * 01729 ELSE 01730 IF( ZGX ) 01731 $ C3 = 'ZGX' 01732 IF( ZXV ) 01733 $ C3 = 'ZXV' 01734 END IF 01735 * 01736 * Reset the random number seed. 01737 * 01738 IF( NEWSD.EQ.0 ) THEN 01739 DO 250 K = 1, 4 01740 ISEED( K ) = IOLDSD( K ) 01741 250 CONTINUE 01742 END IF 01743 * 01744 IF( LSAMEN( 3, C3, 'ZHS' ) .OR. LSAMEN( 3, C3, 'NEP' ) ) THEN 01745 * 01746 * ------------------------------------- 01747 * NEP: Nonsymmetric Eigenvalue Problem 01748 * ------------------------------------- 01749 * Vary the parameters 01750 * NB = block size 01751 * NBMIN = minimum block size 01752 * NX = crossover point 01753 * NS = number of shifts 01754 * MAXB = minimum submatrix size 01755 * 01756 MAXTYP = 21 01757 NTYPES = MIN( MAXTYP, NTYPES ) 01758 CALL ALAREQ( C3, NTYPES, DOTYPE, MAXTYP, NIN, NOUT ) 01759 CALL XLAENV( 1, 1 ) 01760 IF( TSTERR ) 01761 $ CALL ZERRHS( 'ZHSEQR', NOUT ) 01762 DO 270 I = 1, NPARMS 01763 CALL XLAENV( 1, NBVAL( I ) ) 01764 CALL XLAENV( 2, NBMIN( I ) ) 01765 CALL XLAENV( 3, NXVAL( I ) ) 01766 CALL XLAENV(12, MAX( 11, INMIN( I ) ) ) 01767 CALL XLAENV(13, INWIN( I ) ) 01768 CALL XLAENV(14, INIBL( I ) ) 01769 CALL XLAENV(15, ISHFTS( I ) ) 01770 CALL XLAENV(16, IACC22( I ) ) 01771 * 01772 IF( NEWSD.EQ.0 ) THEN 01773 DO 260 K = 1, 4 01774 ISEED( K ) = IOLDSD( K ) 01775 260 CONTINUE 01776 END IF 01777 WRITE( NOUT, FMT = 9961 )C3, NBVAL( I ), NBMIN( I ), 01778 $ NXVAL( I ), MAX( 11, INMIN(I)), 01779 $ INWIN( I ), INIBL( I ), ISHFTS( I ), IACC22( I ) 01780 CALL ZCHKHS( NN, NVAL, MAXTYP, DOTYPE, ISEED, THRESH, NOUT, 01781 $ A( 1, 1 ), NMAX, A( 1, 2 ), A( 1, 3 ), 01782 $ A( 1, 4 ), A( 1, 5 ), NMAX, A( 1, 6 ), 01783 $ A( 1, 7 ), DC( 1, 1 ), DC( 1, 2 ), A( 1, 8 ), 01784 $ A( 1, 9 ), A( 1, 10 ), A( 1, 11 ), A( 1, 12 ), 01785 $ DC( 1, 3 ), WORK, LWORK, RWORK, IWORK, LOGWRK, 01786 $ RESULT, INFO ) 01787 IF( INFO.NE.0 ) 01788 $ WRITE( NOUT, FMT = 9980 )'ZCHKHS', INFO 01789 270 CONTINUE 01790 * 01791 ELSE IF( LSAMEN( 3, C3, 'ZST' ) .OR. LSAMEN( 3, C3, 'SEP' ) ) THEN 01792 * 01793 * ---------------------------------- 01794 * SEP: Symmetric Eigenvalue Problem 01795 * ---------------------------------- 01796 * Vary the parameters 01797 * NB = block size 01798 * NBMIN = minimum block size 01799 * NX = crossover point 01800 * 01801 MAXTYP = 21 01802 NTYPES = MIN( MAXTYP, NTYPES ) 01803 CALL ALAREQ( C3, NTYPES, DOTYPE, MAXTYP, NIN, NOUT ) 01804 CALL XLAENV( 1, 1 ) 01805 CALL XLAENV( 9, 25 ) 01806 IF( TSTERR ) 01807 $ CALL ZERRST( 'ZST', NOUT ) 01808 DO 290 I = 1, NPARMS 01809 CALL XLAENV( 1, NBVAL( I ) ) 01810 CALL XLAENV( 2, NBMIN( I ) ) 01811 CALL XLAENV( 3, NXVAL( I ) ) 01812 * 01813 IF( NEWSD.EQ.0 ) THEN 01814 DO 280 K = 1, 4 01815 ISEED( K ) = IOLDSD( K ) 01816 280 CONTINUE 01817 END IF 01818 WRITE( NOUT, FMT = 9997 )C3, NBVAL( I ), NBMIN( I ), 01819 $ NXVAL( I ) 01820 IF( TSTCHK ) THEN 01821 CALL ZCHKST( NN, NVAL, MAXTYP, DOTYPE, ISEED, THRESH, 01822 $ NOUT, A( 1, 1 ), NMAX, A( 1, 2 ), 01823 $ DR( 1, 1 ), DR( 1, 2 ), DR( 1, 3 ), 01824 $ DR( 1, 4 ), DR( 1, 5 ), DR( 1, 6 ), 01825 $ DR( 1, 7 ), DR( 1, 8 ), DR( 1, 9 ), 01826 $ DR( 1, 10 ), DR( 1, 11 ), A( 1, 3 ), NMAX, 01827 $ A( 1, 4 ), A( 1, 5 ), DC( 1, 1 ), A( 1, 6 ), 01828 $ WORK, LWORK, RWORK, LWORK, IWORK, LIWORK, 01829 $ RESULT, INFO ) 01830 IF( INFO.NE.0 ) 01831 $ WRITE( NOUT, FMT = 9980 )'ZCHKST', INFO 01832 END IF 01833 IF( TSTDRV ) THEN 01834 CALL ZDRVST( NN, NVAL, 18, DOTYPE, ISEED, THRESH, NOUT, 01835 $ A( 1, 1 ), NMAX, DR( 1, 3 ), DR( 1, 4 ), 01836 $ DR( 1, 5 ), DR( 1, 8 ), DR( 1, 9 ), 01837 $ DR( 1, 10 ), A( 1, 2 ), NMAX, A( 1, 3 ), 01838 $ DC( 1, 1 ), A( 1, 4 ), WORK, LWORK, RWORK, 01839 $ LWORK, IWORK, LIWORK, RESULT, INFO ) 01840 IF( INFO.NE.0 ) 01841 $ WRITE( NOUT, FMT = 9980 )'ZDRVST', INFO 01842 END IF 01843 290 CONTINUE 01844 * 01845 ELSE IF( LSAMEN( 3, C3, 'ZSG' ) ) THEN 01846 * 01847 * ---------------------------------------------- 01848 * ZSG: Hermitian Generalized Eigenvalue Problem 01849 * ---------------------------------------------- 01850 * Vary the parameters 01851 * NB = block size 01852 * NBMIN = minimum block size 01853 * NX = crossover point 01854 * 01855 MAXTYP = 21 01856 NTYPES = MIN( MAXTYP, NTYPES ) 01857 CALL ALAREQ( C3, NTYPES, DOTYPE, MAXTYP, NIN, NOUT ) 01858 CALL XLAENV( 9, 25 ) 01859 DO 310 I = 1, NPARMS 01860 CALL XLAENV( 1, NBVAL( I ) ) 01861 CALL XLAENV( 2, NBMIN( I ) ) 01862 CALL XLAENV( 3, NXVAL( I ) ) 01863 * 01864 IF( NEWSD.EQ.0 ) THEN 01865 DO 300 K = 1, 4 01866 ISEED( K ) = IOLDSD( K ) 01867 300 CONTINUE 01868 END IF 01869 WRITE( NOUT, FMT = 9997 )C3, NBVAL( I ), NBMIN( I ), 01870 $ NXVAL( I ) 01871 IF( TSTCHK ) THEN 01872 CALL ZDRVSG( NN, NVAL, MAXTYP, DOTYPE, ISEED, THRESH, 01873 $ NOUT, A( 1, 1 ), NMAX, A( 1, 2 ), NMAX, 01874 $ DR( 1, 3 ), A( 1, 3 ), NMAX, A( 1, 4 ), 01875 $ A( 1, 5 ), A( 1, 6 ), A( 1, 7 ), WORK, 01876 $ LWORK, RWORK, LWORK, IWORK, LIWORK, RESULT, 01877 $ INFO ) 01878 IF( INFO.NE.0 ) 01879 $ WRITE( NOUT, FMT = 9980 )'ZDRVSG', INFO 01880 END IF 01881 310 CONTINUE 01882 * 01883 ELSE IF( LSAMEN( 3, C3, 'ZBD' ) .OR. LSAMEN( 3, C3, 'SVD' ) ) THEN 01884 * 01885 * ---------------------------------- 01886 * SVD: Singular Value Decomposition 01887 * ---------------------------------- 01888 * Vary the parameters 01889 * NB = block size 01890 * NBMIN = minimum block size 01891 * NX = crossover point 01892 * NRHS = number of right hand sides 01893 * 01894 MAXTYP = 16 01895 NTYPES = MIN( MAXTYP, NTYPES ) 01896 CALL ALAREQ( C3, NTYPES, DOTYPE, MAXTYP, NIN, NOUT ) 01897 CALL XLAENV( 9, 25 ) 01898 * 01899 * Test the error exits 01900 * 01901 CALL XLAENV( 1, 1 ) 01902 IF( TSTERR .AND. TSTCHK ) 01903 $ CALL ZERRBD( 'ZBD', NOUT ) 01904 IF( TSTERR .AND. TSTDRV ) 01905 $ CALL ZERRED( 'ZBD', NOUT ) 01906 * 01907 DO 330 I = 1, NPARMS 01908 NRHS = NSVAL( I ) 01909 CALL XLAENV( 1, NBVAL( I ) ) 01910 CALL XLAENV( 2, NBMIN( I ) ) 01911 CALL XLAENV( 3, NXVAL( I ) ) 01912 IF( NEWSD.EQ.0 ) THEN 01913 DO 320 K = 1, 4 01914 ISEED( K ) = IOLDSD( K ) 01915 320 CONTINUE 01916 END IF 01917 WRITE( NOUT, FMT = 9995 )C3, NBVAL( I ), NBMIN( I ), 01918 $ NXVAL( I ), NRHS 01919 IF( TSTCHK ) THEN 01920 CALL ZCHKBD( NN, MVAL, NVAL, MAXTYP, DOTYPE, NRHS, ISEED, 01921 $ THRESH, A( 1, 1 ), NMAX, DR( 1, 1 ), 01922 $ DR( 1, 2 ), DR( 1, 3 ), DR( 1, 4 ), 01923 $ A( 1, 2 ), NMAX, A( 1, 3 ), A( 1, 4 ), 01924 $ A( 1, 5 ), NMAX, A( 1, 6 ), NMAX, A( 1, 7 ), 01925 $ A( 1, 8 ), WORK, LWORK, RWORK, NOUT, INFO ) 01926 IF( INFO.NE.0 ) 01927 $ WRITE( NOUT, FMT = 9980 )'ZCHKBD', INFO 01928 END IF 01929 IF( TSTDRV ) 01930 $ CALL ZDRVBD( NN, MVAL, NVAL, MAXTYP, DOTYPE, ISEED, 01931 $ THRESH, A( 1, 1 ), NMAX, A( 1, 2 ), NMAX, 01932 $ A( 1, 3 ), NMAX, A( 1, 4 ), A( 1, 5 ), 01933 $ A( 1, 6 ), DR( 1, 1 ), DR( 1, 2 ), 01934 $ DR( 1, 3 ), WORK, LWORK, RWORK, IWORK, NOUT, 01935 $ INFO ) 01936 330 CONTINUE 01937 * 01938 ELSE IF( LSAMEN( 3, C3, 'ZEV' ) ) THEN 01939 * 01940 * -------------------------------------------- 01941 * ZEV: Nonsymmetric Eigenvalue Problem Driver 01942 * ZGEEV (eigenvalues and eigenvectors) 01943 * -------------------------------------------- 01944 * 01945 MAXTYP = 21 01946 NTYPES = MIN( MAXTYP, NTYPES ) 01947 IF( NTYPES.LE.0 ) THEN 01948 WRITE( NOUT, FMT = 9990 )C3 01949 ELSE 01950 IF( TSTERR ) 01951 $ CALL ZERRED( C3, NOUT ) 01952 CALL ALAREQ( C3, NTYPES, DOTYPE, MAXTYP, NIN, NOUT ) 01953 CALL ZDRVEV( NN, NVAL, NTYPES, DOTYPE, ISEED, THRESH, NOUT, 01954 $ A( 1, 1 ), NMAX, A( 1, 2 ), DC( 1, 1 ), 01955 $ DC( 1, 2 ), A( 1, 3 ), NMAX, A( 1, 4 ), NMAX, 01956 $ A( 1, 5 ), NMAX, RESULT, WORK, LWORK, RWORK, 01957 $ IWORK, INFO ) 01958 IF( INFO.NE.0 ) 01959 $ WRITE( NOUT, FMT = 9980 )'ZGEEV', INFO 01960 END IF 01961 WRITE( NOUT, FMT = 9973 ) 01962 GO TO 10 01963 * 01964 ELSE IF( LSAMEN( 3, C3, 'ZES' ) ) THEN 01965 * 01966 * -------------------------------------------- 01967 * ZES: Nonsymmetric Eigenvalue Problem Driver 01968 * ZGEES (Schur form) 01969 * -------------------------------------------- 01970 * 01971 MAXTYP = 21 01972 NTYPES = MIN( MAXTYP, NTYPES ) 01973 IF( NTYPES.LE.0 ) THEN 01974 WRITE( NOUT, FMT = 9990 )C3 01975 ELSE 01976 IF( TSTERR ) 01977 $ CALL ZERRED( C3, NOUT ) 01978 CALL ALAREQ( C3, NTYPES, DOTYPE, MAXTYP, NIN, NOUT ) 01979 CALL ZDRVES( NN, NVAL, NTYPES, DOTYPE, ISEED, THRESH, NOUT, 01980 $ A( 1, 1 ), NMAX, A( 1, 2 ), A( 1, 3 ), 01981 $ DC( 1, 1 ), DC( 1, 2 ), A( 1, 4 ), NMAX, 01982 $ RESULT, WORK, LWORK, RWORK, IWORK, LOGWRK, 01983 $ INFO ) 01984 IF( INFO.NE.0 ) 01985 $ WRITE( NOUT, FMT = 9980 )'ZGEES', INFO 01986 END IF 01987 WRITE( NOUT, FMT = 9973 ) 01988 GO TO 10 01989 * 01990 ELSE IF( LSAMEN( 3, C3, 'ZVX' ) ) THEN 01991 * 01992 * -------------------------------------------------------------- 01993 * ZVX: Nonsymmetric Eigenvalue Problem Expert Driver 01994 * ZGEEVX (eigenvalues, eigenvectors and condition numbers) 01995 * -------------------------------------------------------------- 01996 * 01997 MAXTYP = 21 01998 NTYPES = MIN( MAXTYP, NTYPES ) 01999 IF( NTYPES.LT.0 ) THEN 02000 WRITE( NOUT, FMT = 9990 )C3 02001 ELSE 02002 IF( TSTERR ) 02003 $ CALL ZERRED( C3, NOUT ) 02004 CALL ALAREQ( C3, NTYPES, DOTYPE, MAXTYP, NIN, NOUT ) 02005 CALL ZDRVVX( NN, NVAL, NTYPES, DOTYPE, ISEED, THRESH, NIN, 02006 $ NOUT, A( 1, 1 ), NMAX, A( 1, 2 ), DC( 1, 1 ), 02007 $ DC( 1, 2 ), A( 1, 3 ), NMAX, A( 1, 4 ), NMAX, 02008 $ A( 1, 5 ), NMAX, DR( 1, 1 ), DR( 1, 2 ), 02009 $ DR( 1, 3 ), DR( 1, 4 ), DR( 1, 5 ), DR( 1, 6 ), 02010 $ DR( 1, 7 ), DR( 1, 8 ), RESULT, WORK, LWORK, 02011 $ RWORK, INFO ) 02012 IF( INFO.NE.0 ) 02013 $ WRITE( NOUT, FMT = 9980 )'ZGEEVX', INFO 02014 END IF 02015 WRITE( NOUT, FMT = 9973 ) 02016 GO TO 10 02017 * 02018 ELSE IF( LSAMEN( 3, C3, 'ZSX' ) ) THEN 02019 * 02020 * --------------------------------------------------- 02021 * ZSX: Nonsymmetric Eigenvalue Problem Expert Driver 02022 * ZGEESX (Schur form and condition numbers) 02023 * --------------------------------------------------- 02024 * 02025 MAXTYP = 21 02026 NTYPES = MIN( MAXTYP, NTYPES ) 02027 IF( NTYPES.LT.0 ) THEN 02028 WRITE( NOUT, FMT = 9990 )C3 02029 ELSE 02030 IF( TSTERR ) 02031 $ CALL ZERRED( C3, NOUT ) 02032 CALL ALAREQ( C3, NTYPES, DOTYPE, MAXTYP, NIN, NOUT ) 02033 CALL ZDRVSX( NN, NVAL, NTYPES, DOTYPE, ISEED, THRESH, NIN, 02034 $ NOUT, A( 1, 1 ), NMAX, A( 1, 2 ), A( 1, 3 ), 02035 $ DC( 1, 1 ), DC( 1, 2 ), DC( 1, 3 ), A( 1, 4 ), 02036 $ NMAX, A( 1, 5 ), RESULT, WORK, LWORK, RWORK, 02037 $ LOGWRK, INFO ) 02038 IF( INFO.NE.0 ) 02039 $ WRITE( NOUT, FMT = 9980 )'ZGEESX', INFO 02040 END IF 02041 WRITE( NOUT, FMT = 9973 ) 02042 GO TO 10 02043 * 02044 ELSE IF( LSAMEN( 3, C3, 'ZGG' ) ) THEN 02045 * 02046 * ------------------------------------------------- 02047 * ZGG: Generalized Nonsymmetric Eigenvalue Problem 02048 * ------------------------------------------------- 02049 * Vary the parameters 02050 * NB = block size 02051 * NBMIN = minimum block size 02052 * NS = number of shifts 02053 * MAXB = minimum submatrix size 02054 * NBCOL = minimum column dimension for blocks 02055 * 02056 MAXTYP = 26 02057 NTYPES = MIN( MAXTYP, NTYPES ) 02058 CALL ALAREQ( C3, NTYPES, DOTYPE, MAXTYP, NIN, NOUT ) 02059 IF( TSTCHK .AND. TSTERR ) 02060 $ CALL ZERRGG( C3, NOUT ) 02061 DO 350 I = 1, NPARMS 02062 CALL XLAENV( 1, NBVAL( I ) ) 02063 CALL XLAENV( 2, NBMIN( I ) ) 02064 CALL XLAENV( 4, NSVAL( I ) ) 02065 CALL XLAENV( 8, MXBVAL( I ) ) 02066 CALL XLAENV( 5, NBCOL( I ) ) 02067 * 02068 IF( NEWSD.EQ.0 ) THEN 02069 DO 340 K = 1, 4 02070 ISEED( K ) = IOLDSD( K ) 02071 340 CONTINUE 02072 END IF 02073 WRITE( NOUT, FMT = 9996 )C3, NBVAL( I ), NBMIN( I ), 02074 $ NSVAL( I ), MXBVAL( I ), NBCOL( I ) 02075 TSTDIF = .FALSE. 02076 THRSHN = 10.D0 02077 IF( TSTCHK ) THEN 02078 CALL ZCHKGG( NN, NVAL, MAXTYP, DOTYPE, ISEED, THRESH, 02079 $ TSTDIF, THRSHN, NOUT, A( 1, 1 ), NMAX, 02080 $ A( 1, 2 ), A( 1, 3 ), A( 1, 4 ), A( 1, 5 ), 02081 $ A( 1, 6 ), A( 1, 7 ), A( 1, 8 ), A( 1, 9 ), 02082 $ NMAX, A( 1, 10 ), A( 1, 11 ), A( 1, 12 ), 02083 $ DC( 1, 1 ), DC( 1, 2 ), DC( 1, 3 ), 02084 $ DC( 1, 4 ), A( 1, 13 ), A( 1, 14 ), WORK, 02085 $ LWORK, RWORK, LOGWRK, RESULT, INFO ) 02086 IF( INFO.NE.0 ) 02087 $ WRITE( NOUT, FMT = 9980 )'ZCHKGG', INFO 02088 END IF 02089 CALL XLAENV( 1, 1 ) 02090 IF( TSTDRV ) THEN 02091 CALL ZDRVGG( NN, NVAL, MAXTYP, DOTYPE, ISEED, THRESH, 02092 $ THRSHN, NOUT, A( 1, 1 ), NMAX, A( 1, 2 ), 02093 $ A( 1, 3 ), A( 1, 4 ), A( 1, 5 ), A( 1, 6 ), 02094 $ A( 1, 7 ), NMAX, A( 1, 8 ), DC( 1, 1 ), 02095 $ DC( 1, 2 ), DC( 1, 3 ), DC( 1, 4 ), 02096 $ A( 1, 8 ), A( 1, 9 ), WORK, LWORK, RWORK, 02097 $ RESULT, INFO ) 02098 IF( INFO.NE.0 ) 02099 $ WRITE( NOUT, FMT = 9980 )'ZDRVGG', INFO 02100 END IF 02101 350 CONTINUE 02102 * 02103 ELSE IF( LSAMEN( 3, C3, 'ZGS' ) ) THEN 02104 * 02105 * ------------------------------------------------- 02106 * ZGS: Generalized Nonsymmetric Eigenvalue Problem 02107 * ZGGES (Schur form) 02108 * ------------------------------------------------- 02109 * 02110 MAXTYP = 26 02111 NTYPES = MIN( MAXTYP, NTYPES ) 02112 IF( NTYPES.LE.0 ) THEN 02113 WRITE( NOUT, FMT = 9990 )C3 02114 ELSE 02115 IF( TSTERR ) 02116 $ CALL ZERRGG( C3, NOUT ) 02117 CALL ALAREQ( C3, NTYPES, DOTYPE, MAXTYP, NIN, NOUT ) 02118 CALL ZDRGES( NN, NVAL, MAXTYP, DOTYPE, ISEED, THRESH, NOUT, 02119 $ A( 1, 1 ), NMAX, A( 1, 2 ), A( 1, 3 ), 02120 $ A( 1, 4 ), A( 1, 7 ), NMAX, A( 1, 8 ), 02121 $ DC( 1, 1 ), DC( 1, 2 ), WORK, LWORK, RWORK, 02122 $ RESULT, LOGWRK, INFO ) 02123 * 02124 IF( INFO.NE.0 ) 02125 $ WRITE( NOUT, FMT = 9980 )'ZDRGES', INFO 02126 END IF 02127 WRITE( NOUT, FMT = 9973 ) 02128 GO TO 10 02129 * 02130 ELSE IF( ZGX ) THEN 02131 * 02132 * ------------------------------------------------- 02133 * ZGX Generalized Nonsymmetric Eigenvalue Problem 02134 * ZGGESX (Schur form and condition numbers) 02135 * ------------------------------------------------- 02136 * 02137 MAXTYP = 5 02138 NTYPES = MAXTYP 02139 IF( NN.LT.0 ) THEN 02140 WRITE( NOUT, FMT = 9990 )C3 02141 ELSE 02142 IF( TSTERR ) 02143 $ CALL ZERRGG( C3, NOUT ) 02144 CALL ALAREQ( C3, NTYPES, DOTYPE, MAXTYP, NIN, NOUT ) 02145 CALL XLAENV( 5, 2 ) 02146 CALL ZDRGSX( NN, NCMAX, THRESH, NIN, NOUT, A( 1, 1 ), NMAX, 02147 $ A( 1, 2 ), A( 1, 3 ), A( 1, 4 ), A( 1, 5 ), 02148 $ A( 1, 6 ), DC( 1, 1 ), DC( 1, 2 ), C, 02149 $ NCMAX*NCMAX, S, WORK, LWORK, RWORK, IWORK, 02150 $ LIWORK, LOGWRK, INFO ) 02151 IF( INFO.NE.0 ) 02152 $ WRITE( NOUT, FMT = 9980 )'ZDRGSX', INFO 02153 END IF 02154 WRITE( NOUT, FMT = 9973 ) 02155 GO TO 10 02156 * 02157 ELSE IF( LSAMEN( 3, C3, 'ZGV' ) ) THEN 02158 * 02159 * ------------------------------------------------- 02160 * ZGV: Generalized Nonsymmetric Eigenvalue Problem 02161 * ZGGEV (Eigenvalue/vector form) 02162 * ------------------------------------------------- 02163 * 02164 MAXTYP = 26 02165 NTYPES = MIN( MAXTYP, NTYPES ) 02166 IF( NTYPES.LE.0 ) THEN 02167 WRITE( NOUT, FMT = 9990 )C3 02168 ELSE 02169 IF( TSTERR ) 02170 $ CALL ZERRGG( C3, NOUT ) 02171 CALL ALAREQ( C3, NTYPES, DOTYPE, MAXTYP, NIN, NOUT ) 02172 CALL ZDRGEV( NN, NVAL, MAXTYP, DOTYPE, ISEED, THRESH, NOUT, 02173 $ A( 1, 1 ), NMAX, A( 1, 2 ), A( 1, 3 ), 02174 $ A( 1, 4 ), A( 1, 7 ), NMAX, A( 1, 8 ), 02175 $ A( 1, 9 ), NMAX, DC( 1, 1 ), DC( 1, 2 ), 02176 $ DC( 1, 3 ), DC( 1, 4 ), WORK, LWORK, RWORK, 02177 $ RESULT, INFO ) 02178 IF( INFO.NE.0 ) 02179 $ WRITE( NOUT, FMT = 9980 )'ZDRGEV', INFO 02180 END IF 02181 WRITE( NOUT, FMT = 9973 ) 02182 GO TO 10 02183 * 02184 ELSE IF( ZXV ) THEN 02185 * 02186 * ------------------------------------------------- 02187 * ZXV: Generalized Nonsymmetric Eigenvalue Problem 02188 * ZGGEVX (eigenvalue/vector with condition numbers) 02189 * ------------------------------------------------- 02190 * 02191 MAXTYP = 2 02192 NTYPES = MAXTYP 02193 IF( NN.LT.0 ) THEN 02194 WRITE( NOUT, FMT = 9990 )C3 02195 ELSE 02196 IF( TSTERR ) 02197 $ CALL ZERRGG( C3, NOUT ) 02198 CALL ALAREQ( C3, NTYPES, DOTYPE, MAXTYP, NIN, NOUT ) 02199 CALL ZDRGVX( NN, THRESH, NIN, NOUT, A( 1, 1 ), NMAX, 02200 $ A( 1, 2 ), A( 1, 3 ), A( 1, 4 ), DC( 1, 1 ), 02201 $ DC( 1, 2 ), A( 1, 5 ), A( 1, 6 ), IWORK( 1 ), 02202 $ IWORK( 2 ), DR( 1, 1 ), DR( 1, 2 ), DR( 1, 3 ), 02203 $ DR( 1, 4 ), DR( 1, 5 ), DR( 1, 6 ), WORK, 02204 $ LWORK, RWORK, IWORK( 3 ), LIWORK-2, RESULT, 02205 $ LOGWRK, INFO ) 02206 * 02207 IF( INFO.NE.0 ) 02208 $ WRITE( NOUT, FMT = 9980 )'ZDRGVX', INFO 02209 END IF 02210 WRITE( NOUT, FMT = 9973 ) 02211 GO TO 10 02212 * 02213 ELSE IF( LSAMEN( 3, C3, 'ZHB' ) ) THEN 02214 * 02215 * ------------------------------ 02216 * ZHB: Hermitian Band Reduction 02217 * ------------------------------ 02218 * 02219 MAXTYP = 15 02220 NTYPES = MIN( MAXTYP, NTYPES ) 02221 CALL ALAREQ( C3, NTYPES, DOTYPE, MAXTYP, NIN, NOUT ) 02222 IF( TSTERR ) 02223 $ CALL ZERRST( 'ZHB', NOUT ) 02224 CALL ZCHKHB( NN, NVAL, NK, KVAL, MAXTYP, DOTYPE, ISEED, THRESH, 02225 $ NOUT, A( 1, 1 ), NMAX, DR( 1, 1 ), DR( 1, 2 ), 02226 $ A( 1, 2 ), NMAX, WORK, LWORK, RWORK, RESULT, 02227 $ INFO ) 02228 IF( INFO.NE.0 ) 02229 $ WRITE( NOUT, FMT = 9980 )'ZCHKHB', INFO 02230 * 02231 ELSE IF( LSAMEN( 3, C3, 'ZBB' ) ) THEN 02232 * 02233 * ------------------------------ 02234 * ZBB: General Band Reduction 02235 * ------------------------------ 02236 * 02237 MAXTYP = 15 02238 NTYPES = MIN( MAXTYP, NTYPES ) 02239 CALL ALAREQ( C3, NTYPES, DOTYPE, MAXTYP, NIN, NOUT ) 02240 DO 370 I = 1, NPARMS 02241 NRHS = NSVAL( I ) 02242 * 02243 IF( NEWSD.EQ.0 ) THEN 02244 DO 360 K = 1, 4 02245 ISEED( K ) = IOLDSD( K ) 02246 360 CONTINUE 02247 END IF 02248 WRITE( NOUT, FMT = 9966 )C3, NRHS 02249 CALL ZCHKBB( NN, MVAL, NVAL, NK, KVAL, MAXTYP, DOTYPE, NRHS, 02250 $ ISEED, THRESH, NOUT, A( 1, 1 ), NMAX, 02251 $ A( 1, 2 ), 2*NMAX, DR( 1, 1 ), DR( 1, 2 ), 02252 $ A( 1, 4 ), NMAX, A( 1, 5 ), NMAX, A( 1, 6 ), 02253 $ NMAX, A( 1, 7 ), WORK, LWORK, RWORK, RESULT, 02254 $ INFO ) 02255 IF( INFO.NE.0 ) 02256 $ WRITE( NOUT, FMT = 9980 )'ZCHKBB', INFO 02257 370 CONTINUE 02258 * 02259 ELSE IF( LSAMEN( 3, C3, 'GLM' ) ) THEN 02260 * 02261 * ----------------------------------------- 02262 * GLM: Generalized Linear Regression Model 02263 * ----------------------------------------- 02264 * 02265 CALL XLAENV( 1, 1 ) 02266 IF( TSTERR ) 02267 $ CALL ZERRGG( 'GLM', NOUT ) 02268 CALL ZCKGLM( NN, NVAL, MVAL, PVAL, NTYPES, ISEED, THRESH, NMAX, 02269 $ A( 1, 1 ), A( 1, 2 ), B( 1, 1 ), B( 1, 2 ), X, 02270 $ WORK, DR( 1, 1 ), NIN, NOUT, INFO ) 02271 IF( INFO.NE.0 ) 02272 $ WRITE( NOUT, FMT = 9980 )'ZCKGLM', INFO 02273 * 02274 ELSE IF( LSAMEN( 3, C3, 'GQR' ) ) THEN 02275 * 02276 * ------------------------------------------ 02277 * GQR: Generalized QR and RQ factorizations 02278 * ------------------------------------------ 02279 * 02280 CALL XLAENV( 1, 1 ) 02281 IF( TSTERR ) 02282 $ CALL ZERRGG( 'GQR', NOUT ) 02283 CALL ZCKGQR( NN, MVAL, NN, PVAL, NN, NVAL, NTYPES, ISEED, 02284 $ THRESH, NMAX, A( 1, 1 ), A( 1, 2 ), A( 1, 3 ), 02285 $ A( 1, 4 ), TAUA, B( 1, 1 ), B( 1, 2 ), B( 1, 3 ), 02286 $ B( 1, 4 ), B( 1, 5 ), TAUB, WORK, DR( 1, 1 ), NIN, 02287 $ NOUT, INFO ) 02288 IF( INFO.NE.0 ) 02289 $ WRITE( NOUT, FMT = 9980 )'ZCKGQR', INFO 02290 * 02291 ELSE IF( LSAMEN( 3, C3, 'GSV' ) ) THEN 02292 * 02293 * ---------------------------------------------- 02294 * GSV: Generalized Singular Value Decomposition 02295 * ---------------------------------------------- 02296 * 02297 IF( TSTERR ) 02298 $ CALL ZERRGG( 'GSV', NOUT ) 02299 CALL ZCKGSV( NN, MVAL, PVAL, NVAL, NTYPES, ISEED, THRESH, NMAX, 02300 $ A( 1, 1 ), A( 1, 2 ), B( 1, 1 ), B( 1, 2 ), 02301 $ A( 1, 3 ), B( 1, 3 ), A( 1, 4 ), ALPHA, BETA, 02302 $ B( 1, 4 ), IWORK, WORK, DR( 1, 1 ), NIN, NOUT, 02303 $ INFO ) 02304 IF( INFO.NE.0 ) 02305 $ WRITE( NOUT, FMT = 9980 )'ZCKGSV', INFO 02306 * 02307 ELSE IF( LSAMEN( 3, C3, 'CSD' ) ) THEN 02308 * 02309 * ---------------------------------------------- 02310 * CSD: CS Decomposition 02311 * ---------------------------------------------- 02312 * 02313 IF( TSTERR ) 02314 $ CALL ZERRGG( 'CSD', NOUT ) 02315 CALL ZCKCSD( NN, MVAL, PVAL, NVAL, NTYPES, ISEED, THRESH, NMAX, 02316 $ A( 1, 1 ), A( 1, 2 ), A( 1, 3 ), A( 1, 4 ), 02317 $ A( 1, 5 ), A( 1, 6 ), A( 1, 7 ), IWORK, WORK, 02318 $ DR( 1, 1 ), NIN, NOUT, INFO ) 02319 IF( INFO.NE.0 ) 02320 $ WRITE( NOUT, FMT = 9980 )'ZCKCSD', INFO 02321 * 02322 ELSE IF( LSAMEN( 3, C3, 'LSE' ) ) THEN 02323 * 02324 * -------------------------------------- 02325 * LSE: Constrained Linear Least Squares 02326 * -------------------------------------- 02327 * 02328 CALL XLAENV( 1, 1 ) 02329 IF( TSTERR ) 02330 $ CALL ZERRGG( 'LSE', NOUT ) 02331 CALL ZCKLSE( NN, MVAL, PVAL, NVAL, NTYPES, ISEED, THRESH, NMAX, 02332 $ A( 1, 1 ), A( 1, 2 ), B( 1, 1 ), B( 1, 2 ), X, 02333 $ WORK, DR( 1, 1 ), NIN, NOUT, INFO ) 02334 IF( INFO.NE.0 ) 02335 $ WRITE( NOUT, FMT = 9980 )'ZCKLSE', INFO 02336 ELSE 02337 WRITE( NOUT, FMT = * ) 02338 WRITE( NOUT, FMT = * ) 02339 WRITE( NOUT, FMT = 9992 )C3 02340 END IF 02341 IF( .NOT.( ZGX .OR. ZXV ) ) 02342 $ GO TO 190 02343 380 CONTINUE 02344 WRITE( NOUT, FMT = 9994 ) 02345 S2 = DSECND( ) 02346 WRITE( NOUT, FMT = 9993 )S2 - S1 02347 * 02348 9999 FORMAT( / ' Execution not attempted due to input errors' ) 02349 9998 FORMAT( / / 1X, A3, ': NB =', I4, ', NBMIN =', I4, ', NX =', I4, 02350 $ ', NS =', I4, ', MAXB =', I4 ) 02351 9997 FORMAT( / / 1X, A3, ': NB =', I4, ', NBMIN =', I4, ', NX =', I4 ) 02352 9996 FORMAT( / / 1X, A3, ': NB =', I4, ', NBMIN =', I4, ', NS =', I4, 02353 $ ', MAXB =', I4, ', NBCOL =', I4 ) 02354 9995 FORMAT( / / 1X, A3, ': NB =', I4, ', NBMIN =', I4, ', NX =', I4, 02355 $ ', NRHS =', I4 ) 02356 9994 FORMAT( / / ' End of tests' ) 02357 9993 FORMAT( ' Total time used = ', F12.2, ' seconds', / ) 02358 9992 FORMAT( 1X, A3, ': Unrecognized path name' ) 02359 9991 FORMAT( / / ' *** Invalid integer value in column ', I2, 02360 $ ' of input', ' line:', / A79 ) 02361 9990 FORMAT( / / 1X, A3, ' routines were not tested' ) 02362 9989 FORMAT( ' Invalid input value: ', A, '=', I6, '; must be >=', 02363 $ I6 ) 02364 9988 FORMAT( ' Invalid input value: ', A, '=', I6, '; must be <=', 02365 $ I6 ) 02366 9987 FORMAT( ' Tests of the Nonsymmetric Eigenvalue Problem routines' ) 02367 9986 FORMAT( ' Tests of the Hermitian Eigenvalue Problem routines' ) 02368 9985 FORMAT( ' Tests of the Singular Value Decomposition routines' ) 02369 9984 FORMAT( / ' The following parameter values will be used:' ) 02370 9983 FORMAT( 4X, A, 10I6, / 10X, 10I6 ) 02371 9982 FORMAT( / ' Routines pass computational tests if test ratio is ', 02372 $ 'less than', F8.2, / ) 02373 9981 FORMAT( ' Relative machine ', A, ' is taken to be', D16.6 ) 02374 9980 FORMAT( ' *** Error code from ', A, ' = ', I4 ) 02375 9979 FORMAT( / ' Tests of the Nonsymmetric Eigenvalue Problem Driver', 02376 $ / ' ZGEEV (eigenvalues and eigevectors)' ) 02377 9978 FORMAT( / ' Tests of the Nonsymmetric Eigenvalue Problem Driver', 02378 $ / ' ZGEES (Schur form)' ) 02379 9977 FORMAT( / ' Tests of the Nonsymmetric Eigenvalue Problem Expert', 02380 $ ' Driver', / ' ZGEEVX (eigenvalues, eigenvectors and', 02381 $ ' condition numbers)' ) 02382 9976 FORMAT( / ' Tests of the Nonsymmetric Eigenvalue Problem Expert', 02383 $ ' Driver', / ' ZGEESX (Schur form and condition', 02384 $ ' numbers)' ) 02385 9975 FORMAT( / ' Tests of the Generalized Nonsymmetric Eigenvalue ', 02386 $ 'Problem routines' ) 02387 9974 FORMAT( ' Tests of ZHBTRD', / ' (reduction of a Hermitian band ', 02388 $ 'matrix to real tridiagonal form)' ) 02389 9973 FORMAT( / 1X, 71( '-' ) ) 02390 9972 FORMAT( / ' LAPACK VERSION ', I1, '.', I1, '.', I1 ) 02391 9971 FORMAT( / ' Tests of the Generalized Linear Regression Model ', 02392 $ 'routines' ) 02393 9970 FORMAT( / ' Tests of the Generalized QR and RQ routines' ) 02394 9969 FORMAT( / ' Tests of the Generalized Singular Value', 02395 $ ' Decomposition routines' ) 02396 9968 FORMAT( / ' Tests of the Linear Least Squares routines' ) 02397 9967 FORMAT( ' Tests of ZGBBRD', / ' (reduction of a general band ', 02398 $ 'matrix to real bidiagonal form)' ) 02399 9966 FORMAT( / / 1X, A3, ': NRHS =', I4 ) 02400 9965 FORMAT( / ' Tests of the Generalized Nonsymmetric Eigenvalue ', 02401 $ 'Problem Expert Driver ZGGESX' ) 02402 9964 FORMAT( / ' Tests of the Generalized Nonsymmetric Eigenvalue ', 02403 $ 'Problem Driver ZGGES' ) 02404 9963 FORMAT( / ' Tests of the Generalized Nonsymmetric Eigenvalue ', 02405 $ 'Problem Driver ZGGEV' ) 02406 9962 FORMAT( / ' Tests of the Generalized Nonsymmetric Eigenvalue ', 02407 $ 'Problem Expert Driver ZGGEVX' ) 02408 9961 FORMAT( / / 1X, A3, ': NB =', I4, ', NBMIN =', I4, ', NX =', I4, 02409 $ ', INMIN=', I4, 02410 $ ', INWIN =', I4, ', INIBL =', I4, ', ISHFTS =', I4, 02411 $ ', IACC22 =', I4) 02412 9960 FORMAT( / ' Tests of the CS Decomposition routines' ) 02413 * 02414 * End of ZCHKEE 02415 * 02416 END
1.7.3 
