diff options
| author | philippe.theveny <philippe.theveny@8a072113-8704-0410-8d35-dd094bca7971> | 2015-02-25 23:40:47 +0000 |
|---|---|---|
| committer | philippe.theveny <philippe.theveny@8a072113-8704-0410-8d35-dd094bca7971> | 2015-02-25 23:40:47 +0000 |
| commit | 409a29e546c393994eccc8e90a370d13b6f10d4c (patch) | |
| tree | 53634db59c218c99f0e124785ef3ee6f9782edf2 | |
| parent | ba12bdc7bc6c5795eb5e11c49fcd4d9058f1b6da (diff) | |
| download | lapack-409a29e546c393994eccc8e90a370d13b6f10d4c.tar.gz lapack-409a29e546c393994eccc8e90a370d13b6f10d4c.tar.bz2 lapack-409a29e546c393994eccc8e90a370d13b6f10d4c.zip | |
Fix wrong update when the right Schur vector is not computed.
| -rw-r--r-- | SRC/cgghd3.f | 12 | ||||
| -rw-r--r-- | SRC/dgghd3.f | 12 | ||||
| -rw-r--r-- | SRC/sgghd3.f | 12 | ||||
| -rw-r--r-- | SRC/zgghd3.f | 12 |
4 files changed, 20 insertions, 28 deletions
diff --git a/SRC/cgghd3.f b/SRC/cgghd3.f index 347d799f..b20cf44f 100644 --- a/SRC/cgghd3.f +++ b/SRC/cgghd3.f @@ -754,13 +754,11 @@ END DO END DO ELSE - - DO J = ILO, ILO+NNB - DO I = J+2, IHI - A( I, J ) = CZERO - B( I, J ) = CZERO - END DO - END DO +* + CALL CLASET( 'Lower', IHI - JCOL + 1, NNB, CZERO, CZERO, + $ A( JCOL + 2, JCOL ), LDA ) + CALL CLASET( 'Lower', IHI - JCOL + 1, NNB, CZERO, CZERO, + $ B( JCOL + 2, JCOL ), LDB ) END IF * * Apply accumulated unitary matrices to A and B. diff --git a/SRC/dgghd3.f b/SRC/dgghd3.f index 7bed5cc4..ed6e98ad 100644 --- a/SRC/dgghd3.f +++ b/SRC/dgghd3.f @@ -751,13 +751,11 @@ END DO END DO ELSE - - DO J = ILO, ILO+NNB - DO I = J+2, IHI - A( I, J ) = ZERO - B( I, J ) = ZERO - END DO - END DO +* + CALL DLASET( 'Lower', IHI - JCOL + 1, NNB, ZERO, ZERO, + $ A( JCOL + 2, JCOL ), LDA ) + CALL DLASET( 'Lower', IHI - JCOL + 1, NNB, ZERO, ZERO, + $ B( JCOL + 2, JCOL ), LDB ) END IF * * Apply accumulated orthogonal matrices to A and B. diff --git a/SRC/sgghd3.f b/SRC/sgghd3.f index bf91f559..abd56fa4 100644 --- a/SRC/sgghd3.f +++ b/SRC/sgghd3.f @@ -751,13 +751,11 @@ END DO END DO ELSE - - DO J = ILO, ILO+NNB - DO I = J+2, IHI - A( I, J ) = ZERO - B( I, J ) = ZERO - END DO - END DO +* + CALL SLASET( 'Lower', IHI - JCOL + 1, NNB, ZERO, ZERO, + $ A( JCOL + 2, JCOL ), LDA ) + CALL SLASET( 'Lower', IHI - JCOL + 1, NNB, ZERO, ZERO, + $ B( JCOL + 2, JCOL ), LDB ) END IF * * Apply accumulated orthogonal matrices to A and B. diff --git a/SRC/zgghd3.f b/SRC/zgghd3.f index 55952a4b..47e7d448 100644 --- a/SRC/zgghd3.f +++ b/SRC/zgghd3.f @@ -749,13 +749,11 @@ END DO END DO ELSE - - DO J = ILO, ILO+NNB - DO I = J+2, IHI - A( I, J ) = CZERO - B( I, J ) = CZERO - END DO - END DO +* + CALL ZLASET( 'Lower', IHI - JCOL + 1, NNB, CZERO, CZERO, + $ A( JCOL + 2, JCOL ), LDA ) + CALL ZLASET( 'Lower', IHI - JCOL + 1, NNB, CZERO, CZERO, + $ B( JCOL + 2, JCOL ), LDB ) END IF * * Apply accumulated unitary matrices to A and B. |