--- rpl/lapack/lapack/dla_gbamv.f	2010/12/21 13:53:27	1.4
+++ rpl/lapack/lapack/dla_gbamv.f	2011/07/22 07:38:06	1.5
@@ -1,7 +1,7 @@
       SUBROUTINE DLA_GBAMV( TRANS, M, N, KL, KU, ALPHA, AB, LDAB, X,
      $                      INCX, BETA, Y, INCY )
 *
-*     -- LAPACK routine (version 3.2.2)                                 --
+*     -- LAPACK routine (version 3.3.1)                                 --
 *     -- Contributed by James Demmel, Deaglan Halligan, Yozo Hida and --
 *     -- Jason Riedy of Univ. of California Berkeley.                 --
 *     -- June 2010                                                    --
@@ -25,7 +25,7 @@
 *  DLA_GBAMV  performs one of the matrix-vector operations
 *
 *          y := alpha*abs(A)*abs(x) + beta*abs(y),
-*     or   y := alpha*abs(A)'*abs(x) + beta*abs(y),
+*     or   y := alpha*abs(A)**T*abs(x) + beta*abs(y),
 *
 *  where alpha and beta are scalars, x and y are vectors and A is an
 *  m by n matrix.
@@ -47,43 +47,43 @@
 *           follows:
 *
 *             BLAS_NO_TRANS      y := alpha*abs(A)*abs(x) + beta*abs(y)
-*             BLAS_TRANS         y := alpha*abs(A')*abs(x) + beta*abs(y)
-*             BLAS_CONJ_TRANS    y := alpha*abs(A')*abs(x) + beta*abs(y)
+*             BLAS_TRANS         y := alpha*abs(A**T)*abs(x) + beta*abs(y)
+*             BLAS_CONJ_TRANS    y := alpha*abs(A**T)*abs(x) + beta*abs(y)
 *
 *           Unchanged on exit.
 *
-*  M       (input) INTEGER
+*  M        (input) INTEGER
 *           On entry, M specifies the number of rows of the matrix A.
 *           M must be at least zero.
 *           Unchanged on exit.
 *
-*  N       (input) INTEGER
+*  N        (input) INTEGER
 *           On entry, N specifies the number of columns of the matrix A.
 *           N must be at least zero.
 *           Unchanged on exit.
 *
-*  KL      (input) INTEGER
+*  KL       (input) INTEGER
 *           The number of subdiagonals within the band of A.  KL >= 0.
 *
-*  KU      (input) INTEGER
+*  KU       (input) INTEGER
 *           The number of superdiagonals within the band of A.  KU >= 0.
 *
-*  ALPHA  - DOUBLE PRECISION
+*  ALPHA    (input) DOUBLE PRECISION
 *           On entry, ALPHA specifies the scalar alpha.
 *           Unchanged on exit.
 *
-*  A      - DOUBLE PRECISION   array of DIMENSION ( LDA, n )
-*           Before entry, the leading m by n part of the array A must
+*  AB       (input) DOUBLE PRECISION array of DIMENSION ( LDAB, n )
+*           Before entry, the leading m by n part of the array AB must
 *           contain the matrix of coefficients.
 *           Unchanged on exit.
 *
-*  LDA     (input) INTEGER
-*           On entry, LDA specifies the first dimension of A as declared
-*           in the calling (sub) program. LDA must be at least
+*  LDAB     (input) INTEGER
+*           On entry, LDA specifies the first dimension of AB as declared
+*           in the calling (sub) program. LDAB must be at least
 *           max( 1, m ).
 *           Unchanged on exit.
 *
-*  X       (input) DOUBLE PRECISION array, dimension
+*  X        (input) DOUBLE PRECISION array, dimension
 *           ( 1 + ( n - 1 )*abs( INCX ) ) when TRANS = 'N' or 'n'
 *           and at least
 *           ( 1 + ( m - 1 )*abs( INCX ) ) otherwise.
@@ -91,17 +91,17 @@
 *           vector x.
 *           Unchanged on exit.
 *
-*  INCX    (input) INTEGER
+*  INCX     (input) INTEGER
 *           On entry, INCX specifies the increment for the elements of
 *           X. INCX must not be zero.
 *           Unchanged on exit.
 *
-*  BETA   - DOUBLE PRECISION
+*  BETA     (input)  DOUBLE PRECISION
 *           On entry, BETA specifies the scalar beta. When BETA is
 *           supplied as zero then Y need not be set on input.
 *           Unchanged on exit.
 *
-*  Y       (input/output) DOUBLE PRECISION  array, dimension
+*  Y        (input/output) DOUBLE PRECISION  array, dimension
 *           ( 1 + ( m - 1 )*abs( INCY ) ) when TRANS = 'N' or 'n'
 *           and at least
 *           ( 1 + ( n - 1 )*abs( INCY ) ) otherwise.
@@ -109,7 +109,7 @@
 *           must contain the vector y. On exit, Y is overwritten by the
 *           updated vector y.
 *
-*  INCY    (input) INTEGER
+*  INCY     (input) INTEGER
 *           On entry, INCY specifies the increment for the elements of
 *           Y. INCY must not be zero.
 *           Unchanged on exit.
@@ -118,7 +118,7 @@
 *  Level 2 Blas routine.
 *
 *  =====================================================================
-
+*
 *     .. Parameters ..
       DOUBLE PRECISION   ONE, ZERO
       PARAMETER          ( ONE = 1.0D+0, ZERO = 0.0D+0 )