diff options
| author | julie <julielangou@users.noreply.github.com> | 2011-04-02 11:46:46 +0000 |
|---|---|---|
| committer | julie <julielangou@users.noreply.github.com> | 2011-04-02 11:46:46 +0000 |
| commit | d9d50d133bdf00870f37a31faefea58f0926a0a6 (patch) | |
| tree | eedf9e8992fc27e42ce5e71d7e7bf3dc40fc736c /BLAS | |
| parent | f2953573ede24d7f8c01fdb18de48f65f00a9943 (diff) | |
| download | lapack-d9d50d133bdf00870f37a31faefea58f0926a0a6.tar.gz lapack-d9d50d133bdf00870f37a31faefea58f0926a0a6.tar.bz2 lapack-d9d50d133bdf00870f37a31faefea58f0926a0a6.zip | |
First pass for BLAS to homgenize notation for transpose (**T) and conjugate transpose (**H)
Corresponds to bug0024
Diffstat (limited to 'BLAS')
80 files changed, 322 insertions, 322 deletions
diff --git a/BLAS/SRC/cgbmv.f b/BLAS/SRC/cgbmv.f index a3dd2799..6392e1d8 100644 --- a/BLAS/SRC/cgbmv.f +++ b/BLAS/SRC/cgbmv.f @@ -13,9 +13,9 @@ * * 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 @@ * * 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 @@ 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 diff --git a/BLAS/SRC/cgemm.f b/BLAS/SRC/cgemm.f index 9eb7c6a2..24548db1 100644 --- a/BLAS/SRC/cgemm.f +++ b/BLAS/SRC/cgemm.f @@ -17,7 +17,7 @@ * * 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 @@ * * 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 @@ * * 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 @@ 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 @@ 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 @@ 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 @@ 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 @@ 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 @@ 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 @@ 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 @@ 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 diff --git a/BLAS/SRC/cgemv.f b/BLAS/SRC/cgemv.f index b39e8dc5..d9e55f9a 100644 --- a/BLAS/SRC/cgemv.f +++ b/BLAS/SRC/cgemv.f @@ -13,9 +13,9 @@ * * 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 @@ * * 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 @@ 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 diff --git a/BLAS/SRC/cgerc.f b/BLAS/SRC/cgerc.f index 3860871b..ffc578a1 100644 --- a/BLAS/SRC/cgerc.f +++ b/BLAS/SRC/cgerc.f @@ -12,7 +12,7 @@ * * 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. diff --git a/BLAS/SRC/cgeru.f b/BLAS/SRC/cgeru.f index f9715b6b..712aba2a 100644 --- a/BLAS/SRC/cgeru.f +++ b/BLAS/SRC/cgeru.f @@ -12,7 +12,7 @@ * * 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. diff --git a/BLAS/SRC/cher.f b/BLAS/SRC/cher.f index 1ef7154e..b8d5ac01 100644 --- a/BLAS/SRC/cher.f +++ b/BLAS/SRC/cher.f @@ -13,7 +13,7 @@ * * 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. diff --git a/BLAS/SRC/cher2.f b/BLAS/SRC/cher2.f index 861cd94e..3cc1a7d0 100644 --- a/BLAS/SRC/cher2.f +++ b/BLAS/SRC/cher2.f @@ -13,7 +13,7 @@ * * 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. diff --git a/BLAS/SRC/cher2k.f b/BLAS/SRC/cher2k.f index 328fcce1..368c1878 100644 --- a/BLAS/SRC/cher2k.f +++ b/BLAS/SRC/cher2k.f @@ -14,11 +14,11 @@ * * 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 @@ * 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 @@ * 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 @@ 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 diff --git a/BLAS/SRC/cherk.f b/BLAS/SRC/cherk.f index d957e578..b2c1d556 100644 --- a/BLAS/SRC/cherk.f +++ b/BLAS/SRC/cherk.f @@ -13,11 +13,11 @@ * * 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 @@ * 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 @@ * 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 @@ 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 diff --git a/BLAS/SRC/chpr.f b/BLAS/SRC/chpr.f index 11bd5c6e..747b74af 100644 --- a/BLAS/SRC/chpr.f +++ b/BLAS/SRC/chpr.f @@ -13,7 +13,7 @@ * * 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. diff --git a/BLAS/SRC/chpr2.f b/BLAS/SRC/chpr2.f index a0020ef3..71734cee 100644 --- a/BLAS/SRC/chpr2.f +++ b/BLAS/SRC/chpr2.f @@ -13,7 +13,7 @@ * * CHPR2 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, supplied in packed form. diff --git a/BLAS/SRC/csyr2k.f b/BLAS/SRC/csyr2k.f index fa1f6022..ed1a56cd 100644 --- a/BLAS/SRC/csyr2k.f +++ b/BLAS/SRC/csyr2k.f @@ -13,11 +13,11 @@ * * CSYR2K performs one of the symmetric rank 2k operations * -* C := alpha*A*B' + alpha*B*A' + beta*C, +* C := alpha*A*B**T + alpha*B*A**T + beta*C, * * or * -* C := alpha*A'*B + alpha*B'*A + beta*C, +* C := alpha*A**T*B + alpha*B**T*A + beta*C, * * where alpha and beta are scalars, C is an n by n symmetric matrix * and A and B are n by k matrices in the first case and k by n @@ -43,10 +43,10 @@ * On entry, TRANS specifies the operation to be performed as * follows: * -* TRANS = 'N' or 'n' C := alpha*A*B' + alpha*B*A' + +* TRANS = 'N' or 'n' C := alpha*A*B**T + alpha*B*A**T + * beta*C. * -* TRANS = 'T' or 't' C := alpha*A'*B + alpha*B'*A + +* TRANS = 'T' or 't' C := alpha*A**T*B + alpha*B**T*A + * beta*C. * * Unchanged on exit. @@ -231,7 +231,7 @@ * IF (LSAME(TRANS,'N')) THEN * -* Form C := alpha*A*B' + alpha*B*A' + C. +* Form C := alpha*A*B**T + alpha*B*A**T + C. * IF (UPPER) THEN DO 130 J = 1,N @@ -280,7 +280,7 @@ END IF ELSE * -* Form C := alpha*A'*B + alpha*B'*A + C. +* Form C := alpha*A**T*B + alpha*B**T*A + C. * IF (UPPER) THEN DO 210 J = 1,N diff --git a/BLAS/SRC/csyrk.f b/BLAS/SRC/csyrk.f index 3846631d..7b346ff0 100644 --- a/BLAS/SRC/csyrk.f +++ b/BLAS/SRC/csyrk.f @@ -13,11 +13,11 @@ * * CSYRK performs one of the symmetric rank k operations * -* C := alpha*A*A' + beta*C, +* C := alpha*A*A**T + beta*C, * * or * -* C := alpha*A'*A + beta*C, +* C := alpha*A**T*A + beta*C, * * where alpha and beta are scalars, C is an n by n symmetric matrix * and A is an n by k matrix in the first case and a k by n matrix @@ -43,9 +43,9 @@ * On entry, TRANS specifies the operation to be performed as * follows: * -* TRANS = 'N' or 'n' C := alpha*A*A' + beta*C. +* TRANS = 'N' or 'n' C := alpha*A*A**T + beta*C. * -* TRANS = 'T' or 't' C := alpha*A'*A + beta*C. +* TRANS = 'T' or 't' C := alpha*A**T*A + beta*C. * * Unchanged on exit. * @@ -212,7 +212,7 @@ * IF (LSAME(TRANS,'N')) THEN * -* Form C := alpha*A*A' + beta*C. +* Form C := alpha*A*A**T + beta*C. * IF (UPPER) THEN DO 130 J = 1,N @@ -257,7 +257,7 @@ END IF ELSE * -* Form C := alpha*A'*A + beta*C. +* Form C := alpha*A**T*A + beta*C. * IF (UPPER) THEN DO 210 J = 1,N diff --git a/BLAS/SRC/ctbmv.f b/BLAS/SRC/ctbmv.f index a9f4d896..37e0e7f8 100644 --- a/BLAS/SRC/ctbmv.f +++ b/BLAS/SRC/ctbmv.f @@ -12,7 +12,7 @@ * * CTBMV performs one of the matrix-vector operations * -* x := A*x, or x := A'*x, or x := conjg( A' )*x, +* x := A*x, or x := A**T*x, or x := A**H*x, * * where x is an n element vector and A is an n by n unit, or non-unit, * upper or lower triangular band matrix, with ( k + 1 ) diagonals. @@ -36,9 +36,9 @@ * * TRANS = 'N' or 'n' x := A*x. * -* TRANS = 'T' or 't' x := A'*x. +* TRANS = 'T' or 't' x := A**T*x. * -* TRANS = 'C' or 'c' x := conjg( A' )*x. +* TRANS = 'C' or 'c' x := A**H*x. * * Unchanged on exit. * @@ -269,7 +269,7 @@ END IF ELSE * -* Form x := A'*x or x := conjg( A' )*x. +* Form x := A**T*x or x := A**H*x. * IF (LSAME(UPLO,'U')) THEN KPLUS1 = K + 1 diff --git a/BLAS/SRC/ctbsv.f b/BLAS/SRC/ctbsv.f index 853b9d75..bdbd9a3f 100644 --- a/BLAS/SRC/ctbsv.f +++ b/BLAS/SRC/ctbsv.f @@ -12,7 +12,7 @@ * * CTBSV solves one of the systems of equations * -* A*x = b, or A'*x = b, or conjg( A' )*x = b, +* A*x = b, or A**T*x = b, or A**H*x = b, * * where b and x are n element vectors and A is an n by n unit, or * non-unit, upper or lower triangular band matrix, with ( k + 1 ) @@ -40,9 +40,9 @@ * * TRANS = 'N' or 'n' A*x = b. * -* TRANS = 'T' or 't' A'*x = b. +* TRANS = 'T' or 't' A**T*x = b. * -* TRANS = 'C' or 'c' conjg( A' )*x = b. +* TRANS = 'C' or 'c' A**H*x = b. * * Unchanged on exit. * @@ -272,7 +272,7 @@ END IF ELSE * -* Form x := inv( A' )*x or x := inv( conjg( A') )*x. +* Form x := inv( A**T )*x or x := inv( A**H )*x. * IF (LSAME(UPLO,'U')) THEN KPLUS1 = K + 1 diff --git a/BLAS/SRC/ctpmv.f b/BLAS/SRC/ctpmv.f index 03660f11..53b1b165 100644 --- a/BLAS/SRC/ctpmv.f +++ b/BLAS/SRC/ctpmv.f @@ -12,7 +12,7 @@ * * CTPMV performs one of the matrix-vector operations * -* x := A*x, or x := A'*x, or x := conjg( A' )*x, +* x := A*x, or x := A**T*x, or x := A**H*x, * * where x is an n element vector and A is an n by n unit, or non-unit, * upper or lower triangular matrix, supplied in packed form. @@ -36,9 +36,9 @@ * * TRANS = 'N' or 'n' x := A*x. * -* TRANS = 'T' or 't' x := A'*x. +* TRANS = 'T' or 't' x := A**T*x. * -* TRANS = 'C' or 'c' x := conjg( A' )*x. +* TRANS = 'C' or 'c' x := A**H*x. * * Unchanged on exit. * @@ -228,7 +228,7 @@ END IF ELSE * -* Form x := A'*x or x := conjg( A' )*x. +* Form x := A**T*x or x := A**H*x. * IF (LSAME(UPLO,'U')) THEN KK = (N* (N+1))/2 diff --git a/BLAS/SRC/ctpsv.f b/BLAS/SRC/ctpsv.f index 1804797e..3a7da00c 100644 --- a/BLAS/SRC/ctpsv.f +++ b/BLAS/SRC/ctpsv.f @@ -12,7 +12,7 @@ * * CTPSV solves one of the systems of equations * -* A*x = b, or A'*x = b, or conjg( A' )*x = b, +* A*x = b, or A**T*x = b, or A**H*x = b, * * where b and x are n element vectors and A is an n by n unit, or * non-unit, upper or lower triangular matrix, supplied in packed form. @@ -39,9 +39,9 @@ * * TRANS = 'N' or 'n' A*x = b. * -* TRANS = 'T' or 't' A'*x = b. +* TRANS = 'T' or 't' A**T*x = b. * -* TRANS = 'C' or 'c' conjg( A' )*x = b. +* TRANS = 'C' or 'c' A**H*x = b. * * Unchanged on exit. * @@ -229,7 +229,7 @@ END IF ELSE * -* Form x := inv( A' )*x or x := inv( conjg( A' ) )*x. +* Form x := inv( A**T )*x or x := inv( A**H )*x. * IF (LSAME(UPLO,'U')) THEN KK = 1 diff --git a/BLAS/SRC/ctrmm.f b/BLAS/SRC/ctrmm.f index 4f599b29..e05a623e 100644 --- a/BLAS/SRC/ctrmm.f +++ b/BLAS/SRC/ctrmm.f @@ -18,7 +18,7 @@ * where alpha is a scalar, B is an m by n matrix, A is a unit, or * non-unit, upper or lower triangular matrix and op( A ) is one of * -* op( A ) = A or op( A ) = A' or op( A ) = conjg( A' ). +* op( A ) = A or op( A ) = A**T or op( A ) = A**H. * * Arguments * ========== @@ -49,9 +49,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 ) = conjg( A' ). +* TRANSA = 'C' or 'c' op( A ) = A**H. * * Unchanged on exit. * @@ -237,7 +237,7 @@ END IF ELSE * -* Form B := alpha*A'*B or B := alpha*conjg( A' )*B. +* Form B := alpha*A**T*B or B := alpha*A**H*B. * IF (UPPER) THEN DO 120 J = 1,N @@ -317,7 +317,7 @@ END IF ELSE * -* Form B := alpha*B*A' or B := alpha*B*conjg( A' ). +* Form B := alpha*B*A**T or B := alpha*B*A**H. * IF (UPPER) THEN DO 280 K = 1,N diff --git a/BLAS/SRC/ctrmv.f b/BLAS/SRC/ctrmv.f index 631df7ac..ac0afb73 100644 --- a/BLAS/SRC/ctrmv.f +++ b/BLAS/SRC/ctrmv.f @@ -12,7 +12,7 @@ * * CTRMV performs one of the matrix-vector operations * -* x := A*x, or x := A'*x, or x := conjg( A' )*x, +* x := A*x, or x := A**T*x, or x := A**H*x, * * where x is an n element vector and A is an n by n unit, or non-unit, * upper or lower triangular matrix. @@ -36,9 +36,9 @@ * * TRANS = 'N' or 'n' x := A*x. * -* TRANS = 'T' or 't' x := A'*x. +* TRANS = 'T' or 't' x := A**T*x. * -* TRANS = 'C' or 'c' x := conjg( A' )*x. +* TRANS = 'C' or 'c' x := A**H*x. * * Unchanged on exit. * @@ -223,7 +223,7 @@ END IF ELSE * -* Form x := A'*x or x := conjg( A' )*x. +* Form x := A**T*x or x := A**H*x. * IF (LSAME(UPLO,'U')) THEN IF (INCX.EQ.1) THEN diff --git a/BLAS/SRC/ctrsm.f b/BLAS/SRC/ctrsm.f index 96ff19d5..2aeb8727 100644 --- a/BLAS/SRC/ctrsm.f +++ b/BLAS/SRC/ctrsm.f @@ -18,7 +18,7 @@ * where alpha is a scalar, X and B are m by n matrices, A is a unit, or * non-unit, upper or lower triangular matrix and op( A ) is one of * -* op( A ) = A or op( A ) = A' or op( A ) = conjg( A' ). +* op( A ) = A or op( A ) = A**T or op( A ) = A**H. * * The matrix X is overwritten on B. * @@ -51,9 +51,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 ) = conjg( A' ). +* TRANSA = 'C' or 'c' op( A ) = A**H. * * Unchanged on exit. * @@ -245,8 +245,8 @@ END IF ELSE * -* Form B := alpha*inv( A' )*B -* or B := alpha*inv( conjg( A' ) )*B. +* Form B := alpha*inv( A**T )*B +* or B := alpha*inv( A**H )*B. * IF (UPPER) THEN DO 140 J = 1,N @@ -336,8 +336,8 @@ END IF ELSE * -* Form B := alpha*B*inv( A' ) -* or B := alpha*B*inv( conjg( A' ) ). +* Form B := alpha*B*inv( A**T ) +* or B := alpha*B*inv( A**H ). * IF (UPPER) THEN DO 330 K = N,1,-1 diff --git a/BLAS/SRC/ctrsv.f b/BLAS/SRC/ctrsv.f index 97e150a6..93e7e44f 100644 --- a/BLAS/SRC/ctrsv.f +++ b/BLAS/SRC/ctrsv.f @@ -12,7 +12,7 @@ * * CTRSV solves one of the systems of equations * -* A*x = b, or A'*x = b, or conjg( A' )*x = b, +* A*x = b, or A**T*x = b, or A**H*x = b, * * where b and x are n element vectors and A is an n by n unit, or * non-unit, upper or lower triangular matrix. @@ -39,9 +39,9 @@ * * TRANS = 'N' or 'n' A*x = b. * -* TRANS = 'T' or 't' A'*x = b. +* TRANS = 'T' or 't' A**T*x = b. * -* TRANS = 'C' or 'c' conjg( A' )*x = b. +* TRANS = 'C' or 'c' A**H*x = b. * * Unchanged on exit. * @@ -224,7 +224,7 @@ END IF ELSE * -* Form x := inv( A' )*x or x := inv( conjg( A' ) )*x. +* Form x := inv( A**T )*x or x := inv( A**H )*x. * IF (LSAME(UPLO,'U')) THEN IF (INCX.EQ.1) THEN diff --git a/BLAS/SRC/dgbmv.f b/BLAS/SRC/dgbmv.f index 630878b7..1fd5a589 100644 --- a/BLAS/SRC/dgbmv.f +++ b/BLAS/SRC/dgbmv.f @@ -13,7 +13,7 @@ * * DGBMV performs one of the matrix-vector operations * -* y := alpha*A*x + beta*y, or y := alpha*A'*x + beta*y, +* y := alpha*A*x + beta*y, or y := alpha*A**T*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. @@ -27,9 +27,9 @@ * * 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*A'*x + beta*y. +* TRANS = 'C' or 'c' y := alpha*A**T*x + beta*y. * * Unchanged on exit. * @@ -266,7 +266,7 @@ END IF ELSE * -* Form y := alpha*A'*x + y. +* Form y := alpha*A**T*x + y. * JY = KY IF (INCX.EQ.1) THEN diff --git a/BLAS/SRC/dgemm.f b/BLAS/SRC/dgemm.f index 4afccb16..7ac8c46b 100644 --- a/BLAS/SRC/dgemm.f +++ b/BLAS/SRC/dgemm.f @@ -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 diff --git a/BLAS/SRC/dgemv.f b/BLAS/SRC/dgemv.f index 064bdb9f..a4125941 100644 --- a/BLAS/SRC/dgemv.f +++ b/BLAS/SRC/dgemv.f @@ -13,7 +13,7 @@ * * DGEMV performs one of the matrix-vector operations * -* y := alpha*A*x + beta*y, or y := alpha*A'*x + beta*y, +* y := alpha*A*x + beta*y, or y := alpha*A**T*x + beta*y, * * where alpha and beta are scalars, x and y are vectors and A is an * m by n matrix. @@ -27,9 +27,9 @@ * * 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*A'*x + beta*y. +* TRANS = 'C' or 'c' y := alpha*A**T*x + beta*y. * * Unchanged on exit. * @@ -232,7 +232,7 @@ END IF ELSE * -* Form y := alpha*A'*x + y. +* Form y := alpha*A**T*x + y. * JY = KY IF (INCX.EQ.1) THEN diff --git a/BLAS/SRC/dger.f b/BLAS/SRC/dger.f index 0a9dcf2a..1d95257e 100644 --- a/BLAS/SRC/dger.f +++ b/BLAS/SRC/dger.f @@ -12,7 +12,7 @@ * * DGER 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. diff --git a/BLAS/SRC/dspr.f b/BLAS/SRC/dspr.f index 538e4f76..64567fc7 100644 --- a/BLAS/SRC/dspr.f +++ b/BLAS/SRC/dspr.f @@ -13,7 +13,7 @@ * * DSPR performs the symmetric rank 1 operation * -* A := alpha*x*x' + A, +* A := alpha*x*x**T + A, * * where alpha is a real scalar, x is an n element vector and A is an * n by n symmetric matrix, supplied in packed form. diff --git a/BLAS/SRC/dspr2.f b/BLAS/SRC/dspr2.f index 6f6b54a8..22900ad0 100644 --- a/BLAS/SRC/dspr2.f +++ b/BLAS/SRC/dspr2.f @@ -13,7 +13,7 @@ * * DSPR2 performs the symmetric rank 2 operation * -* A := alpha*x*y' + alpha*y*x' + A, +* A := alpha*x*y**T + alpha*y*x**T + A, * * where alpha is a scalar, x and y are n element vectors and A is an * n by n symmetric matrix, supplied in packed form. diff --git a/BLAS/SRC/dsyr.f b/BLAS/SRC/dsyr.f index b418d3c3..5cc4bcb5 100644 --- a/BLAS/SRC/dsyr.f +++ b/BLAS/SRC/dsyr.f @@ -13,7 +13,7 @@ * * DSYR performs the symmetric rank 1 operation * -* A := alpha*x*x' + A, +* A := alpha*x*x**T + A, * * where alpha is a real scalar, x is an n element vector and A is an * n by n symmetric matrix. diff --git a/BLAS/SRC/dsyr2.f b/BLAS/SRC/dsyr2.f index 387969d4..b2b2f3d7 100644 --- a/BLAS/SRC/dsyr2.f +++ b/BLAS/SRC/dsyr2.f @@ -13,7 +13,7 @@ * * DSYR2 performs the symmetric rank 2 operation * -* A := alpha*x*y' + alpha*y*x' + A, +* A := alpha*x*y**T + alpha*y*x**T + A, * * where alpha is a scalar, x and y are n element vectors and A is an n * by n symmetric matrix. diff --git a/BLAS/SRC/dsyr2k.f b/BLAS/SRC/dsyr2k.f index 588e71c1..7289b055 100644 --- a/BLAS/SRC/dsyr2k.f +++ b/BLAS/SRC/dsyr2k.f @@ -13,11 +13,11 @@ * * DSYR2K performs one of the symmetric rank 2k operations * -* C := alpha*A*B' + alpha*B*A' + beta*C, +* C := alpha*A*B**T + alpha*B*A**T + beta*C, * * or * -* C := alpha*A'*B + alpha*B'*A + beta*C, +* C := alpha*A**T*B + alpha*B**T*A + beta*C, * * where alpha and beta are scalars, C is an n by n symmetric matrix * and A and B are n by k matrices in the first case and k by n @@ -43,13 +43,13 @@ * On entry, TRANS specifies the operation to be performed as * follows: * -* TRANS = 'N' or 'n' C := alpha*A*B' + alpha*B*A' + +* TRANS = 'N' or 'n' C := alpha*A*B**T + alpha*B*A**T + * beta*C. * -* TRANS = 'T' or 't' C := alpha*A'*B + alpha*B'*A + +* TRANS = 'T' or 't' C := alpha*A**T*B + alpha*B**T*A + * beta*C. * -* TRANS = 'C' or 'c' C := alpha*A'*B + alpha*B'*A + +* TRANS = 'C' or 'c' C := alpha*A**T*B + alpha*B**T*A + * beta*C. * * Unchanged on exit. @@ -234,7 +234,7 @@ * IF (LSAME(TRANS,'N')) THEN * -* Form C := alpha*A*B' + alpha*B*A' + C. +* Form C := alpha*A*B**T + alpha*B*A**T + C. * IF (UPPER) THEN DO 130 J = 1,N @@ -283,7 +283,7 @@ END IF ELSE * -* Form C := alpha*A'*B + alpha*B'*A + C. +* Form C := alpha*A**T*B + alpha*B**T*A + C. * IF (UPPER) THEN DO 210 J = 1,N diff --git a/BLAS/SRC/dsyrk.f b/BLAS/SRC/dsyrk.f index 05e1f5b0..00fb33ae 100644 --- a/BLAS/SRC/dsyrk.f +++ b/BLAS/SRC/dsyrk.f @@ -13,11 +13,11 @@ * * DSYRK performs one of the symmetric rank k operations * -* C := alpha*A*A' + beta*C, +* C := alpha*A*A**T + beta*C, * * or * -* C := alpha*A'*A + beta*C, +* C := alpha*A**T*A + beta*C, * * where alpha and beta are scalars, C is an n by n symmetric matrix * and A is an n by k matrix in the first case and a k by n matrix @@ -43,11 +43,11 @@ * On entry, TRANS specifies the operation to be performed as * follows: * -* TRANS = 'N' or 'n' C := alpha*A*A' + beta*C. +* TRANS = 'N' or 'n' C := alpha*A*A**T + beta*C. * -* TRANS = 'T' or 't' C := alpha*A'*A + beta*C. +* TRANS = 'T' or 't' C := alpha*A**T*A + beta*C. * -* TRANS = 'C' or 'c' C := alpha*A'*A + beta*C. +* TRANS = 'C' or 'c' C := alpha*A**T*A + beta*C. * * Unchanged on exit. * @@ -213,7 +213,7 @@ * IF (LSAME(TRANS,'N')) THEN * -* Form C := alpha*A*A' + beta*C. +* Form C := alpha*A*A**T + beta*C. * IF (UPPER) THEN DO 130 J = 1,N @@ -258,7 +258,7 @@ END IF ELSE * -* Form C := alpha*A'*A + beta*C. +* Form C := alpha*A**T*A + beta*C. * IF (UPPER) THEN DO 210 J = 1,N diff --git a/BLAS/SRC/dtbmv.f b/BLAS/SRC/dtbmv.f index cfcdb15a..4f8778ad 100644 --- a/BLAS/SRC/dtbmv.f +++ b/BLAS/SRC/dtbmv.f @@ -12,7 +12,7 @@ * * DTBMV performs one of the matrix-vector operations * -* x := A*x, or x := A'*x, +* x := A*x, or x := A**T*x, * * where x is an n element vector and A is an n by n unit, or non-unit, * upper or lower triangular band matrix, with ( k + 1 ) diagonals. @@ -36,9 +36,9 @@ * * TRANS = 'N' or 'n' x := A*x. * -* TRANS = 'T' or 't' x := A'*x. +* TRANS = 'T' or 't' x := A**T*x. * -* TRANS = 'C' or 'c' x := A'*x. +* TRANS = 'C' or 'c' x := A**T*x. * * Unchanged on exit. * @@ -268,7 +268,7 @@ END IF ELSE * -* Form x := A'*x. +* Form x := A**T*x. * IF (LSAME(UPLO,'U')) THEN KPLUS1 = K + 1 diff --git a/BLAS/SRC/dtbsv.f b/BLAS/SRC/dtbsv.f index cfeb0b82..621c1ec7 100644 --- a/BLAS/SRC/dtbsv.f +++ b/BLAS/SRC/dtbsv.f @@ -12,7 +12,7 @@ * * DTBSV solves one of the systems of equations * -* A*x = b, or A'*x = b, +* A*x = b, or A**T*x = b, * * where b and x are n element vectors and A is an n by n unit, or * non-unit, upper or lower triangular band matrix, with ( k + 1 ) @@ -40,9 +40,9 @@ * * TRANS = 'N' or 'n' A*x = b. * -* TRANS = 'T' or 't' A'*x = b. +* TRANS = 'T' or 't' A**T*x = b. * -* TRANS = 'C' or 'c' A'*x = b. +* TRANS = 'C' or 'c' A**T*x = b. * * Unchanged on exit. * @@ -271,7 +271,7 @@ END IF ELSE * -* Form x := inv( A')*x. +* Form x := inv( A**T)*x. * IF (LSAME(UPLO,'U')) THEN KPLUS1 = K + 1 diff --git a/BLAS/SRC/dtpmv.f b/BLAS/SRC/dtpmv.f index b25a4058..1def763e 100644 --- a/BLAS/SRC/dtpmv.f +++ b/BLAS/SRC/dtpmv.f @@ -12,7 +12,7 @@ * * DTPMV performs one of the matrix-vector operations * -* x := A*x, or x := A'*x, +* x := A*x, or x := A**T*x, * * where x is an n element vector and A is an n by n unit, or non-unit, * upper or lower triangular matrix, supplied in packed form. @@ -36,9 +36,9 @@ * * TRANS = 'N' or 'n' x := A*x. * -* TRANS = 'T' or 't' x := A'*x. +* TRANS = 'T' or 't' x := A**T*x. * -* TRANS = 'C' or 'c' x := A'*x. +* TRANS = 'C' or 'c' x := A**T*x. * * Unchanged on exit. * @@ -224,7 +224,7 @@ END IF ELSE * -* Form x := A'*x. +* Form x := A**T*x. * IF (LSAME(UPLO,'U')) THEN KK = (N* (N+1))/2 diff --git a/BLAS/SRC/dtpsv.f b/BLAS/SRC/dtpsv.f index c7e58d32..49b65237 100644 --- a/BLAS/SRC/dtpsv.f +++ b/BLAS/SRC/dtpsv.f @@ -12,7 +12,7 @@ * * DTPSV solves one of the systems of equations * -* A*x = b, or A'*x = b, +* A*x = b, or A**T*x = b, * * where b and x are n element vectors and A is an n by n unit, or * non-unit, upper or lower triangular matrix, supplied in packed form. @@ -39,9 +39,9 @@ * * TRANS = 'N' or 'n' A*x = b. * -* TRANS = 'T' or 't' A'*x = b. +* TRANS = 'T' or 't' A**T*x = b. * -* TRANS = 'C' or 'c' A'*x = b. +* TRANS = 'C' or 'c' A**T*x = b. * * Unchanged on exit. * @@ -225,7 +225,7 @@ END IF ELSE * -* Form x := inv( A' )*x. +* Form x := inv( A**T )*x. * IF (LSAME(UPLO,'U')) THEN KK = 1 diff --git a/BLAS/SRC/dtrmm.f b/BLAS/SRC/dtrmm.f index 931748db..fc03769f 100644 --- a/BLAS/SRC/dtrmm.f +++ b/BLAS/SRC/dtrmm.f @@ -18,7 +18,7 @@ * where alpha is a scalar, B is an m by n matrix, A is a unit, or * non-unit, upper or lower triangular matrix and op( A ) is one of * -* op( A ) = A or op( A ) = A'. +* op( A ) = A or op( A ) = A**T. * * Arguments * ========== @@ -49,9 +49,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. * @@ -234,7 +234,7 @@ END IF ELSE * -* Form B := alpha*A'*B. +* Form B := alpha*A**T*B. * IF (UPPER) THEN DO 110 J = 1,N @@ -300,7 +300,7 @@ END IF ELSE * -* Form B := alpha*B*A'. +* Form B := alpha*B*A**T. * IF (UPPER) THEN DO 260 K = 1,N diff --git a/BLAS/SRC/dtrmv.f b/BLAS/SRC/dtrmv.f index 0381e1d2..5356cbbc 100644 --- a/BLAS/SRC/dtrmv.f +++ b/BLAS/SRC/dtrmv.f @@ -12,7 +12,7 @@ * * DTRMV performs one of the matrix-vector operations * -* x := A*x, or x := A'*x, +* x := A*x, or x := A**T*x, * * where x is an n element vector and A is an n by n unit, or non-unit, * upper or lower triangular matrix. @@ -36,9 +36,9 @@ * * TRANS = 'N' or 'n' x := A*x. * -* TRANS = 'T' or 't' x := A'*x. +* TRANS = 'T' or 't' x := A**T*x. * -* TRANS = 'C' or 'c' x := A'*x. +* TRANS = 'C' or 'c' x := A**T*x. * * Unchanged on exit. * @@ -222,7 +222,7 @@ END IF ELSE * -* Form x := A'*x. +* Form x := A**T*x. * IF (LSAME(UPLO,'U')) THEN IF (INCX.EQ.1) THEN diff --git a/BLAS/SRC/dtrsm.f b/BLAS/SRC/dtrsm.f index 88811104..dec5c8d0 100644 --- a/BLAS/SRC/dtrsm.f +++ b/BLAS/SRC/dtrsm.f @@ -18,7 +18,7 @@ * where alpha is a scalar, X and B are m by n matrices, A is a unit, or * non-unit, upper or lower triangular matrix and op( A ) is one of * -* op( A ) = A or op( A ) = A'. +* op( A ) = A or op( A ) = A**T. * * The matrix X is overwritten on B. * @@ -51,9 +51,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. * @@ -243,7 +243,7 @@ END IF ELSE * -* Form B := alpha*inv( A' )*B. +* Form B := alpha*inv( A**T )*B. * IF (UPPER) THEN DO 130 J = 1,N @@ -319,7 +319,7 @@ END IF ELSE * -* Form B := alpha*B*inv( A' ). +* Form B := alpha*B*inv( A**T ). * IF (UPPER) THEN DO 310 K = N,1,-1 diff --git a/BLAS/SRC/dtrsv.f b/BLAS/SRC/dtrsv.f index 847fd391..f901e699 100644 --- a/BLAS/SRC/dtrsv.f +++ b/BLAS/SRC/dtrsv.f @@ -12,7 +12,7 @@ * * DTRSV solves one of the systems of equations * -* A*x = b, or A'*x = b, +* A*x = b, or A**T*x = b, * * where b and x are n element vectors and A is an n by n unit, or * non-unit, upper or lower triangular matrix. @@ -39,9 +39,9 @@ * * TRANS = 'N' or 'n' A*x = b. * -* TRANS = 'T' or 't' A'*x = b. +* TRANS = 'T' or 't' A**T*x = b. * -* TRANS = 'C' or 'c' A'*x = b. +* TRANS = 'C' or 'c' A**T*x = b. * * Unchanged on exit. * @@ -221,7 +221,7 @@ END IF ELSE * -* Form x := inv( A' )*x. +* Form x := inv( A**T )*x. * IF (LSAME(UPLO,'U')) THEN IF (INCX.EQ.1) THEN diff --git a/BLAS/SRC/dznrm2.f b/BLAS/SRC/dznrm2.f index d5e76c54..71b5674b 100644 --- a/BLAS/SRC/dznrm2.f +++ b/BLAS/SRC/dznrm2.f @@ -12,7 +12,7 @@ * DZNRM2 returns the euclidean norm of a vector via the function * name, so that * -* DZNRM2 := sqrt( conjg( x' )*x ) +* DZNRM2 := sqrt( x**H*x ) * * Further Details * =============== diff --git a/BLAS/SRC/scnrm2.f b/BLAS/SRC/scnrm2.f index f8288026..1738c315 100644 --- a/BLAS/SRC/scnrm2.f +++ b/BLAS/SRC/scnrm2.f @@ -12,7 +12,7 @@ * SCNRM2 returns the euclidean norm of a vector via the function * name, so that * -* SCNRM2 := sqrt( conjg( x' )*x ) +* SCNRM2 := sqrt( x**H*x ) * * Further Details * =============== diff --git a/BLAS/SRC/sgbmv.f b/BLAS/SRC/sgbmv.f index 1136f702..9c81731a 100644 --- a/BLAS/SRC/sgbmv.f +++ b/BLAS/SRC/sgbmv.f @@ -13,7 +13,7 @@ * * SGBMV performs one of the matrix-vector operations * -* y := alpha*A*x + beta*y, or y := alpha*A'*x + beta*y, +* y := alpha*A*x + beta*y, or y := alpha*A**T*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. @@ -27,9 +27,9 @@ * * 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*A'*x + beta*y. +* TRANS = 'C' or 'c' y := alpha*A**T*x + beta*y. * * Unchanged on exit. * @@ -266,7 +266,7 @@ END IF ELSE * -* Form y := alpha*A'*x + y. +* Form y := alpha*A**T*x + y. * JY = KY IF (INCX.EQ.1) THEN diff --git a/BLAS/SRC/sgemm.f b/BLAS/SRC/sgemm.f index caafaba0..36c00d03 100644 --- a/BLAS/SRC/sgemm.f +++ b/BLAS/SRC/sgemm.f @@ -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 diff --git a/BLAS/SRC/sgemv.f b/BLAS/SRC/sgemv.f index 90eb6394..afae2698 100644 --- a/BLAS/SRC/sgemv.f +++ b/BLAS/SRC/sgemv.f @@ -13,7 +13,7 @@ * * SGEMV performs one of the matrix-vector operations * -* y := alpha*A*x + beta*y, or y := alpha*A'*x + beta*y, +* y := alpha*A*x + beta*y, or y := alpha*A**T*x + beta*y, * * where alpha and beta are scalars, x and y are vectors and A is an * m by n matrix. @@ -27,9 +27,9 @@ * * 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*A'*x + beta*y. +* TRANS = 'C' or 'c' y := alpha*A**T*x + beta*y. * * Unchanged on exit. * @@ -232,7 +232,7 @@ END IF ELSE * -* Form y := alpha*A'*x + y. +* Form y := alpha*A**T*x + y. * JY = KY IF (INCX.EQ.1) THEN diff --git a/BLAS/SRC/sger.f b/BLAS/SRC/sger.f index 3ddd10c7..828b7454 100644 --- a/BLAS/SRC/sger.f +++ b/BLAS/SRC/sger.f @@ -12,7 +12,7 @@ * * SGER 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. diff --git a/BLAS/SRC/sspr.f b/BLAS/SRC/sspr.f index bae92612..6c2a21c6 100644 --- a/BLAS/SRC/sspr.f +++ b/BLAS/SRC/sspr.f @@ -13,7 +13,7 @@ * * SSPR performs the symmetric rank 1 operation * -* A := alpha*x*x' + A, +* A := alpha*x*x**T + A, * * where alpha is a real scalar, x is an n element vector and A is an * n by n symmetric matrix, supplied in packed form. diff --git a/BLAS/SRC/sspr2.f b/BLAS/SRC/sspr2.f index cd27c734..51831071 100644 --- a/BLAS/SRC/sspr2.f +++ b/BLAS/SRC/sspr2.f @@ -13,7 +13,7 @@ * * SSPR2 performs the symmetric rank 2 operation * -* A := alpha*x*y' + alpha*y*x' + A, +* A := alpha*x*y**T + alpha*y*x**T + A, * * where alpha is a scalar, x and y are n element vectors and A is an * n by n symmetric matrix, supplied in packed form. diff --git a/BLAS/SRC/ssyr.f b/BLAS/SRC/ssyr.f index 6c7a8770..1a546009 100644 --- a/BLAS/SRC/ssyr.f +++ b/BLAS/SRC/ssyr.f @@ -13,7 +13,7 @@ * * SSYR performs the symmetric rank 1 operation * -* A := alpha*x*x' + A, +* A := alpha*x*x**T + A, * * where alpha is a real scalar, x is an n element vector and A is an * n by n symmetric matrix. diff --git a/BLAS/SRC/ssyr2.f b/BLAS/SRC/ssyr2.f index cc472ccd..545cfe55 100644 --- a/BLAS/SRC/ssyr2.f +++ b/BLAS/SRC/ssyr2.f @@ -13,7 +13,7 @@ * * SSYR2 performs the symmetric rank 2 operation * -* A := alpha*x*y' + alpha*y*x' + A, +* A := alpha*x*y**T + alpha*y*x**T + A, * * where alpha is a scalar, x and y are n element vectors and A is an n * by n symmetric matrix. diff --git a/BLAS/SRC/ssyr2k.f b/BLAS/SRC/ssyr2k.f index 7f0f9875..c9f44aa8 100644 --- a/BLAS/SRC/ssyr2k.f +++ b/BLAS/SRC/ssyr2k.f @@ -13,11 +13,11 @@ * * SSYR2K performs one of the symmetric rank 2k operations * -* C := alpha*A*B' + alpha*B*A' + beta*C, +* C := alpha*A*B**T + alpha*B*A**T + beta*C, * * or * -* C := alpha*A'*B + alpha*B'*A + beta*C, +* C := alpha*A**T*B + alpha*B**T*A + beta*C, * * where alpha and beta are scalars, C is an n by n symmetric matrix * and A and B are n by k matrices in the first case and k by n @@ -43,13 +43,13 @@ * On entry, TRANS specifies the operation to be performed as * follows: * -* TRANS = 'N' or 'n' C := alpha*A*B' + alpha*B*A' + +* TRANS = 'N' or 'n' C := alpha*A*B**T + alpha*B*A**T + * beta*C. * -* TRANS = 'T' or 't' C := alpha*A'*B + alpha*B'*A + +* TRANS = 'T' or 't' C := alpha*A**T*B + alpha*B**T*A + * beta*C. * -* TRANS = 'C' or 'c' C := alpha*A'*B + alpha*B'*A + +* TRANS = 'C' or 'c' C := alpha*A**T*B + alpha*B**T*A + * beta*C. * * Unchanged on exit. @@ -234,7 +234,7 @@ * IF (LSAME(TRANS,'N')) THEN * -* Form C := alpha*A*B' + alpha*B*A' + C. +* Form C := alpha*A*B**T + alpha*B*A**T + C. * IF (UPPER) THEN DO 130 J = 1,N @@ -283,7 +283,7 @@ END IF ELSE * -* Form C := alpha*A'*B + alpha*B'*A + C. +* Form C := alpha*A**T*B + alpha*B**T*A + C. * IF (UPPER) THEN DO 210 J = 1,N diff --git a/BLAS/SRC/ssyrk.f b/BLAS/SRC/ssyrk.f index a56ba40c..33b7296f 100644 --- a/BLAS/SRC/ssyrk.f +++ b/BLAS/SRC/ssyrk.f @@ -13,11 +13,11 @@ * * SSYRK performs one of the symmetric rank k operations * -* C := alpha*A*A' + beta*C, +* C := alpha*A*A**T + beta*C, * * or * -* C := alpha*A'*A + beta*C, +* C := alpha*A**T*A + beta*C, * * where alpha and beta are scalars, C is an n by n symmetric matrix * and A is an n by k matrix in the first case and a k by n matrix @@ -43,11 +43,11 @@ * On entry, TRANS specifies the operation to be performed as * follows: * -* TRANS = 'N' or 'n' C := alpha*A*A' + beta*C. +* TRANS = 'N' or 'n' C := alpha*A*A**T + beta*C. * -* TRANS = 'T' or 't' C := alpha*A'*A + beta*C. +* TRANS = 'T' or 't' C := alpha*A**T*A + beta*C. * -* TRANS = 'C' or 'c' C := alpha*A'*A + beta*C. +* TRANS = 'C' or 'c' C := alpha*A**T*A + beta*C. * * Unchanged on exit. * @@ -213,7 +213,7 @@ * IF (LSAME(TRANS,'N')) THEN * -* Form C := alpha*A*A' + beta*C. +* Form C := alpha*A*A**T + beta*C. * IF (UPPER) THEN DO 130 J = 1,N @@ -258,7 +258,7 @@ END IF ELSE * -* Form C := alpha*A'*A + beta*C. +* Form C := alpha*A**T*A + beta*C. * IF (UPPER) THEN DO 210 J = 1,N diff --git a/BLAS/SRC/stbmv.f b/BLAS/SRC/stbmv.f index 6e57a9d5..7446b4e6 100644 --- a/BLAS/SRC/stbmv.f +++ b/BLAS/SRC/stbmv.f @@ -12,7 +12,7 @@ * * STBMV performs one of the matrix-vector operations * -* x := A*x, or x := A'*x, +* x := A*x, or x := A**T*x, * * where x is an n element vector and A is an n by n unit, or non-unit, * upper or lower triangular band matrix, with ( k + 1 ) diagonals. @@ -36,9 +36,9 @@ * * TRANS = 'N' or 'n' x := A*x. * -* TRANS = 'T' or 't' x := A'*x. +* TRANS = 'T' or 't' x := A**T*x. * -* TRANS = 'C' or 'c' x := A'*x. +* TRANS = 'C' or 'c' x := A**T*x. * * Unchanged on exit. * @@ -268,7 +268,7 @@ END IF ELSE * -* Form x := A'*x. +* Form x := A**T*x. * IF (LSAME(UPLO,'U')) THEN KPLUS1 = K + 1 diff --git a/BLAS/SRC/stbsv.f b/BLAS/SRC/stbsv.f index b846be85..d6d247e2 100644 --- a/BLAS/SRC/stbsv.f +++ b/BLAS/SRC/stbsv.f @@ -12,7 +12,7 @@ * * STBSV solves one of the systems of equations * -* A*x = b, or A'*x = b, +* A*x = b, or A**T*x = b, * * where b and x are n element vectors and A is an n by n unit, or * non-unit, upper or lower triangular band matrix, with ( k + 1 ) @@ -40,9 +40,9 @@ * * TRANS = 'N' or 'n' A*x = b. * -* TRANS = 'T' or 't' A'*x = b. +* TRANS = 'T' or 't' A**T*x = b. * -* TRANS = 'C' or 'c' A'*x = b. +* TRANS = 'C' or 'c' A**T*x = b. * * Unchanged on exit. * @@ -271,7 +271,7 @@ END IF ELSE * -* Form x := inv( A')*x. +* Form x := inv( A**T)*x. * IF (LSAME(UPLO,'U')) THEN KPLUS1 = K + 1 diff --git a/BLAS/SRC/stpmv.f b/BLAS/SRC/stpmv.f index 3c02c275..573288e6 100644 --- a/BLAS/SRC/stpmv.f +++ b/BLAS/SRC/stpmv.f @@ -12,7 +12,7 @@ * * STPMV performs one of the matrix-vector operations * -* x := A*x, or x := A'*x, +* x := A*x, or x := A**T*x, * * where x is an n element vector and A is an n by n unit, or non-unit, * upper or lower triangular matrix, supplied in packed form. @@ -36,9 +36,9 @@ * * TRANS = 'N' or 'n' x := A*x. * -* TRANS = 'T' or 't' x := A'*x. +* TRANS = 'T' or 't' x := A**T*x. * -* TRANS = 'C' or 'c' x := A'*x. +* TRANS = 'C' or 'c' x := A**T*x. * * Unchanged on exit. * @@ -224,7 +224,7 @@ END IF ELSE * -* Form x := A'*x. +* Form x := A**T*x. * IF (LSAME(UPLO,'U')) THEN KK = (N* (N+1))/2 diff --git a/BLAS/SRC/stpsv.f b/BLAS/SRC/stpsv.f index 7d95efbd..aca8a7fc 100644 --- a/BLAS/SRC/stpsv.f +++ b/BLAS/SRC/stpsv.f @@ -12,7 +12,7 @@ * * STPSV solves one of the systems of equations * -* A*x = b, or A'*x = b, +* A*x = b, or A**T*x = b, * * where b and x are n element vectors and A is an n by n unit, or * non-unit, upper or lower triangular matrix, supplied in packed form. @@ -39,9 +39,9 @@ * * TRANS = 'N' or 'n' A*x = b. * -* TRANS = 'T' or 't' A'*x = b. +* TRANS = 'T' or 't' A**T*x = b. * -* TRANS = 'C' or 'c' A'*x = b. +* TRANS = 'C' or 'c' A**T*x = b. * * Unchanged on exit. * @@ -225,7 +225,7 @@ END IF ELSE * -* Form x := inv( A' )*x. +* Form x := inv( A**T )*x. * IF (LSAME(UPLO,'U')) THEN KK = 1 diff --git a/BLAS/SRC/strmm.f b/BLAS/SRC/strmm.f index 97fcf6c4..9d76a16d 100644 --- a/BLAS/SRC/strmm.f +++ b/BLAS/SRC/strmm.f @@ -18,7 +18,7 @@ * where alpha is a scalar, B is an m by n matrix, A is a unit, or * non-unit, upper or lower triangular matrix and op( A ) is one of * -* op( A ) = A or op( A ) = A'. +* op( A ) = A or op( A ) = A**T. * * Arguments * ========== @@ -49,9 +49,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. * @@ -234,7 +234,7 @@ END IF ELSE * -* Form B := alpha*A'*B. +* Form B := alpha*A**T*B. * IF (UPPER) THEN DO 110 J = 1,N @@ -300,7 +300,7 @@ END IF ELSE * -* Form B := alpha*B*A'. +* Form B := alpha*B*A**T. * IF (UPPER) THEN DO 260 K = 1,N diff --git a/BLAS/SRC/strmv.f b/BLAS/SRC/strmv.f index 3f1f1ae3..f20cfa2d 100644 --- a/BLAS/SRC/strmv.f +++ b/BLAS/SRC/strmv.f @@ -12,7 +12,7 @@ * * STRMV performs one of the matrix-vector operations * -* x := A*x, or x := A'*x, +* x := A*x, or x := A**T*x, * * where x is an n element vector and A is an n by n unit, or non-unit, * upper or lower triangular matrix. @@ -36,9 +36,9 @@ * * TRANS = 'N' or 'n' x := A*x. * -* TRANS = 'T' or 't' x := A'*x. +* TRANS = 'T' or 't' x := A**T*x. * -* TRANS = 'C' or 'c' x := A'*x. +* TRANS = 'C' or 'c' x := A**T*x. * * Unchanged on exit. * @@ -222,7 +222,7 @@ END IF ELSE * -* Form x := A'*x. +* Form x := A**T*x. * IF (LSAME(UPLO,'U')) THEN IF (INCX.EQ.1) THEN diff --git a/BLAS/SRC/strsm.f b/BLAS/SRC/strsm.f index 4572784e..88880696 100644 --- a/BLAS/SRC/strsm.f +++ b/BLAS/SRC/strsm.f @@ -18,7 +18,7 @@ * where alpha is a scalar, X and B are m by n matrices, A is a unit, or * non-unit, upper or lower triangular matrix and op( A ) is one of * -* op( A ) = A or op( A ) = A'. +* op( A ) = A or op( A ) = A**T. * * The matrix X is overwritten on B. * @@ -51,9 +51,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. * @@ -243,7 +243,7 @@ END IF ELSE * -* Form B := alpha*inv( A' )*B. +* Form B := alpha*inv( A**T )*B. * IF (UPPER) THEN DO 130 J = 1,N @@ -319,7 +319,7 @@ END IF ELSE * -* Form B := alpha*B*inv( A' ). +* Form B := alpha*B*inv( A**T ). * IF (UPPER) THEN DO 310 K = N,1,-1 diff --git a/BLAS/SRC/strsv.f b/BLAS/SRC/strsv.f index 0f820fe4..02420fbf 100644 --- a/BLAS/SRC/strsv.f +++ b/BLAS/SRC/strsv.f @@ -12,7 +12,7 @@ * * STRSV solves one of the systems of equations * -* A*x = b, or A'*x = b, +* A*x = b, or A**T*x = b, * * where b and x are n element vectors and A is an n by n unit, or * non-unit, upper or lower triangular matrix. @@ -39,9 +39,9 @@ * * TRANS = 'N' or 'n' A*x = b. * -* TRANS = 'T' or 't' A'*x = b. +* TRANS = 'T' or 't' A**T*x = b. * -* TRANS = 'C' or 'c' A'*x = b. +* TRANS = 'C' or 'c' A**T*x = b. * * Unchanged on exit. * @@ -223,7 +223,7 @@ END IF ELSE * -* Form x := inv( A' )*x. +* Form x := inv( A**T )*x. * IF (LSAME(UPLO,'U')) THEN IF (INCX.EQ.1) THEN diff --git a/BLAS/SRC/zgbmv.f b/BLAS/SRC/zgbmv.f index 35ca8044..04a29d46 100644 --- a/BLAS/SRC/zgbmv.f +++ b/BLAS/SRC/zgbmv.f @@ -13,9 +13,9 @@ * * ZGBMV 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 @@ * * 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 @@ 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 diff --git a/BLAS/SRC/zgemm.f b/BLAS/SRC/zgemm.f index 198e57f5..f4f98981 100644 --- a/BLAS/SRC/zgemm.f +++ b/BLAS/SRC/zgemm.f @@ -17,7 +17,7 @@ * * 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 @@ * * 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 @@ * * 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 @@ 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 @@ 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 @@ 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 @@ 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 @@ 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 @@ 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 @@ 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 @@ 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 diff --git a/BLAS/SRC/zgemv.f b/BLAS/SRC/zgemv.f index 143be030..bb2ae4fc 100644 --- a/BLAS/SRC/zgemv.f +++ b/BLAS/SRC/zgemv.f @@ -13,9 +13,9 @@ * * ZGEMV 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 @@ * * 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 @@ 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 diff --git a/BLAS/SRC/zgerc.f b/BLAS/SRC/zgerc.f index 62fed2c3..06fdf597 100644 --- a/BLAS/SRC/zgerc.f +++ b/BLAS/SRC/zgerc.f @@ -12,7 +12,7 @@ * * ZGERC 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. diff --git a/BLAS/SRC/zgeru.f b/BLAS/SRC/zgeru.f index 3dd3bd46..dd43b7c5 100644 --- a/BLAS/SRC/zgeru.f +++ b/BLAS/SRC/zgeru.f @@ -12,7 +12,7 @@ * * ZGERU 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. diff --git a/BLAS/SRC/zher.f b/BLAS/SRC/zher.f index 848a3702..2f82e019 100644 --- a/BLAS/SRC/zher.f +++ b/BLAS/SRC/zher.f @@ -13,7 +13,7 @@ * * ZHER 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. diff --git a/BLAS/SRC/zher2.f b/BLAS/SRC/zher2.f index c0b6a306..a2cc1164 100644 --- a/BLAS/SRC/zher2.f +++ b/BLAS/SRC/zher2.f @@ -13,7 +13,7 @@ * * ZHER2 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. diff --git a/BLAS/SRC/zher2k.f b/BLAS/SRC/zher2k.f index 0dfb32bf..e20f2040 100644 --- a/BLAS/SRC/zher2k.f +++ b/BLAS/SRC/zher2k.f @@ -14,11 +14,11 @@ * * ZHER2K 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 @@ * 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 @@ * 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 @@ 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 diff --git a/BLAS/SRC/zherk.f b/BLAS/SRC/zherk.f index ed51aa18..d0684707 100644 --- a/BLAS/SRC/zherk.f +++ b/BLAS/SRC/zherk.f @@ -13,11 +13,11 @@ * * ZHERK 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 @@ * 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 @@ * 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 @@ 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 diff --git a/BLAS/SRC/zhpr.f b/BLAS/SRC/zhpr.f index 40efbc7d..577d52cc 100644 --- a/BLAS/SRC/zhpr.f +++ b/BLAS/SRC/zhpr.f @@ -13,7 +13,7 @@ * * ZHPR 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. diff --git a/BLAS/SRC/zhpr2.f b/BLAS/SRC/zhpr2.f index 99977462..024fc88f 100644 --- a/BLAS/SRC/zhpr2.f +++ b/BLAS/SRC/zhpr2.f @@ -13,7 +13,7 @@ * * ZHPR2 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, supplied in packed form. diff --git a/BLAS/SRC/zsyr2k.f b/BLAS/SRC/zsyr2k.f index 94ed2e36..808b7ede 100644 --- a/BLAS/SRC/zsyr2k.f +++ b/BLAS/SRC/zsyr2k.f @@ -13,11 +13,11 @@ * * ZSYR2K performs one of the symmetric rank 2k operations * -* C := alpha*A*B' + alpha*B*A' + beta*C, +* C := alpha*A*B**T + alpha*B*A**T + beta*C, * * or * -* C := alpha*A'*B + alpha*B'*A + beta*C, +* C := alpha*A**T*B + alpha*B**T*A + beta*C, * * where alpha and beta are scalars, C is an n by n symmetric matrix * and A and B are n by k matrices in the first case and k by n @@ -43,10 +43,10 @@ * On entry, TRANS specifies the operation to be performed as * follows: * -* TRANS = 'N' or 'n' C := alpha*A*B' + alpha*B*A' + +* TRANS = 'N' or 'n' C := alpha*A*B**T + alpha*B*A**T + * beta*C. * -* TRANS = 'T' or 't' C := alpha*A'*B + alpha*B'*A + +* TRANS = 'T' or 't' C := alpha*A**T*B + alpha*B**T*A + * beta*C. * * Unchanged on exit. @@ -231,7 +231,7 @@ * IF (LSAME(TRANS,'N')) THEN * -* Form C := alpha*A*B' + alpha*B*A' + C. +* Form C := alpha*A*B**T + alpha*B*A**T + C. * IF (UPPER) THEN DO 130 J = 1,N @@ -280,7 +280,7 @@ END IF ELSE * -* Form C := alpha*A'*B + alpha*B'*A + C. +* Form C := alpha*A**T*B + alpha*B**T*A + C. * IF (UPPER) THEN DO 210 J = 1,N diff --git a/BLAS/SRC/zsyrk.f b/BLAS/SRC/zsyrk.f index 486a1744..ecc1285d 100644 --- a/BLAS/SRC/zsyrk.f +++ b/BLAS/SRC/zsyrk.f @@ -13,11 +13,11 @@ * * ZSYRK performs one of the symmetric rank k operations * -* C := alpha*A*A' + beta*C, +* C := alpha*A*A**T + beta*C, * * or * -* C := alpha*A'*A + beta*C, +* C := alpha*A**T*A + beta*C, * * where alpha and beta are scalars, C is an n by n symmetric matrix * and A is an n by k matrix in the first case and a k by n matrix @@ -43,9 +43,9 @@ * On entry, TRANS specifies the operation to be performed as * follows: * -* TRANS = 'N' or 'n' C := alpha*A*A' + beta*C. +* TRANS = 'N' or 'n' C := alpha*A*A**T + beta*C. * -* TRANS = 'T' or 't' C := alpha*A'*A + beta*C. +* TRANS = 'T' or 't' C := alpha*A**T*A + beta*C. * * Unchanged on exit. * @@ -212,7 +212,7 @@ * IF (LSAME(TRANS,'N')) THEN * -* Form C := alpha*A*A' + beta*C. +* Form C := alpha*A*A**T + beta*C. * IF (UPPER) THEN DO 130 J = 1,N @@ -257,7 +257,7 @@ END IF ELSE * -* Form C := alpha*A'*A + beta*C. +* Form C := alpha*A**T*A + beta*C. * IF (UPPER) THEN DO 210 J = 1,N diff --git a/BLAS/SRC/ztbmv.f b/BLAS/SRC/ztbmv.f index 23d5f4c4..c7fdb8f3 100644 --- a/BLAS/SRC/ztbmv.f +++ b/BLAS/SRC/ztbmv.f @@ -12,7 +12,7 @@ * * ZTBMV performs one of the matrix-vector operations * -* x := A*x, or x := A'*x, or x := conjg( A' )*x, +* x := A*x, or x := A**T*x, or x := A**H*x, * * where x is an n element vector and A is an n by n unit, or non-unit, * upper or lower triangular band matrix, with ( k + 1 ) diagonals. @@ -36,9 +36,9 @@ * * TRANS = 'N' or 'n' x := A*x. * -* TRANS = 'T' or 't' x := A'*x. +* TRANS = 'T' or 't' x := A**T*x. * -* TRANS = 'C' or 'c' x := conjg( A' )*x. +* TRANS = 'C' or 'c' x := A**H*x. * * Unchanged on exit. * @@ -269,7 +269,7 @@ END IF ELSE * -* Form x := A'*x or x := conjg( A' )*x. +* Form x := A**T*x or x := A**H*x. * IF (LSAME(UPLO,'U')) THEN KPLUS1 = K + 1 diff --git a/BLAS/SRC/ztbsv.f b/BLAS/SRC/ztbsv.f index 42b234a7..9f3bf510 100644 --- a/BLAS/SRC/ztbsv.f +++ b/BLAS/SRC/ztbsv.f @@ -12,7 +12,7 @@ * * ZTBSV solves one of the systems of equations * -* A*x = b, or A'*x = b, or conjg( A' )*x = b, +* A*x = b, or A**T*x = b, or A**H*x = b, * * where b and x are n element vectors and A is an n by n unit, or * non-unit, upper or lower triangular band matrix, with ( k + 1 ) @@ -40,9 +40,9 @@ * * TRANS = 'N' or 'n' A*x = b. * -* TRANS = 'T' or 't' A'*x = b. +* TRANS = 'T' or 't' A**T*x = b. * -* TRANS = 'C' or 'c' conjg( A' )*x = b. +* TRANS = 'C' or 'c' A**H*x = b. * * Unchanged on exit. * @@ -272,7 +272,7 @@ END IF ELSE * -* Form x := inv( A' )*x or x := inv( conjg( A') )*x. +* Form x := inv( A**T )*x or x := inv( A**H )*x. * IF (LSAME(UPLO,'U')) THEN KPLUS1 = K + 1 diff --git a/BLAS/SRC/ztpmv.f b/BLAS/SRC/ztpmv.f index 5832692a..bf74ec04 100644 --- a/BLAS/SRC/ztpmv.f +++ b/BLAS/SRC/ztpmv.f @@ -12,7 +12,7 @@ * * ZTPMV performs one of the matrix-vector operations * -* x := A*x, or x := A'*x, or x := conjg( A' )*x, +* x := A*x, or x := A**T*x, or x := A**H*x, * * where x is an n element vector and A is an n by n unit, or non-unit, * upper or lower triangular matrix, supplied in packed form. @@ -36,9 +36,9 @@ * * TRANS = 'N' or 'n' x := A*x. * -* TRANS = 'T' or 't' x := A'*x. +* TRANS = 'T' or 't' x := A**T*x. * -* TRANS = 'C' or 'c' x := conjg( A' )*x. +* TRANS = 'C' or 'c' x := A**H*x. * * Unchanged on exit. * @@ -228,7 +228,7 @@ END IF ELSE * -* Form x := A'*x or x := conjg( A' )*x. +* Form x := A**T*x or x := A**H*x. * IF (LSAME(UPLO,'U')) THEN KK = (N* (N+1))/2 diff --git a/BLAS/SRC/ztpsv.f b/BLAS/SRC/ztpsv.f index b56e1d8c..d8aae2ed 100644 --- a/BLAS/SRC/ztpsv.f +++ b/BLAS/SRC/ztpsv.f @@ -12,7 +12,7 @@ * * ZTPSV solves one of the systems of equations * -* A*x = b, or A'*x = b, or conjg( A' )*x = b, +* A*x = b, or A**T*x = b, or A**H*x = b, * * where b and x are n element vectors and A is an n by n unit, or * non-unit, upper or lower triangular matrix, supplied in packed form. @@ -39,9 +39,9 @@ * * TRANS = 'N' or 'n' A*x = b. * -* TRANS = 'T' or 't' A'*x = b. +* TRANS = 'T' or 't' A**T*x = b. * -* TRANS = 'C' or 'c' conjg( A' )*x = b. +* TRANS = 'C' or 'c' A**H*x = b. * * Unchanged on exit. * @@ -229,7 +229,7 @@ END IF ELSE * -* Form x := inv( A' )*x or x := inv( conjg( A' ) )*x. +* Form x := inv( A**T )*x or x := inv( A**H )*x. * IF (LSAME(UPLO,'U')) THEN KK = 1 diff --git a/BLAS/SRC/ztrmm.f b/BLAS/SRC/ztrmm.f index 56f3360f..9ea6b9b9 100644 --- a/BLAS/SRC/ztrmm.f +++ b/BLAS/SRC/ztrmm.f @@ -18,7 +18,7 @@ * where alpha is a scalar, B is an m by n matrix, A is a unit, or * non-unit, upper or lower triangular matrix and op( A ) is one of * -* op( A ) = A or op( A ) = A' or op( A ) = conjg( A' ). +* op( A ) = A or op( A ) = A**T or op( A ) = A**H. * * Arguments * ========== @@ -49,9 +49,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 ) = conjg( A' ). +* TRANSA = 'C' or 'c' op( A ) = A**H. * * Unchanged on exit. * @@ -237,7 +237,7 @@ END IF ELSE * -* Form B := alpha*A'*B or B := alpha*conjg( A' )*B. +* Form B := alpha*A**T*B or B := alpha*A**H*B. * IF (UPPER) THEN DO 120 J = 1,N @@ -317,7 +317,7 @@ END IF ELSE * -* Form B := alpha*B*A' or B := alpha*B*conjg( A' ). +* Form B := alpha*B*A**T or B := alpha*B*A**H. * IF (UPPER) THEN DO 280 K = 1,N diff --git a/BLAS/SRC/ztrmv.f b/BLAS/SRC/ztrmv.f index 8a9e7ea5..50708e1a 100644 --- a/BLAS/SRC/ztrmv.f +++ b/BLAS/SRC/ztrmv.f @@ -12,7 +12,7 @@ * * ZTRMV performs one of the matrix-vector operations * -* x := A*x, or x := A'*x, or x := conjg( A' )*x, +* x := A*x, or x := A**T*x, or x := A**H*x, * * where x is an n element vector and A is an n by n unit, or non-unit, * upper or lower triangular matrix. @@ -36,9 +36,9 @@ * * TRANS = 'N' or 'n' x := A*x. * -* TRANS = 'T' or 't' x := A'*x. +* TRANS = 'T' or 't' x := A**T*x. * -* TRANS = 'C' or 'c' x := conjg( A' )*x. +* TRANS = 'C' or 'c' x := A**H*x. * * Unchanged on exit. * @@ -223,7 +223,7 @@ END IF ELSE * -* Form x := A'*x or x := conjg( A' )*x. +* Form x := A**T*x or x := A**H*x. * IF (LSAME(UPLO,'U')) THEN IF (INCX.EQ.1) THEN diff --git a/BLAS/SRC/ztrsm.f b/BLAS/SRC/ztrsm.f index 7d81de51..e8bee0eb 100644 --- a/BLAS/SRC/ztrsm.f +++ b/BLAS/SRC/ztrsm.f @@ -18,7 +18,7 @@ * where alpha is a scalar, X and B are m by n matrices, A is a unit, or * non-unit, upper or lower triangular matrix and op( A ) is one of * -* op( A ) = A or op( A ) = A' or op( A ) = conjg( A' ). +* op( A ) = A or op( A ) = A**T or op( A ) = A**H. * * The matrix X is overwritten on B. * @@ -51,9 +51,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 ) = conjg( A' ). +* TRANSA = 'C' or 'c' op( A ) = A**H. * * Unchanged on exit. * @@ -245,8 +245,8 @@ END IF ELSE * -* Form B := alpha*inv( A' )*B -* or B := alpha*inv( conjg( A' ) )*B. +* Form B := alpha*inv( A**T )*B +* or B := alpha*inv( A**H )*B. * IF (UPPER) THEN DO 140 J = 1,N @@ -336,8 +336,8 @@ END IF ELSE * -* Form B := alpha*B*inv( A' ) -* or B := alpha*B*inv( conjg( A' ) ). +* Form B := alpha*B*inv( A**T ) +* or B := alpha*B*inv( A**H ). * IF (UPPER) THEN DO 330 K = N,1,-1 diff --git a/BLAS/SRC/ztrsv.f b/BLAS/SRC/ztrsv.f index c48d9945..d5d54837 100644 --- a/BLAS/SRC/ztrsv.f +++ b/BLAS/SRC/ztrsv.f @@ -12,7 +12,7 @@ * * ZTRSV solves one of the systems of equations * -* A*x = b, or A'*x = b, or conjg( A' )*x = b, +* A*x = b, or A**T*x = b, or A**H*x = b, * * where b and x are n element vectors and A is an n by n unit, or * non-unit, upper or lower triangular matrix. @@ -39,9 +39,9 @@ * * TRANS = 'N' or 'n' A*x = b. * -* TRANS = 'T' or 't' A'*x = b. +* TRANS = 'T' or 't' A**T*x = b. * -* TRANS = 'C' or 'c' conjg( A' )*x = b. +* TRANS = 'C' or 'c' A**H*x = b. * * Unchanged on exit. * @@ -224,7 +224,7 @@ END IF ELSE * -* Form x := inv( A' )*x or x := inv( conjg( A' ) )*x. +* Form x := inv( A**T )*x or x := inv( A**H )*x. * IF (LSAME(UPLO,'U')) THEN IF (INCX.EQ.1) THEN |