Follow-up: fix for new routines to handle LDA

See bug0067

diff --git a/SRC/csytrs2.f b/SRC/csytrs2.f
index d3bb8ca537de523960f482b753f99798e761d449..d90ff4787c7a580706229b47276b327446df6b25 100644 (file)
--- 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 d00e815ee875e418588251747075aac751c8e7d5..3e2e67483b80b57fab81fba8aaea7fe617ce470d 100644 (file)
--- 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 885373dae4bcac17e23f2b2460f67e4ef25c8683..ff212ea12b8145d39d5ef0c57536f7bf291b463e 100644 (file)
--- 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 2c530921d7d02bec00168c2879fef30008b7d445..31441e585c986e97875003b46c5cbe067bc8478f 100644 (file)
--- 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) )) ]
 *