don't subtract N_DEFLATED from istop twice

thijssteel · thijssteel · commit 65850558e010 · 2023-02-10T08:51:57.000+01:00
diff --git a/SRC/claqz0.f b/SRC/claqz0.f
@@ -648,7 +648,7 @@ RECURSIVE SUBROUTINE CLAQZ0( WANTS, WANTQ, WANTZ, N, ILO, IHI, A,
 
          NS = MIN( NSHIFTS, ISTOP-ISTART2 )
          NS = MIN( NS, N_UNDEFLATED )
-         SHIFTPOS = ISTOP-N_DEFLATED-N_UNDEFLATED+1
+         SHIFTPOS = ISTOP-N_UNDEFLATED+1
 
          IF ( MOD( LD, 6 ) .EQ. 0 ) THEN
 * 
diff --git a/SRC/dlaqz0.f b/SRC/dlaqz0.f
@@ -682,7 +682,7 @@ RECURSIVE SUBROUTINE DLAQZ0( WANTS, WANTQ, WANTZ, N, ILO, IHI, A,
 
          NS = MIN( NSHIFTS, ISTOP-ISTART2 )
          NS = MIN( NS, N_UNDEFLATED )
-         SHIFTPOS = ISTOP-N_DEFLATED-N_UNDEFLATED+1
+         SHIFTPOS = ISTOP-N_UNDEFLATED+1
 *
 *        Shuffle shifts to put double shifts in front
 *        This ensures that we don't split up a double shift
diff --git a/SRC/slaqz0.f b/SRC/slaqz0.f
@@ -678,7 +678,7 @@ RECURSIVE SUBROUTINE SLAQZ0( WANTS, WANTQ, WANTZ, N, ILO, IHI, A,
 
          NS = MIN( NSHIFTS, ISTOP-ISTART2 )
          NS = MIN( NS, N_UNDEFLATED )
-         SHIFTPOS = ISTOP-N_DEFLATED-N_UNDEFLATED+1
+         SHIFTPOS = ISTOP-N_UNDEFLATED+1
 *
 *        Shuffle shifts to put double shifts in front
 *        This ensures that we don't split up a double shift
diff --git a/SRC/zlaqz0.f b/SRC/zlaqz0.f
@@ -649,7 +649,7 @@ RECURSIVE SUBROUTINE ZLAQZ0( WANTS, WANTQ, WANTZ, N, ILO, IHI, A,
 
          NS = MIN( NSHIFTS, ISTOP-ISTART2 )
          NS = MIN( NS, N_UNDEFLATED )
-         SHIFTPOS = ISTOP-N_DEFLATED-N_UNDEFLATED+1
+         SHIFTPOS = ISTOP-N_UNDEFLATED+1
 
          IF ( MOD( LD, 6 ) .EQ. 0 ) THEN
 * 

Original file line number	Diff line number	Diff line change
`@@ -648,7 +648,7 @@ RECURSIVE SUBROUTINE CLAQZ0( WANTS, WANTQ, WANTZ, N, ILO, IHI, A,`
`648`	`648`
`649`	`649`	`NS = MIN( NSHIFTS, ISTOP-ISTART2 )`
`650`	`650`	`NS = MIN( NS, N_UNDEFLATED )`
`651`		`- SHIFTPOS = ISTOP-N_DEFLATED-N_UNDEFLATED+1`
	`651`	`+ SHIFTPOS = ISTOP-N_UNDEFLATED+1`
`652`	`652`
`653`	`653`	`IF ( MOD( LD, 6 ) .EQ. 0 ) THEN`
`654`	`654`	`*`
Original file line number	Diff line number	Diff line change
`@@ -682,7 +682,7 @@ RECURSIVE SUBROUTINE DLAQZ0( WANTS, WANTQ, WANTZ, N, ILO, IHI, A,`
`682`	`682`
`683`	`683`	`NS = MIN( NSHIFTS, ISTOP-ISTART2 )`
`684`	`684`	`NS = MIN( NS, N_UNDEFLATED )`
`685`		`- SHIFTPOS = ISTOP-N_DEFLATED-N_UNDEFLATED+1`
	`685`	`+ SHIFTPOS = ISTOP-N_UNDEFLATED+1`
`686`	`686`	`*`
`687`	`687`	`* Shuffle shifts to put double shifts in front`
`688`	`688`	`* This ensures that we don't split up a double shift`
Original file line number	Diff line number	Diff line change
`@@ -678,7 +678,7 @@ RECURSIVE SUBROUTINE SLAQZ0( WANTS, WANTQ, WANTZ, N, ILO, IHI, A,`
`678`	`678`
`679`	`679`	`NS = MIN( NSHIFTS, ISTOP-ISTART2 )`
`680`	`680`	`NS = MIN( NS, N_UNDEFLATED )`
`681`		`- SHIFTPOS = ISTOP-N_DEFLATED-N_UNDEFLATED+1`
	`681`	`+ SHIFTPOS = ISTOP-N_UNDEFLATED+1`
`682`	`682`	`*`
`683`	`683`	`* Shuffle shifts to put double shifts in front`
`684`	`684`	`* This ensures that we don't split up a double shift`
Original file line number	Diff line number	Diff line change
`@@ -649,7 +649,7 @@ RECURSIVE SUBROUTINE ZLAQZ0( WANTS, WANTQ, WANTZ, N, ILO, IHI, A,`
`649`	`649`
`650`	`650`	`NS = MIN( NSHIFTS, ISTOP-ISTART2 )`
`651`	`651`	`NS = MIN( NS, N_UNDEFLATED )`
`652`		`- SHIFTPOS = ISTOP-N_DEFLATED-N_UNDEFLATED+1`
	`652`	`+ SHIFTPOS = ISTOP-N_UNDEFLATED+1`
`653`	`653`
`654`	`654`	`IF ( MOD( LD, 6 ) .EQ. 0 ) THEN`
`655`	`655`	`*`