diff options
| author | langou <langou@users.noreply.github.com> | 2010-11-14 15:46:03 +0000 |
|---|---|---|
| committer | langou <langou@users.noreply.github.com> | 2010-11-14 15:46:03 +0000 |
| commit | f5029122585e96e8f8ad69c71e04c0bd53993ca6 (patch) | |
| tree | 46ef2541e63f67fe83ca577344976a9205c9bd17 /SRC/zsyconv.f | |
| parent | b0f66356a7156ef35bf9628cd5831eab0c10a6d1 (diff) | |
| download | lapack-f5029122585e96e8f8ad69c71e04c0bd53993ca6.tar.gz lapack-f5029122585e96e8f8ad69c71e04c0bd53993ca6.tar.bz2 lapack-f5029122585e96e8f8ad69c71e04c0bd53993ca6.zip | |
correct indentation
Diffstat (limited to 'SRC/zsyconv.f')
| -rw-r--r-- | SRC/zsyconv.f | 138 |
1 files changed, 69 insertions, 69 deletions
diff --git a/SRC/zsyconv.f b/SRC/zsyconv.f index 26af333d..ae01e61c 100644 --- a/SRC/zsyconv.f +++ b/SRC/zsyconv.f @@ -20,8 +20,8 @@ * Purpose * ======= * -* ZSYCONV convert A given by TRF into L and D and vice-versa. -* Get Non-diag elements of D (returned in workspace) and +* ZSYCONV converts A given by ZHETRF into L and D or vice-versa. +* Get nondiagonal elements of D (returned in workspace) and * apply or reverse permutation done in TRF. * * Arguments @@ -110,13 +110,14 @@ * IF( UPPER ) THEN * -* A is UPPER +* A is UPPER * -* Convert A (A is upper) + IF ( CONVERT ) THEN * -* Convert VALUE +* Convert A (A is upper) +* +* Convert VALUE * - IF ( CONVERT ) THEN I=N WORK(1)=ZERO DO WHILE ( I .GT. 1 ) @@ -130,39 +131,38 @@ I=I-1 END DO * -* Convert PERMUTATIONS +* Convert PERMUTATIONS * - I=N - DO WHILE ( I .GE. 1 ) - IF( IPIV(I) .GT. 0) THEN - IP=IPIV(I) - IF( I .LT. N) THEN - DO 12 J= I+1,N - TEMP=A(IP,J) - A(IP,J)=A(I,J) - A(I,J)=TEMP - 12 CONTINUE + I=N + DO WHILE ( I .GE. 1 ) + IF( IPIV(I) .GT. 0) THEN + IP=IPIV(I) + IF( I .LT. N) THEN + DO 12 J= I+1,N + TEMP=A(IP,J) + A(IP,J)=A(I,J) + A(I,J)=TEMP + 12 CONTINUE + ENDIF + ELSE + IP=-IPIV(I) + IF( I .LT. N) THEN + DO 13 J= I+1,N + TEMP=A(IP,J) + A(IP,J)=A(I-1,J) + A(I-1,J)=TEMP + 13 CONTINUE + ENDIF + I=I-1 ENDIF - ELSE - IP=-IPIV(I) - IF( I .LT. N) THEN - DO 13 J= I+1,N - TEMP=A(IP,J) - A(IP,J)=A(I-1,J) - A(I-1,J)=TEMP - 13 CONTINUE - ENDIF - I=I-1 - ENDIF - I=I-1 - END DO - - ELSE + I=I-1 + END DO * -* Revert A (A is upper) + ELSE * +* Revert A (A is upper) * -* Revert PERMUTATIONS +* Revert PERMUTATIONS * I=1 DO WHILE ( I .LE. N ) @@ -189,7 +189,7 @@ I=I+1 END DO * -* Revert VALUE +* Revert VALUE * I=N DO WHILE ( I .GT. 1 ) @@ -200,16 +200,16 @@ I=I-1 END DO END IF +* ELSE * -* A is LOWER +* A is LOWER * IF ( CONVERT ) THEN * -* Convert A (A is lower) -* +* Convert A (A is lower) * -* Convert VALUE +* Convert VALUE * I=1 WORK(N)=ZERO @@ -224,38 +224,38 @@ I=I+1 END DO * -* Convert PERMUTATIONS -* - I=1 - DO WHILE ( I .LE. N ) - IF( IPIV(I) .GT. 0 ) THEN - IP=IPIV(I) - IF (I .GT. 1) THEN - DO 22 J= 1,I-1 - TEMP=A(IP,J) - A(IP,J)=A(I,J) - A(I,J)=TEMP - 22 CONTINUE +* Convert PERMUTATIONS +* + I=1 + DO WHILE ( I .LE. N ) + IF( IPIV(I) .GT. 0 ) THEN + IP=IPIV(I) + IF (I .GT. 1) THEN + DO 22 J= 1,I-1 + TEMP=A(IP,J) + A(IP,J)=A(I,J) + A(I,J)=TEMP + 22 CONTINUE + ENDIF + ELSE + IP=-IPIV(I) + IF (I .GT. 1) THEN + DO 23 J= 1,I-1 + TEMP=A(IP,J) + A(IP,J)=A(I+1,J) + A(I+1,J)=TEMP + 23 CONTINUE + ENDIF + I=I+1 ENDIF - ELSE - IP=-IPIV(I) - IF (I .GT. 1) THEN - DO 23 J= 1,I-1 - TEMP=A(IP,J) - A(IP,J)=A(I+1,J) - A(I+1,J)=TEMP - 23 CONTINUE - ENDIF - I=I+1 - ENDIF - I=I+1 - END DO - ELSE + I=I+1 + END DO * -* Revert A (A is lower) + ELSE * +* Revert A (A is lower) * -* Revert PERMUTATIONS +* Revert PERMUTATIONS * I=N DO WHILE ( I .GE. 1 ) @@ -282,7 +282,7 @@ I=I-1 END DO * -* Revert VALUE +* Revert VALUE * I=1 DO WHILE ( I .LE. N-1 ) @@ -294,7 +294,7 @@ END DO END IF END IF - +* RETURN * * End of ZSYCONV |