Fix bugs

Author: weslleyspereira

SRC/clahqr.f

@@ -220,7 +220,7 @@ SUBROUTINE CLAHQR( WANTT, WANTZ, N, ILO, IHI, H, LDH, W, ILOZ,
       REAL               DAT1
       PARAMETER          ( DAT1 = 3.0e0 / 4.0e0 )
       INTEGER            KEXSH
-      PARAMETER          ( KEXSH = 6 )
+      PARAMETER          ( KEXSH = 10 )
 *     ..
 *     .. Local Scalars ..
       COMPLEX            CDUM, H11, H11S, H22, SC, SUM, T, T1, TEMP, U,
@@ -320,7 +320,7 @@ SUBROUTINE CLAHQR( WANTT, WANTZ, N, ILO, IHI, H, LDH, W, ILOZ,
 *
 *     KDEFL counts the number of iterations since a deflation
 *
-      KDEFL = -2
+      KDEFL = 0
 *
 *     The main loop begins here. I is the loop index and decreases from
 *     IHI to ILO in steps of 1. Each iteration of the loop works

SRC/dlahqr.f

@@ -229,7 +229,7 @@ SUBROUTINE DLAHQR( WANTT, WANTZ, N, ILO, IHI, H, LDH, WR, WI,
       DOUBLE PRECISION   DAT1, DAT2
       PARAMETER          ( DAT1 = 3.0d0 / 4.0d0, DAT2 = -0.4375d0 )
       INTEGER            KEXSH
-      PARAMETER          ( KEXSH = 6 )
+      PARAMETER          ( KEXSH = 10 )
 *     ..
 *     .. Local Scalars ..
       DOUBLE PRECISION   AA, AB, BA, BB, CS, DET, H11, H12, H21, H21S,
@@ -300,7 +300,7 @@ SUBROUTINE DLAHQR( WANTT, WANTZ, N, ILO, IHI, H, LDH, WR, WI,
 *
 *     KDEFL counts the number of iterations since a deflation
 *
-      KDEFL = -2
+      KDEFL = 0
 *
 *     The main loop begins here. I is the loop index and decreases from
 *     IHI to ILO in steps of 1 or 2. Each iteration of the loop works

SRC/slahqr.f

@@ -229,7 +229,7 @@ SUBROUTINE SLAHQR( WANTT, WANTZ, N, ILO, IHI, H, LDH, WR, WI,
       REAL               DAT1, DAT2
       PARAMETER          ( DAT1 = 3.0e0 / 4.0e0, DAT2 = -0.4375e0 )
       INTEGER            KEXSH
-      PARAMETER          ( KEXSH = 6 )
+      PARAMETER          ( KEXSH = 10 )
 *     ..
 *     .. Local Scalars ..
       REAL               AA, AB, BA, BB, CS, DET, H11, H12, H21, H21S,
@@ -300,7 +300,7 @@ SUBROUTINE SLAHQR( WANTT, WANTZ, N, ILO, IHI, H, LDH, WR, WI,
 *
 *     KDEFL counts the number of iterations since a deflation
 *
-      KDEFL = -2
+      KDEFL = 0
 *
 *     The main loop begins here. I is the loop index and decreases from
 *     IHI to ILO in steps of 1 or 2. Each iteration of the loop works
@@ -390,6 +390,7 @@ SUBROUTINE SLAHQR( WANTT, WANTZ, N, ILO, IHI, H, LDH, WR, WI,
             H12 = DAT2*S
             H21 = S
             H22 = H11
+         ELSE
 *
 *           Prepare to use Francis' double shift
 *           (i.e. 2nd degree generalized Rayleigh quotient)

SRC/zlahqr.f

@@ -220,7 +220,7 @@ SUBROUTINE ZLAHQR( WANTT, WANTZ, N, ILO, IHI, H, LDH, W, ILOZ,
       DOUBLE PRECISION   DAT1
       PARAMETER          ( DAT1 = 3.0d0 / 4.0d0 )
       INTEGER            KEXSH
-      PARAMETER          ( KEXSH = 6 )
+      PARAMETER          ( KEXSH = 10 )
 *     ..
 *     .. Local Scalars ..
       COMPLEX*16         CDUM, H11, H11S, H22, SC, SUM, T, T1, TEMP, U,
@@ -320,7 +320,7 @@ SUBROUTINE ZLAHQR( WANTT, WANTZ, N, ILO, IHI, H, LDH, W, ILOZ,
 *
 *     KDEFL counts the number of iterations since a deflation
 *
-      KDEFL = -2
+      KDEFL = 0
 *
 *     The main loop begins here. I is the loop index and decreases from
 *     IHI to ILO in steps of 1. Each iteration of the loop works