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/zsytrs2.f | |
parent | 890ff5ee9c7fd351957abd45b6c64990b4c7759d (diff) |
Follow-up: fix for new routines to handle LDA
See bug0067
Diffstat (limited to 'SRC/zsytrs2.f')
-rw-r--r-- | SRC/zsytrs2.f | 12 |
1 files changed, 6 insertions, 6 deletions
diff --git a/SRC/zsytrs2.f b/SRC/zsytrs2.f index 2c530921..31441e58 100644 --- a/SRC/zsytrs2.f +++ b/SRC/zsytrs2.f @@ -141,14 +141,14 @@ * * Compute (U \P' * B) -> B [ (U \P' * B) ] * - CALL ZTRSM('L','U','N','U',N,NRHS,ONE,A,N,B,N) + CALL ZTRSM('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 ZSCAL( NRHS, ONE / A( I, I ), B( I, 1 ), N ) + CALL ZSCAL( 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) @@ -169,7 +169,7 @@ * * Compute (U' \ B) -> B [ U' \ (D \ (U \P' * B) ) ] * - CALL ZTRSM('L','U','T','U',N,NRHS,ONE,A,N,B,N) + CALL ZTRSM('L','U','T','U',N,NRHS,ONE,A,LDA,B,LDB) * * P * B [ P * (U' \ (D \ (U \P' * B) )) ] * @@ -218,14 +218,14 @@ * * Compute (L \P' * B) -> B [ (L \P' * B) ] * - CALL ZTRSM('L','L','N','U',N,NRHS,ONE,A,N,B,N) + CALL ZTRSM('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 ZSCAL( NRHS, ONE / A( I, I ), B( I, 1 ), N ) + CALL ZSCAL( NRHS, ONE / A( I, I ), B( I, 1 ), LDB ) ELSE AKM1K = WORK(I) AKM1 = A( I, I ) / AKM1K @@ -244,7 +244,7 @@ * * Compute (L' \ B) -> B [ L' \ (D \ (L \P' * B) ) ] * - CALL ZTRSM('L','L','T','U',N,NRHS,ONE,A,N,B,N) + CALL ZTRSM('L','L','T','U',N,NRHS,ONE,A,LDA,B,LDB) * * P * B [ P * (L' \ (D \ (L \P' * B) )) ] * |