First pass for BLAS to homgenize notation for transpose (**T) and conjugate transpose (**H)

Corresponds to bug0024

Author: julielangou

BLAS/SRC/cgbmv.f

@@ -13,9 +13,9 @@ SUBROUTINE CGBMV(TRANS,M,N,KL,KU,ALPHA,A,LDA,X,INCX,BETA,Y,INCY)
 *
 *  CGBMV  performs one of the matrix-vector operations
 *
-*     y := alpha*A*x + beta*y,   or   y := alpha*A'*x + beta*y,   or
+*     y := alpha*A*x + beta*y,   or   y := alpha*A**T*x + beta*y,   or
 *
-*     y := alpha*conjg( A' )*x + beta*y,
+*     y := alpha*A**H*x + beta*y,
 *
 *  where alpha and beta are scalars, x and y are vectors and A is an
 *  m by n band matrix, with kl sub-diagonals and ku super-diagonals.
@@ -29,9 +29,9 @@ SUBROUTINE CGBMV(TRANS,M,N,KL,KU,ALPHA,A,LDA,X,INCX,BETA,Y,INCY)
 *
 *              TRANS = 'N' or 'n'   y := alpha*A*x + beta*y.
 *
-*              TRANS = 'T' or 't'   y := alpha*A'*x + beta*y.
+*              TRANS = 'T' or 't'   y := alpha*A**T*x + beta*y.
 *
-*              TRANS = 'C' or 'c'   y := alpha*conjg( A' )*x + beta*y.
+*              TRANS = 'C' or 'c'   y := alpha*A**H*x + beta*y.
 *
 *           Unchanged on exit.
 *
@@ -274,7 +274,7 @@ SUBROUTINE CGBMV(TRANS,M,N,KL,KU,ALPHA,A,LDA,X,INCX,BETA,Y,INCY)
           END IF
       ELSE
 *
-*        Form  y := alpha*A'*x + y  or  y := alpha*conjg( A' )*x + y.
+*        Form  y := alpha*A**T*x + y  or  y := alpha*A**H*x + y.
 *
           JY = KY
           IF (INCX.EQ.1) THEN

BLAS/SRC/cgemm.f

@@ -17,7 +17,7 @@ SUBROUTINE CGEMM(TRANSA,TRANSB,M,N,K,ALPHA,A,LDA,B,LDB,BETA,C,LDC)
 *
 *  where  op( X ) is one of
 *
-*     op( X ) = X   or   op( X ) = X'   or   op( X ) = conjg( X' ),
+*     op( X ) = X   or   op( X ) = X**T   or   op( X ) = X**H,
 *
 *  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 @@ SUBROUTINE CGEMM(TRANSA,TRANSB,M,N,K,ALPHA,A,LDA,B,LDB,BETA,C,LDC)
 *
 *              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 ) = conjg( A' ).
+*              TRANSA = 'C' or 'c',  op( A ) = A**H.
 *
 *           Unchanged on exit.
 *
@@ -43,9 +43,9 @@ SUBROUTINE CGEMM(TRANSA,TRANSB,M,N,K,ALPHA,A,LDA,B,LDB,BETA,C,LDC)
 *
 *              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 ) = conjg( B' ).
+*              TRANSB = 'C' or 'c',  op( B ) = B**H.
 *
 *           Unchanged on exit.
 *
@@ -255,7 +255,7 @@ SUBROUTINE CGEMM(TRANSA,TRANSB,M,N,K,ALPHA,A,LDA,B,LDB,BETA,C,LDC)
    90         CONTINUE
           ELSE IF (CONJA) THEN
 *
-*           Form  C := alpha*conjg( A' )*B + beta*C.
+*           Form  C := alpha*A**H*B + beta*C.
 *
               DO 120 J = 1,N
                   DO 110 I = 1,M
@@ -272,7 +272,7 @@ SUBROUTINE CGEMM(TRANSA,TRANSB,M,N,K,ALPHA,A,LDA,B,LDB,BETA,C,LDC)
   120         CONTINUE
           ELSE
 *
-*           Form  C := alpha*A'*B + beta*C
+*           Form  C := alpha*A**T*B + beta*C
 *
               DO 150 J = 1,N
                   DO 140 I = 1,M
@@ -291,7 +291,7 @@ SUBROUTINE CGEMM(TRANSA,TRANSB,M,N,K,ALPHA,A,LDA,B,LDB,BETA,C,LDC)
       ELSE IF (NOTA) THEN
           IF (CONJB) THEN
 *
-*           Form  C := alpha*A*conjg( B' ) + beta*C.
+*           Form  C := alpha*A*B**H + beta*C.
 *
               DO 200 J = 1,N
                   IF (BETA.EQ.ZERO) THEN
@@ -314,7 +314,7 @@ SUBROUTINE CGEMM(TRANSA,TRANSB,M,N,K,ALPHA,A,LDA,B,LDB,BETA,C,LDC)
   200         CONTINUE
           ELSE
 *
-*           Form  C := alpha*A*B'          + beta*C
+*           Form  C := alpha*A*B**T          + beta*C
 *
               DO 250 J = 1,N
                   IF (BETA.EQ.ZERO) THEN
@@ -339,7 +339,7 @@ SUBROUTINE CGEMM(TRANSA,TRANSB,M,N,K,ALPHA,A,LDA,B,LDB,BETA,C,LDC)
       ELSE IF (CONJA) THEN
           IF (CONJB) THEN
 *
-*           Form  C := alpha*conjg( A' )*conjg( B' ) + beta*C.
+*           Form  C := alpha*A**H*B**H + beta*C.
 *
               DO 280 J = 1,N
                   DO 270 I = 1,M
@@ -356,7 +356,7 @@ SUBROUTINE CGEMM(TRANSA,TRANSB,M,N,K,ALPHA,A,LDA,B,LDB,BETA,C,LDC)
   280         CONTINUE
           ELSE
 *
-*           Form  C := alpha*conjg( A' )*B' + beta*C
+*           Form  C := alpha*A**H*B**T + beta*C
 *
               DO 310 J = 1,N
                   DO 300 I = 1,M
@@ -375,7 +375,7 @@ SUBROUTINE CGEMM(TRANSA,TRANSB,M,N,K,ALPHA,A,LDA,B,LDB,BETA,C,LDC)
       ELSE
           IF (CONJB) THEN
 *
-*           Form  C := alpha*A'*conjg( B' ) + beta*C
+*           Form  C := alpha*A**T*B**H + beta*C
 *
               DO 340 J = 1,N
                   DO 330 I = 1,M
@@ -392,7 +392,7 @@ SUBROUTINE CGEMM(TRANSA,TRANSB,M,N,K,ALPHA,A,LDA,B,LDB,BETA,C,LDC)
   340         CONTINUE
           ELSE
 *
-*           Form  C := alpha*A'*B' + beta*C
+*           Form  C := alpha*A**T*B**T + beta*C
 *
               DO 370 J = 1,N
                   DO 360 I = 1,M

BLAS/SRC/cgemv.f

@@ -13,9 +13,9 @@ SUBROUTINE CGEMV(TRANS,M,N,ALPHA,A,LDA,X,INCX,BETA,Y,INCY)
 *
 *  CGEMV performs one of the matrix-vector operations
 *
-*     y := alpha*A*x + beta*y,   or   y := alpha*A'*x + beta*y,   or
+*     y := alpha*A*x + beta*y,   or   y := alpha*A**T*x + beta*y,   or
 *
-*     y := alpha*conjg( A' )*x + beta*y,
+*     y := alpha*A**H*x + beta*y,
 *
 *  where alpha and beta are scalars, x and y are vectors and A is an
 *  m by n matrix.
@@ -29,9 +29,9 @@ SUBROUTINE CGEMV(TRANS,M,N,ALPHA,A,LDA,X,INCX,BETA,Y,INCY)
 *
 *              TRANS = 'N' or 'n'   y := alpha*A*x + beta*y.
 *
-*              TRANS = 'T' or 't'   y := alpha*A'*x + beta*y.
+*              TRANS = 'T' or 't'   y := alpha*A**T*x + beta*y.
 *
-*              TRANS = 'C' or 'c'   y := alpha*conjg( A' )*x + beta*y.
+*              TRANS = 'C' or 'c'   y := alpha*A**H*x + beta*y.
 *
 *           Unchanged on exit.
 *
@@ -239,7 +239,7 @@ SUBROUTINE CGEMV(TRANS,M,N,ALPHA,A,LDA,X,INCX,BETA,Y,INCY)
           END IF
       ELSE
 *
-*        Form  y := alpha*A'*x + y  or  y := alpha*conjg( A' )*x + y.
+*        Form  y := alpha*A**T*x + y  or  y := alpha*A**H*x + y.
 *
           JY = KY
           IF (INCX.EQ.1) THEN

BLAS/SRC/cgerc.f

@@ -12,7 +12,7 @@ SUBROUTINE CGERC(M,N,ALPHA,X,INCX,Y,INCY,A,LDA)
 *
 *  CGERC  performs the rank 1 operation
 *
-*     A := alpha*x*conjg( y' ) + A,
+*     A := alpha*x*y**H + A,
 *
 *  where alpha is a scalar, x is an m element vector, y is an n element
 *  vector and A is an m by n matrix.

BLAS/SRC/cgeru.f

@@ -12,7 +12,7 @@ SUBROUTINE CGERU(M,N,ALPHA,X,INCX,Y,INCY,A,LDA)
 *
 *  CGERU  performs the rank 1 operation
 *
-*     A := alpha*x*y' + A,
+*     A := alpha*x*y**T + A,
 *
 *  where alpha is a scalar, x is an m element vector, y is an n element
 *  vector and A is an m by n matrix.

BLAS/SRC/cher.f

@@ -13,7 +13,7 @@ SUBROUTINE CHER(UPLO,N,ALPHA,X,INCX,A,LDA)
 *
 *  CHER   performs the hermitian rank 1 operation
 *
-*     A := alpha*x*conjg( x' ) + A,
+*     A := alpha*x*x**H + A,
 *
 *  where alpha is a real scalar, x is an n element vector and A is an
 *  n by n hermitian matrix.

BLAS/SRC/cher2.f

@@ -13,7 +13,7 @@ SUBROUTINE CHER2(UPLO,N,ALPHA,X,INCX,Y,INCY,A,LDA)
 *
 *  CHER2  performs the hermitian rank 2 operation
 *
-*     A := alpha*x*conjg( y' ) + conjg( alpha )*y*conjg( x' ) + A,
+*     A := alpha*x*y**H + conjg( alpha )*y*x**H + A,
 *
 *  where alpha is a scalar, x and y are n element vectors and A is an n
 *  by n hermitian matrix.

BLAS/SRC/cher2k.f

@@ -14,11 +14,11 @@ SUBROUTINE CHER2K(UPLO,TRANS,N,K,ALPHA,A,LDA,B,LDB,BETA,C,LDC)
 *
 *  CHER2K  performs one of the hermitian rank 2k operations
 *
-*     C := alpha*A*conjg( B' ) + conjg( alpha )*B*conjg( A' ) + beta*C,
+*     C := alpha*A*B**H + conjg( alpha )*B*A**H + beta*C,
 *
 *  or
 *
-*     C := alpha*conjg( A' )*B + conjg( alpha )*conjg( B' )*A + beta*C,
+*     C := alpha*A**H*B + conjg( alpha )*B**H*A + beta*C,
 *
 *  where  alpha and beta  are scalars with  beta  real,  C is an  n by n
 *  hermitian matrix and  A and B  are  n by k matrices in the first case
@@ -44,12 +44,12 @@ SUBROUTINE CHER2K(UPLO,TRANS,N,K,ALPHA,A,LDA,B,LDB,BETA,C,LDC)
 *           On entry,  TRANS  specifies the operation to be performed as
 *           follows:
 *
-*              TRANS = 'N' or 'n'    C := alpha*A*conjg( B' )          +
-*                                         conjg( alpha )*B*conjg( A' ) +
+*              TRANS = 'N' or 'n'    C := alpha*A*B**H          +
+*                                         conjg( alpha )*B*A**H +
 *                                         beta*C.
 *
-*              TRANS = 'C' or 'c'    C := alpha*conjg( A' )*B          +
-*                                         conjg( alpha )*conjg( B' )*A +
+*              TRANS = 'C' or 'c'    C := alpha*A**H*B          +
+*                                         conjg( alpha )*B**H*A +
 *                                         beta*C.
 *
 *           Unchanged on exit.
@@ -242,7 +242,7 @@ SUBROUTINE CHER2K(UPLO,TRANS,N,K,ALPHA,A,LDA,B,LDB,BETA,C,LDC)
 *
       IF (LSAME(TRANS,'N')) THEN
 *
-*        Form  C := alpha*A*conjg( B' ) + conjg( alpha )*B*conjg( A' ) +
+*        Form  C := alpha*A*B**H + conjg( alpha )*B*A**H +
 *                   C.
 *
           IF (UPPER) THEN
@@ -302,7 +302,7 @@ SUBROUTINE CHER2K(UPLO,TRANS,N,K,ALPHA,A,LDA,B,LDB,BETA,C,LDC)
           END IF
       ELSE
 *
-*        Form  C := alpha*conjg( A' )*B + conjg( alpha )*conjg( B' )*A +
+*        Form  C := alpha*A**H*B + conjg( alpha )*B**H*A +
 *                   C.
 *
           IF (UPPER) THEN

BLAS/SRC/cherk.f

@@ -13,11 +13,11 @@ SUBROUTINE CHERK(UPLO,TRANS,N,K,ALPHA,A,LDA,BETA,C,LDC)
 *
 *  CHERK  performs one of the hermitian rank k operations
 *
-*     C := alpha*A*conjg( A' ) + beta*C,
+*     C := alpha*A*A**H + beta*C,
 *
 *  or
 *
-*     C := alpha*conjg( A' )*A + beta*C,
+*     C := alpha*A**H*A + beta*C,
 *
 *  where  alpha and beta  are  real scalars,  C is an  n by n  hermitian
 *  matrix and  A  is an  n by k  matrix in the  first case and a  k by n
@@ -43,9 +43,9 @@ SUBROUTINE CHERK(UPLO,TRANS,N,K,ALPHA,A,LDA,BETA,C,LDC)
 *           On entry,  TRANS  specifies the operation to be performed as
 *           follows:
 *
-*              TRANS = 'N' or 'n'   C := alpha*A*conjg( A' ) + beta*C.
+*              TRANS = 'N' or 'n'   C := alpha*A*A**H + beta*C.
 *
-*              TRANS = 'C' or 'c'   C := alpha*conjg( A' )*A + beta*C.
+*              TRANS = 'C' or 'c'   C := alpha*A**H*A + beta*C.
 *
 *           Unchanged on exit.
 *
@@ -219,7 +219,7 @@ SUBROUTINE CHERK(UPLO,TRANS,N,K,ALPHA,A,LDA,BETA,C,LDC)
 *
       IF (LSAME(TRANS,'N')) THEN
 *
-*        Form  C := alpha*A*conjg( A' ) + beta*C.
+*        Form  C := alpha*A*A**H + beta*C.
 *
           IF (UPPER) THEN
               DO 130 J = 1,N
@@ -272,7 +272,7 @@ SUBROUTINE CHERK(UPLO,TRANS,N,K,ALPHA,A,LDA,BETA,C,LDC)
           END IF
       ELSE
 *
-*        Form  C := alpha*conjg( A' )*A + beta*C.
+*        Form  C := alpha*A**H*A + beta*C.
 *
           IF (UPPER) THEN
               DO 220 J = 1,N

BLAS/SRC/chpr.f

@@ -13,7 +13,7 @@ SUBROUTINE CHPR(UPLO,N,ALPHA,X,INCX,AP)
 *
 *  CHPR    performs the hermitian rank 1 operation
 *
-*     A := alpha*x*conjg( x' ) + A,
+*     A := alpha*x*x**H + A,
 *
 *  where alpha is a real scalar, x is an n element vector and A is an
 *  n by n hermitian matrix, supplied in packed form.