From 8785529faec28defac27df15fb1a976f9c6f33a5 Mon Sep 17 00:00:00 2001 From: Hao Xu Date: Fri, 19 Jul 2024 15:24:42 -0400 Subject: [PATCH] Fix d/sSTEBZ --- SRC/dstebz.f | 6 +++--- SRC/sstebz.f | 8 ++++---- 2 files changed, 7 insertions(+), 7 deletions(-) diff --git a/SRC/dstebz.f b/SRC/dstebz.f index 9588b0be20..4dd69600f9 100644 --- a/SRC/dstebz.f +++ b/SRC/dstebz.f @@ -597,7 +597,7 @@ SUBROUTINE DSTEBZ( RANGE, ORDER, N, VL, VU, IL, IU, ABSTOL, D, END IF * IF( IRANGE.GT.1 ) THEN - IF( GU.LT.WL ) THEN + IF( GU.LE.WL ) THEN NWL = NWL + IN NWU = NWU + IN GO TO 70 @@ -701,7 +701,7 @@ SUBROUTINE DSTEBZ( RANGE, ORDER, N, VL, VU, IL, IU, ABSTOL, D, WKILL = W( JE ) END IF 90 CONTINUE - IBLOCK( IW ) = 0 + IF( IW .NE. 0) IBLOCK( IW ) = 0 100 CONTINUE END IF IF( IDISCU.GT.0 ) THEN @@ -716,7 +716,7 @@ SUBROUTINE DSTEBZ( RANGE, ORDER, N, VL, VU, IL, IU, ABSTOL, D, WKILL = W( JE ) END IF 110 CONTINUE - IBLOCK( IW ) = 0 + IF( IW .NE. 0) IBLOCK( IW ) = 0 120 CONTINUE END IF IM = 0 diff --git a/SRC/sstebz.f b/SRC/sstebz.f index b40c88b945..4c53e6124c 100644 --- a/SRC/sstebz.f +++ b/SRC/sstebz.f @@ -293,7 +293,7 @@ SUBROUTINE SSTEBZ( RANGE, ORDER, N, VL, VU, IL, IU, ABSTOL, D, PARAMETER ( ZERO = 0.0E0, ONE = 1.0E0, TWO = 2.0E0, $ HALF = 1.0E0 / TWO ) REAL FUDGE, RELFAC - PARAMETER ( FUDGE = 2.1E0, RELFAC = 2.0E0 ) + PARAMETER ( FUDGE = 2.3E0, RELFAC = 2.0E0 ) * .. * .. Local Scalars .. LOGICAL NCNVRG, TOOFEW @@ -596,7 +596,7 @@ SUBROUTINE SSTEBZ( RANGE, ORDER, N, VL, VU, IL, IU, ABSTOL, D, END IF * IF( IRANGE.GT.1 ) THEN - IF( GU.LT.WL ) THEN + IF( GU.LE.WL ) THEN NWL = NWL + IN NWU = NWU + IN GO TO 70 @@ -700,7 +700,7 @@ SUBROUTINE SSTEBZ( RANGE, ORDER, N, VL, VU, IL, IU, ABSTOL, D, WKILL = W( JE ) END IF 90 CONTINUE - IBLOCK( IW ) = 0 + IF( IW .NE. 0) IBLOCK( IW ) = 0 100 CONTINUE END IF IF( IDISCU.GT.0 ) THEN @@ -715,7 +715,7 @@ SUBROUTINE SSTEBZ( RANGE, ORDER, N, VL, VU, IL, IU, ABSTOL, D, WKILL = W( JE ) END IF 110 CONTINUE - IBLOCK( IW ) = 0 + IF( IW .NE. 0) IBLOCK( IW ) = 0 120 CONTINUE END IF IM = 0