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/zcgesv.f | |
parent | 86579dddaf0142000a291a7fdac0850278b56317 (diff) |
Formating code, use of * for comment, and follow other LAPACK style
Diffstat (limited to 'SRC/zcgesv.f')
-rw-r--r-- | SRC/zcgesv.f | 28 |
1 files changed, 14 insertions, 14 deletions
diff --git a/SRC/zcgesv.f b/SRC/zcgesv.f index ebf98aae..71611f90 100644 --- a/SRC/zcgesv.f +++ b/SRC/zcgesv.f @@ -1,5 +1,5 @@ SUBROUTINE ZCGESV( N, NRHS, A, LDA, IPIV, B, LDB, X, LDX, WORK, - + SWORK, RWORK, ITER, INFO ) + $ SWORK, RWORK, ITER, INFO ) * * -- LAPACK PROTOTYPE driver routine (version 3.2.2) -- * -- LAPACK is a software package provided by Univ. of Tennessee, -- @@ -15,7 +15,7 @@ DOUBLE PRECISION RWORK( * ) COMPLEX SWORK( * ) COMPLEX*16 A( LDA, * ), B( LDB, * ), WORK( N, * ), - + X( LDX, * ) + $ X( LDX, * ) * .. * * Purpose @@ -125,7 +125,7 @@ * factor U is exactly singular, so the solution * could not be computed. * -* ========= +* ===================================================================== * * .. Parameters .. LOGICAL DOITREF @@ -139,7 +139,7 @@ * COMPLEX*16 NEGONE, ONE PARAMETER ( NEGONE = ( -1.0D+00, 0.0D+00 ), - + ONE = ( 1.0D+00, 0.0D+00 ) ) + $ ONE = ( 1.0D+00, 0.0D+00 ) ) * * .. Local Scalars .. INTEGER I, IITER, PTSA, PTSX @@ -148,7 +148,7 @@ * * .. External Subroutines .. EXTERNAL CGETRS, CGETRF, CLAG2Z, XERBLA, ZAXPY, ZGEMM, - + ZLACPY, ZLAG2C + $ ZLACPY, ZLAG2C * .. * .. External Functions .. INTEGER IZAMAX @@ -190,7 +190,7 @@ * Quick return if (N.EQ.0). * IF( N.EQ.0 ) - + RETURN + $ RETURN * * Skip single precision iterative refinement if a priori slower * than double precision factorization. @@ -243,7 +243,7 @@ * Solve the system SA*SX = SB. * CALL CGETRS( 'No transpose', N, NRHS, SWORK( PTSA ), N, IPIV, - + SWORK( PTSX ), N, INFO ) + $ SWORK( PTSX ), N, INFO ) * * Convert SX back to double precision * @@ -254,7 +254,7 @@ CALL ZLACPY( 'All', N, NRHS, B, LDB, WORK, N ) * CALL ZGEMM( 'No Transpose', 'No Transpose', N, NRHS, N, NEGONE, A, - + LDA, X, LDX, ONE, WORK, N ) + $ LDA, X, LDX, ONE, WORK, N ) * * Check whether the NRHS normwise backward errors satisfy the * stopping criterion. If yes, set ITER=0 and return. @@ -263,7 +263,7 @@ XNRM = CABS1( X( IZAMAX( N, X( 1, I ), 1 ), I ) ) RNRM = CABS1( WORK( IZAMAX( N, WORK( 1, I ), 1 ), I ) ) IF( RNRM.GT.XNRM*CTE ) - + GO TO 10 + $ GO TO 10 END DO * * If we are here, the NRHS normwise backward errors satisfy the @@ -289,7 +289,7 @@ * Solve the system SA*SX = SR. * CALL CGETRS( 'No transpose', N, NRHS, SWORK( PTSA ), N, IPIV, - + SWORK( PTSX ), N, INFO ) + $ SWORK( PTSX ), N, INFO ) * * Convert SX back to double precision and update the current * iterate. @@ -305,7 +305,7 @@ CALL ZLACPY( 'All', N, NRHS, B, LDB, WORK, N ) * CALL ZGEMM( 'No Transpose', 'No Transpose', N, NRHS, N, NEGONE, - + A, LDA, X, LDX, ONE, WORK, N ) + $ A, LDA, X, LDX, ONE, WORK, N ) * * Check whether the NRHS normwise backward errors satisfy the * stopping criterion. If yes, set ITER=IITER>0 and return. @@ -314,7 +314,7 @@ XNRM = CABS1( X( IZAMAX( N, X( 1, I ), 1 ), I ) ) RNRM = CABS1( WORK( IZAMAX( N, WORK( 1, I ), 1 ), I ) ) IF( RNRM.GT.XNRM*CTE ) - + GO TO 20 + $ GO TO 20 END DO * * If we are here, the NRHS normwise backward errors satisfy the @@ -343,11 +343,11 @@ CALL ZGETRF( N, N, A, LDA, IPIV, INFO ) * IF( INFO.NE.0 ) - + RETURN + $ RETURN * CALL ZLACPY( 'All', N, NRHS, B, LDB, X, LDX ) CALL ZGETRS( 'No transpose', N, NRHS, A, LDA, IPIV, X, LDX, - + INFO ) + $ INFO ) * RETURN * |