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 * 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 *