--- rpl/lapack/lapack/dlansf.f	2010/12/21 13:53:30	1.5
+++ rpl/lapack/lapack/dlansf.f	2011/07/22 07:38:07	1.6
@@ -1,9 +1,9 @@
       DOUBLE PRECISION FUNCTION DLANSF( NORM, TRANSR, UPLO, N, A, WORK )
 *
-*  -- LAPACK routine (version 3.3.0)                                    --
+*  -- LAPACK routine (version 3.3.1)                                    --
 *
 *  -- Contributed by Fred Gustavson of the IBM Watson Research Center --
-*     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..--
@@ -193,19 +193,19 @@
 *
       NOE = 1
       IF( MOD( N, 2 ).EQ.0 )
-     +   NOE = 0
+     $   NOE = 0
 *
 *     set ifm = 0 when form='T or 't' and 1 otherwise
 *
       IFM = 1
       IF( LSAME( TRANSR, 'T' ) )
-     +   IFM = 0
+     $   IFM = 0
 *
 *     set ilu = 0 when uplo='U or 'u' and 1 otherwise
 *
       ILU = 1
       IF( LSAME( UPLO, 'U' ) )
-     +   ILU = 0
+     $   ILU = 0
 *
 *     set lda = (n+1)/2 when ifm = 0
 *     set lda = n when ifm = 1 and noe = 1
@@ -265,7 +265,7 @@
             END IF
          END IF
       ELSE IF( ( LSAME( NORM, 'I' ) ) .OR. ( LSAME( NORM, 'O' ) ) .OR.
-     +         ( NORM.EQ.'1' ) ) THEN
+     $         ( NORM.EQ.'1' ) ) THEN
 *
 *        Find normI(A) ( = norm1(A), since A is symmetric).
 *
@@ -289,7 +289,7 @@
 *                    -> A(j+k,j+k)
                      WORK( J+K ) = S + AA
                      IF( I.EQ.K+K )
-     +                  GO TO 10
+     $                  GO TO 10
                      I = I + 1
                      AA = ABS( A( I+J*LDA ) )
 *                    -> A(j,j)
@@ -724,7 +724,7 @@
             ELSE
 *              A is xpose
                IF( ILU.EQ.0 ) THEN
-*                 A' is upper
+*                 A**T is upper
                   DO J = 1, K - 2
                      CALL DLASSQ( J, A( 0+( K+J )*LDA ), 1, SCALE, S )
 *                    U at A(0,k)
@@ -735,7 +735,7 @@
                   END DO
                   DO J = 0, K - 2
                      CALL DLASSQ( K-J-1, A( J+1+( J+K-1 )*LDA ), 1,
-     +                            SCALE, S )
+     $                            SCALE, S )
 *                    L at A(0,k-1)
                   END DO
                   S = S + S
@@ -745,7 +745,7 @@
                   CALL DLASSQ( K, A( 0+( K-1 )*LDA ), LDA+1, SCALE, S )
 *                 tri L at A(0,k-1)
                ELSE
-*                 A' is lower
+*                 A**T is lower
                   DO J = 1, K - 1
                      CALL DLASSQ( J, A( 0+J*LDA ), 1, SCALE, S )
 *                    U at A(0,0)
@@ -806,7 +806,7 @@
             ELSE
 *              A is xpose
                IF( ILU.EQ.0 ) THEN
-*                 A' is upper
+*                 A**T is upper
                   DO J = 1, K - 1
                      CALL DLASSQ( J, A( 0+( K+1+J )*LDA ), 1, SCALE, S )
 *                    U at A(0,k+1)
@@ -817,7 +817,7 @@
                   END DO
                   DO J = 0, K - 2
                      CALL DLASSQ( K-J-1, A( J+1+( J+K )*LDA ), 1, SCALE,
-     +                            S )
+     $                            S )
 *                    L at A(0,k)
                   END DO
                   S = S + S
@@ -827,7 +827,7 @@
                   CALL DLASSQ( K, A( 0+K*LDA ), LDA+1, SCALE, S )
 *                 tri L at A(0,k)
                ELSE
-*                 A' is lower
+*                 A**T is lower
                   DO J = 1, K - 1
                      CALL DLASSQ( J, A( 0+( J+1 )*LDA ), 1, SCALE, S )
 *                    U at A(0,1)