Bypass the calculation of *kteqr

sh-zheng · Mar 16, 2024 · bd6ff6c · bd6ff6c
1 parent 99fa055
commit bd6ff6c
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 16 deletions.
diff --git 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
@@ -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