diff options
author | julie <julielangou@users.noreply.github.com> | 2015-01-30 19:50:52 +0000 |
---|---|---|
committer | julie <julielangou@users.noreply.github.com> | 2015-01-30 19:50:52 +0000 |
commit | 024f62998124bd38a60f525b410ac43dc28989c6 (patch) | |
tree | 726feb5392231ec7ce9d1e5c704c684e6763d141 | |
parent | c6bcd83265f7b6d44e05ed2336d664ee8622f70a (diff) | |
download | lapack-024f62998124bd38a60f525b410ac43dc28989c6.tar.gz lapack-024f62998124bd38a60f525b410ac43dc28989c6.tar.bz2 lapack-024f62998124bd38a60f525b410ac43dc28989c6.zip |
Ed Anderson Improvement - Jan 26 2015
STGEX2 can call Level 3 BLAS and I did this in LAPACK3E.
I replaced STGEX2 with the LAPACK3E version (when I say STGEX2, I mean DTGEX2 too)
and also fixed an incorrect leading dimension in DTGEX2.
-rw-r--r-- | SRC/dtgex2.f | 6 | ||||
-rw-r--r-- | SRC/stgex2.f | 4 |
2 files changed, 3 insertions, 7 deletions
diff --git a/SRC/dtgex2.f b/SRC/dtgex2.f index a5813390..fa58ce48 100644 --- a/SRC/dtgex2.f +++ b/SRC/dtgex2.f @@ -595,9 +595,7 @@ * * Standardize existing 2-by-2 blocks. * - DO 50 I = 1, M*M - WORK(I) = ZERO - 50 CONTINUE + CALL DLASET( 'Full', M, M, ZERO, ZERO, WORK, M ) WORK( 1 ) = ONE T( 1, 1 ) = ONE IDUM = LWORK - M*M - 2 @@ -668,7 +666,7 @@ $ A( J1, I ), LDA, ZERO, WORK, M ) CALL DLACPY( 'Full', M, N-I+1, WORK, M, A( J1, I ), LDA ) CALL DGEMM( 'T', 'N', M, N-I+1, M, ONE, LI, LDST, - $ B( J1, I ), LDA, ZERO, WORK, M ) + $ B( J1, I ), LDB, ZERO, WORK, M ) CALL DLACPY( 'Full', M, N-I+1, WORK, M, B( J1, I ), LDB ) END IF I = J1 - 1 diff --git a/SRC/stgex2.f b/SRC/stgex2.f index 7630b119..1cbf1a45 100644 --- a/SRC/stgex2.f +++ b/SRC/stgex2.f @@ -595,9 +595,7 @@ * * Standardize existing 2-by-2 blocks. * - DO 50 I = 1, M*M - WORK(I) = ZERO - 50 CONTINUE + CALL SLASET( 'Full', M, M, ZERO, ZERO, WORK, M ) WORK( 1 ) = ONE T( 1, 1 ) = ONE IDUM = LWORK - M*M - 2 |