From 8e844f51fe7acc5ed08d2b7fbac3aa890cb04236 Mon Sep 17 00:00:00 2001 From: julie Date: Thu, 27 Jan 2011 14:21:35 +0000 Subject: [PATCH] Follow-up: fix for new routines to handle LDA See bug0067 --- SRC/csytrs2.f | 12 ++++++------ SRC/dsytrs2.f | 12 ++++++------ SRC/ssytrs2.f | 12 ++++++------ SRC/zsytrs2.f | 12 ++++++------ 4 files changed, 24 insertions(+), 24 deletions(-) diff --git a/SRC/csytrs2.f b/SRC/csytrs2.f index d3bb8ca..d90ff47 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 d00e815..3e2e674 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 885373d..ff212ea 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 2c53092..31441e5 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) )) ] * -- 2.7.4