summaryrefslogtreecommitdiff
path: root/TESTING/LIN/schksy_aa.f
diff options
context:
space:
mode:
authorlangou <julien.langou@ucdenver.edu>2017-06-12 05:16:16 +0000
committerGitHub <noreply@github.com>2017-06-12 05:16:16 +0000
commita26a82ff4a30a9240767621e8f0fe30390c55446 (patch)
tree9cb06867074c60a318b468c6c13d7c7044eaef21 /TESTING/LIN/schksy_aa.f
parent5c392680ef6362c234c2f2ab3dadd9a9148ffb8f (diff)
parentfd5cae5d069bf8544067a105e6c46d7c02598865 (diff)
downloadlapack-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.f69
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.