Break out of potentially infinite rescaling loop in LAPACK xLARGV/xLARTG/xLARTGP

Reference-LAPACK issue 411
2020-06-07 14:30:20 +02:00
parent 0f9a935a5a
commit 522aaf53bf
8 changed files with 8 additions and 8 deletions
--- a/lapack-netlib/SRC/dlartgp.f
+++ b/lapack-netlib/SRC/dlartgp.f
@@ -161,7 +161,7 @@
            F1 = F1*SAFMN2
            G1 = G1*SAFMN2
            SCALE = MAX( ABS( F1 ), ABS( G1 ) )
-            IF( SCALE.GE.SAFMX2 )
+            IF( SCALE.GE.SAFMX2 .AND. COUNT .LT. 20 )
     $         GO TO 10
            R = SQRT( F1**2+G1**2 )
            CS = F1 / R