diff --git a/SRC/dstebz.f b/SRC/dstebz.f index 9588b0be2..4dd69600f 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 b40c88b94..4c53e6124 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