diff options
| author | lawrence.mulholland <lawrence.mulholland@8a072113-8704-0410-8d35-dd094bca7971> | 2015-10-13 12:10:47 +0000 |
|---|---|---|
| committer | lawrence.mulholland <lawrence.mulholland@8a072113-8704-0410-8d35-dd094bca7971> | 2015-10-13 12:10:47 +0000 |
| commit | 2491c9367a9c5e827a935679a4755387b348e159 (patch) | |
| tree | f72b0ed01a40ed2f6251e60862ca73a230d9375a | |
| parent | 292c6d6dc87df038746bb9d705ccc3d9fa9987ec (diff) | |
| download | lapack-2491c9367a9c5e827a935679a4755387b348e159.tar.gz lapack-2491c9367a9c5e827a935679a4755387b348e159.tar.bz2 lapack-2491c9367a9c5e827a935679a4755387b348e159.zip | |
bugs in ?GGHD3 : B,LDA-->B,LDB and zeroing of too large a submatrix
| -rw-r--r-- | SRC/cgghd3.f | 6 | ||||
| -rw-r--r-- | SRC/dgghd3.f | 6 | ||||
| -rw-r--r-- | SRC/sgghd3.f | 6 | ||||
| -rw-r--r-- | SRC/zgghd3.f | 6 |
4 files changed, 12 insertions, 12 deletions
diff --git a/SRC/cgghd3.f b/SRC/cgghd3.f index b20cf44f..112b41a1 100644 --- 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 ed6e98ad..812df3f2 100644 --- 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 abd56fa4..3c58aea7 100644 --- 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 47e7d448..9d6e3628 100644 --- 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 * |