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 /SRC/sgghd3.f | |
| 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
Diffstat (limited to 'SRC/sgghd3.f')
| -rw-r--r-- | SRC/sgghd3.f | 6 |
1 files changed, 3 insertions, 3 deletions
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 * |