--- rpl/lapack/lapack/dsposv.f	2010/12/21 13:53:38	1.5
+++ rpl/lapack/lapack/dsposv.f	2011/07/22 07:38:11	1.6
@@ -1,9 +1,9 @@
       SUBROUTINE DSPOSV( UPLO, N, NRHS, A, LDA, B, LDB, X, LDX, WORK,
-     +                   SWORK, ITER, INFO )
+     $                   SWORK, ITER, INFO )
 *
-*  -- LAPACK PROTOTYPE driver routine (version 3.3.0) --
+*  -- LAPACK PROTOTYPE driver routine (version 3.3.1) --
 *     Univ. of Tennessee, Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..
-*     November 2010
+*  -- April 2011                                                      --
 *
 *     ..
 *     .. Scalar Arguments ..
@@ -13,7 +13,7 @@
 *     .. Array Arguments ..
       REAL               SWORK( * )
       DOUBLE PRECISION   A( LDA, * ), B( LDB, * ), WORK( N, * ),
-     +                   X( LDX, * )
+     $                   X( LDX, * )
 *     ..
 *
 *  Purpose
@@ -126,7 +126,7 @@
 *                factorization could not be completed, and the solution
 *                has not been computed.
 *
-*  =========
+*  =====================================================================
 *
 *     .. Parameters ..
       LOGICAL            DOITREF
@@ -147,7 +147,7 @@
 *
 *     .. External Subroutines ..
       EXTERNAL           DAXPY, DSYMM, DLACPY, DLAT2S, DLAG2S, SLAG2D,
-     +                   SPOTRF, SPOTRS, XERBLA
+     $                   SPOTRF, SPOTRS, XERBLA
 *     ..
 *     .. External Functions ..
       INTEGER            IDAMAX
@@ -186,7 +186,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.
@@ -239,7 +239,7 @@
 *     Solve the system SA*SX = SB.
 *
       CALL SPOTRS( UPLO, N, NRHS, SWORK( PTSA ), N, SWORK( PTSX ), N,
-     +             INFO )
+     $             INFO )
 *
 *     Convert SX back to double precision
 *
@@ -250,7 +250,7 @@
       CALL DLACPY( 'All', N, NRHS, B, LDB, WORK, N )
 *
       CALL DSYMM( '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.
@@ -259,7 +259,7 @@
          XNRM = ABS( X( IDAMAX( N, X( 1, I ), 1 ), I ) )
          RNRM = ABS( WORK( IDAMAX( 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
@@ -285,7 +285,7 @@
 *        Solve the system SA*SX = SR.
 *
          CALL SPOTRS( UPLO, N, NRHS, SWORK( PTSA ), N, SWORK( PTSX ), N,
-     +                INFO )
+     $                INFO )
 *
 *        Convert SX back to double precision and update the current
 *        iterate.
@@ -301,7 +301,7 @@
          CALL DLACPY( 'All', N, NRHS, B, LDB, WORK, N )
 *
          CALL DSYMM( '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.
@@ -310,7 +310,7 @@
             XNRM = ABS( X( IDAMAX( N, X( 1, I ), 1 ), I ) )
             RNRM = ABS( WORK( IDAMAX( 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
@@ -339,7 +339,7 @@
       CALL DPOTRF( UPLO, N, A, LDA, INFO )
 *
       IF( INFO.NE.0 )
-     +   RETURN
+     $   RETURN
 *
       CALL DLACPY( 'All', N, NRHS, B, LDB, X, LDX )
       CALL DPOTRS( UPLO, N, NRHS, A, LDA, X, LDX, INFO )