 |
LAPACK 3.3.1
Linear Algebra PACKage
|
|
LAPACK 3.3.1
Linear Algebra PACKage
|
00001 PROGRAM DCHKAA 00002 * 00003 * -- LAPACK test routine (version 3.1.1) -- 00004 * Univ. of Tennessee, Univ. of California Berkeley and NAG Ltd.. 00005 * January 2007 00006 * 00007 * Purpose 00008 * ======= 00009 * 00010 * DCHKAA is the main test program for the DOUBLE PRECISION LAPACK 00011 * linear equation routines 00012 * 00013 * The program must be driven by a short data file. The first 14 records 00014 * specify problem dimensions and program options using list-directed 00015 * input. The remaining lines specify the LAPACK test paths and the 00016 * number of matrix types to use in testing. An annotated example of a 00017 * data file can be obtained by deleting the first 3 characters from the 00018 * following 36 lines: 00019 * Data file for testing DOUBLE PRECISION LAPACK linear eqn. routines 00020 * 7 Number of values of M 00021 * 0 1 2 3 5 10 16 Values of M (row dimension) 00022 * 7 Number of values of N 00023 * 0 1 2 3 5 10 16 Values of N (column dimension) 00024 * 1 Number of values of NRHS 00025 * 2 Values of NRHS (number of right hand sides) 00026 * 5 Number of values of NB 00027 * 1 3 3 3 20 Values of NB (the blocksize) 00028 * 1 0 5 9 1 Values of NX (crossover point) 00029 * 3 Number of values of RANK 00030 * 30 50 90 Values of rank (as a % of N) 00031 * 20.0 Threshold value of test ratio 00032 * T Put T to test the LAPACK routines 00033 * T Put T to test the driver routines 00034 * T Put T to test the error exits 00035 * DGE 11 List types on next line if 0 < NTYPES < 11 00036 * DGB 8 List types on next line if 0 < NTYPES < 8 00037 * DGT 12 List types on next line if 0 < NTYPES < 12 00038 * DPO 9 List types on next line if 0 < NTYPES < 9 00039 * DPS 9 List types on next line if 0 < NTYPES < 9 00040 * DPP 9 List types on next line if 0 < NTYPES < 9 00041 * DPB 8 List types on next line if 0 < NTYPES < 8 00042 * DPT 12 List types on next line if 0 < NTYPES < 12 00043 * DSY 10 List types on next line if 0 < NTYPES < 10 00044 * DSP 10 List types on next line if 0 < NTYPES < 10 00045 * DTR 18 List types on next line if 0 < NTYPES < 18 00046 * DTP 18 List types on next line if 0 < NTYPES < 18 00047 * DTB 17 List types on next line if 0 < NTYPES < 17 00048 * DQR 8 List types on next line if 0 < NTYPES < 8 00049 * DRQ 8 List types on next line if 0 < NTYPES < 8 00050 * DLQ 8 List types on next line if 0 < NTYPES < 8 00051 * DQL 8 List types on next line if 0 < NTYPES < 8 00052 * DQP 6 List types on next line if 0 < NTYPES < 6 00053 * DTZ 3 List types on next line if 0 < NTYPES < 3 00054 * DLS 6 List types on next line if 0 < NTYPES < 6 00055 * DEQ 00056 * 00057 * Internal Parameters 00058 * =================== 00059 * 00060 * NMAX INTEGER 00061 * The maximum allowable value for N 00062 * 00063 * MAXIN INTEGER 00064 * The number of different values that can be used for each of 00065 * M, N, NRHS, NB, and NX 00066 * 00067 * MAXRHS INTEGER 00068 * The maximum number of right hand sides 00069 * 00070 * NIN INTEGER 00071 * The unit number for input 00072 * 00073 * NOUT INTEGER 00074 * The unit number for output 00075 * 00076 * ===================================================================== 00077 * 00078 * .. Parameters .. 00079 INTEGER NMAX 00080 PARAMETER ( NMAX = 132 ) 00081 INTEGER MAXIN 00082 PARAMETER ( MAXIN = 12 ) 00083 INTEGER MAXRHS 00084 PARAMETER ( MAXRHS = 16 ) 00085 INTEGER MATMAX 00086 PARAMETER ( MATMAX = 30 ) 00087 INTEGER NIN, NOUT 00088 PARAMETER ( NIN = 5, NOUT = 6 ) 00089 INTEGER KDMAX 00090 PARAMETER ( KDMAX = NMAX+( NMAX+1 ) / 4 ) 00091 * .. 00092 * .. Local Scalars .. 00093 LOGICAL FATAL, TSTCHK, TSTDRV, TSTERR 00094 CHARACTER C1 00095 CHARACTER*2 C2 00096 CHARACTER*3 PATH 00097 CHARACTER*10 INTSTR 00098 CHARACTER*72 ALINE 00099 INTEGER I, IC, J, K, LA, LAFAC, LDA, NB, NM, NMATS, NN, 00100 $ NNB, NNB2, NNS, NRHS, NTYPES, NRANK, 00101 $ VERS_MAJOR, VERS_MINOR, VERS_PATCH 00102 DOUBLE PRECISION EPS, S1, S2, THREQ, THRESH 00103 * .. 00104 * .. Local Arrays .. 00105 LOGICAL DOTYPE( MATMAX ) 00106 INTEGER IWORK( 25*NMAX ), MVAL( MAXIN ), 00107 $ NBVAL( MAXIN ), NBVAL2( MAXIN ), 00108 $ NSVAL( MAXIN ), NVAL( MAXIN ), NXVAL( MAXIN ), 00109 $ RANKVAL( MAXIN ), PIV( NMAX ) 00110 DOUBLE PRECISION A( ( KDMAX+1 )*NMAX, 7 ), B( NMAX*MAXRHS, 4 ), 00111 $ RWORK( 5*NMAX+2*MAXRHS ), S( 2*NMAX ), 00112 $ WORK( NMAX, NMAX+MAXRHS+30 ) 00113 * .. 00114 * .. External Functions .. 00115 LOGICAL LSAME, LSAMEN 00116 DOUBLE PRECISION DLAMCH, DSECND 00117 EXTERNAL LSAME, LSAMEN, DLAMCH, DSECND 00118 * .. 00119 * .. External Subroutines .. 00120 EXTERNAL ALAREQ, DCHKEQ, DCHKGB, DCHKGE, DCHKGT, DCHKLQ, 00121 $ DCHKPB, DCHKPO, DCHKPS, DCHKPP, DCHKPT, DCHKQ3, 00122 $ DCHKQL, DCHKQP, DCHKQR, DCHKRQ, DCHKSP, DCHKSY, 00123 $ DCHKTB, DCHKTP, DCHKTR, DCHKTZ, DDRVGB, DDRVGE, 00124 $ DDRVGT, DDRVLS, DDRVPB, DDRVPO, DDRVPP, DDRVPT, 00125 $ DDRVSP, DDRVSY, ILAVER 00126 * .. 00127 * .. Scalars in Common .. 00128 LOGICAL LERR, OK 00129 CHARACTER*32 SRNAMT 00130 INTEGER INFOT, NUNIT 00131 * .. 00132 * .. Arrays in Common .. 00133 INTEGER IPARMS( 100 ) 00134 * .. 00135 * .. Common blocks .. 00136 COMMON / INFOC / INFOT, NUNIT, OK, LERR 00137 COMMON / SRNAMC / SRNAMT 00138 COMMON / CLAENV / IPARMS 00139 * .. 00140 * .. Data statements .. 00141 DATA THREQ / 2.0D0 / , INTSTR / '0123456789' / 00142 * .. 00143 * .. Executable Statements .. 00144 * 00145 S1 = DSECND( ) 00146 LDA = NMAX 00147 FATAL = .FALSE. 00148 * 00149 * Read a dummy line. 00150 * 00151 READ( NIN, FMT = * ) 00152 * 00153 * Report values of parameters. 00154 * 00155 CALL ILAVER( VERS_MAJOR, VERS_MINOR, VERS_PATCH ) 00156 WRITE( NOUT, FMT = 9994 ) VERS_MAJOR, VERS_MINOR, VERS_PATCH 00157 * 00158 * Read the values of M 00159 * 00160 READ( NIN, FMT = * )NM 00161 IF( NM.LT.1 ) THEN 00162 WRITE( NOUT, FMT = 9996 )' NM ', NM, 1 00163 NM = 0 00164 FATAL = .TRUE. 00165 ELSE IF( NM.GT.MAXIN ) THEN 00166 WRITE( NOUT, FMT = 9995 )' NM ', NM, MAXIN 00167 NM = 0 00168 FATAL = .TRUE. 00169 END IF 00170 READ( NIN, FMT = * )( MVAL( I ), I = 1, NM ) 00171 DO 10 I = 1, NM 00172 IF( MVAL( I ).LT.0 ) THEN 00173 WRITE( NOUT, FMT = 9996 )' M ', MVAL( I ), 0 00174 FATAL = .TRUE. 00175 ELSE IF( MVAL( I ).GT.NMAX ) THEN 00176 WRITE( NOUT, FMT = 9995 )' M ', MVAL( I ), NMAX 00177 FATAL = .TRUE. 00178 END IF 00179 10 CONTINUE 00180 IF( NM.GT.0 ) 00181 $ WRITE( NOUT, FMT = 9993 )'M ', ( MVAL( I ), I = 1, NM ) 00182 * 00183 * Read the values of N 00184 * 00185 READ( NIN, FMT = * )NN 00186 IF( NN.LT.1 ) THEN 00187 WRITE( NOUT, FMT = 9996 )' NN ', NN, 1 00188 NN = 0 00189 FATAL = .TRUE. 00190 ELSE IF( NN.GT.MAXIN ) THEN 00191 WRITE( NOUT, FMT = 9995 )' NN ', NN, MAXIN 00192 NN = 0 00193 FATAL = .TRUE. 00194 END IF 00195 READ( NIN, FMT = * )( NVAL( I ), I = 1, NN ) 00196 DO 20 I = 1, NN 00197 IF( NVAL( I ).LT.0 ) THEN 00198 WRITE( NOUT, FMT = 9996 )' N ', NVAL( I ), 0 00199 FATAL = .TRUE. 00200 ELSE IF( NVAL( I ).GT.NMAX ) THEN 00201 WRITE( NOUT, FMT = 9995 )' N ', NVAL( I ), NMAX 00202 FATAL = .TRUE. 00203 END IF 00204 20 CONTINUE 00205 IF( NN.GT.0 ) 00206 $ WRITE( NOUT, FMT = 9993 )'N ', ( NVAL( I ), I = 1, NN ) 00207 * 00208 * Read the values of NRHS 00209 * 00210 READ( NIN, FMT = * )NNS 00211 IF( NNS.LT.1 ) THEN 00212 WRITE( NOUT, FMT = 9996 )' NNS', NNS, 1 00213 NNS = 0 00214 FATAL = .TRUE. 00215 ELSE IF( NNS.GT.MAXIN ) THEN 00216 WRITE( NOUT, FMT = 9995 )' NNS', NNS, MAXIN 00217 NNS = 0 00218 FATAL = .TRUE. 00219 END IF 00220 READ( NIN, FMT = * )( NSVAL( I ), I = 1, NNS ) 00221 DO 30 I = 1, NNS 00222 IF( NSVAL( I ).LT.0 ) THEN 00223 WRITE( NOUT, FMT = 9996 )'NRHS', NSVAL( I ), 0 00224 FATAL = .TRUE. 00225 ELSE IF( NSVAL( I ).GT.MAXRHS ) THEN 00226 WRITE( NOUT, FMT = 9995 )'NRHS', NSVAL( I ), MAXRHS 00227 FATAL = .TRUE. 00228 END IF 00229 30 CONTINUE 00230 IF( NNS.GT.0 ) 00231 $ WRITE( NOUT, FMT = 9993 )'NRHS', ( NSVAL( I ), I = 1, NNS ) 00232 * 00233 * Read the values of NB 00234 * 00235 READ( NIN, FMT = * )NNB 00236 IF( NNB.LT.1 ) THEN 00237 WRITE( NOUT, FMT = 9996 )'NNB ', NNB, 1 00238 NNB = 0 00239 FATAL = .TRUE. 00240 ELSE IF( NNB.GT.MAXIN ) THEN 00241 WRITE( NOUT, FMT = 9995 )'NNB ', NNB, MAXIN 00242 NNB = 0 00243 FATAL = .TRUE. 00244 END IF 00245 READ( NIN, FMT = * )( NBVAL( I ), I = 1, NNB ) 00246 DO 40 I = 1, NNB 00247 IF( NBVAL( I ).LT.0 ) THEN 00248 WRITE( NOUT, FMT = 9996 )' NB ', NBVAL( I ), 0 00249 FATAL = .TRUE. 00250 END IF 00251 40 CONTINUE 00252 IF( NNB.GT.0 ) 00253 $ WRITE( NOUT, FMT = 9993 )'NB ', ( NBVAL( I ), I = 1, NNB ) 00254 * 00255 * Set NBVAL2 to be the set of unique values of NB 00256 * 00257 NNB2 = 0 00258 DO 60 I = 1, NNB 00259 NB = NBVAL( I ) 00260 DO 50 J = 1, NNB2 00261 IF( NB.EQ.NBVAL2( J ) ) 00262 $ GO TO 60 00263 50 CONTINUE 00264 NNB2 = NNB2 + 1 00265 NBVAL2( NNB2 ) = NB 00266 60 CONTINUE 00267 * 00268 * Read the values of NX 00269 * 00270 READ( NIN, FMT = * )( NXVAL( I ), I = 1, NNB ) 00271 DO 70 I = 1, NNB 00272 IF( NXVAL( I ).LT.0 ) THEN 00273 WRITE( NOUT, FMT = 9996 )' NX ', NXVAL( I ), 0 00274 FATAL = .TRUE. 00275 END IF 00276 70 CONTINUE 00277 IF( NNB.GT.0 ) 00278 $ WRITE( NOUT, FMT = 9993 )'NX ', ( NXVAL( I ), I = 1, NNB ) 00279 * 00280 * Read the values of RANKVAL 00281 * 00282 READ( NIN, FMT = * )NRANK 00283 IF( NN.LT.1 ) THEN 00284 WRITE( NOUT, FMT = 9996 )' NRANK ', NRANK, 1 00285 NRANK = 0 00286 FATAL = .TRUE. 00287 ELSE IF( NN.GT.MAXIN ) THEN 00288 WRITE( NOUT, FMT = 9995 )' NRANK ', NRANK, MAXIN 00289 NRANK = 0 00290 FATAL = .TRUE. 00291 END IF 00292 READ( NIN, FMT = * )( RANKVAL( I ), I = 1, NRANK ) 00293 DO I = 1, NRANK 00294 IF( RANKVAL( I ).LT.0 ) THEN 00295 WRITE( NOUT, FMT = 9996 )' RANK ', RANKVAL( I ), 0 00296 FATAL = .TRUE. 00297 ELSE IF( RANKVAL( I ).GT.100 ) THEN 00298 WRITE( NOUT, FMT = 9995 )' RANK ', RANKVAL( I ), 100 00299 FATAL = .TRUE. 00300 END IF 00301 END DO 00302 IF( NRANK.GT.0 ) 00303 $ WRITE( NOUT, FMT = 9993 )'RANK % OF N', 00304 $ ( RANKVAL( I ), I = 1, NRANK ) 00305 * 00306 * Read the threshold value for the test ratios. 00307 * 00308 READ( NIN, FMT = * )THRESH 00309 WRITE( NOUT, FMT = 9992 )THRESH 00310 * 00311 * Read the flag that indicates whether to test the LAPACK routines. 00312 * 00313 READ( NIN, FMT = * )TSTCHK 00314 * 00315 * Read the flag that indicates whether to test the driver routines. 00316 * 00317 READ( NIN, FMT = * )TSTDRV 00318 * 00319 * Read the flag that indicates whether to test the error exits. 00320 * 00321 READ( NIN, FMT = * )TSTERR 00322 * 00323 IF( FATAL ) THEN 00324 WRITE( NOUT, FMT = 9999 ) 00325 STOP 00326 END IF 00327 * 00328 * Calculate and print the machine dependent constants. 00329 * 00330 EPS = DLAMCH( 'Underflow threshold' ) 00331 WRITE( NOUT, FMT = 9991 )'underflow', EPS 00332 EPS = DLAMCH( 'Overflow threshold' ) 00333 WRITE( NOUT, FMT = 9991 )'overflow ', EPS 00334 EPS = DLAMCH( 'Epsilon' ) 00335 WRITE( NOUT, FMT = 9991 )'precision', EPS 00336 WRITE( NOUT, FMT = * ) 00337 * 00338 80 CONTINUE 00339 * 00340 * Read a test path and the number of matrix types to use. 00341 * 00342 READ( NIN, FMT = '(A72)', END = 140 )ALINE 00343 PATH = ALINE( 1: 3 ) 00344 NMATS = MATMAX 00345 I = 3 00346 90 CONTINUE 00347 I = I + 1 00348 IF( I.GT.72 ) THEN 00349 NMATS = MATMAX 00350 GO TO 130 00351 END IF 00352 IF( ALINE( I: I ).EQ.' ' ) 00353 $ GO TO 90 00354 NMATS = 0 00355 100 CONTINUE 00356 C1 = ALINE( I: I ) 00357 DO 110 K = 1, 10 00358 IF( C1.EQ.INTSTR( K: K ) ) THEN 00359 IC = K - 1 00360 GO TO 120 00361 END IF 00362 110 CONTINUE 00363 GO TO 130 00364 120 CONTINUE 00365 NMATS = NMATS*10 + IC 00366 I = I + 1 00367 IF( I.GT.72 ) 00368 $ GO TO 130 00369 GO TO 100 00370 130 CONTINUE 00371 C1 = PATH( 1: 1 ) 00372 C2 = PATH( 2: 3 ) 00373 NRHS = NSVAL( 1 ) 00374 * 00375 * Check first character for correct precision. 00376 * 00377 IF( .NOT.LSAME( C1, 'Double precision' ) ) THEN 00378 WRITE( NOUT, FMT = 9990 )PATH 00379 * 00380 ELSE IF( NMATS.LE.0 ) THEN 00381 * 00382 * Check for a positive number of tests requested. 00383 * 00384 WRITE( NOUT, FMT = 9989 )PATH 00385 * 00386 ELSE IF( LSAMEN( 2, C2, 'GE' ) ) THEN 00387 * 00388 * GE: general matrices 00389 * 00390 NTYPES = 11 00391 CALL ALAREQ( PATH, NMATS, DOTYPE, NTYPES, NIN, NOUT ) 00392 * 00393 IF( TSTCHK ) THEN 00394 CALL DCHKGE( DOTYPE, NM, MVAL, NN, NVAL, NNB2, NBVAL2, NNS, 00395 $ NSVAL, THRESH, TSTERR, LDA, A( 1, 1 ), 00396 $ A( 1, 2 ), A( 1, 3 ), B( 1, 1 ), B( 1, 2 ), 00397 $ B( 1, 3 ), WORK, RWORK, IWORK, NOUT ) 00398 ELSE 00399 WRITE( NOUT, FMT = 9989 )PATH 00400 END IF 00401 * 00402 IF( TSTDRV ) THEN 00403 CALL DDRVGE( DOTYPE, NN, NVAL, NRHS, THRESH, TSTERR, LDA, 00404 $ A( 1, 1 ), A( 1, 2 ), A( 1, 3 ), B( 1, 1 ), 00405 $ B( 1, 2 ), B( 1, 3 ), B( 1, 4 ), S, WORK, 00406 $ RWORK, IWORK, NOUT ) 00407 ELSE 00408 WRITE( NOUT, FMT = 9988 )PATH 00409 END IF 00410 * 00411 ELSE IF( LSAMEN( 2, C2, 'GB' ) ) THEN 00412 * 00413 * GB: general banded matrices 00414 * 00415 LA = ( 2*KDMAX+1 )*NMAX 00416 LAFAC = ( 3*KDMAX+1 )*NMAX 00417 NTYPES = 8 00418 CALL ALAREQ( PATH, NMATS, DOTYPE, NTYPES, NIN, NOUT ) 00419 * 00420 IF( TSTCHK ) THEN 00421 CALL DCHKGB( DOTYPE, NM, MVAL, NN, NVAL, NNB2, NBVAL2, NNS, 00422 $ NSVAL, THRESH, TSTERR, A( 1, 1 ), LA, 00423 $ A( 1, 3 ), LAFAC, B( 1, 1 ), B( 1, 2 ), 00424 $ B( 1, 3 ), WORK, RWORK, IWORK, NOUT ) 00425 ELSE 00426 WRITE( NOUT, FMT = 9989 )PATH 00427 END IF 00428 * 00429 IF( TSTDRV ) THEN 00430 CALL DDRVGB( DOTYPE, NN, NVAL, NRHS, THRESH, TSTERR, 00431 $ A( 1, 1 ), LA, A( 1, 3 ), LAFAC, A( 1, 6 ), 00432 $ B( 1, 1 ), B( 1, 2 ), B( 1, 3 ), B( 1, 4 ), S, 00433 $ WORK, RWORK, IWORK, NOUT ) 00434 ELSE 00435 WRITE( NOUT, FMT = 9988 )PATH 00436 END IF 00437 * 00438 ELSE IF( LSAMEN( 2, C2, 'GT' ) ) THEN 00439 * 00440 * GT: general tridiagonal matrices 00441 * 00442 NTYPES = 12 00443 CALL ALAREQ( PATH, NMATS, DOTYPE, NTYPES, NIN, NOUT ) 00444 * 00445 IF( TSTCHK ) THEN 00446 CALL DCHKGT( DOTYPE, NN, NVAL, NNS, NSVAL, THRESH, TSTERR, 00447 $ A( 1, 1 ), A( 1, 2 ), B( 1, 1 ), B( 1, 2 ), 00448 $ B( 1, 3 ), WORK, RWORK, IWORK, NOUT ) 00449 ELSE 00450 WRITE( NOUT, FMT = 9989 )PATH 00451 END IF 00452 * 00453 IF( TSTDRV ) THEN 00454 CALL DDRVGT( DOTYPE, NN, NVAL, NRHS, THRESH, TSTERR, 00455 $ A( 1, 1 ), A( 1, 2 ), B( 1, 1 ), B( 1, 2 ), 00456 $ B( 1, 3 ), WORK, RWORK, IWORK, NOUT ) 00457 ELSE 00458 WRITE( NOUT, FMT = 9988 )PATH 00459 END IF 00460 * 00461 ELSE IF( LSAMEN( 2, C2, 'PO' ) ) THEN 00462 * 00463 * PO: positive definite matrices 00464 * 00465 NTYPES = 9 00466 CALL ALAREQ( PATH, NMATS, DOTYPE, NTYPES, NIN, NOUT ) 00467 * 00468 IF( TSTCHK ) THEN 00469 CALL DCHKPO( DOTYPE, NN, NVAL, NNB2, NBVAL2, NNS, NSVAL, 00470 $ THRESH, TSTERR, LDA, A( 1, 1 ), A( 1, 2 ), 00471 $ A( 1, 3 ), B( 1, 1 ), B( 1, 2 ), B( 1, 3 ), 00472 $ WORK, RWORK, IWORK, NOUT ) 00473 ELSE 00474 WRITE( NOUT, FMT = 9989 )PATH 00475 END IF 00476 * 00477 IF( TSTDRV ) THEN 00478 CALL DDRVPO( DOTYPE, NN, NVAL, NRHS, THRESH, TSTERR, LDA, 00479 $ A( 1, 1 ), A( 1, 2 ), A( 1, 3 ), B( 1, 1 ), 00480 $ B( 1, 2 ), B( 1, 3 ), B( 1, 4 ), S, WORK, 00481 $ RWORK, IWORK, NOUT ) 00482 ELSE 00483 WRITE( NOUT, FMT = 9988 )PATH 00484 END IF 00485 * 00486 ELSE IF( LSAMEN( 2, C2, 'PS' ) ) THEN 00487 * 00488 * PS: positive semi-definite matrices 00489 * 00490 NTYPES = 9 00491 * 00492 CALL ALAREQ( PATH, NMATS, DOTYPE, NTYPES, NIN, NOUT ) 00493 * 00494 IF( TSTCHK ) THEN 00495 CALL DCHKPS( DOTYPE, NN, NVAL, NNB2, NBVAL2, NRANK, 00496 $ RANKVAL, THRESH, TSTERR, LDA, A( 1, 1 ), 00497 $ A( 1, 2 ), A( 1, 3 ), PIV, WORK, RWORK, 00498 $ NOUT ) 00499 ELSE 00500 WRITE( NOUT, FMT = 9989 )PATH 00501 END IF 00502 * 00503 ELSE IF( LSAMEN( 2, C2, 'PP' ) ) THEN 00504 * 00505 * PP: positive definite packed matrices 00506 * 00507 NTYPES = 9 00508 CALL ALAREQ( PATH, NMATS, DOTYPE, NTYPES, NIN, NOUT ) 00509 * 00510 IF( TSTCHK ) THEN 00511 CALL DCHKPP( DOTYPE, NN, NVAL, NNS, NSVAL, THRESH, TSTERR, 00512 $ LDA, A( 1, 1 ), A( 1, 2 ), A( 1, 3 ), 00513 $ B( 1, 1 ), B( 1, 2 ), B( 1, 3 ), WORK, RWORK, 00514 $ IWORK, NOUT ) 00515 ELSE 00516 WRITE( NOUT, FMT = 9989 )PATH 00517 END IF 00518 * 00519 IF( TSTDRV ) THEN 00520 CALL DDRVPP( DOTYPE, NN, NVAL, NRHS, THRESH, TSTERR, LDA, 00521 $ A( 1, 1 ), A( 1, 2 ), A( 1, 3 ), B( 1, 1 ), 00522 $ B( 1, 2 ), B( 1, 3 ), B( 1, 4 ), S, WORK, 00523 $ RWORK, IWORK, NOUT ) 00524 ELSE 00525 WRITE( NOUT, FMT = 9988 )PATH 00526 END IF 00527 * 00528 ELSE IF( LSAMEN( 2, C2, 'PB' ) ) THEN 00529 * 00530 * PB: positive definite banded matrices 00531 * 00532 NTYPES = 8 00533 CALL ALAREQ( PATH, NMATS, DOTYPE, NTYPES, NIN, NOUT ) 00534 * 00535 IF( TSTCHK ) THEN 00536 CALL DCHKPB( DOTYPE, NN, NVAL, NNB2, NBVAL2, NNS, NSVAL, 00537 $ THRESH, TSTERR, LDA, A( 1, 1 ), A( 1, 2 ), 00538 $ A( 1, 3 ), B( 1, 1 ), B( 1, 2 ), B( 1, 3 ), 00539 $ WORK, RWORK, IWORK, NOUT ) 00540 ELSE 00541 WRITE( NOUT, FMT = 9989 )PATH 00542 END IF 00543 * 00544 IF( TSTDRV ) THEN 00545 CALL DDRVPB( DOTYPE, NN, NVAL, NRHS, THRESH, TSTERR, LDA, 00546 $ A( 1, 1 ), A( 1, 2 ), A( 1, 3 ), B( 1, 1 ), 00547 $ B( 1, 2 ), B( 1, 3 ), B( 1, 4 ), S, WORK, 00548 $ RWORK, IWORK, NOUT ) 00549 ELSE 00550 WRITE( NOUT, FMT = 9988 )PATH 00551 END IF 00552 * 00553 ELSE IF( LSAMEN( 2, C2, 'PT' ) ) THEN 00554 * 00555 * PT: positive definite tridiagonal matrices 00556 * 00557 NTYPES = 12 00558 CALL ALAREQ( PATH, NMATS, DOTYPE, NTYPES, NIN, NOUT ) 00559 * 00560 IF( TSTCHK ) THEN 00561 CALL DCHKPT( DOTYPE, NN, NVAL, NNS, NSVAL, THRESH, TSTERR, 00562 $ A( 1, 1 ), A( 1, 2 ), A( 1, 3 ), B( 1, 1 ), 00563 $ B( 1, 2 ), B( 1, 3 ), WORK, RWORK, NOUT ) 00564 ELSE 00565 WRITE( NOUT, FMT = 9989 )PATH 00566 END IF 00567 * 00568 IF( TSTDRV ) THEN 00569 CALL DDRVPT( DOTYPE, NN, NVAL, NRHS, THRESH, TSTERR, 00570 $ A( 1, 1 ), A( 1, 2 ), A( 1, 3 ), B( 1, 1 ), 00571 $ B( 1, 2 ), B( 1, 3 ), WORK, RWORK, NOUT ) 00572 ELSE 00573 WRITE( NOUT, FMT = 9988 )PATH 00574 END IF 00575 * 00576 ELSE IF( LSAMEN( 2, C2, 'SY' ) ) THEN 00577 * 00578 * SY: symmetric indefinite matrices 00579 * 00580 NTYPES = 10 00581 CALL ALAREQ( PATH, NMATS, DOTYPE, NTYPES, NIN, NOUT ) 00582 * 00583 IF( TSTCHK ) THEN 00584 CALL DCHKSY( DOTYPE, NN, NVAL, NNB2, NBVAL2, NNS, NSVAL, 00585 $ THRESH, TSTERR, LDA, A( 1, 1 ), A( 1, 2 ), 00586 $ A( 1, 3 ), B( 1, 1 ), B( 1, 2 ), B( 1, 3 ), 00587 $ WORK, RWORK, IWORK, NOUT ) 00588 ELSE 00589 WRITE( NOUT, FMT = 9989 )PATH 00590 END IF 00591 * 00592 IF( TSTDRV ) THEN 00593 CALL DDRVSY( DOTYPE, NN, NVAL, NRHS, THRESH, TSTERR, LDA, 00594 $ A( 1, 1 ), A( 1, 2 ), A( 1, 3 ), B( 1, 1 ), 00595 $ B( 1, 2 ), B( 1, 3 ), WORK, RWORK, IWORK, 00596 $ NOUT ) 00597 ELSE 00598 WRITE( NOUT, FMT = 9988 )PATH 00599 END IF 00600 * 00601 ELSE IF( LSAMEN( 2, C2, 'SP' ) ) THEN 00602 * 00603 * SP: symmetric indefinite packed matrices 00604 * 00605 NTYPES = 10 00606 CALL ALAREQ( PATH, NMATS, DOTYPE, NTYPES, NIN, NOUT ) 00607 * 00608 IF( TSTCHK ) THEN 00609 CALL DCHKSP( DOTYPE, NN, NVAL, NNS, NSVAL, THRESH, TSTERR, 00610 $ LDA, A( 1, 1 ), A( 1, 2 ), A( 1, 3 ), 00611 $ B( 1, 1 ), B( 1, 2 ), B( 1, 3 ), WORK, RWORK, 00612 $ IWORK, NOUT ) 00613 ELSE 00614 WRITE( NOUT, FMT = 9989 )PATH 00615 END IF 00616 * 00617 IF( TSTDRV ) THEN 00618 CALL DDRVSP( DOTYPE, NN, NVAL, NRHS, THRESH, TSTERR, LDA, 00619 $ A( 1, 1 ), A( 1, 2 ), A( 1, 3 ), B( 1, 1 ), 00620 $ B( 1, 2 ), B( 1, 3 ), WORK, RWORK, IWORK, 00621 $ NOUT ) 00622 ELSE 00623 WRITE( NOUT, FMT = 9988 )PATH 00624 END IF 00625 * 00626 ELSE IF( LSAMEN( 2, C2, 'TR' ) ) THEN 00627 * 00628 * TR: triangular matrices 00629 * 00630 NTYPES = 18 00631 CALL ALAREQ( PATH, NMATS, DOTYPE, NTYPES, NIN, NOUT ) 00632 * 00633 IF( TSTCHK ) THEN 00634 CALL DCHKTR( DOTYPE, NN, NVAL, NNB2, NBVAL2, NNS, NSVAL, 00635 $ THRESH, TSTERR, LDA, A( 1, 1 ), A( 1, 2 ), 00636 $ B( 1, 1 ), B( 1, 2 ), B( 1, 3 ), WORK, RWORK, 00637 $ IWORK, NOUT ) 00638 ELSE 00639 WRITE( NOUT, FMT = 9989 )PATH 00640 END IF 00641 * 00642 ELSE IF( LSAMEN( 2, C2, 'TP' ) ) THEN 00643 * 00644 * TP: triangular packed matrices 00645 * 00646 NTYPES = 18 00647 CALL ALAREQ( PATH, NMATS, DOTYPE, NTYPES, NIN, NOUT ) 00648 * 00649 IF( TSTCHK ) THEN 00650 CALL DCHKTP( DOTYPE, NN, NVAL, NNS, NSVAL, THRESH, TSTERR, 00651 $ LDA, A( 1, 1 ), A( 1, 2 ), B( 1, 1 ), 00652 $ B( 1, 2 ), B( 1, 3 ), WORK, RWORK, IWORK, 00653 $ NOUT ) 00654 ELSE 00655 WRITE( NOUT, FMT = 9989 )PATH 00656 END IF 00657 * 00658 ELSE IF( LSAMEN( 2, C2, 'TB' ) ) THEN 00659 * 00660 * TB: triangular banded matrices 00661 * 00662 NTYPES = 17 00663 CALL ALAREQ( PATH, NMATS, DOTYPE, NTYPES, NIN, NOUT ) 00664 * 00665 IF( TSTCHK ) THEN 00666 CALL DCHKTB( DOTYPE, NN, NVAL, NNS, NSVAL, THRESH, TSTERR, 00667 $ LDA, A( 1, 1 ), A( 1, 2 ), B( 1, 1 ), 00668 $ B( 1, 2 ), B( 1, 3 ), WORK, RWORK, IWORK, 00669 $ NOUT ) 00670 ELSE 00671 WRITE( NOUT, FMT = 9989 )PATH 00672 END IF 00673 * 00674 ELSE IF( LSAMEN( 2, C2, 'QR' ) ) THEN 00675 * 00676 * QR: QR factorization 00677 * 00678 NTYPES = 8 00679 CALL ALAREQ( PATH, NMATS, DOTYPE, NTYPES, NIN, NOUT ) 00680 * 00681 IF( TSTCHK ) THEN 00682 CALL DCHKQR( DOTYPE, NM, MVAL, NN, NVAL, NNB, NBVAL, NXVAL, 00683 $ NRHS, THRESH, TSTERR, NMAX, A( 1, 1 ), 00684 $ A( 1, 2 ), A( 1, 3 ), A( 1, 4 ), A( 1, 5 ), 00685 $ B( 1, 1 ), B( 1, 2 ), B( 1, 3 ), B( 1, 4 ), 00686 $ WORK, RWORK, IWORK, NOUT ) 00687 ELSE 00688 WRITE( NOUT, FMT = 9989 )PATH 00689 END IF 00690 * 00691 ELSE IF( LSAMEN( 2, C2, 'LQ' ) ) THEN 00692 * 00693 * LQ: LQ factorization 00694 * 00695 NTYPES = 8 00696 CALL ALAREQ( PATH, NMATS, DOTYPE, NTYPES, NIN, NOUT ) 00697 * 00698 IF( TSTCHK ) THEN 00699 CALL DCHKLQ( DOTYPE, NM, MVAL, NN, NVAL, NNB, NBVAL, NXVAL, 00700 $ NRHS, THRESH, TSTERR, NMAX, A( 1, 1 ), 00701 $ A( 1, 2 ), A( 1, 3 ), A( 1, 4 ), A( 1, 5 ), 00702 $ B( 1, 1 ), B( 1, 2 ), B( 1, 3 ), B( 1, 4 ), 00703 $ WORK, RWORK, IWORK, NOUT ) 00704 ELSE 00705 WRITE( NOUT, FMT = 9989 )PATH 00706 END IF 00707 * 00708 ELSE IF( LSAMEN( 2, C2, 'QL' ) ) THEN 00709 * 00710 * QL: QL factorization 00711 * 00712 NTYPES = 8 00713 CALL ALAREQ( PATH, NMATS, DOTYPE, NTYPES, NIN, NOUT ) 00714 * 00715 IF( TSTCHK ) THEN 00716 CALL DCHKQL( DOTYPE, NM, MVAL, NN, NVAL, NNB, NBVAL, NXVAL, 00717 $ NRHS, THRESH, TSTERR, NMAX, A( 1, 1 ), 00718 $ A( 1, 2 ), A( 1, 3 ), A( 1, 4 ), A( 1, 5 ), 00719 $ B( 1, 1 ), B( 1, 2 ), B( 1, 3 ), B( 1, 4 ), 00720 $ WORK, RWORK, IWORK, NOUT ) 00721 ELSE 00722 WRITE( NOUT, FMT = 9989 )PATH 00723 END IF 00724 * 00725 ELSE IF( LSAMEN( 2, C2, 'RQ' ) ) THEN 00726 * 00727 * RQ: RQ factorization 00728 * 00729 NTYPES = 8 00730 CALL ALAREQ( PATH, NMATS, DOTYPE, NTYPES, NIN, NOUT ) 00731 * 00732 IF( TSTCHK ) THEN 00733 CALL DCHKRQ( DOTYPE, NM, MVAL, NN, NVAL, NNB, NBVAL, NXVAL, 00734 $ NRHS, THRESH, TSTERR, NMAX, A( 1, 1 ), 00735 $ A( 1, 2 ), A( 1, 3 ), A( 1, 4 ), A( 1, 5 ), 00736 $ B( 1, 1 ), B( 1, 2 ), B( 1, 3 ), B( 1, 4 ), 00737 $ WORK, RWORK, IWORK, NOUT ) 00738 ELSE 00739 WRITE( NOUT, FMT = 9989 )PATH 00740 END IF 00741 * 00742 ELSE IF( LSAMEN( 2, C2, 'QP' ) ) THEN 00743 * 00744 * QP: QR factorization with pivoting 00745 * 00746 NTYPES = 6 00747 CALL ALAREQ( PATH, NMATS, DOTYPE, NTYPES, NIN, NOUT ) 00748 * 00749 IF( TSTCHK ) THEN 00750 CALL DCHKQP( DOTYPE, NM, MVAL, NN, NVAL, THRESH, TSTERR, 00751 $ A( 1, 1 ), A( 1, 2 ), B( 1, 1 ), B( 1, 2 ), 00752 $ B( 1, 3 ), WORK, IWORK, NOUT ) 00753 CALL DCHKQ3( DOTYPE, NM, MVAL, NN, NVAL, NNB, NBVAL, NXVAL, 00754 $ THRESH, A( 1, 1 ), A( 1, 2 ), B( 1, 1 ), 00755 $ B( 1, 2 ), B( 1, 3 ), WORK, IWORK, NOUT ) 00756 ELSE 00757 WRITE( NOUT, FMT = 9989 )PATH 00758 END IF 00759 * 00760 ELSE IF( LSAMEN( 2, C2, 'TZ' ) ) THEN 00761 * 00762 * TZ: Trapezoidal matrix 00763 * 00764 NTYPES = 3 00765 CALL ALAREQ( PATH, NMATS, DOTYPE, NTYPES, NIN, NOUT ) 00766 * 00767 IF( TSTCHK ) THEN 00768 CALL DCHKTZ( DOTYPE, NM, MVAL, NN, NVAL, THRESH, TSTERR, 00769 $ A( 1, 1 ), A( 1, 2 ), B( 1, 1 ), B( 1, 2 ), 00770 $ B( 1, 3 ), WORK, NOUT ) 00771 ELSE 00772 WRITE( NOUT, FMT = 9989 )PATH 00773 END IF 00774 * 00775 ELSE IF( LSAMEN( 2, C2, 'LS' ) ) THEN 00776 * 00777 * LS: Least squares drivers 00778 * 00779 NTYPES = 6 00780 CALL ALAREQ( PATH, NMATS, DOTYPE, NTYPES, NIN, NOUT ) 00781 * 00782 IF( TSTDRV ) THEN 00783 CALL DDRVLS( DOTYPE, NM, MVAL, NN, NVAL, NNS, NSVAL, NNB, 00784 $ NBVAL, NXVAL, THRESH, TSTERR, A( 1, 1 ), 00785 $ A( 1, 2 ), B( 1, 1 ), B( 1, 2 ), B( 1, 3 ), 00786 $ RWORK, RWORK( NMAX+1 ), WORK, IWORK, NOUT ) 00787 ELSE 00788 WRITE( NOUT, FMT = 9988 )PATH 00789 END IF 00790 * 00791 ELSE IF( LSAMEN( 2, C2, 'EQ' ) ) THEN 00792 * 00793 * EQ: Equilibration routines for general and positive definite 00794 * matrices (THREQ should be between 2 and 10) 00795 * 00796 IF( TSTCHK ) THEN 00797 CALL DCHKEQ( THREQ, NOUT ) 00798 ELSE 00799 WRITE( NOUT, FMT = 9989 )PATH 00800 END IF 00801 * 00802 ELSE 00803 * 00804 WRITE( NOUT, FMT = 9990 )PATH 00805 END IF 00806 * 00807 * Go back to get another input line. 00808 * 00809 GO TO 80 00810 * 00811 * Branch to this line when the last record is read. 00812 * 00813 140 CONTINUE 00814 CLOSE ( NIN ) 00815 S2 = DSECND( ) 00816 WRITE( NOUT, FMT = 9998 ) 00817 WRITE( NOUT, FMT = 9997 )S2 - S1 00818 * 00819 9999 FORMAT( / ' Execution not attempted due to input errors' ) 00820 9998 FORMAT( / ' End of tests' ) 00821 9997 FORMAT( ' Total time used = ', F12.2, ' seconds', / ) 00822 9996 FORMAT( ' Invalid input value: ', A4, '=', I6, '; must be >=', 00823 $ I6 ) 00824 9995 FORMAT( ' Invalid input value: ', A4, '=', I6, '; must be <=', 00825 $ I6 ) 00826 9994 FORMAT( ' Tests of the DOUBLE PRECISION LAPACK routines ', 00827 $ / ' LAPACK VERSION ', I1, '.', I1, '.', I1, 00828 $ / / ' The following parameter values will be used:' ) 00829 9993 FORMAT( 4X, A4, ': ', 10I6, / 11X, 10I6 ) 00830 9992 FORMAT( / ' Routines pass computational tests if test ratio is ', 00831 $ 'less than', F8.2, / ) 00832 9991 FORMAT( ' Relative machine ', A, ' is taken to be', D16.6 ) 00833 9990 FORMAT( / 1X, A3, ': Unrecognized path name' ) 00834 9989 FORMAT( / 1X, A3, ' routines were not tested' ) 00835 9988 FORMAT( / 1X, A3, ' driver routines were not tested' ) 00836 * 00837 * End of DCHKAA 00838 * 00839 END
1.7.3 
