--- rpl/lapack/lapack/zcposv.f	2010/12/21 13:53:42	1.5
+++ rpl/lapack/lapack/zcposv.f	2011/07/22 07:38:13	1.6
@@ -1,9 +1,9 @@
       SUBROUTINE ZCPOSV( UPLO, N, NRHS, A, LDA, B, LDB, X, LDX, WORK,
-     +                   SWORK, RWORK, ITER, INFO )
+     $                   SWORK, RWORK, ITER, INFO )
 *
-*  -- LAPACK PROTOTYPE driver routine (version 3.3.0)                 --
+*  -- LAPACK PROTOTYPE driver routine (version 3.3.1)                 --
 *
-*     November 2010
+*  -- April 2011                                                      --
 *
 *  -- LAPACK is a software package provided by Univ. of Tennessee,    --
 *  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
@@ -16,7 +16,7 @@
       DOUBLE PRECISION   RWORK( * )
       COMPLEX            SWORK( * )
       COMPLEX*16         A( LDA, * ), B( LDB, * ), WORK( N, * ),
-     +                   X( LDX, * )
+     $                   X( LDX, * )
 *     ..
 *
 *  Purpose
@@ -134,7 +134,7 @@
 *                factorization could not be completed, and the solution
 *                has not been computed.
 *
-*  =========
+*  =====================================================================
 *
 *     .. Parameters ..
       LOGICAL            DOITREF
@@ -148,7 +148,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
@@ -157,7 +157,7 @@
 *
 *     .. External Subroutines ..
       EXTERNAL           ZAXPY, ZHEMM, ZLACPY, ZLAT2C, ZLAG2C, CLAG2Z,
-     +                   CPOTRF, CPOTRS, XERBLA
+     $                   CPOTRF, CPOTRS, XERBLA
 *     ..
 *     .. External Functions ..
       INTEGER            IZAMAX
@@ -201,7 +201,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.
@@ -254,7 +254,7 @@
 *     Solve the system SA*SX = SB.
 *
       CALL CPOTRS( UPLO, N, NRHS, SWORK( PTSA ), N, SWORK( PTSX ), N,
-     +             INFO )
+     $             INFO )
 *
 *     Convert SX back to COMPLEX*16
 *
@@ -265,7 +265,7 @@
       CALL ZLACPY( 'All', N, NRHS, B, LDB, WORK, N )
 *
       CALL ZHEMM( 'Left', UPLO, N, NRHS, NEGONE, A, LDA, X, LDX, ONE,
-     +            WORK, N )
+     $            WORK, N )
 *
 *     Check whether the NRHS normwise backward errors satisfy the
 *     stopping criterion. If yes, set ITER=0 and return.
@@ -274,7 +274,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
@@ -300,7 +300,7 @@
 *        Solve the system SA*SX = SR.
 *
          CALL CPOTRS( UPLO, N, NRHS, SWORK( PTSA ), N, SWORK( PTSX ), N,
-     +                INFO )
+     $                INFO )
 *
 *        Convert SX back to double precision and update the current
 *        iterate.
@@ -316,7 +316,7 @@
          CALL ZLACPY( 'All', N, NRHS, B, LDB, WORK, N )
 *
          CALL ZHEMM( 'L', UPLO, N, NRHS, NEGONE, A, LDA, X, LDX, ONE,
-     +               WORK, N )
+     $               WORK, N )
 *
 *        Check whether the NRHS normwise backward errors satisfy the
 *        stopping criterion. If yes, set ITER=IITER>0 and return.
@@ -325,7 +325,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
@@ -354,7 +354,7 @@
       CALL ZPOTRF( UPLO, N, A, LDA, INFO )
 *
       IF( INFO.NE.0 )
-     +   RETURN
+     $   RETURN
 *
       CALL ZLACPY( 'All', N, NRHS, B, LDB, X, LDX )
       CALL ZPOTRS( UPLO, N, NRHS, A, LDA, X, LDX, INFO )