Skip to content

Commit 99acde8

Browse files
author
Ichitaro Yamazaki
committed
skip "singularity" check in sytrf_aa tester
1 parent 3ba9cdf commit 99acde8

6 files changed

Lines changed: 229 additions & 210 deletions

File tree

‎TESTING/LIN/cchkhe_aa.f‎

Lines changed: 41 additions & 38 deletions
Original file line numberDiff line numberDiff line change
@@ -434,22 +434,22 @@ SUBROUTINE CCHKHE_AA( DOTYPE, NN, NVAL, NNB, NBVAL, NNS, NSVAL,
434434
* Adjust the expected value of INFO to account for
435435
* pivoting.
436436
*
437-
IF( IZERO.GT.0 ) THEN
438-
J = 1
439-
K = IZERO
440-
100 CONTINUE
441-
IF( J.EQ.K ) THEN
442-
K = IWORK( J )
443-
ELSE IF( IWORK( J ).EQ.K ) THEN
444-
K = J
445-
END IF
446-
IF( J.LT.K ) THEN
447-
J = J + 1
448-
GO TO 100
449-
END IF
450-
ELSE
437+
c IF( IZERO.GT.0 ) THEN
438+
c J = 1
439+
c K = IZERO
440+
c 100 CONTINUE
441+
c IF( J.EQ.K ) THEN
442+
c K = IWORK( J )
443+
c ELSE IF( IWORK( J ).EQ.K ) THEN
444+
c K = J
445+
c END IF
446+
c IF( J.LT.K ) THEN
447+
c J = J + 1
448+
c GO TO 100
449+
c END IF
450+
c ELSE
451451
K = 0
452-
END IF
452+
c END IF
453453
*
454454
* Check error code from CHETRF and handle error.
455455
*
@@ -512,30 +512,33 @@ SUBROUTINE CCHKHE_AA( DOTYPE, NN, NVAL, NNB, NBVAL, NNS, NSVAL,
512512
* Check error code from CHETRS and handle error.
513513
*
514514
IF( INFO.NE.0 ) THEN
515-
CALL ALAERH( PATH, 'CHETRS_AA', INFO, 0,
516-
$ UPLO, N, N, -1, -1, NRHS, IMAT,
517-
$ NFAIL, NERRS, NOUT )
518-
END IF
519-
*
520-
CALL CLACPY( 'Full', N, NRHS, B, LDA, WORK, LDA )
521-
*
522-
* Compute the residual for the solution
523-
*
524-
CALL CPOT02( UPLO, N, NRHS, A, LDA, X, LDA, WORK,
525-
$ LDA, RWORK, RESULT( 2 ) )
526-
*
527-
* Print information about the tests that did not pass
528-
* the threshold.
529-
*
530-
DO 120 K = 2, 2
531-
IF( RESULT( K ).GE.THRESH ) THEN
532-
IF( NFAIL.EQ.0 .AND. NERRS.EQ.0 )
533-
$ CALL ALAHD( NOUT, PATH )
534-
WRITE( NOUT, FMT = 9998 )UPLO, N, NRHS,
535-
$ IMAT, K, RESULT( K )
536-
NFAIL = NFAIL + 1
515+
IF( IZERO.EQ.0 ) THEN
516+
CALL ALAERH( PATH, 'CHETRS_AA', INFO, 0,
517+
$ UPLO, N, N, -1, -1, NRHS, IMAT,
518+
$ NFAIL, NERRS, NOUT )
537519
END IF
538-
120 CONTINUE
520+
ELSE
521+
CALL CLACPY( 'Full', N, NRHS, B, LDA, WORK, LDA
522+
$ )
523+
*
524+
* Compute the residual for the solution
525+
*
526+
CALL CPOT02( UPLO, N, NRHS, A, LDA, X, LDA,
527+
$ WORK, LDA, RWORK, RESULT( 2 ) )
528+
*
529+
* Print information about the tests that did not pass
530+
* the threshold.
531+
*
532+
DO 120 K = 2, 2
533+
IF( RESULT( K ).GE.THRESH ) THEN
534+
IF( NFAIL.EQ.0 .AND. NERRS.EQ.0 )
535+
$ CALL ALAHD( NOUT, PATH )
536+
WRITE( NOUT, FMT = 9998 )UPLO, N, NRHS,
537+
$ IMAT, K, RESULT( K )
538+
NFAIL = NFAIL + 1
539+
END IF
540+
120 CONTINUE
541+
END IF
539542
NRUN = NRUN + 1
540543
*
541544
* End do for each value of NRHS in NSVAL.

‎TESTING/LIN/cchksy_aa.f‎

Lines changed: 38 additions & 35 deletions
Original file line numberDiff line numberDiff line change
@@ -435,22 +435,22 @@ SUBROUTINE CCHKSY_AA( DOTYPE, NN, NVAL, NNB, NBVAL, NNS, NSVAL,
435435
* Adjust the expected value of INFO to account for
436436
* pivoting.
437437
*
438-
IF( IZERO.GT.0 ) THEN
439-
J = 1
440-
K = IZERO
441-
100 CONTINUE
442-
IF( J.EQ.K ) THEN
443-
K = IWORK( J )
444-
ELSE IF( IWORK( J ).EQ.K ) THEN
445-
K = J
446-
END IF
447-
IF( J.LT.K ) THEN
448-
J = J + 1
449-
GO TO 100
450-
END IF
451-
ELSE
438+
c IF( IZERO.GT.0 ) THEN
439+
c J = 1
440+
c K = IZERO
441+
c 100 CONTINUE
442+
c IF( J.EQ.K ) THEN
443+
c K = IWORK( J )
444+
c ELSE IF( IWORK( J ).EQ.K ) THEN
445+
c K = J
446+
c END IF
447+
c IF( J.LT.K ) THEN
448+
c J = J + 1
449+
c GO TO 100
450+
c END IF
451+
c ELSE
452452
K = 0
453-
END IF
453+
c END IF
454454
*
455455
* Check error code from CSYTRF and handle error.
456456
*
@@ -514,31 +514,34 @@ SUBROUTINE CCHKSY_AA( DOTYPE, NN, NVAL, NNB, NBVAL, NNS, NSVAL,
514514
* Check error code from CSYTRS and handle error.
515515
*
516516
IF( INFO.NE.0 ) THEN
517-
CALL ALAERH( PATH, 'CSYTRS_AA', INFO, 0,
518-
$ UPLO, N, N, -1, -1, NRHS, IMAT,
519-
$ NFAIL, NERRS, NOUT )
520-
END IF
521-
*
522-
CALL CLACPY( 'Full', N, NRHS, B, LDA, WORK, LDA )
517+
IF( IZERO.EQ.0 ) THEN
518+
CALL ALAERH( PATH, 'CSYTRS_AA', INFO, 0,
519+
$ UPLO, N, N, -1, -1, NRHS, IMAT,
520+
$ NFAIL, NERRS, NOUT )
521+
END IF
522+
ELSE
523+
CALL CLACPY( 'Full', N, NRHS, B, LDA, WORK, LDA
524+
$ )
523525
*
524-
* Compute the residual for the solution
526+
* Compute the residual for the solution
525527
*
526-
CALL CSYT02( UPLO, N, NRHS, A, LDA, X, LDA, WORK,
527-
$ LDA, RWORK, RESULT( 2 ) )
528+
CALL CSYT02( UPLO, N, NRHS, A, LDA, X, LDA,
529+
$ WORK, LDA, RWORK, RESULT( 2 ) )
528530
*
529531
*
530-
* Print information about the tests that did not pass
531-
* the threshold.
532+
* Print information about the tests that did not pass
533+
* the threshold.
532534
*
533-
DO 120 K = 2, 2
534-
IF( RESULT( K ).GE.THRESH ) THEN
535-
IF( NFAIL.EQ.0 .AND. NERRS.EQ.0 )
536-
$ CALL ALAHD( NOUT, PATH )
537-
WRITE( NOUT, FMT = 9998 )UPLO, N, NRHS,
538-
$ IMAT, K, RESULT( K )
539-
NFAIL = NFAIL + 1
540-
END IF
541-
120 CONTINUE
535+
DO 120 K = 2, 2
536+
IF( RESULT( K ).GE.THRESH ) THEN
537+
IF( NFAIL.EQ.0 .AND. NERRS.EQ.0 )
538+
$ CALL ALAHD( NOUT, PATH )
539+
WRITE( NOUT, FMT = 9998 )UPLO, N, NRHS,
540+
$ IMAT, K, RESULT( K )
541+
NFAIL = NFAIL + 1
542+
END IF
543+
120 CONTINUE
544+
END IF
542545
NRUN = NRUN + 1
543546
*
544547
* End do for each value of NRHS in NSVAL.

‎TESTING/LIN/dchksy_aa.f‎

Lines changed: 38 additions & 35 deletions
Original file line numberDiff line numberDiff line change
@@ -432,22 +432,22 @@ SUBROUTINE DCHKSY_AA( DOTYPE, NN, NVAL, NNB, NBVAL, NNS, NSVAL,
432432
* Adjust the expected value of INFO to account for
433433
* pivoting.
434434
*
435-
IF( IZERO.GT.0 ) THEN
436-
J = 1
437-
K = IZERO
438-
100 CONTINUE
439-
IF( J.EQ.K ) THEN
440-
K = IWORK( J )
441-
ELSE IF( IWORK( J ).EQ.K ) THEN
442-
K = J
443-
END IF
444-
IF( J.LT.K ) THEN
445-
J = J + 1
446-
GO TO 100
447-
END IF
448-
ELSE
435+
c IF( IZERO.GT.0 ) THEN
436+
c J = 1
437+
c K = IZERO
438+
c 100 CONTINUE
439+
c IF( J.EQ.K ) THEN
440+
c K = IWORK( J )
441+
c ELSE IF( IWORK( J ).EQ.K ) THEN
442+
c K = J
443+
c END IF
444+
c IF( J.LT.K ) THEN
445+
c J = J + 1
446+
c GO TO 100
447+
c END IF
448+
c ELSE
449449
K = 0
450-
END IF
450+
c END IF
451451
*
452452
* Check error code from DSYTRF and handle error.
453453
*
@@ -511,31 +511,34 @@ SUBROUTINE DCHKSY_AA( DOTYPE, NN, NVAL, NNB, NBVAL, NNS, NSVAL,
511511
* Check error code from DSYTRS and handle error.
512512
*
513513
IF( INFO.NE.0 ) THEN
514-
CALL ALAERH( PATH, 'DSYTRS_AA', INFO, 0,
515-
$ UPLO, N, N, -1, -1, NRHS, IMAT,
516-
$ NFAIL, NERRS, NOUT )
517-
END IF
518-
*
519-
CALL DLACPY( 'Full', N, NRHS, B, LDA, WORK, LDA )
514+
IF( IZERO.EQ.0 ) THEN
515+
CALL ALAERH( PATH, 'DSYTRS_AA', INFO, 0,
516+
$ UPLO, N, N, -1, -1, NRHS, IMAT,
517+
$ NFAIL, NERRS, NOUT )
518+
END IF
519+
ELSE
520+
CALL DLACPY( 'Full', N, NRHS, B, LDA, WORK, LDA
521+
$ )
520522
*
521-
* Compute the residual for the solution
523+
* Compute the residual for the solution
522524
*
523-
CALL DPOT02( UPLO, N, NRHS, A, LDA, X, LDA, WORK,
524-
$ LDA, RWORK, RESULT( 2 ) )
525+
CALL DPOT02( UPLO, N, NRHS, A, LDA, X, LDA,
526+
$ WORK, LDA, RWORK, RESULT( 2 ) )
525527
*
526528
*
527-
* Print information about the tests that did not pass
528-
* the threshold.
529+
* Print information about the tests that did not pass
530+
* the threshold.
529531
*
530-
DO 120 K = 2, 2
531-
IF( RESULT( K ).GE.THRESH ) THEN
532-
IF( NFAIL.EQ.0 .AND. NERRS.EQ.0 )
533-
$ CALL ALAHD( NOUT, PATH )
534-
WRITE( NOUT, FMT = 9998 )UPLO, N, NRHS,
535-
$ IMAT, K, RESULT( K )
536-
NFAIL = NFAIL + 1
537-
END IF
538-
120 CONTINUE
532+
DO 120 K = 2, 2
533+
IF( RESULT( K ).GE.THRESH ) THEN
534+
IF( NFAIL.EQ.0 .AND. NERRS.EQ.0 )
535+
$ CALL ALAHD( NOUT, PATH )
536+
WRITE( NOUT, FMT = 9998 )UPLO, N, NRHS,
537+
$ IMAT, K, RESULT( K )
538+
NFAIL = NFAIL + 1
539+
END IF
540+
120 CONTINUE
541+
END IF
539542
NRUN = NRUN + 1
540543
*
541544
* End do for each value of NRHS in NSVAL.

‎TESTING/LIN/schksy_aa.f‎

Lines changed: 36 additions & 33 deletions
Original file line numberDiff line numberDiff line change
@@ -433,22 +433,22 @@ SUBROUTINE SCHKSY_AA( DOTYPE, NN, NVAL, NNB, NBVAL, NNS, NSVAL,
433433
* Adjust the expected value of INFO to account for
434434
* pivoting.
435435
*
436-
IF( IZERO.GT.0 ) THEN
437-
J = 1
438-
K = IZERO
439-
100 CONTINUE
440-
IF( J.EQ.K ) THEN
441-
K = IWORK( J )
442-
ELSE IF( IWORK( J ).EQ.K ) THEN
443-
K = J
444-
END IF
445-
IF( J.LT.K ) THEN
446-
J = J + 1
447-
GO TO 100
448-
END IF
449-
ELSE
436+
c IF( IZERO.GT.0 ) THEN
437+
c J = 1
438+
c K = IZERO
439+
c 100 CONTINUE
440+
c IF( J.EQ.K ) THEN
441+
c K = IWORK( J )
442+
c ELSE IF( IWORK( J ).EQ.K ) THEN
443+
c K = J
444+
c END IF
445+
c IF( J.LT.K ) THEN
446+
c J = J + 1
447+
c GO TO 100
448+
c END IF
449+
c ELSE
450450
K = 0
451-
END IF
451+
c END IF
452452
*
453453
* Check error code from SSYTRF and handle error.
454454
*
@@ -512,31 +512,34 @@ SUBROUTINE SCHKSY_AA( DOTYPE, NN, NVAL, NNB, NBVAL, NNS, NSVAL,
512512
* Check error code from SSYTRS and handle error.
513513
*
514514
IF( INFO.NE.0 ) THEN
515-
CALL ALAERH( PATH, 'SSYTRS_AA', INFO, 0,
516-
$ UPLO, N, N, -1, -1, NRHS, IMAT,
517-
$ NFAIL, NERRS, NOUT )
518-
END IF
519-
*
520-
CALL SLACPY( 'Full', N, NRHS, B, LDA, WORK, LDA )
515+
IF( IZERO.EQ.0 ) THEN
516+
CALL ALAERH( PATH, 'SSYTRS_AA', INFO, 0,
517+
$ UPLO, N, N, -1, -1, NRHS, IMAT,
518+
$ NFAIL, NERRS, NOUT )
519+
END IF
520+
ELSE
521+
CALL SLACPY( 'Full', N, NRHS, B, LDA, WORK, LDA
522+
$ )
521523
*
522-
* Compute the residual for the solution
524+
* Compute the residual for the solution
523525
*
524-
CALL SPOT02( UPLO, N, NRHS, A, LDA, X, LDA, WORK,
525-
$ LDA, RWORK, RESULT( 2 ) )
526+
CALL SPOT02( UPLO, N, NRHS, A, LDA, X, LDA,
527+
$ WORK, LDA, RWORK, RESULT( 2 ) )
526528
*
527529
*
528530
* Print information about the tests that did not pass
529531
* the threshold.
530532
*
531-
DO 120 K = 2, 2
532-
IF( RESULT( K ).GE.THRESH ) THEN
533-
IF( NFAIL.EQ.0 .AND. NERRS.EQ.0 )
534-
$ CALL ALAHD( NOUT, PATH )
535-
WRITE( NOUT, FMT = 9998 )UPLO, N, NRHS,
536-
$ IMAT, K, RESULT( K )
537-
NFAIL = NFAIL + 1
538-
END IF
539-
120 CONTINUE
533+
DO 120 K = 2, 2
534+
IF( RESULT( K ).GE.THRESH ) THEN
535+
IF( NFAIL.EQ.0 .AND. NERRS.EQ.0 )
536+
$ CALL ALAHD( NOUT, PATH )
537+
WRITE( NOUT, FMT = 9998 )UPLO, N, NRHS,
538+
$ IMAT, K, RESULT( K )
539+
NFAIL = NFAIL + 1
540+
END IF
541+
120 CONTINUE
542+
END IF
540543
NRUN = NRUN + 1
541544
*
542545
* End do for each value of NRHS in NSVAL.

0 commit comments

Comments
 (0)