diff --git a/lapack-netlib/SRC/dlaqz0.c b/lapack-netlib/SRC/dlaqz0.c
index 074c125982..e001a77d6b 100644
--- a/lapack-netlib/SRC/dlaqz0.c
+++ b/lapack-netlib/SRC/dlaqz0.c
@@ -720,7 +720,7 @@ static integer c__1 = 1;
 /*        Shuffle shifts to put double shifts in front */
 /*        This ensures that we don't split up a double shift */
 
-	i__2 = shiftpos + n_undeflated__ - 1;
+	i__2 = shiftpos + n_undeflated__ - 3;
 	for (i__ = shiftpos; i__ <= i__2; i__ += 2) {
 	    if (alphai[i__] != -alphai[i__ + 1]) {
 
diff --git a/lapack-netlib/SRC/dlaqz0.f b/lapack-netlib/SRC/dlaqz0.f
index 67da04824a..f91ac447b2 100644
--- a/lapack-netlib/SRC/dlaqz0.f
+++ b/lapack-netlib/SRC/dlaqz0.f
@@ -691,7 +691,7 @@ RECURSIVE SUBROUTINE DLAQZ0( WANTS, WANTQ, WANTZ, N, ILO, IHI,
 *        Shuffle shifts to put double shifts in front
 *        This ensures that we don't split up a double shift
 *
-         DO I = SHIFTPOS, SHIFTPOS+N_UNDEFLATED-1, 2
+         DO I = SHIFTPOS, SHIFTPOS+N_UNDEFLATED-3, 2
             IF( ALPHAI( I ).NE.-ALPHAI( I+1 ) ) THEN
 *
                SWAP = ALPHAR( I )
diff --git a/lapack-netlib/SRC/slaqz0.c b/lapack-netlib/SRC/slaqz0.c
index 733a4d17eb..7ed8c38eb1 100644
--- a/lapack-netlib/SRC/slaqz0.c
+++ b/lapack-netlib/SRC/slaqz0.c
@@ -956,7 +956,7 @@ static integer c__1 = 1;
 /*        Shuffle shifts to put double shifts in front */
 /*        This ensures that we don't split up a double shift */
 
-	i__2 = shiftpos + n_undeflated__ - 1;
+	i__2 = shiftpos + n_undeflated__ - 3;
 	for (i__ = shiftpos; i__ <= i__2; i__ += 2) {
 	    if (alphai[i__] != -alphai[i__ + 1]) {
 
diff --git a/lapack-netlib/SRC/slaqz0.f b/lapack-netlib/SRC/slaqz0.f
index 2e1a263bd2..504bf488eb 100644
--- a/lapack-netlib/SRC/slaqz0.f
+++ b/lapack-netlib/SRC/slaqz0.f
@@ -689,7 +689,7 @@ RECURSIVE SUBROUTINE SLAQZ0( WANTS, WANTQ, WANTZ, N, ILO, IHI,
 *        Shuffle shifts to put double shifts in front
 *        This ensures that we don't split up a double shift
 *
-         DO I = SHIFTPOS, SHIFTPOS+N_UNDEFLATED-1, 2
+         DO I = SHIFTPOS, SHIFTPOS+N_UNDEFLATED-3, 2
             IF( ALPHAI( I ).NE.-ALPHAI( I+1 ) ) THEN
 *
                SWAP = ALPHAR( I )