diff options
author | julie <julielangou@users.noreply.github.com> | 2011-01-27 14:21:35 +0000 |
---|---|---|
committer | julie <julielangou@users.noreply.github.com> | 2011-01-27 14:21:35 +0000 |
commit | 8e844f51fe7acc5ed08d2b7fbac3aa890cb04236 (patch) | |
tree | 14cf0f78ca1ac0694fbe780f0698579cb7299540 /SRC/dsytrs2.f | |
parent | 890ff5ee9c7fd351957abd45b6c64990b4c7759d (diff) | |
download | lapack-8e844f51fe7acc5ed08d2b7fbac3aa890cb04236.tar.gz lapack-8e844f51fe7acc5ed08d2b7fbac3aa890cb04236.tar.bz2 lapack-8e844f51fe7acc5ed08d2b7fbac3aa890cb04236.zip |
Follow-up: fix for new routines to handle LDA
See bug0067
Diffstat (limited to 'SRC/dsytrs2.f')
-rw-r--r-- | SRC/dsytrs2.f | 12 |
1 files changed, 6 insertions, 6 deletions
diff --git a/SRC/dsytrs2.f b/SRC/dsytrs2.f index d00e815e..3e2e6748 100644 --- a/SRC/dsytrs2.f +++ b/SRC/dsytrs2.f @@ -140,14 +140,14 @@ * * Compute (U \P' * B) -> B [ (U \P' * B) ] * - CALL DTRSM('L','U','N','U',N,NRHS,ONE,A,N,B,N) + CALL DTRSM('L','U','N','U',N,NRHS,ONE,A,LDA,B,LDB) * * Compute D \ B -> B [ D \ (U \P' * B) ] * I=N DO WHILE ( I .GE. 1 ) IF( IPIV(I) .GT. 0 ) THEN - CALL DSCAL( NRHS, ONE / A( I, I ), B( I, 1 ), N ) + CALL DSCAL( NRHS, ONE / A( I, I ), B( I, 1 ), LDB ) ELSEIF ( I .GT. 1) THEN IF ( IPIV(I-1) .EQ. IPIV(I) ) THEN AKM1K = WORK(I) @@ -168,7 +168,7 @@ * * Compute (U' \ B) -> B [ U' \ (D \ (U \P' * B) ) ] * - CALL DTRSM('L','U','T','U',N,NRHS,ONE,A,N,B,N) + CALL DTRSM('L','U','T','U',N,NRHS,ONE,A,LDA,B,LDB) * * P * B [ P * (U' \ (D \ (U \P' * B) )) ] * @@ -217,14 +217,14 @@ * * Compute (L \P' * B) -> B [ (L \P' * B) ] * - CALL DTRSM('L','L','N','U',N,NRHS,ONE,A,N,B,N) + CALL DTRSM('L','L','N','U',N,NRHS,ONE,A,LDA,B,LDB) * * Compute D \ B -> B [ D \ (L \P' * B) ] * I=1 DO WHILE ( I .LE. N ) IF( IPIV(I) .GT. 0 ) THEN - CALL DSCAL( NRHS, ONE / A( I, I ), B( I, 1 ), N ) + CALL DSCAL( NRHS, ONE / A( I, I ), B( I, 1 ), LDB ) ELSE AKM1K = WORK(I) AKM1 = A( I, I ) / AKM1K @@ -243,7 +243,7 @@ * * Compute (L' \ B) -> B [ L' \ (D \ (L \P' * B) ) ] * - CALL DTRSM('L','L','T','U',N,NRHS,ONE,A,N,B,N) + CALL DTRSM('L','L','T','U',N,NRHS,ONE,A,LDA,B,LDB) * * P * B [ P * (L' \ (D \ (L \P' * B) )) ] * |