Follow-up: fix for new routines to handle LDA

See bug0067

diff --git a/SRC/csytrs2.f b/SRC/csytrs2.f
index d3bb8ca..d90ff47 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 d00e815..3e2e674 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 885373d..ff212ea 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 2c53092..31441e5 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) )) ]
 *