diff options
author | langou <julien.langou@ucdenver.edu> | 2017-06-12 05:16:16 +0000 |
---|---|---|
committer | GitHub <noreply@github.com> | 2017-06-12 05:16:16 +0000 |
commit | a26a82ff4a30a9240767621e8f0fe30390c55446 (patch) | |
tree | 9cb06867074c60a318b468c6c13d7c7044eaef21 /TESTING/LIN/schksy_aa.f | |
parent | 5c392680ef6362c234c2f2ab3dadd9a9148ffb8f (diff) | |
parent | fd5cae5d069bf8544067a105e6c46d7c02598865 (diff) | |
download | lapack-a26a82ff4a30a9240767621e8f0fe30390c55446.tar.gz lapack-a26a82ff4a30a9240767621e8f0fe30390c55446.tar.bz2 lapack-a26a82ff4a30a9240767621e8f0fe30390c55446.zip |
Merge pull request #166 from iyamazaki/master
fixing the definition of "info" in Aasen's
Diffstat (limited to 'TESTING/LIN/schksy_aa.f')
-rw-r--r-- | TESTING/LIN/schksy_aa.f | 69 |
1 files changed, 36 insertions, 33 deletions
diff --git a/TESTING/LIN/schksy_aa.f b/TESTING/LIN/schksy_aa.f index 226f675f..a23b9cb3 100644 --- a/TESTING/LIN/schksy_aa.f +++ b/TESTING/LIN/schksy_aa.f @@ -433,22 +433,22 @@ * Adjust the expected value of INFO to account for * pivoting. * - IF( IZERO.GT.0 ) THEN - J = 1 - K = IZERO - 100 CONTINUE - IF( J.EQ.K ) THEN - K = IWORK( J ) - ELSE IF( IWORK( J ).EQ.K ) THEN - K = J - END IF - IF( J.LT.K ) THEN - J = J + 1 - GO TO 100 - END IF - ELSE +c IF( IZERO.GT.0 ) THEN +c J = 1 +c K = IZERO +c 100 CONTINUE +c IF( J.EQ.K ) THEN +c K = IWORK( J ) +c ELSE IF( IWORK( J ).EQ.K ) THEN +c K = J +c END IF +c IF( J.LT.K ) THEN +c J = J + 1 +c GO TO 100 +c END IF +c ELSE K = 0 - END IF +c END IF * * Check error code from SSYTRF and handle error. * @@ -512,31 +512,34 @@ * Check error code from SSYTRS and handle error. * IF( INFO.NE.0 ) THEN - CALL ALAERH( PATH, 'SSYTRS_AA', INFO, 0, - $ UPLO, N, N, -1, -1, NRHS, IMAT, - $ NFAIL, NERRS, NOUT ) - END IF -* - CALL SLACPY( 'Full', N, NRHS, B, LDA, WORK, LDA ) + IF( IZERO.EQ.0 ) THEN + CALL ALAERH( PATH, 'SSYTRS_AA', INFO, 0, + $ UPLO, N, N, -1, -1, NRHS, IMAT, + $ NFAIL, NERRS, NOUT ) + END IF + ELSE + CALL SLACPY( 'Full', N, NRHS, B, LDA, WORK, LDA + $ ) * -* Compute the residual for the solution +* Compute the residual for the solution * - CALL SPOT02( UPLO, N, NRHS, A, LDA, X, LDA, WORK, - $ LDA, RWORK, RESULT( 2 ) ) + CALL SPOT02( UPLO, N, NRHS, A, LDA, X, LDA, + $ WORK, LDA, RWORK, RESULT( 2 ) ) * * * Print information about the tests that did not pass * the threshold. * - DO 120 K = 2, 2 - IF( RESULT( K ).GE.THRESH ) THEN - IF( NFAIL.EQ.0 .AND. NERRS.EQ.0 ) - $ CALL ALAHD( NOUT, PATH ) - WRITE( NOUT, FMT = 9998 )UPLO, N, NRHS, - $ IMAT, K, RESULT( K ) - NFAIL = NFAIL + 1 - END IF - 120 CONTINUE + DO 120 K = 2, 2 + IF( RESULT( K ).GE.THRESH ) THEN + IF( NFAIL.EQ.0 .AND. NERRS.EQ.0 ) + $ CALL ALAHD( NOUT, PATH ) + WRITE( NOUT, FMT = 9998 )UPLO, N, NRHS, + $ IMAT, K, RESULT( K ) + NFAIL = NFAIL + 1 + END IF + 120 CONTINUE + END IF NRUN = NRUN + 1 * * End do for each value of NRHS in NSVAL. |