--- rpl/lapack/lapack/zcgesv.f	2010/12/21 13:53:42	1.8
+++ rpl/lapack/lapack/zcgesv.f	2011/07/22 07:38:13	1.9
@@ -1,10 +1,10 @@
       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 PROTOTYPE driver routine (version 3.3.1) --
 *  -- LAPACK is a software package provided by Univ. of Tennessee,    --
 *  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
-*     January 2007
+*  -- April 2011                                                      --
 *
 *     ..
 *     .. Scalar Arguments ..
@@ -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
 *