Skip to content

Commit

Permalink
Merge pull request #1024 from cenewcombe/ddrvst_bug
Browse files Browse the repository at this point in the history
Fix infinite loop when an error occurs in tests ddrvst and sdrvst
  • Loading branch information
langou authored Jun 13, 2024
2 parents 163c34b + 7f30ba8 commit 8b468db
Show file tree
Hide file tree
Showing 3 changed files with 12 additions and 8 deletions.
8 changes: 5 additions & 3 deletions TESTING/EIG/ddrvst.f
Original file line number Diff line number Diff line change
Expand Up @@ -2772,7 +2772,7 @@ SUBROUTINE DDRVST( NSIZES, NN, NTYPES, DOTYPE, ISEED, THRESH,
RESULT( NTEST ) = ULPINV
RESULT( NTEST+1 ) = ULPINV
RESULT( NTEST+2 ) = ULPINV
GO TO 700
GO TO 1750
END IF
END IF
*
Expand All @@ -2797,13 +2797,13 @@ SUBROUTINE DDRVST( NSIZES, NN, NTYPES, DOTYPE, ISEED, THRESH,
RETURN
ELSE
RESULT( NTEST ) = ULPINV
GO TO 700
GO TO 1750
END IF
END IF
*
IF( M3.EQ.0 .AND. N.GT.0 ) THEN
RESULT( NTEST ) = ULPINV
GO TO 700
GO TO 1750
END IF
*
* Do test 78 (or +54)
Expand All @@ -2819,6 +2819,8 @@ SUBROUTINE DDRVST( NSIZES, NN, NTYPES, DOTYPE, ISEED, THRESH,
$ MAX( UNFL, TEMP3*ULP )
*
CALL DLACPY( ' ', N, N, V, LDU, A, LDA )
*
1750 CONTINUE
*
1720 CONTINUE
*
Expand Down
4 changes: 2 additions & 2 deletions TESTING/EIG/dlahd2.f
Original file line number Diff line number Diff line change
Expand Up @@ -534,8 +534,8 @@ SUBROUTINE DLAHD2( IOUNIT, PATH )
$ / ' 2: norm( I - Q'' Q ) / ( m ulp )',
$ / ' 3: norm( I - PT PT'' ) / ( n ulp )',
$ / ' 4: norm( Y - Q'' C ) / ( norm(Y) max(m,nrhs) ulp )' )
9968 FORMAT( / ' Tests performed: See sdrvst.f' )
9967 FORMAT( / ' Tests performed: See cdrvst.f' )
9968 FORMAT( / ' Tests performed: See ddrvst.f' )
9967 FORMAT( / ' Tests performed: See zdrvst.f' )
*
* End of DLAHD2
*
Expand Down
8 changes: 5 additions & 3 deletions TESTING/EIG/sdrvst.f
Original file line number Diff line number Diff line change
Expand Up @@ -2772,7 +2772,7 @@ SUBROUTINE SDRVST( NSIZES, NN, NTYPES, DOTYPE, ISEED, THRESH,
RESULT( NTEST ) = ULPINV
RESULT( NTEST+1 ) = ULPINV
RESULT( NTEST+2 ) = ULPINV
GO TO 700
GO TO 1750
END IF
END IF
*
Expand All @@ -2797,13 +2797,13 @@ SUBROUTINE SDRVST( NSIZES, NN, NTYPES, DOTYPE, ISEED, THRESH,
RETURN
ELSE
RESULT( NTEST ) = ULPINV
GO TO 700
GO TO 1750
END IF
END IF
*
IF( M3.EQ.0 .AND. N.GT.0 ) THEN
RESULT( NTEST ) = ULPINV
GO TO 700
GO TO 1750
END IF
*
* Do test 78 (or +54)
Expand All @@ -2819,6 +2819,8 @@ SUBROUTINE SDRVST( NSIZES, NN, NTYPES, DOTYPE, ISEED, THRESH,
$ MAX( UNFL, TEMP3*ULP )
*
CALL SLACPY( ' ', N, N, V, LDU, A, LDA )
*
1750 CONTINUE
*
1720 CONTINUE
*
Expand Down

0 comments on commit 8b468db

Please sign in to comment.