diff options
Diffstat (limited to 'SRC/ssytri2x.f')
| -rw-r--r-- | SRC/ssytri2x.f | 16 |
1 files changed, 8 insertions, 8 deletions
diff --git a/SRC/ssytri2x.f b/SRC/ssytri2x.f index f628aa82..29168dba 100644 --- a/SRC/ssytri2x.f +++ b/SRC/ssytri2x.f @@ -308,15 +308,15 @@ DO WHILE ( I .LE. N ) IF( IPIV(I) .GT. 0 ) THEN IP=IPIV(I) - IF (I .LT. IP) CALL SSYSWAPR( UPLO, N, A, I ,IP ) - IF (I .GT. IP) CALL SSYSWAPR( UPLO, N, A, IP ,I ) + IF (I .LT. IP) CALL SSYSWAPR( UPLO, N, A, LDA, I ,IP ) + IF (I .GT. IP) CALL SSYSWAPR( UPLO, N, A, LDA, IP ,I ) ELSE IP=-IPIV(I) I=I+1 IF ( (I-1) .LT. IP) - $ CALL SSYSWAPR( UPLO, N, A, I-1 ,IP ) + $ CALL SSYSWAPR( UPLO, N, A, LDA, I-1 ,IP ) IF ( (I-1) .GT. IP) - $ CALL SSYSWAPR( UPLO, N, A, IP ,I-1 ) + $ CALL SSYSWAPR( UPLO, N, A, LDA, IP ,I-1 ) ENDIF I=I+1 END DO @@ -487,12 +487,12 @@ DO WHILE ( I .GE. 1 ) IF( IPIV(I) .GT. 0 ) THEN IP=IPIV(I) - IF (I .LT. IP) CALL SSYSWAPR( UPLO, N, A, I ,IP ) - IF (I .GT. IP) CALL SSYSWAPR( UPLO, N, A, IP ,I ) + IF (I .LT. IP) CALL SSYSWAPR( UPLO, N, A, LDA, I ,IP ) + IF (I .GT. IP) CALL SSYSWAPR( UPLO, N, A, LDA, IP ,I ) ELSE IP=-IPIV(I) - IF ( I .LT. IP) CALL SSYSWAPR( UPLO, N, A, I ,IP ) - IF ( I .GT. IP) CALL SSYSWAPR( UPLO, N, A, IP ,I ) + IF ( I .LT. IP) CALL SSYSWAPR( UPLO, N, A, LDA, I ,IP ) + IF ( I .GT. IP) CALL SSYSWAPR( UPLO, N, A, LDA, IP ,I ) I=I-1 ENDIF I=I-1 |