--- rpl/lapack/lapack/zpttrs.f	2010/12/21 13:53:55	1.7
+++ rpl/lapack/lapack/zpttrs.f	2011/07/22 07:38:20	1.8
@@ -1,9 +1,9 @@
       SUBROUTINE ZPTTRS( UPLO, N, NRHS, D, E, 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
@@ -19,7 +19,7 @@
 *
 *  ZPTTRS solves a tridiagonal system of the form
 *     A * X = B
-*  using the factorization A = U'*D*U or A = L*D*L' computed by ZPTTRF.
+*  using the factorization A = U**H *D* U or A = L*D*L**H computed by ZPTTRF.
 *  D is a diagonal matrix specified in the vector D, U (or L) is a unit
 *  bidiagonal matrix whose superdiagonal (subdiagonal) is specified in
 *  the vector E, and X and B are N by NRHS matrices.
@@ -31,8 +31,8 @@
 *          Specifies the form of the factorization and whether the
 *          vector E is the superdiagonal of the upper bidiagonal factor
 *          U or the subdiagonal of the lower bidiagonal factor L.
-*          = 'U':  A = U'*D*U, E is the superdiagonal of U
-*          = 'L':  A = L*D*L', E is the subdiagonal of L
+*          = 'U':  A = U**H *D*U, E is the superdiagonal of U
+*          = 'L':  A = L*D*L**H, E is the subdiagonal of L
 *
 *  N       (input) INTEGER
 *          The order of the tridiagonal matrix A.  N >= 0.
@@ -43,13 +43,13 @@
 *
 *  D       (input) DOUBLE PRECISION array, dimension (N)
 *          The n diagonal elements of the diagonal matrix D from the
-*          factorization A = U'*D*U or A = L*D*L'.
+*          factorization A = U**H *D*U or A = L*D*L**H.
 *
 *  E       (input) COMPLEX*16 array, dimension (N-1)
 *          If UPLO = 'U', the (n-1) superdiagonal elements of the unit
-*          bidiagonal factor U from the factorization A = U'*D*U.
+*          bidiagonal factor U from the factorization A = U**H*D*U.
 *          If UPLO = 'L', the (n-1) subdiagonal elements of the unit
-*          bidiagonal factor L from the factorization A = L*D*L'.
+*          bidiagonal factor L from the factorization A = L*D*L**H.
 *
 *  B       (input/output) DOUBLE PRECISION array, dimension (LDB,NRHS)
 *          On entry, the right hand side vectors B for the system of