Revert "Avoid memory corruption with non-finite inputs in ?GEEV"

martin-frbg · web-flow · commit caf674f0d261 · 2025-06-26T16:41:21.000+02:00
diff --git a/SRC/cgeev.f b/SRC/cgeev.f
@@ -496,7 +496,7 @@ SUBROUTINE CGEEV( JOBVL, JOBVR, N, A, LDA, W, VL, LDVL, VR,
 *     Undo scaling if necessary
 *
    50 CONTINUE
-      IF( SCALEA .AND. INFO.GE.0) THEN
+      IF( SCALEA ) THEN
          CALL CLASCL( 'G', 0, 0, CSCALE, ANRM, N-INFO, 1,
      $                W( INFO+1 ),
      $                MAX( N-INFO, 1 ), IERR )
diff --git a/SRC/dgeev.f b/SRC/dgeev.f
@@ -517,7 +517,7 @@ SUBROUTINE DGEEV( JOBVL, JOBVR, N, A, LDA, WR, WI, VL, LDVL,
 *     Undo scaling if necessary
 *
    50 CONTINUE
-      IF( SCALEA .AND. INFO.GE.0 ) THEN
+      IF( SCALEA ) THEN
          CALL DLASCL( 'G', 0, 0, CSCALE, ANRM, N-INFO, 1,
      $                WR( INFO+1 ),
      $                MAX( N-INFO, 1 ), IERR )
diff --git a/SRC/sgeev.f b/SRC/sgeev.f
@@ -519,7 +519,7 @@ SUBROUTINE SGEEV( JOBVL, JOBVR, N, A, LDA, WR, WI, VL, LDVL,
 *     Undo scaling if necessary
 *
    50 CONTINUE
-      IF( SCALEA .AND. INFO.GE.0 ) THEN
+      IF( SCALEA ) THEN
          CALL SLASCL( 'G', 0, 0, CSCALE, ANRM, N-INFO, 1,
      $                WR( INFO+1 ),
      $                MAX( N-INFO, 1 ), IERR )
diff --git a/SRC/zgeev.f b/SRC/zgeev.f
@@ -493,7 +493,7 @@ SUBROUTINE ZGEEV( JOBVL, JOBVR, N, A, LDA, W, VL, LDVL, VR,
 *     Undo scaling if necessary
 *
    50 CONTINUE
-      IF( SCALEA .AND.INFO.GE.0) THEN
+      IF( SCALEA ) THEN
          CALL ZLASCL( 'G', 0, 0, CSCALE, ANRM, N-INFO, 1,
      $                W( INFO+1 ),
      $                MAX( N-INFO, 1 ), IERR )

Original file line number	Diff line number	Diff line change
`@@ -496,7 +496,7 @@ SUBROUTINE CGEEV( JOBVL, JOBVR, N, A, LDA, W, VL, LDVL, VR,`
`496`	`496`	`* Undo scaling if necessary`
`497`	`497`	`*`
`498`	`498`	`50 CONTINUE`
`499`		`- IF( SCALEA .AND. INFO.GE.0) THEN`
	`499`	`+ IF( SCALEA ) THEN`
`500`	`500`	`CALL CLASCL( 'G', 0, 0, CSCALE, ANRM, N-INFO, 1,`
`501`	`501`	`$ W( INFO+1 ),`
`502`	`502`	`$ MAX( N-INFO, 1 ), IERR )`
Original file line number	Diff line number	Diff line change
`@@ -517,7 +517,7 @@ SUBROUTINE DGEEV( JOBVL, JOBVR, N, A, LDA, WR, WI, VL, LDVL,`
`517`	`517`	`* Undo scaling if necessary`
`518`	`518`	`*`
`519`	`519`	`50 CONTINUE`
`520`		`- IF( SCALEA .AND. INFO.GE.0 ) THEN`
	`520`	`+ IF( SCALEA ) THEN`
`521`	`521`	`CALL DLASCL( 'G', 0, 0, CSCALE, ANRM, N-INFO, 1,`
`522`	`522`	`$ WR( INFO+1 ),`
`523`	`523`	`$ MAX( N-INFO, 1 ), IERR )`
Original file line number	Diff line number	Diff line change
`@@ -519,7 +519,7 @@ SUBROUTINE SGEEV( JOBVL, JOBVR, N, A, LDA, WR, WI, VL, LDVL,`
`519`	`519`	`* Undo scaling if necessary`
`520`	`520`	`*`
`521`	`521`	`50 CONTINUE`
`522`		`- IF( SCALEA .AND. INFO.GE.0 ) THEN`
	`522`	`+ IF( SCALEA ) THEN`
`523`	`523`	`CALL SLASCL( 'G', 0, 0, CSCALE, ANRM, N-INFO, 1,`
`524`	`524`	`$ WR( INFO+1 ),`
`525`	`525`	`$ MAX( N-INFO, 1 ), IERR )`
Original file line number	Diff line number	Diff line change
`@@ -493,7 +493,7 @@ SUBROUTINE ZGEEV( JOBVL, JOBVR, N, A, LDA, W, VL, LDVL, VR,`
`493`	`493`	`* Undo scaling if necessary`
`494`	`494`	`*`
`495`	`495`	`50 CONTINUE`
`496`		`- IF( SCALEA .AND.INFO.GE.0) THEN`
	`496`	`+ IF( SCALEA ) THEN`
`497`	`497`	`CALL ZLASCL( 'G', 0, 0, CSCALE, ANRM, N-INFO, 1,`
`498`	`498`	`$ W( INFO+1 ),`
`499`	`499`	`$ MAX( N-INFO, 1 ), IERR )`