unify comments: expression for residual

This resolves discussion in PR Reference-LAPACK#582.

Author: jip

TESTING/LIN/ctbt02.f

@@ -30,10 +30,10 @@
 *>
 *> CTBT02 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
-*> number of right hand sides of
-*>    norm(op(A)*X - B) / ( norm(op(A)) * norm(X) * EPS ),
-*> where op(A) = A, A**T, or A**H, and EPS is the machine epsilon.
+*> triangular band matrix. The test ratio is the maximum over
+*>    norm(b - op(A)*x) / ( ||op(A)||_1 * norm(x) * EPS ),
+*> where op(A) = A, A**T, or A**H, b is the column of B, x is the
+*> solution vector, and EPS is the machine epsilon.
 *> \endverbatim
 *
 *  Arguments:

TESTING/LIN/ctpt02.f

@@ -30,9 +30,10 @@
 *> CTPT02 computes the residual for the computed solution to a
 *> 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(op(A)*X - B) / ( norm(op(A)) * norm(X) * EPS ),
-*> where op(A) = A, A**T, or A**H, and EPS is the machine epsilon.
+*> the maximum over
+*>    norm(b - op(A)*x) / ( ||op(A)||_1 * norm(x) * EPS ),
+*> where op(A) = A, A**T, or A**H, b is the column of B, x is the
+*> solution vector, and EPS is the machine epsilon.
 *> \endverbatim
 *
 *  Arguments:

TESTING/LIN/ctrt02.f

@@ -30,10 +30,10 @@
 *>
 *> CTRT02 computes the residual for the computed solution to a
 *> 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(op(A)*X - B) / ( norm(op(A)) * norm(X) * EPS ),
-*> where op(A) = A, A**T, or A**H, and EPS is the machine epsilon.
+*> triangular matrix. The test ratio is the maximum over
+*>    norm(b - op(A)*x) / ( ||op(A)||_1 * norm(x) * EPS ),
+*> where op(A) = A, A**T, or A**H, b is the column of B, x is the
+*> solution vector, and EPS is the machine epsilon.
 *> \endverbatim
 *
 *  Arguments:

TESTING/LIN/dtbt02.f

@@ -29,10 +29,10 @@
 *>
 *> DTBT02 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
-*> number of right hand sides of
-*>    norm(op(A)*X - B) / ( norm(op(A)) * norm(X) * EPS ),
-*> where op(A) = A, A**T, or A**H and EPS is the machine epsilon.
+*> triangular band matrix. The test ratio is the maximum over
+*>    norm(b - op(A)*x) / ( ||op(A)||_1 * norm(x) * EPS ),
+*> where op(A) = A or A**T, b is the column of B, x is the solution
+*> vector, and EPS is the machine epsilon.
 *> The norm used is the 1-norm.
 *> \endverbatim
 *

TESTING/LIN/dtpt02.f

@@ -29,9 +29,10 @@
 *> DTPT02 computes the residual for the computed solution to a
 *> 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(op(A)*X - B) / ( norm(op(A)) * norm(X) * EPS ),
-*> where op(A) = A or A**T, and EPS is the machine epsilon.
+*> the maximum over
+*>    norm(b - op(A)*x) / ( ||op(A)||_1 * norm(x) * EPS ),
+*> where op(A) = A or A**T, b is the column of B, x is the solution
+*> vector, and EPS is the machine epsilon.
 *> The norm used is the 1-norm.
 *> \endverbatim
 *

TESTING/LIN/dtrt02.f

@@ -29,10 +29,10 @@
 *>
 *> DTRT02 computes the residual for the computed solution to a
 *> 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(op(A)*X - B) / ( norm(op(A)) * norm(X) * EPS ),
-*> where op(A) = A or A**T, and EPS is the machine epsilon.
+*> triangular matrix. The test ratio is the maximum over
+*>    norm(b - op(A)*x) / ( ||op(A)||_1 * norm(x) * EPS ),
+*> where op(A) = A or A**T, b is the column of B, x is the solution
+*> vector, and EPS is the machine epsilon.
 *> The norm used is the 1-norm.
 *> \endverbatim
 *

TESTING/LIN/stbt02.f

@@ -29,10 +29,10 @@
 *>
 *> 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
-*> number of right hand sides of
-*>    norm(op(A)*X - B) / ( norm(op(A)) * norm(X) * EPS ),
-*> where op(A) = A or A**T, and EPS is the machine epsilon.
+*> triangular band matrix. The test ratio is the maximum over
+*>    norm(b - op(A)*x) / ( ||op(A)||_1 * norm(x) * EPS ),
+*> where op(A) = A or A**T, b is the column of B, x is the solution
+*> vector, and EPS is the machine epsilon.
 *> The norm used is the 1-norm.
 *> \endverbatim
 *

TESTING/LIN/stpt02.f

@@ -29,9 +29,10 @@
 *> STPT02 computes the residual for the computed solution to a
 *> 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(op(A)*X - B) / ( norm(op(A)) * norm(X) * EPS ),
-*> where op(A) = A or A**T, and EPS is the machine epsilon.
+*> the maximum over
+*>    norm(b - op(A)*x) / ( ||op(A)||_1 * norm(x) * EPS ),
+*> where op(A) = A or A**T, b is the column of B, x is the solution
+*> vector, and EPS is the machine epsilon.
 *> The norm used is the 1-norm.
 *> \endverbatim
 *

TESTING/LIN/strt02.f

@@ -29,10 +29,10 @@
 *>
 *> STRT02 computes the residual for the computed solution to a
 *> 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(op(A)*X - B) / ( norm(op(A)) * norm(X) * EPS ),
-*> where op(A) = A or A**T, and EPS is the machine epsilon.
+*> triangular matrix. The test ratio is the maximum over
+*>    norm(b - op(A)*x) / ( ||op(A)||_1 * norm(x) * EPS ),
+*> where op(A) = A or A**T, b is the column of B, x is the solution
+*> vector, and EPS is the machine epsilon.
 *> The norm used is the 1-norm.
 *> \endverbatim
 *

TESTING/LIN/ztbt02.f

@@ -30,10 +30,10 @@
 *>
 *> ZTBT02 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
-*> number of right hand sides of
-*>    norm(op(A)*X - B) / ( norm(op(A)) * norm(X) * EPS ),
-*> where op(A) = A, A**T, or A**H, and EPS is the machine epsilon.
+*> triangular band matrix. The test ratio is the maximum over
+*>    norm(b - op(A)*x) / ( ||op(A)||_1 * norm(x) * EPS ),
+*> where op(A) = A, A**T, or A**H, b is the column of B, x is the
+*> solution vector, and EPS is the machine epsilon.
 *> \endverbatim
 *
 *  Arguments: