unify comments: homogenize notation for transpose (**T) and conjugate transpose (**H)

This commit completes changes made in:
* 457afcf 2021-06-03

Author: jip

TESTING/LIN/ctbt02.f

@@ -29,13 +29,11 @@
 *> \verbatim
 *>
 *> CTBT02 computes the residual for the computed solution to a
-*> triangular system of linear equations  A*x = b,  A**T *x = b,  or
-*> A**H *x = b  when A is a triangular band matrix.  Here A**T denotes
-*> the transpose of A, A**H denotes the conjugate transpose of A, and
-*> x and b are N by NRHS matrices.  The test ratio is the maximum over
-*> the number of right hand sides of
+*> triangular system of linear equations op(A)*X = B, when A is a
+*> triangular band matrix. The test ratio is the maximum over the
+*> number of right hand sides of
 *>    norm(b - op(A)*x) / ( norm(op(A)) * norm(x) * EPS ),
-*> where op(A) denotes A, A**T, or A**H, and EPS is the machine epsilon.
+*> where op(A) = A, A**T, or A**H, and EPS is the machine epsilon.
 *> \endverbatim
 *
 *  Arguments:

TESTING/LIN/ctpt02.f

@@ -28,14 +28,11 @@
 *> \verbatim
 *>
 *> CTPT02 computes the residual for the computed solution to a
-*> triangular system of linear equations  A*x = b,  A**T *x = b,  or
-*> A**H *x = b, when the triangular matrix A is stored in packed format.
-*> Here A**T denotes the transpose of A, A**H denotes the conjugate
-*> transpose of A, and x and b are N by NRHS matrices.  The test ratio
-*> is the maximum over the number of right hand sides of
+*> triangular system of linear equations op(A)*X = B, when the
+*> triangular matrix A is stored in packed format. The test ratio is
 *> the maximum over the number of right hand sides of
 *>    norm(b - op(A)*x) / ( norm(op(A)) * norm(x) * EPS ),
-*> where op(A) denotes A, A**T, or A**H, and EPS is the machine epsilon.
+*> where op(A) = A, A**T, or A**H, and EPS is the machine epsilon.
 *> \endverbatim
 *
 *  Arguments:
@@ -191,7 +188,7 @@ SUBROUTINE CTPT02( UPLO, TRANS, DIAG, N, NRHS, AP, X, LDX, B, LDB,
          RETURN
       END IF
 *
-*     Compute the 1-norm of A or A**H.
+*     Compute the 1-norm of op(A).
 *
       IF( LSAME( TRANS, 'N' ) ) THEN
          ANORM = CLANTP( '1', UPLO, DIAG, N, AP, RWORK )

TESTING/LIN/ctrt02.f

@@ -29,13 +29,11 @@
 *> \verbatim
 *>
 *> CTRT02 computes the residual for the computed solution to a
-*> triangular system of linear equations  A*x = b,  A**T *x = b,
-*> or A**H *x = b.  Here A is a triangular matrix, A**T is the transpose
-*> of A, A**H is the conjugate transpose of A, and x and b are N by NRHS
-*> matrices.  The test ratio is the maximum over the number of right
-*> hand sides of
+*> triangular system of linear equations op(A)*X = B, where A is a
+*> triangular matrix. The test ratio is the maximum over the number of
+*> right hand sides of
 *>    norm(b - op(A)*x) / ( norm(op(A)) * norm(x) * EPS ),
-*> where op(A) denotes A, A**T, or A**H, and EPS is the machine epsilon.
+*> where op(A) = A, A**T, or A**H, and EPS is the machine epsilon.
 *> \endverbatim
 *
 *  Arguments:
@@ -200,7 +198,7 @@ SUBROUTINE CTRT02( UPLO, TRANS, DIAG, N, NRHS, A, LDA, X, LDX, B,
          RETURN
       END IF
 *
-*     Compute the 1-norm of A or A**H.
+*     Compute the 1-norm of op(A).
 *
       IF( LSAME( TRANS, 'N' ) ) THEN
          ANORM = CLANTR( '1', UPLO, DIAG, N, N, A, LDA, RWORK )

TESTING/LIN/dtbt02.f

@@ -28,12 +28,11 @@
 *> \verbatim
 *>
 *> DTBT02 computes the residual for the computed solution to a
-*> triangular system of linear equations  A*x = b  or  A' *x = b when
-*> A is a triangular band matrix.  Here A' is the transpose of A and
-*> x and b are N by NRHS matrices.  The test ratio is the maximum over
-*> the number of right hand sides of
+*> triangular system of linear equations op(A)*X = B, when A is a
+*> triangular band matrix. The test ratio is the maximum over the
+*> number of right hand sides of
 *>    norm(b - op(A)*x) / ( norm(op(A)) * norm(x) * EPS ),
-*> where op(A) denotes A or A' and EPS is the machine epsilon.
+*> where op(A) = A, A**T, or A**H and EPS is the machine epsilon.
 *> The norm used is the 1-norm.
 *> \endverbatim
 *

TESTING/LIN/dtpt02.f

@@ -27,12 +27,11 @@
 *> \verbatim
 *>
 *> DTPT02 computes the residual for the computed solution to a
-*> triangular system of linear equations  A*x = b  or  A'*x = b  when
-*> the triangular matrix A is stored in packed format.  Here A' is the
-*> transpose of A and x and b are N by NRHS matrices.  The test ratio is
+*> triangular system of linear equations op(A)*X = B, when the
+*> triangular matrix A is stored in packed format. The test ratio is
 *> the maximum over the number of right hand sides of
 *>    norm(b - op(A)*x) / ( norm(op(A)) * norm(x) * EPS ),
-*> where op(A) denotes A or A' and EPS is the machine epsilon.
+*> where op(A) = A or A**T, and EPS is the machine epsilon.
 *> The norm used is the 1-norm.
 *> \endverbatim
 *
@@ -51,9 +50,9 @@
 *> \verbatim
 *>          TRANS is CHARACTER*1
 *>          Specifies the operation applied to A.
-*>          = 'N':  A *x = b  (No transpose)
-*>          = 'T':  A'*x = b  (Transpose)
-*>          = 'C':  A'*x = b  (Conjugate transpose = Transpose)
+*>          = 'N':  A    * X = B  (No transpose)
+*>          = 'T':  A**T * X = B  (Transpose)
+*>          = 'C':  A**H * X = B  (Conjugate transpose = Transpose)
 *> \endverbatim
 *>
 *> \param[in] DIAG
@@ -183,7 +182,7 @@ SUBROUTINE DTPT02( UPLO, TRANS, DIAG, N, NRHS, AP, X, LDX, B, LDB,
          RETURN
       END IF
 *
-*     Compute the 1-norm of A or A'.
+*     Compute the 1-norm of op(A).
 *
       IF( LSAME( TRANS, 'N' ) ) THEN
          ANORM = DLANTP( '1', UPLO, DIAG, N, AP, WORK )

TESTING/LIN/dtrt02.f

@@ -28,12 +28,11 @@
 *> \verbatim
 *>
 *> DTRT02 computes the residual for the computed solution to a
-*> triangular system of linear equations  A*x = b  or  A'*x = b.
-*> Here A is a triangular matrix, A' is the transpose of A, and x and b
-*> are N by NRHS matrices.  The test ratio is the maximum over the
-*> number of right hand sides of
+*> triangular system of linear equations op(A)*X = B, where A is a
+*> triangular matrix. The test ratio is the maximum over the number of
+*> right hand sides of
 *>    norm(b - op(A)*x) / ( norm(op(A)) * norm(x) * EPS ),
-*> where op(A) denotes A or A' and EPS is the machine epsilon.
+*> where op(A) = A or A**T, and EPS is the machine epsilon.
 *> The norm used is the 1-norm.
 *> \endverbatim
 *
@@ -52,9 +51,9 @@
 *> \verbatim
 *>          TRANS is CHARACTER*1
 *>          Specifies the operation applied to A.
-*>          = 'N':  A *x = b  (No transpose)
-*>          = 'T':  A'*x = b  (Transpose)
-*>          = 'C':  A'*x = b  (Conjugate transpose = Transpose)
+*>          = 'N':  A    * X = B  (No transpose)
+*>          = 'T':  A**T * X = B  (Transpose)
+*>          = 'C':  A**H * X = B  (Conjugate transpose = Transpose)
 *> \endverbatim
 *>
 *> \param[in] DIAG
@@ -193,7 +192,7 @@ SUBROUTINE DTRT02( UPLO, TRANS, DIAG, N, NRHS, A, LDA, X, LDX, B,
          RETURN
       END IF
 *
-*     Compute the 1-norm of A or A'.
+*     Compute the 1-norm of op(A).
 *
       IF( LSAME( TRANS, 'N' ) ) THEN
          ANORM = DLANTR( '1', UPLO, DIAG, N, N, A, LDA, WORK )

TESTING/LIN/stbt02.f

@@ -28,11 +28,11 @@
 *> \verbatim
 *>
 *> STBT02 computes the residual for the computed solution to a
-*> triangular system of linear equations op(A)*X = B when A is a
-*> triangular band matrix.  The test ratio is the maximum over the
+*> triangular system of linear equations op(A)*X = B, when A is a
+*> triangular band matrix. The test ratio is the maximum over the
 *> number of right hand sides of
 *>    norm(b - op(A)*x) / ( norm(op(A)) * norm(x) * EPS ),
-*> where op(A) denotes A or A' and EPS is the machine epsilon.
+*> where op(A) = A or A**T, and EPS is the machine epsilon.
 *> The norm used is the 1-norm.
 *> \endverbatim
 *

TESTING/LIN/stpt02.f

@@ -27,12 +27,11 @@
 *> \verbatim
 *>
 *> STPT02 computes the residual for the computed solution to a
-*> triangular system of linear equations  A*x = b  or  A'*x = b  when
-*> the triangular matrix A is stored in packed format.  Here A' is the
-*> transpose of A and x and b are N by NRHS matrices.  The test ratio is
+*> triangular system of linear equations op(A)*X = B, when the
+*> triangular matrix A is stored in packed format. The test ratio is
 *> the maximum over the number of right hand sides of
 *>    norm(b - op(A)*x) / ( norm(op(A)) * norm(x) * EPS ),
-*> where op(A) denotes A or A' and EPS is the machine epsilon.
+*> where op(A) = A or A**T, and EPS is the machine epsilon.
 *> The norm used is the 1-norm.
 *> \endverbatim
 *
@@ -51,9 +50,9 @@
 *> \verbatim
 *>          TRANS is CHARACTER*1
 *>          Specifies the operation applied to A.
-*>          = 'N':  A *x = b  (No transpose)
-*>          = 'T':  A'*x = b  (Transpose)
-*>          = 'C':  A'*x = b  (Conjugate transpose = Transpose)
+*>          = 'N':  A    * X = B  (No transpose)
+*>          = 'T':  A**T * X = B  (Transpose)
+*>          = 'C':  A**H * X = B  (Conjugate transpose = Transpose)
 *> \endverbatim
 *>
 *> \param[in] DIAG
@@ -183,7 +182,7 @@ SUBROUTINE STPT02( UPLO, TRANS, DIAG, N, NRHS, AP, X, LDX, B, LDB,
          RETURN
       END IF
 *
-*     Compute the 1-norm of A or A'.
+*     Compute the 1-norm of op(A).
 *
       IF( LSAME( TRANS, 'N' ) ) THEN
          ANORM = SLANTP( '1', UPLO, DIAG, N, AP, WORK )

TESTING/LIN/strt02.f

@@ -28,12 +28,11 @@
 *> \verbatim
 *>
 *> STRT02 computes the residual for the computed solution to a
-*> triangular system of linear equations  A*x = b  or  A'*x = b.
-*> Here A is a triangular matrix, A' is the transpose of A, and x and b
-*> are N by NRHS matrices.  The test ratio is the maximum over the
-*> number of right hand sides of
+*> triangular system of linear equations op(A)*X = B, where A is a
+*> triangular matrix. The test ratio is the maximum over the number of
+*> right hand sides of
 *>    norm(b - op(A)*x) / ( norm(op(A)) * norm(x) * EPS ),
-*> where op(A) denotes A or A' and EPS is the machine epsilon.
+*> where op(A) = A or A**T, and EPS is the machine epsilon.
 *> The norm used is the 1-norm.
 *> \endverbatim
 *
@@ -52,9 +51,9 @@
 *> \verbatim
 *>          TRANS is CHARACTER*1
 *>          Specifies the operation applied to A.
-*>          = 'N':  A *x = b  (No transpose)
-*>          = 'T':  A'*x = b  (Transpose)
-*>          = 'C':  A'*x = b  (Conjugate transpose = Transpose)
+*>          = 'N':  A    * X = B  (No transpose)
+*>          = 'T':  A**T * X = B  (Transpose)
+*>          = 'C':  A**H * X = B  (Conjugate transpose = Transpose)
 *> \endverbatim
 *>
 *> \param[in] DIAG
@@ -193,7 +192,7 @@ SUBROUTINE STRT02( UPLO, TRANS, DIAG, N, NRHS, A, LDA, X, LDX, B,
          RETURN
       END IF
 *
-*     Compute the 1-norm of A or A'.
+*     Compute the 1-norm of op(A).
 *
       IF( LSAME( TRANS, 'N' ) ) THEN
          ANORM = SLANTR( '1', UPLO, DIAG, N, N, A, LDA, WORK )

TESTING/LIN/ztbt02.f

@@ -29,13 +29,11 @@
 *> \verbatim
 *>
 *> ZTBT02 computes the residual for the computed solution to a
-*> triangular system of linear equations  A*x = b,  A**T *x = b,  or
-*> A**H *x = b  when A is a triangular band matrix.  Here A**T denotes
-*> the transpose of A, A**H denotes the conjugate transpose of A, and
-*> x and b are N by NRHS matrices.  The test ratio is the maximum over
-*> the number of right hand sides of
+*> triangular system of linear equations op(A)*X = B, when A is a
+*> triangular band matrix. The test ratio is the maximum over the
+*> number of right hand sides of
 *>    norm(b - op(A)*x) / ( norm(op(A)) * norm(x) * EPS ),
-*> where op(A) denotes A, A**T, or A**H, and EPS is the machine epsilon.
+*> where op(A) = A, A**T, or A**H, and EPS is the machine epsilon.
 *> \endverbatim
 *
 *  Arguments: