--- rpl/lapack/blas/dgemm.f	2010/12/21 13:51:24	1.6
+++ rpl/lapack/blas/dgemm.f	2011/07/22 07:38:01	1.7
@@ -17,7 +17,7 @@
 *
 *  where  op( X ) is one of
 *
-*     op( X ) = X   or   op( X ) = X',
+*     op( X ) = X   or   op( X ) = X**T,
 *
 *  alpha and beta are scalars, and A, B and C are matrices, with op( A )
 *  an m by k matrix,  op( B )  a  k by n matrix and  C an m by n matrix.
@@ -31,9 +31,9 @@
 *
 *              TRANSA = 'N' or 'n',  op( A ) = A.
 *
-*              TRANSA = 'T' or 't',  op( A ) = A'.
+*              TRANSA = 'T' or 't',  op( A ) = A**T.
 *
-*              TRANSA = 'C' or 'c',  op( A ) = A'.
+*              TRANSA = 'C' or 'c',  op( A ) = A**T.
 *
 *           Unchanged on exit.
 *
@@ -43,9 +43,9 @@
 *
 *              TRANSB = 'N' or 'n',  op( B ) = B.
 *
-*              TRANSB = 'T' or 't',  op( B ) = B'.
+*              TRANSB = 'T' or 't',  op( B ) = B**T.
 *
-*              TRANSB = 'C' or 'c',  op( B ) = B'.
+*              TRANSB = 'C' or 'c',  op( B ) = B**T.
 *
 *           Unchanged on exit.
 *
@@ -249,7 +249,7 @@
    90         CONTINUE
           ELSE
 *
-*           Form  C := alpha*A'*B + beta*C
+*           Form  C := alpha*A**T*B + beta*C
 *
               DO 120 J = 1,N
                   DO 110 I = 1,M
@@ -268,7 +268,7 @@
       ELSE
           IF (NOTA) THEN
 *
-*           Form  C := alpha*A*B' + beta*C
+*           Form  C := alpha*A*B**T + beta*C
 *
               DO 170 J = 1,N
                   IF (BETA.EQ.ZERO) THEN
@@ -291,7 +291,7 @@
   170         CONTINUE
           ELSE
 *
-*           Form  C := alpha*A'*B' + beta*C
+*           Form  C := alpha*A**T*B**T + beta*C
 *
               DO 200 J = 1,N
                   DO 190 I = 1,M