 |
LAPACK 3.3.1
Linear Algebra PACKage
|
|
LAPACK 3.3.1
Linear Algebra PACKage
|
00001 SUBROUTINE DDRVGE( DOTYPE, NN, NVAL, NRHS, THRESH, TSTERR, NMAX, 00002 $ A, AFAC, ASAV, B, BSAV, X, XACT, S, WORK, 00003 $ RWORK, IWORK, NOUT ) 00004 * 00005 * -- LAPACK test routine (version 3.1) -- 00006 * Univ. of Tennessee, Univ. of California Berkeley and NAG Ltd.. 00007 * November 2006 00008 * 00009 * .. Scalar Arguments .. 00010 LOGICAL TSTERR 00011 INTEGER NMAX, NN, NOUT, NRHS 00012 DOUBLE PRECISION THRESH 00013 * .. 00014 * .. Array Arguments .. 00015 LOGICAL DOTYPE( * ) 00016 INTEGER IWORK( * ), NVAL( * ) 00017 DOUBLE PRECISION A( * ), AFAC( * ), ASAV( * ), B( * ), 00018 $ BSAV( * ), RWORK( * ), S( * ), WORK( * ), 00019 $ X( * ), XACT( * ) 00020 * .. 00021 * 00022 * Purpose 00023 * ======= 00024 * 00025 * DDRVGE tests the driver routines DGESV and -SVX. 00026 * 00027 * Arguments 00028 * ========= 00029 * 00030 * DOTYPE (input) LOGICAL array, dimension (NTYPES) 00031 * The matrix types to be used for testing. Matrices of type j 00032 * (for 1 <= j <= NTYPES) are used for testing if DOTYPE(j) = 00033 * .TRUE.; if DOTYPE(j) = .FALSE., then type j is not used. 00034 * 00035 * NN (input) INTEGER 00036 * The number of values of N contained in the vector NVAL. 00037 * 00038 * NVAL (input) INTEGER array, dimension (NN) 00039 * The values of the matrix column dimension N. 00040 * 00041 * NRHS (input) INTEGER 00042 * The number of right hand side vectors to be generated for 00043 * each linear system. 00044 * 00045 * THRESH (input) DOUBLE PRECISION 00046 * The threshold value for the test ratios. A result is 00047 * included in the output file if RESULT >= THRESH. To have 00048 * every test ratio printed, use THRESH = 0. 00049 * 00050 * TSTERR (input) LOGICAL 00051 * Flag that indicates whether error exits are to be tested. 00052 * 00053 * NMAX (input) INTEGER 00054 * The maximum value permitted for N, used in dimensioning the 00055 * work arrays. 00056 * 00057 * A (workspace) DOUBLE PRECISION array, dimension (NMAX*NMAX) 00058 * 00059 * AFAC (workspace) DOUBLE PRECISION array, dimension (NMAX*NMAX) 00060 * 00061 * ASAV (workspace) DOUBLE PRECISION array, dimension (NMAX*NMAX) 00062 * 00063 * B (workspace) DOUBLE PRECISION array, dimension (NMAX*NRHS) 00064 * 00065 * BSAV (workspace) DOUBLE PRECISION array, dimension (NMAX*NRHS) 00066 * 00067 * X (workspace) DOUBLE PRECISION array, dimension (NMAX*NRHS) 00068 * 00069 * XACT (workspace) DOUBLE PRECISION array, dimension (NMAX*NRHS) 00070 * 00071 * S (workspace) DOUBLE PRECISION array, dimension (2*NMAX) 00072 * 00073 * WORK (workspace) DOUBLE PRECISION array, dimension 00074 * (NMAX*max(3,NRHS)) 00075 * 00076 * RWORK (workspace) DOUBLE PRECISION array, dimension (2*NRHS+NMAX) 00077 * 00078 * IWORK (workspace) INTEGER array, dimension (2*NMAX) 00079 * 00080 * NOUT (input) INTEGER 00081 * The unit number for output. 00082 * 00083 * ===================================================================== 00084 * 00085 * .. Parameters .. 00086 DOUBLE PRECISION ONE, ZERO 00087 PARAMETER ( ONE = 1.0D+0, ZERO = 0.0D+0 ) 00088 INTEGER NTYPES 00089 PARAMETER ( NTYPES = 11 ) 00090 INTEGER NTESTS 00091 PARAMETER ( NTESTS = 7 ) 00092 INTEGER NTRAN 00093 PARAMETER ( NTRAN = 3 ) 00094 * .. 00095 * .. Local Scalars .. 00096 LOGICAL EQUIL, NOFACT, PREFAC, TRFCON, ZEROT 00097 CHARACTER DIST, EQUED, FACT, TRANS, TYPE, XTYPE 00098 CHARACTER*3 PATH 00099 INTEGER I, IEQUED, IFACT, IMAT, IN, INFO, IOFF, ITRAN, 00100 $ IZERO, K, K1, KL, KU, LDA, LWORK, MODE, N, NB, 00101 $ NBMIN, NERRS, NFACT, NFAIL, NIMAT, NRUN, NT 00102 DOUBLE PRECISION AINVNM, AMAX, ANORM, ANORMI, ANORMO, CNDNUM, 00103 $ COLCND, RCOND, RCONDC, RCONDI, RCONDO, ROLDC, 00104 $ ROLDI, ROLDO, ROWCND, RPVGRW 00105 * .. 00106 * .. Local Arrays .. 00107 CHARACTER EQUEDS( 4 ), FACTS( 3 ), TRANSS( NTRAN ) 00108 INTEGER ISEED( 4 ), ISEEDY( 4 ) 00109 DOUBLE PRECISION RESULT( NTESTS ) 00110 * .. 00111 * .. External Functions .. 00112 LOGICAL LSAME 00113 DOUBLE PRECISION DGET06, DLAMCH, DLANGE, DLANTR 00114 EXTERNAL LSAME, DGET06, DLAMCH, DLANGE, DLANTR 00115 * .. 00116 * .. External Subroutines .. 00117 EXTERNAL ALADHD, ALAERH, ALASVM, DERRVX, DGEEQU, DGESV, 00118 $ DGESVX, DGET01, DGET02, DGET04, DGET07, DGETRF, 00119 $ DGETRI, DLACPY, DLAQGE, DLARHS, DLASET, DLATB4, 00120 $ DLATMS, XLAENV 00121 * .. 00122 * .. Intrinsic Functions .. 00123 INTRINSIC ABS, MAX 00124 * .. 00125 * .. Scalars in Common .. 00126 LOGICAL LERR, OK 00127 CHARACTER*32 SRNAMT 00128 INTEGER INFOT, NUNIT 00129 * .. 00130 * .. Common blocks .. 00131 COMMON / INFOC / INFOT, NUNIT, OK, LERR 00132 COMMON / SRNAMC / SRNAMT 00133 * .. 00134 * .. Data statements .. 00135 DATA ISEEDY / 1988, 1989, 1990, 1991 / 00136 DATA TRANSS / 'N', 'T', 'C' / 00137 DATA FACTS / 'F', 'N', 'E' / 00138 DATA EQUEDS / 'N', 'R', 'C', 'B' / 00139 * .. 00140 * .. Executable Statements .. 00141 * 00142 * Initialize constants and the random number seed. 00143 * 00144 PATH( 1: 1 ) = 'Double precision' 00145 PATH( 2: 3 ) = 'GE' 00146 NRUN = 0 00147 NFAIL = 0 00148 NERRS = 0 00149 DO 10 I = 1, 4 00150 ISEED( I ) = ISEEDY( I ) 00151 10 CONTINUE 00152 * 00153 * Test the error exits 00154 * 00155 IF( TSTERR ) 00156 $ CALL DERRVX( PATH, NOUT ) 00157 INFOT = 0 00158 * 00159 * Set the block size and minimum block size for testing. 00160 * 00161 NB = 1 00162 NBMIN = 2 00163 CALL XLAENV( 1, NB ) 00164 CALL XLAENV( 2, NBMIN ) 00165 * 00166 * Do for each value of N in NVAL 00167 * 00168 DO 90 IN = 1, NN 00169 N = NVAL( IN ) 00170 LDA = MAX( N, 1 ) 00171 XTYPE = 'N' 00172 NIMAT = NTYPES 00173 IF( N.LE.0 ) 00174 $ NIMAT = 1 00175 * 00176 DO 80 IMAT = 1, NIMAT 00177 * 00178 * Do the tests only if DOTYPE( IMAT ) is true. 00179 * 00180 IF( .NOT.DOTYPE( IMAT ) ) 00181 $ GO TO 80 00182 * 00183 * Skip types 5, 6, or 7 if the matrix size is too small. 00184 * 00185 ZEROT = IMAT.GE.5 .AND. IMAT.LE.7 00186 IF( ZEROT .AND. N.LT.IMAT-4 ) 00187 $ GO TO 80 00188 * 00189 * Set up parameters with DLATB4 and generate a test matrix 00190 * with DLATMS. 00191 * 00192 CALL DLATB4( PATH, IMAT, N, N, TYPE, KL, KU, ANORM, MODE, 00193 $ CNDNUM, DIST ) 00194 RCONDC = ONE / CNDNUM 00195 * 00196 SRNAMT = 'DLATMS' 00197 CALL DLATMS( N, N, DIST, ISEED, TYPE, RWORK, MODE, CNDNUM, 00198 $ ANORM, KL, KU, 'No packing', A, LDA, WORK, 00199 $ INFO ) 00200 * 00201 * Check error code from DLATMS. 00202 * 00203 IF( INFO.NE.0 ) THEN 00204 CALL ALAERH( PATH, 'DLATMS', INFO, 0, ' ', N, N, -1, -1, 00205 $ -1, IMAT, NFAIL, NERRS, NOUT ) 00206 GO TO 80 00207 END IF 00208 * 00209 * For types 5-7, zero one or more columns of the matrix to 00210 * test that INFO is returned correctly. 00211 * 00212 IF( ZEROT ) THEN 00213 IF( IMAT.EQ.5 ) THEN 00214 IZERO = 1 00215 ELSE IF( IMAT.EQ.6 ) THEN 00216 IZERO = N 00217 ELSE 00218 IZERO = N / 2 + 1 00219 END IF 00220 IOFF = ( IZERO-1 )*LDA 00221 IF( IMAT.LT.7 ) THEN 00222 DO 20 I = 1, N 00223 A( IOFF+I ) = ZERO 00224 20 CONTINUE 00225 ELSE 00226 CALL DLASET( 'Full', N, N-IZERO+1, ZERO, ZERO, 00227 $ A( IOFF+1 ), LDA ) 00228 END IF 00229 ELSE 00230 IZERO = 0 00231 END IF 00232 * 00233 * Save a copy of the matrix A in ASAV. 00234 * 00235 CALL DLACPY( 'Full', N, N, A, LDA, ASAV, LDA ) 00236 * 00237 DO 70 IEQUED = 1, 4 00238 EQUED = EQUEDS( IEQUED ) 00239 IF( IEQUED.EQ.1 ) THEN 00240 NFACT = 3 00241 ELSE 00242 NFACT = 1 00243 END IF 00244 * 00245 DO 60 IFACT = 1, NFACT 00246 FACT = FACTS( IFACT ) 00247 PREFAC = LSAME( FACT, 'F' ) 00248 NOFACT = LSAME( FACT, 'N' ) 00249 EQUIL = LSAME( FACT, 'E' ) 00250 * 00251 IF( ZEROT ) THEN 00252 IF( PREFAC ) 00253 $ GO TO 60 00254 RCONDO = ZERO 00255 RCONDI = ZERO 00256 * 00257 ELSE IF( .NOT.NOFACT ) THEN 00258 * 00259 * Compute the condition number for comparison with 00260 * the value returned by DGESVX (FACT = 'N' reuses 00261 * the condition number from the previous iteration 00262 * with FACT = 'F'). 00263 * 00264 CALL DLACPY( 'Full', N, N, ASAV, LDA, AFAC, LDA ) 00265 IF( EQUIL .OR. IEQUED.GT.1 ) THEN 00266 * 00267 * Compute row and column scale factors to 00268 * equilibrate the matrix A. 00269 * 00270 CALL DGEEQU( N, N, AFAC, LDA, S, S( N+1 ), 00271 $ ROWCND, COLCND, AMAX, INFO ) 00272 IF( INFO.EQ.0 .AND. N.GT.0 ) THEN 00273 IF( LSAME( EQUED, 'R' ) ) THEN 00274 ROWCND = ZERO 00275 COLCND = ONE 00276 ELSE IF( LSAME( EQUED, 'C' ) ) THEN 00277 ROWCND = ONE 00278 COLCND = ZERO 00279 ELSE IF( LSAME( EQUED, 'B' ) ) THEN 00280 ROWCND = ZERO 00281 COLCND = ZERO 00282 END IF 00283 * 00284 * Equilibrate the matrix. 00285 * 00286 CALL DLAQGE( N, N, AFAC, LDA, S, S( N+1 ), 00287 $ ROWCND, COLCND, AMAX, EQUED ) 00288 END IF 00289 END IF 00290 * 00291 * Save the condition number of the non-equilibrated 00292 * system for use in DGET04. 00293 * 00294 IF( EQUIL ) THEN 00295 ROLDO = RCONDO 00296 ROLDI = RCONDI 00297 END IF 00298 * 00299 * Compute the 1-norm and infinity-norm of A. 00300 * 00301 ANORMO = DLANGE( '1', N, N, AFAC, LDA, RWORK ) 00302 ANORMI = DLANGE( 'I', N, N, AFAC, LDA, RWORK ) 00303 * 00304 * Factor the matrix A. 00305 * 00306 CALL DGETRF( N, N, AFAC, LDA, IWORK, INFO ) 00307 * 00308 * Form the inverse of A. 00309 * 00310 CALL DLACPY( 'Full', N, N, AFAC, LDA, A, LDA ) 00311 LWORK = NMAX*MAX( 3, NRHS ) 00312 CALL DGETRI( N, A, LDA, IWORK, WORK, LWORK, INFO ) 00313 * 00314 * Compute the 1-norm condition number of A. 00315 * 00316 AINVNM = DLANGE( '1', N, N, A, LDA, RWORK ) 00317 IF( ANORMO.LE.ZERO .OR. AINVNM.LE.ZERO ) THEN 00318 RCONDO = ONE 00319 ELSE 00320 RCONDO = ( ONE / ANORMO ) / AINVNM 00321 END IF 00322 * 00323 * Compute the infinity-norm condition number of A. 00324 * 00325 AINVNM = DLANGE( 'I', N, N, A, LDA, RWORK ) 00326 IF( ANORMI.LE.ZERO .OR. AINVNM.LE.ZERO ) THEN 00327 RCONDI = ONE 00328 ELSE 00329 RCONDI = ( ONE / ANORMI ) / AINVNM 00330 END IF 00331 END IF 00332 * 00333 DO 50 ITRAN = 1, NTRAN 00334 * 00335 * Do for each value of TRANS. 00336 * 00337 TRANS = TRANSS( ITRAN ) 00338 IF( ITRAN.EQ.1 ) THEN 00339 RCONDC = RCONDO 00340 ELSE 00341 RCONDC = RCONDI 00342 END IF 00343 * 00344 * Restore the matrix A. 00345 * 00346 CALL DLACPY( 'Full', N, N, ASAV, LDA, A, LDA ) 00347 * 00348 * Form an exact solution and set the right hand side. 00349 * 00350 SRNAMT = 'DLARHS' 00351 CALL DLARHS( PATH, XTYPE, 'Full', TRANS, N, N, KL, 00352 $ KU, NRHS, A, LDA, XACT, LDA, B, LDA, 00353 $ ISEED, INFO ) 00354 XTYPE = 'C' 00355 CALL DLACPY( 'Full', N, NRHS, B, LDA, BSAV, LDA ) 00356 * 00357 IF( NOFACT .AND. ITRAN.EQ.1 ) THEN 00358 * 00359 * --- Test DGESV --- 00360 * 00361 * Compute the LU factorization of the matrix and 00362 * solve the system. 00363 * 00364 CALL DLACPY( 'Full', N, N, A, LDA, AFAC, LDA ) 00365 CALL DLACPY( 'Full', N, NRHS, B, LDA, X, LDA ) 00366 * 00367 SRNAMT = 'DGESV ' 00368 CALL DGESV( N, NRHS, AFAC, LDA, IWORK, X, LDA, 00369 $ INFO ) 00370 * 00371 * Check error code from DGESV . 00372 * 00373 IF( INFO.NE.IZERO ) 00374 $ CALL ALAERH( PATH, 'DGESV ', INFO, IZERO, 00375 $ ' ', N, N, -1, -1, NRHS, IMAT, 00376 $ NFAIL, NERRS, NOUT ) 00377 * 00378 * Reconstruct matrix from factors and compute 00379 * residual. 00380 * 00381 CALL DGET01( N, N, A, LDA, AFAC, LDA, IWORK, 00382 $ RWORK, RESULT( 1 ) ) 00383 NT = 1 00384 IF( IZERO.EQ.0 ) THEN 00385 * 00386 * Compute residual of the computed solution. 00387 * 00388 CALL DLACPY( 'Full', N, NRHS, B, LDA, WORK, 00389 $ LDA ) 00390 CALL DGET02( 'No transpose', N, N, NRHS, A, 00391 $ LDA, X, LDA, WORK, LDA, RWORK, 00392 $ RESULT( 2 ) ) 00393 * 00394 * Check solution from generated exact solution. 00395 * 00396 CALL DGET04( N, NRHS, X, LDA, XACT, LDA, 00397 $ RCONDC, RESULT( 3 ) ) 00398 NT = 3 00399 END IF 00400 * 00401 * Print information about the tests that did not 00402 * pass the threshold. 00403 * 00404 DO 30 K = 1, NT 00405 IF( RESULT( K ).GE.THRESH ) THEN 00406 IF( NFAIL.EQ.0 .AND. NERRS.EQ.0 ) 00407 $ CALL ALADHD( NOUT, PATH ) 00408 WRITE( NOUT, FMT = 9999 )'DGESV ', N, 00409 $ IMAT, K, RESULT( K ) 00410 NFAIL = NFAIL + 1 00411 END IF 00412 30 CONTINUE 00413 NRUN = NRUN + NT 00414 END IF 00415 * 00416 * --- Test DGESVX --- 00417 * 00418 IF( .NOT.PREFAC ) 00419 $ CALL DLASET( 'Full', N, N, ZERO, ZERO, AFAC, 00420 $ LDA ) 00421 CALL DLASET( 'Full', N, NRHS, ZERO, ZERO, X, LDA ) 00422 IF( IEQUED.GT.1 .AND. N.GT.0 ) THEN 00423 * 00424 * Equilibrate the matrix if FACT = 'F' and 00425 * EQUED = 'R', 'C', or 'B'. 00426 * 00427 CALL DLAQGE( N, N, A, LDA, S, S( N+1 ), ROWCND, 00428 $ COLCND, AMAX, EQUED ) 00429 END IF 00430 * 00431 * Solve the system and compute the condition number 00432 * and error bounds using DGESVX. 00433 * 00434 SRNAMT = 'DGESVX' 00435 CALL DGESVX( FACT, TRANS, N, NRHS, A, LDA, AFAC, 00436 $ LDA, IWORK, EQUED, S, S( N+1 ), B, 00437 $ LDA, X, LDA, RCOND, RWORK, 00438 $ RWORK( NRHS+1 ), WORK, IWORK( N+1 ), 00439 $ INFO ) 00440 * 00441 * Check the error code from DGESVX. 00442 * 00443 IF( INFO.NE.IZERO ) 00444 $ CALL ALAERH( PATH, 'DGESVX', INFO, IZERO, 00445 $ FACT // TRANS, N, N, -1, -1, NRHS, 00446 $ IMAT, NFAIL, NERRS, NOUT ) 00447 * 00448 * Compare WORK(1) from DGESVX with the computed 00449 * reciprocal pivot growth factor RPVGRW 00450 * 00451 IF( INFO.NE.0 ) THEN 00452 RPVGRW = DLANTR( 'M', 'U', 'N', INFO, INFO, 00453 $ AFAC, LDA, WORK ) 00454 IF( RPVGRW.EQ.ZERO ) THEN 00455 RPVGRW = ONE 00456 ELSE 00457 RPVGRW = DLANGE( 'M', N, INFO, A, LDA, 00458 $ WORK ) / RPVGRW 00459 END IF 00460 ELSE 00461 RPVGRW = DLANTR( 'M', 'U', 'N', N, N, AFAC, LDA, 00462 $ WORK ) 00463 IF( RPVGRW.EQ.ZERO ) THEN 00464 RPVGRW = ONE 00465 ELSE 00466 RPVGRW = DLANGE( 'M', N, N, A, LDA, WORK ) / 00467 $ RPVGRW 00468 END IF 00469 END IF 00470 RESULT( 7 ) = ABS( RPVGRW-WORK( 1 ) ) / 00471 $ MAX( WORK( 1 ), RPVGRW ) / 00472 $ DLAMCH( 'E' ) 00473 * 00474 IF( .NOT.PREFAC ) THEN 00475 * 00476 * Reconstruct matrix from factors and compute 00477 * residual. 00478 * 00479 CALL DGET01( N, N, A, LDA, AFAC, LDA, IWORK, 00480 $ RWORK( 2*NRHS+1 ), RESULT( 1 ) ) 00481 K1 = 1 00482 ELSE 00483 K1 = 2 00484 END IF 00485 * 00486 IF( INFO.EQ.0 ) THEN 00487 TRFCON = .FALSE. 00488 * 00489 * Compute residual of the computed solution. 00490 * 00491 CALL DLACPY( 'Full', N, NRHS, BSAV, LDA, WORK, 00492 $ LDA ) 00493 CALL DGET02( TRANS, N, N, NRHS, ASAV, LDA, X, 00494 $ LDA, WORK, LDA, RWORK( 2*NRHS+1 ), 00495 $ RESULT( 2 ) ) 00496 * 00497 * Check solution from generated exact solution. 00498 * 00499 IF( NOFACT .OR. ( PREFAC .AND. LSAME( EQUED, 00500 $ 'N' ) ) ) THEN 00501 CALL DGET04( N, NRHS, X, LDA, XACT, LDA, 00502 $ RCONDC, RESULT( 3 ) ) 00503 ELSE 00504 IF( ITRAN.EQ.1 ) THEN 00505 ROLDC = ROLDO 00506 ELSE 00507 ROLDC = ROLDI 00508 END IF 00509 CALL DGET04( N, NRHS, X, LDA, XACT, LDA, 00510 $ ROLDC, RESULT( 3 ) ) 00511 END IF 00512 * 00513 * Check the error bounds from iterative 00514 * refinement. 00515 * 00516 CALL DGET07( TRANS, N, NRHS, ASAV, LDA, B, LDA, 00517 $ X, LDA, XACT, LDA, RWORK, .TRUE., 00518 $ RWORK( NRHS+1 ), RESULT( 4 ) ) 00519 ELSE 00520 TRFCON = .TRUE. 00521 END IF 00522 * 00523 * Compare RCOND from DGESVX with the computed value 00524 * in RCONDC. 00525 * 00526 RESULT( 6 ) = DGET06( RCOND, RCONDC ) 00527 * 00528 * Print information about the tests that did not pass 00529 * the threshold. 00530 * 00531 IF( .NOT.TRFCON ) THEN 00532 DO 40 K = K1, NTESTS 00533 IF( RESULT( K ).GE.THRESH ) THEN 00534 IF( NFAIL.EQ.0 .AND. NERRS.EQ.0 ) 00535 $ CALL ALADHD( NOUT, PATH ) 00536 IF( PREFAC ) THEN 00537 WRITE( NOUT, FMT = 9997 )'DGESVX', 00538 $ FACT, TRANS, N, EQUED, IMAT, K, 00539 $ RESULT( K ) 00540 ELSE 00541 WRITE( NOUT, FMT = 9998 )'DGESVX', 00542 $ FACT, TRANS, N, IMAT, K, RESULT( K ) 00543 END IF 00544 NFAIL = NFAIL + 1 00545 END IF 00546 40 CONTINUE 00547 NRUN = NRUN + 7 - K1 00548 ELSE 00549 IF( RESULT( 1 ).GE.THRESH .AND. .NOT.PREFAC ) 00550 $ THEN 00551 IF( NFAIL.EQ.0 .AND. NERRS.EQ.0 ) 00552 $ CALL ALADHD( NOUT, PATH ) 00553 IF( PREFAC ) THEN 00554 WRITE( NOUT, FMT = 9997 )'DGESVX', FACT, 00555 $ TRANS, N, EQUED, IMAT, 1, RESULT( 1 ) 00556 ELSE 00557 WRITE( NOUT, FMT = 9998 )'DGESVX', FACT, 00558 $ TRANS, N, IMAT, 1, RESULT( 1 ) 00559 END IF 00560 NFAIL = NFAIL + 1 00561 NRUN = NRUN + 1 00562 END IF 00563 IF( RESULT( 6 ).GE.THRESH ) THEN 00564 IF( NFAIL.EQ.0 .AND. NERRS.EQ.0 ) 00565 $ CALL ALADHD( NOUT, PATH ) 00566 IF( PREFAC ) THEN 00567 WRITE( NOUT, FMT = 9997 )'DGESVX', FACT, 00568 $ TRANS, N, EQUED, IMAT, 6, RESULT( 6 ) 00569 ELSE 00570 WRITE( NOUT, FMT = 9998 )'DGESVX', FACT, 00571 $ TRANS, N, IMAT, 6, RESULT( 6 ) 00572 END IF 00573 NFAIL = NFAIL + 1 00574 NRUN = NRUN + 1 00575 END IF 00576 IF( RESULT( 7 ).GE.THRESH ) THEN 00577 IF( NFAIL.EQ.0 .AND. NERRS.EQ.0 ) 00578 $ CALL ALADHD( NOUT, PATH ) 00579 IF( PREFAC ) THEN 00580 WRITE( NOUT, FMT = 9997 )'DGESVX', FACT, 00581 $ TRANS, N, EQUED, IMAT, 7, RESULT( 7 ) 00582 ELSE 00583 WRITE( NOUT, FMT = 9998 )'DGESVX', FACT, 00584 $ TRANS, N, IMAT, 7, RESULT( 7 ) 00585 END IF 00586 NFAIL = NFAIL + 1 00587 NRUN = NRUN + 1 00588 END IF 00589 * 00590 END IF 00591 * 00592 50 CONTINUE 00593 60 CONTINUE 00594 70 CONTINUE 00595 80 CONTINUE 00596 90 CONTINUE 00597 * 00598 * Print a summary of the results. 00599 * 00600 CALL ALASVM( PATH, NOUT, NFAIL, NRUN, NERRS ) 00601 * 00602 9999 FORMAT( 1X, A, ', N =', I5, ', type ', I2, ', test(', I2, ') =', 00603 $ G12.5 ) 00604 9998 FORMAT( 1X, A, ', FACT=''', A1, ''', TRANS=''', A1, ''', N=', I5, 00605 $ ', type ', I2, ', test(', I1, ')=', G12.5 ) 00606 9997 FORMAT( 1X, A, ', FACT=''', A1, ''', TRANS=''', A1, ''', N=', I5, 00607 $ ', EQUED=''', A1, ''', type ', I2, ', test(', I1, ')=', 00608 $ G12.5 ) 00609 RETURN 00610 * 00611 * End of DDRVGE 00612 * 00613 END
1.7.3 
