From cf2962bdb5a0d4ad9607fee0b0e4ac4a08e11225 Mon Sep 17 00:00:00 2001 From: Martin Kroeker Date: Tue, 18 Jun 2024 11:15:44 +0200 Subject: [PATCH 1/4] fix possible infinite loop on error (Reference-LAPACK PR 1024) --- lapack-netlib/TESTING/EIG/ddrvst.f | 10 ++++++---- 1 file changed, 6 insertions(+), 4 deletions(-) diff --git a/lapack-netlib/TESTING/EIG/ddrvst.f b/lapack-netlib/TESTING/EIG/ddrvst.f index 805fd8271..b58930300 100644 --- a/lapack-netlib/TESTING/EIG/ddrvst.f +++ b/lapack-netlib/TESTING/EIG/ddrvst.f @@ -2772,7 +2772,7 @@ c LIWEDC = 12 RESULT( NTEST ) = ULPINV RESULT( NTEST+1 ) = ULPINV RESULT( NTEST+2 ) = ULPINV - GO TO 700 + GO TO 1750 END IF END IF * @@ -2797,13 +2797,13 @@ c LIWEDC = 12 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) @@ -2819,7 +2819,9 @@ c LIWEDC = 12 $ MAX( UNFL, TEMP3*ULP ) * CALL DLACPY( ' ', N, N, V, LDU, A, LDA ) -* + + 1750 CONTINUE + 1720 CONTINUE * * End of Loop -- Check for RESULT(j) > THRESH From bf521a2ced44f658edca74152401fabb9d894a0d Mon Sep 17 00:00:00 2001 From: Martin Kroeker Date: Tue, 18 Jun 2024 11:18:04 +0200 Subject: [PATCH 2/4] fix possible infinite loop on error (Reference-LAPACK PR 1024) --- lapack-netlib/TESTING/EIG/sdrvst.f | 8 +++++--- 1 file changed, 5 insertions(+), 3 deletions(-) diff --git a/lapack-netlib/TESTING/EIG/sdrvst.f b/lapack-netlib/TESTING/EIG/sdrvst.f index be6d33cee..2d02f54ed 100644 --- a/lapack-netlib/TESTING/EIG/sdrvst.f +++ b/lapack-netlib/TESTING/EIG/sdrvst.f @@ -2772,7 +2772,7 @@ c LIWEDC = 12 RESULT( NTEST ) = ULPINV RESULT( NTEST+1 ) = ULPINV RESULT( NTEST+2 ) = ULPINV - GO TO 700 + GO TO 1750 END IF END IF * @@ -2797,13 +2797,13 @@ c LIWEDC = 12 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) @@ -2819,6 +2819,8 @@ c LIWEDC = 12 $ MAX( UNFL, TEMP3*ULP ) * CALL SLACPY( ' ', N, N, V, LDU, A, LDA ) +* + 1750 CONTINUE * 1720 CONTINUE * From a9817b421235e3a1cba2297149460c7a4da00d24 Mon Sep 17 00:00:00 2001 From: Martin Kroeker Date: Tue, 18 Jun 2024 11:20:22 +0200 Subject: [PATCH 3/4] fix reference in format (Reference-LAPACK PR 1024) --- lapack-netlib/TESTING/EIG/dlahd2.f | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/lapack-netlib/TESTING/EIG/dlahd2.f b/lapack-netlib/TESTING/EIG/dlahd2.f index a6c65635f..cccbd2537 100644 --- a/lapack-netlib/TESTING/EIG/dlahd2.f +++ b/lapack-netlib/TESTING/EIG/dlahd2.f @@ -534,8 +534,8 @@ $ / ' 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 * From 215279662ed1c9ffd6552d51e351ea8cf15b8dc4 Mon Sep 17 00:00:00 2001 From: Martin Kroeker Date: Tue, 18 Jun 2024 11:21:33 +0200 Subject: [PATCH 4/4] fix possible infinite loop on error (Reference-LAPACK PR 1024) --- lapack-netlib/TESTING/EIG/ddrvst.f | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/lapack-netlib/TESTING/EIG/ddrvst.f b/lapack-netlib/TESTING/EIG/ddrvst.f index b58930300..311684238 100644 --- a/lapack-netlib/TESTING/EIG/ddrvst.f +++ b/lapack-netlib/TESTING/EIG/ddrvst.f @@ -2819,9 +2819,9 @@ c LIWEDC = 12 $ MAX( UNFL, TEMP3*ULP ) * CALL DLACPY( ' ', N, N, V, LDU, A, LDA ) - +* 1750 CONTINUE - +* 1720 CONTINUE * * End of Loop -- Check for RESULT(j) > THRESH