--- rpl/lapack/lapack/dsygs2.f	2010/12/21 13:53:39	1.7
+++ rpl/lapack/lapack/dsygs2.f	2011/07/22 07:38:11	1.8
@@ -1,9 +1,9 @@
       SUBROUTINE DSYGS2( ITYPE, UPLO, N, A, LDA, B, LDB, INFO )
 *
-*  -- LAPACK routine (version 3.2) --
+*  -- LAPACK 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 2006
+*  -- April 2011                                                      --
 *
 *     .. Scalar Arguments ..
       CHARACTER          UPLO
@@ -20,19 +20,19 @@
 *  to standard form.
 *
 *  If ITYPE = 1, the problem is A*x = lambda*B*x,
-*  and A is overwritten by inv(U')*A*inv(U) or inv(L)*A*inv(L')
+*  and A is overwritten by inv(U**T)*A*inv(U) or inv(L)*A*inv(L**T)
 *
 *  If ITYPE = 2 or 3, the problem is A*B*x = lambda*x or
-*  B*A*x = lambda*x, and A is overwritten by U*A*U` or L'*A*L.
+*  B*A*x = lambda*x, and A is overwritten by U*A*U**T or L**T *A*L.
 *
-*  B must have been previously factorized as U'*U or L*L' by DPOTRF.
+*  B must have been previously factorized as U**T *U or L*L**T by DPOTRF.
 *
 *  Arguments
 *  =========
 *
 *  ITYPE   (input) INTEGER
-*          = 1: compute inv(U')*A*inv(U) or inv(L)*A*inv(L');
-*          = 2 or 3: compute U*A*U' or L'*A*L.
+*          = 1: compute inv(U**T)*A*inv(U) or inv(L)*A*inv(L**T);
+*          = 2 or 3: compute U*A*U**T or L**T *A*L.
 *
 *  UPLO    (input) CHARACTER*1
 *          Specifies whether the upper or lower triangular part of the
@@ -115,7 +115,7 @@
       IF( ITYPE.EQ.1 ) THEN
          IF( UPPER ) THEN
 *
-*           Compute inv(U')*A*inv(U)
+*           Compute inv(U**T)*A*inv(U)
 *
             DO 10 K = 1, N
 *
@@ -140,7 +140,7 @@
    10       CONTINUE
          ELSE
 *
-*           Compute inv(L)*A*inv(L')
+*           Compute inv(L)*A*inv(L**T)
 *
             DO 20 K = 1, N
 *
@@ -165,7 +165,7 @@
       ELSE
          IF( UPPER ) THEN
 *
-*           Compute U*A*U'
+*           Compute U*A*U**T
 *
             DO 30 K = 1, N
 *
@@ -185,7 +185,7 @@
    30       CONTINUE
          ELSE
 *
-*           Compute L'*A*L
+*           Compute L**T *A*L
 *
             DO 40 K = 1, N
 *