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 | |
parent | 890ff5ee9c7fd351957abd45b6c64990b4c7759d (diff) |
Follow-up: fix for new routines to handle LDA
See bug0067
Diffstat (limited to 'SRC')
-rw-r--r-- | SRC/csytrs2.f | 12 | ||||
-rw-r--r-- | SRC/dsytrs2.f | 12 | ||||
-rw-r--r-- | SRC/ssytrs2.f | 12 | ||||
-rw-r--r-- | SRC/zsytrs2.f | 12 |
4 files changed, 24 insertions, 24 deletions
diff --git a/SRC/csytrs2.f b/SRC/csytrs2.f index d3bb8ca5..d90ff478 100644 --- a/SRC/csytrs2.f +++ b/SRC/csytrs2.f @@ -141,14 +141,14 @@ * * Compute (U \P' * B) -> B [ (U \P' * B) ] * - CALL CTRSM('L','U','N','U',N,NRHS,ONE,A,N,B,N) + CALL CTRSM('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 CSCAL( NRHS, ONE / A( I, I ), B( I, 1 ), N ) + CALL CSCAL( 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 CTRSM('L','U','T','U',N,NRHS,ONE,A,N,B,N) + CALL CTRSM('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 CTRSM('L','L','N','U',N,NRHS,ONE,A,N,B,N) + CALL CTRSM('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 CSCAL( NRHS, ONE / A( I, I ), B( I, 1 ), N ) + CALL CSCAL( 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 CTRSM('L','L','T','U',N,NRHS,ONE,A,N,B,N) + CALL CTRSM('L','L','T','U',N,NRHS,ONE,A,LDA,B,LDB) * * P * B [ P * (L' \ (D \ (L \P' * B) )) ] * 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) )) ] * diff --git a/SRC/ssytrs2.f b/SRC/ssytrs2.f index 885373da..ff212ea1 100644 --- a/SRC/ssytrs2.f +++ b/SRC/ssytrs2.f @@ -141,14 +141,14 @@ * * Compute (U \P' * B) -> B [ (U \P' * B) ] * - CALL STRSM('L','U','N','U',N,NRHS,ONE,A,N,B,N) + CALL STRSM('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 SSCAL( NRHS, ONE / A( I, I ), B( I, 1 ), N ) + CALL SSCAL( 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 STRSM('L','U','T','U',N,NRHS,ONE,A,N,B,N) + CALL STRSM('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 STRSM('L','L','N','U',N,NRHS,ONE,A,N,B,N) + CALL STRSM('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 SSCAL( NRHS, ONE / A( I, I ), B( I, 1 ), N ) + CALL SSCAL( 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 STRSM('L','L','T','U',N,NRHS,ONE,A,N,B,N) + CALL STRSM('L','L','T','U',N,NRHS,ONE,A,LDA,B,LDB) * * P * B [ P * (L' \ (D \ (L \P' * B) )) ] * 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) )) ] * |