diff options
| author | julie <julielangou@users.noreply.github.com> | 2011-03-10 19:09:21 +0000 |
|---|---|---|
| committer | julie <julielangou@users.noreply.github.com> | 2011-03-10 19:09:21 +0000 |
| commit | 1923c69f7e9857d19b67f9e9befd1b1f37c644ef (patch) | |
| tree | a0729e3779db5c74c0dcda417a0d6eff5c8b3f95 /SRC/ctftri.f | |
| parent | 86579dddaf0142000a291a7fdac0850278b56317 (diff) | |
| download | lapack-1923c69f7e9857d19b67f9e9befd1b1f37c644ef.tar.gz lapack-1923c69f7e9857d19b67f9e9befd1b1f37c644ef.tar.bz2 lapack-1923c69f7e9857d19b67f9e9befd1b1f37c644ef.zip | |
Formating code, use of * for comment, and follow other LAPACK style
Diffstat (limited to 'SRC/ctftri.f')
| -rw-r--r-- | SRC/ctftri.f | 84 |
1 files changed, 42 insertions, 42 deletions
diff --git a/SRC/ctftri.f b/SRC/ctftri.f index 0afa056b..a5978124 100644 --- a/SRC/ctftri.f +++ b/SRC/ctftri.f @@ -201,7 +201,7 @@ ELSE IF( .NOT.LOWER .AND. .NOT.LSAME( UPLO, 'U' ) ) THEN INFO = -2 ELSE IF( .NOT.LSAME( DIAG, 'N' ) .AND. .NOT.LSAME( DIAG, 'U' ) ) - + THEN + $ THEN INFO = -3 ELSE IF( N.LT.0 ) THEN INFO = -4 @@ -214,7 +214,7 @@ * Quick return if possible * IF( N.EQ.0 ) - + RETURN + $ RETURN * * If N is odd, set NISODD = .TRUE. * If N is even, set K = N/2 and NISODD = .FALSE. @@ -255,16 +255,16 @@ * CALL CTRTRI( 'L', DIAG, N1, A( 0 ), N, INFO ) IF( INFO.GT.0 ) - + RETURN + $ RETURN CALL CTRMM( 'R', 'L', 'N', DIAG, N2, N1, -CONE, A( 0 ), - + N, A( N1 ), N ) + $ N, A( N1 ), N ) CALL CTRTRI( 'U', DIAG, N2, A( N ), N, INFO ) IF( INFO.GT.0 ) - + INFO = INFO + N1 + $ INFO = INFO + N1 IF( INFO.GT.0 ) - + RETURN + $ RETURN CALL CTRMM( 'L', 'U', 'C', DIAG, N2, N1, CONE, A( N ), N, - + A( N1 ), N ) + $ A( N1 ), N ) * ELSE * @@ -274,16 +274,16 @@ * CALL CTRTRI( 'L', DIAG, N1, A( N2 ), N, INFO ) IF( INFO.GT.0 ) - + RETURN + $ RETURN CALL CTRMM( 'L', 'L', 'C', DIAG, N1, N2, -CONE, A( N2 ), - + N, A( 0 ), N ) + $ N, A( 0 ), N ) CALL CTRTRI( 'U', DIAG, N2, A( N1 ), N, INFO ) IF( INFO.GT.0 ) - + INFO = INFO + N1 + $ INFO = INFO + N1 IF( INFO.GT.0 ) - + RETURN + $ RETURN CALL CTRMM( 'R', 'U', 'N', DIAG, N1, N2, CONE, A( N1 ), - + N, A( 0 ), N ) + $ N, A( 0 ), N ) * END IF * @@ -298,16 +298,16 @@ * CALL CTRTRI( 'U', DIAG, N1, A( 0 ), N1, INFO ) IF( INFO.GT.0 ) - + RETURN + $ RETURN CALL CTRMM( 'L', 'U', 'N', DIAG, N1, N2, -CONE, A( 0 ), - + N1, A( N1*N1 ), N1 ) + $ N1, A( N1*N1 ), N1 ) CALL CTRTRI( 'L', DIAG, N2, A( 1 ), N1, INFO ) IF( INFO.GT.0 ) - + INFO = INFO + N1 + $ INFO = INFO + N1 IF( INFO.GT.0 ) - + RETURN + $ RETURN CALL CTRMM( 'R', 'L', 'C', DIAG, N1, N2, CONE, A( 1 ), - + N1, A( N1*N1 ), N1 ) + $ N1, A( N1*N1 ), N1 ) * ELSE * @@ -316,16 +316,16 @@ * CALL CTRTRI( 'U', DIAG, N1, A( N2*N2 ), N2, INFO ) IF( INFO.GT.0 ) - + RETURN + $ RETURN CALL CTRMM( 'R', 'U', 'C', DIAG, N2, N1, -CONE, - + A( N2*N2 ), N2, A( 0 ), N2 ) + $ A( N2*N2 ), N2, A( 0 ), N2 ) CALL CTRTRI( 'L', DIAG, N2, A( N1*N2 ), N2, INFO ) IF( INFO.GT.0 ) - + INFO = INFO + N1 + $ INFO = INFO + N1 IF( INFO.GT.0 ) - + RETURN + $ RETURN CALL CTRMM( 'L', 'L', 'N', DIAG, N2, N1, CONE, - + A( N1*N2 ), N2, A( 0 ), N2 ) + $ A( N1*N2 ), N2, A( 0 ), N2 ) END IF * END IF @@ -346,16 +346,16 @@ * CALL CTRTRI( 'L', DIAG, K, A( 1 ), N+1, INFO ) IF( INFO.GT.0 ) - + RETURN + $ RETURN CALL CTRMM( 'R', 'L', 'N', DIAG, K, K, -CONE, A( 1 ), - + N+1, A( K+1 ), N+1 ) + $ N+1, A( K+1 ), N+1 ) CALL CTRTRI( 'U', DIAG, K, A( 0 ), N+1, INFO ) IF( INFO.GT.0 ) - + INFO = INFO + K + $ INFO = INFO + K IF( INFO.GT.0 ) - + RETURN + $ RETURN CALL CTRMM( 'L', 'U', 'C', DIAG, K, K, CONE, A( 0 ), N+1, - + A( K+1 ), N+1 ) + $ A( K+1 ), N+1 ) * ELSE * @@ -365,16 +365,16 @@ * CALL CTRTRI( 'L', DIAG, K, A( K+1 ), N+1, INFO ) IF( INFO.GT.0 ) - + RETURN + $ RETURN CALL CTRMM( 'L', 'L', 'C', DIAG, K, K, -CONE, A( K+1 ), - + N+1, A( 0 ), N+1 ) + $ N+1, A( 0 ), N+1 ) CALL CTRTRI( 'U', DIAG, K, A( K ), N+1, INFO ) IF( INFO.GT.0 ) - + INFO = INFO + K + $ INFO = INFO + K IF( INFO.GT.0 ) - + RETURN + $ RETURN CALL CTRMM( 'R', 'U', 'N', DIAG, K, K, CONE, A( K ), N+1, - + A( 0 ), N+1 ) + $ A( 0 ), N+1 ) END IF ELSE * @@ -388,16 +388,16 @@ * CALL CTRTRI( 'U', DIAG, K, A( K ), K, INFO ) IF( INFO.GT.0 ) - + RETURN + $ RETURN CALL CTRMM( 'L', 'U', 'N', DIAG, K, K, -CONE, A( K ), K, - + A( K*( K+1 ) ), K ) + $ A( K*( K+1 ) ), K ) CALL CTRTRI( 'L', DIAG, K, A( 0 ), K, INFO ) IF( INFO.GT.0 ) - + INFO = INFO + K + $ INFO = INFO + K IF( INFO.GT.0 ) - + RETURN + $ RETURN CALL CTRMM( 'R', 'L', 'C', DIAG, K, K, CONE, A( 0 ), K, - + A( K*( K+1 ) ), K ) + $ A( K*( K+1 ) ), K ) ELSE * * SRPA for UPPER, TRANSPOSE and N is even (see paper) @@ -406,16 +406,16 @@ * CALL CTRTRI( 'U', DIAG, K, A( K*( K+1 ) ), K, INFO ) IF( INFO.GT.0 ) - + RETURN + $ RETURN CALL CTRMM( 'R', 'U', 'C', DIAG, K, K, -CONE, - + A( K*( K+1 ) ), K, A( 0 ), K ) + $ A( K*( K+1 ) ), K, A( 0 ), K ) CALL CTRTRI( 'L', DIAG, K, A( K*K ), K, INFO ) IF( INFO.GT.0 ) - + INFO = INFO + K + $ INFO = INFO + K IF( INFO.GT.0 ) - + RETURN + $ RETURN CALL CTRMM( 'L', 'L', 'N', DIAG, K, K, CONE, A( K*K ), K, - + A( 0 ), K ) + $ A( 0 ), K ) END IF END IF END IF |