bugs in ?GGHD3 : B,LDA-->B,LDB and zeroing of too large a submatrix

diff --git a/SRC/cgghd3.f b/SRC/cgghd3.f
index b20cf44..112b41a 100644 (file)
--- a/SRC/cgghd3.f
+++ b/SRC/cgghd3.f
@@ -755,9 +755,9 @@
                END DO
             ELSE
 *
-               CALL CLASET( 'Lower', IHI - JCOL + 1, NNB, CZERO, CZERO,
+               CALL CLASET( 'Lower', IHI - JCOL - 1, NNB, CZERO, CZERO,
      $                      A( JCOL + 2, JCOL ), LDA )
-               CALL CLASET( 'Lower', IHI - JCOL + 1, NNB, CZERO, CZERO,
+               CALL CLASET( 'Lower', IHI - JCOL - 1, NNB, CZERO, CZERO,
      $                      B( JCOL + 2, JCOL ), LDB )
             END IF
 *
@@ -810,7 +810,7 @@
 *
                      CALL CUNM22( 'Right', 'No Transpose', TOP, 2*NNB,
      $                            NNB, NNB, WORK( PPWO ), 2*NNB,
-     $                            B( 1, J ), LDA, WORK( PW ),
+     $                            B( 1, J ), LDB, WORK( PW ),
      $                            LWORK-PW+1, IERR )
                   ELSE
 *

diff --git a/SRC/dgghd3.f b/SRC/dgghd3.f
index ed6e98a..812df3f 100644 (file)
--- a/SRC/dgghd3.f
+++ b/SRC/dgghd3.f
@@ -752,9 +752,9 @@
                END DO
             ELSE
 *
-               CALL DLASET( 'Lower', IHI - JCOL + 1, NNB, ZERO, ZERO,
+               CALL DLASET( 'Lower', IHI - JCOL - 1, NNB, ZERO, ZERO,
      $                      A( JCOL + 2, JCOL ), LDA )
-               CALL DLASET( 'Lower', IHI - JCOL + 1, NNB, ZERO, ZERO,
+               CALL DLASET( 'Lower', IHI - JCOL - 1, NNB, ZERO, ZERO,
      $                      B( JCOL + 2, JCOL ), LDB )
             END IF
 *
@@ -807,7 +807,7 @@
 *
                      CALL DORM22( 'Right', 'No Transpose', TOP, 2*NNB,
      $                            NNB, NNB, WORK( PPWO ), 2*NNB,
-     $                            B( 1, J ), LDA, WORK( PW ),
+     $                            B( 1, J ), LDB, WORK( PW ),
      $                            LWORK-PW+1, IERR )
                   ELSE
 *

diff --git a/SRC/sgghd3.f b/SRC/sgghd3.f
index abd56fa..3c58aea 100644 (file)
--- a/SRC/sgghd3.f
+++ b/SRC/sgghd3.f
@@ -752,9 +752,9 @@
                END DO
             ELSE
 *
-               CALL SLASET( 'Lower', IHI - JCOL + 1, NNB, ZERO, ZERO,
+               CALL SLASET( 'Lower', IHI - JCOL - 1, NNB, ZERO, ZERO,
      $                      A( JCOL + 2, JCOL ), LDA )
-               CALL SLASET( 'Lower', IHI - JCOL + 1, NNB, ZERO, ZERO,
+               CALL SLASET( 'Lower', IHI - JCOL - 1, NNB, ZERO, ZERO,
      $                      B( JCOL + 2, JCOL ), LDB )
             END IF
 *
@@ -807,7 +807,7 @@
 *
                      CALL SORM22( 'Right', 'No Transpose', TOP, 2*NNB,
      $                            NNB, NNB, WORK( PPWO ), 2*NNB,
-     $                            B( 1, J ), LDA, WORK( PW ),
+     $                            B( 1, J ), LDB, WORK( PW ),
      $                            LWORK-PW+1, IERR )
                   ELSE
 *

diff --git a/SRC/zgghd3.f b/SRC/zgghd3.f
index 47e7d44..9d6e362 100644 (file)
--- a/SRC/zgghd3.f
+++ b/SRC/zgghd3.f
@@ -750,9 +750,9 @@
                END DO
             ELSE
 *
-               CALL ZLASET( 'Lower', IHI - JCOL + 1, NNB, CZERO, CZERO,
+               CALL ZLASET( 'Lower', IHI - JCOL - 1, NNB, CZERO, CZERO,
      $                      A( JCOL + 2, JCOL ), LDA )
-               CALL ZLASET( 'Lower', IHI - JCOL + 1, NNB, CZERO, CZERO,
+               CALL ZLASET( 'Lower', IHI - JCOL - 1, NNB, CZERO, CZERO,
      $                      B( JCOL + 2, JCOL ), LDB )
             END IF
 *
@@ -805,7 +805,7 @@
 *
                      CALL ZUNM22( 'Right', 'No Transpose', TOP, 2*NNB,
      $                            NNB, NNB, WORK( PPWO ), 2*NNB,
-     $                            B( 1, J ), LDA, WORK( PW ),
+     $                            B( 1, J ), LDB, WORK( PW ),
      $                            LWORK-PW+1, IERR )
                   ELSE
 *