diff --git a/SRC/dkteqr.f b/SRC/dkteqr.f index 08dee3472..933bfd6eb 100644 --- a/SRC/dkteqr.f +++ b/SRC/dkteqr.f @@ -374,11 +374,8 @@ SUBROUTINE DKTEQR( COMPZ, N, D, E, Z, LDZ, WORK, INFO ) S = SIGN(DLAPY2( P, R ), P) * IF(S.EQ.ZERO .OR. R/S.EQ.ZERO) THEN - VA = -ONE - VB = ONE - VC = ONE - VD = ZERO E(M-1) = ZERO + GO TO 40 ELSEIF(E(M-2).EQ.E(M-1)) THEN E(M-2) = SQRT2*E(M-2) E(M-1) = ZERO @@ -651,11 +648,8 @@ SUBROUTINE DKTEQR( COMPZ, N, D, E, Z, LDZ, WORK, INFO ) S = SIGN(DLAPY2( P, R ), P) * IF(S.EQ.ZERO .OR. R/S.EQ.ZERO) THEN - VA = -ONE - VB = ONE - VC = ONE - VD = ZERO E(M) = ZERO + GO TO 90 ELSEIF(E(M).EQ.E(M+1)) THEN E(M) = SQRT2*E(M) E(M+1) = ZERO diff --git a/SRC/skteqr.f b/SRC/skteqr.f index e7d1a44b6..5210705ac 100644 --- a/SRC/skteqr.f +++ b/SRC/skteqr.f @@ -374,11 +374,8 @@ SUBROUTINE SKTEQR( COMPZ, N, D, E, Z, LDZ, WORK, INFO ) S = SIGN(SLAPY2( P, R ), P) * IF(S.EQ.ZERO .OR. R/S.EQ.ZERO) THEN - VA = -ONE - VB = ONE - VC = ONE - VD = ZERO E(M-1) = ZERO + GO TO 40 ELSEIF(E(M-2).EQ.E(M-1)) THEN E(M-2) = SQRT2*E(M-2) E(M-1) = ZERO @@ -651,11 +648,8 @@ SUBROUTINE SKTEQR( COMPZ, N, D, E, Z, LDZ, WORK, INFO ) S = SIGN(SLAPY2( P, R ), P) * IF(S.EQ.ZERO .OR. R/S.EQ.ZERO) THEN - VA = -ONE - VB = ONE - VC = ONE - VD = ZERO E(M) = ZERO + GO TO 90 ELSEIF(E(M).EQ.E(M+1)) THEN E(M) = SQRT2*E(M) E(M+1) = ZERO