dgesdd: Handle norm nan value

Closes Reference-LAPACK#469
jschueller · Jan 8, 2021 · af872c2 · af872c2
1 parent 62aeb5a
commit af872c2
Showing 1 changed file with 4 additions and 0 deletions.
diff --git a/SRC/dgesdd.f b/SRC/dgesdd.f
@@ -599,6 +599,10 @@ SUBROUTINE DGESDD( JOBZ, M, N, A, LDA, S, U, LDU, VT, LDVT,
 *     Scale A if max element outside range [SMLNUM,BIGNUM]
 *
       ANRM = DLANGE( 'M', M, N, A, LDA, DUM )
+      IF( ANRM.NE.ANRM ) THEN
+          INFO = 1
+          RETURN
+      END IF
       ISCL = 0
       IF( ANRM.GT.ZERO .AND. ANRM.LT.SMLNUM ) THEN
          ISCL = 1