From f2bcf9267210d21af18d27d5b4e36b869ef19498 Mon Sep 17 00:00:00 2001 From: Martin Kroeker Date: Sun, 21 Sep 2025 17:26:02 +0200 Subject: [PATCH 1/2] Rescale Infinity input early to avoid infinite loop later (Reference-LAPACK PR 1158) --- lapack-netlib/SRC/crscl.f | 6 ++++++ 1 file changed, 6 insertions(+) diff --git a/lapack-netlib/SRC/crscl.f b/lapack-netlib/SRC/crscl.f index 22919cd62c..99a0e73ca5 100644 --- a/lapack-netlib/SRC/crscl.f +++ b/lapack-netlib/SRC/crscl.f @@ -114,6 +114,7 @@ SUBROUTINE CRSCL( N, A, X, INCX ) * .. * .. Intrinsic Functions .. INTRINSIC ABS + INTRINSIC HUGE * .. * .. Executable Statements .. * @@ -121,6 +122,11 @@ SUBROUTINE CRSCL( N, A, X, INCX ) * IF( N.LE.0 ) $ RETURN +* + IF( SA.GT.HUGE(SA) .OR. SA.LT.-HUGE(SA) ) THEN + CALL CSSCAL( N, SA, SX, INCX ) + RETURN + END IF * * Get machine parameters * From 420427efe21ad2f373a85e2dddc4e8254087d2a5 Mon Sep 17 00:00:00 2001 From: Martin Kroeker Date: Sun, 21 Sep 2025 17:28:01 +0200 Subject: [PATCH 2/2] Rescale Infinity input early to avoid infinite loop later (Reference-LAPACK PR 1158) --- lapack-netlib/SRC/zrscl.f | 6 ++++++ 1 file changed, 6 insertions(+) diff --git a/lapack-netlib/SRC/zrscl.f b/lapack-netlib/SRC/zrscl.f index 970f6de752..1eb36b6e4f 100644 --- a/lapack-netlib/SRC/zrscl.f +++ b/lapack-netlib/SRC/zrscl.f @@ -113,6 +113,7 @@ SUBROUTINE ZRSCL( N, A, X, INCX ) * .. * .. Intrinsic Functions .. INTRINSIC ABS + INTRINSIC HUGE * .. * .. Executable Statements .. * @@ -120,6 +121,11 @@ SUBROUTINE ZRSCL( N, A, X, INCX ) * IF( N.LE.0 ) $ RETURN +* + IF( SA.GT.HUGE(SA) .OR. SA.LT.-HUGE(SA) ) THEN + CALL ZDSCAL( N, SA, SX, INCX ) + RETURN + END IF * * Get machine parameters *