--- rpl/lapack/lapack/dsytrs2.f	2010/12/21 13:53:40	1.2
+++ rpl/lapack/lapack/dsytrs2.f	2011/07/22 07:38:12	1.3
@@ -1,10 +1,10 @@
       SUBROUTINE DSYTRS2( UPLO, N, NRHS, A, LDA, IPIV, B, LDB, 
      $                    WORK, INFO )
 *
-*  -- LAPACK PROTOTYPE routine (version 3.3.0) --
+*  -- LAPACK PROTOTYPE 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..--
-*     November 2010
+*  -- April 2011                                                      --
 *
 *  -- Written by Julie Langou of the Univ. of TN    --
 *
@@ -116,9 +116,9 @@
 *
       IF( UPPER ) THEN
 *
-*        Solve A*X = B, where A = U*D*U'.
+*        Solve A*X = B, where A = U*D*U**T.
 *
-*       P' * B  
+*       P**T * B  
         K=N
         DO WHILE ( K .GE. 1 )
          IF( IPIV( K ).GT.0 ) THEN
@@ -138,16 +138,16 @@
          END IF
         END DO
 *
-*  Compute (U \P' * B) -> B    [ (U \P' * B) ]
+*  Compute (U \P**T * B) -> B    [ (U \P**T * B) ]
 *
-        CALL DTRSM('L','U','N','U',N,NRHS,ONE,A,N,B,N)
+        CALL DTRSM('L','U','N','U',N,NRHS,ONE,A,LDA,B,LDB)
 *
-*  Compute D \ B -> B   [ D \ (U \P' * B) ]
+*  Compute D \ B -> B   [ D \ (U \P**T * B) ]
 *       
          I=N
          DO WHILE ( I .GE. 1 )
             IF( IPIV(I) .GT. 0 ) THEN
-              CALL DSCAL( NRHS, ONE / A( I, I ), B( I, 1 ), N )
+              CALL DSCAL( NRHS, ONE / A( I, I ), B( I, 1 ), LDB )
             ELSEIF ( I .GT. 1) THEN
                IF ( IPIV(I-1) .EQ. IPIV(I) ) THEN
                   AKM1K = WORK(I)
@@ -166,11 +166,11 @@
             I = I - 1
          END DO
 *
-*      Compute (U' \ B) -> B   [ U' \ (D \ (U \P' * B) ) ]
+*      Compute (U**T \ B) -> B   [ U**T \ (D \ (U \P**T * B) ) ]
 *
-         CALL DTRSM('L','U','T','U',N,NRHS,ONE,A,N,B,N)
+         CALL DTRSM('L','U','T','U',N,NRHS,ONE,A,LDA,B,LDB)
 *
-*       P * B  [ P * (U' \ (D \ (U \P' * B) )) ]
+*       P * B  [ P * (U**T \ (D \ (U \P**T * B) )) ]
 *
         K=1
         DO WHILE ( K .LE. N )
@@ -193,9 +193,9 @@
 *
       ELSE
 *
-*        Solve A*X = B, where A = L*D*L'.
+*        Solve A*X = B, where A = L*D*L**T.
 *
-*       P' * B  
+*       P**T * B  
         K=1
         DO WHILE ( K .LE. N )
          IF( IPIV( K ).GT.0 ) THEN
@@ -215,16 +215,16 @@
          ENDIF
         END DO
 *
-*  Compute (L \P' * B) -> B    [ (L \P' * B) ]
+*  Compute (L \P**T * B) -> B    [ (L \P**T * B) ]
 *
-        CALL DTRSM('L','L','N','U',N,NRHS,ONE,A,N,B,N)
+        CALL DTRSM('L','L','N','U',N,NRHS,ONE,A,LDA,B,LDB)
 *
-*  Compute D \ B -> B   [ D \ (L \P' * B) ]
+*  Compute D \ B -> B   [ D \ (L \P**T * B) ]
 *       
          I=1
          DO WHILE ( I .LE. N )
             IF( IPIV(I) .GT. 0 ) THEN
-              CALL DSCAL( NRHS, ONE / A( I, I ), B( I, 1 ), N )
+              CALL DSCAL( NRHS, ONE / A( I, I ), B( I, 1 ), LDB )
             ELSE
                   AKM1K = WORK(I)
                   AKM1 = A( I, I ) / AKM1K
@@ -241,11 +241,11 @@
             I = I + 1
          END DO
 *
-*  Compute (L' \ B) -> B   [ L' \ (D \ (L \P' * B) ) ]
+*  Compute (L**T \ B) -> B   [ L**T \ (D \ (L \P**T * B) ) ]
 * 
-        CALL DTRSM('L','L','T','U',N,NRHS,ONE,A,N,B,N)
+        CALL DTRSM('L','L','T','U',N,NRHS,ONE,A,LDA,B,LDB)
 *
-*       P * B  [ P * (L' \ (D \ (L \P' * B) )) ]
+*       P * B  [ P * (L**T \ (D \ (L \P**T * B) )) ]
 *
         K=N
         DO WHILE ( K .GE. 1 )