diff --git a/SRC/cla_gbrfsx_extended.f b/SRC/cla_gbrfsx_extended.f
index 612598eb1c..a54be9c5fd 100644
--- a/SRC/cla_gbrfsx_extended.f
+++ b/SRC/cla_gbrfsx_extended.f
@@ -54,7 +54,7 @@
 *> In addition to normwise error bound, the code provides maximum
 *> componentwise error bound if possible. See comments for ERR_BNDS_NORM
 *> and ERR_BNDS_COMP for details of the error bounds. Note that this
-*> subroutine is only resonsible for setting the second fields of
+*> subroutine is only responsible for setting the second fields of
 *> ERR_BNDS_NORM and ERR_BNDS_COMP.
 *> \endverbatim
 *
diff --git a/SRC/cla_gerfsx_extended.f b/SRC/cla_gerfsx_extended.f
index 03beefe354..4b10311017 100644
--- a/SRC/cla_gerfsx_extended.f
+++ b/SRC/cla_gerfsx_extended.f
@@ -54,7 +54,7 @@
 *> In addition to normwise error bound, the code provides maximum
 *> componentwise error bound if possible. See comments for ERRS_N
 *> and ERRS_C for details of the error bounds. Note that this
-*> subroutine is only resonsible for setting the second fields of
+*> subroutine is only responsible for setting the second fields of
 *> ERRS_N and ERRS_C.
 *> \endverbatim
 *
diff --git a/SRC/cla_herfsx_extended.f b/SRC/cla_herfsx_extended.f
index 96ff805e91..6d007ef580 100644
--- a/SRC/cla_herfsx_extended.f
+++ b/SRC/cla_herfsx_extended.f
@@ -55,7 +55,7 @@
 *> In addition to normwise error bound, the code provides maximum
 *> componentwise error bound if possible. See comments for ERR_BNDS_NORM
 *> and ERR_BNDS_COMP for details of the error bounds. Note that this
-*> subroutine is only resonsible for setting the second fields of
+*> subroutine is only responsible for setting the second fields of
 *> ERR_BNDS_NORM and ERR_BNDS_COMP.
 *> \endverbatim
 *
diff --git a/SRC/cla_porfsx_extended.f b/SRC/cla_porfsx_extended.f
index 49c5324c53..9ced9b1b93 100644
--- a/SRC/cla_porfsx_extended.f
+++ b/SRC/cla_porfsx_extended.f
@@ -54,7 +54,7 @@
 *> In addition to normwise error bound, the code provides maximum
 *> componentwise error bound if possible. See comments for ERR_BNDS_NORM
 *> and ERR_BNDS_COMP for details of the error bounds. Note that this
-*> subroutine is only resonsible for setting the second fields of
+*> subroutine is only responsible for setting the second fields of
 *> ERR_BNDS_NORM and ERR_BNDS_COMP.
 *> \endverbatim
 *
diff --git a/SRC/cla_syrfsx_extended.f b/SRC/cla_syrfsx_extended.f
index 498759ebfe..4fe538a989 100644
--- a/SRC/cla_syrfsx_extended.f
+++ b/SRC/cla_syrfsx_extended.f
@@ -55,7 +55,7 @@
 *> In addition to normwise error bound, the code provides maximum
 *> componentwise error bound if possible. See comments for ERR_BNDS_NORM
 *> and ERR_BNDS_COMP for details of the error bounds. Note that this
-*> subroutine is only resonsible for setting the second fields of
+*> subroutine is only responsible for setting the second fields of
 *> ERR_BNDS_NORM and ERR_BNDS_COMP.
 *> \endverbatim
 *
diff --git a/SRC/dla_gbrfsx_extended.f b/SRC/dla_gbrfsx_extended.f
index ecb0b9f09f..5454b1bfcd 100644
--- a/SRC/dla_gbrfsx_extended.f
+++ b/SRC/dla_gbrfsx_extended.f
@@ -55,7 +55,7 @@
 *> In addition to normwise error bound, the code provides maximum
 *> componentwise error bound if possible. See comments for ERR_BNDS_NORM
 *> and ERR_BNDS_COMP for details of the error bounds. Note that this
-*> subroutine is only resonsible for setting the second fields of
+*> subroutine is only responsible for setting the second fields of
 *> ERR_BNDS_NORM and ERR_BNDS_COMP.
 *> \endverbatim
 *
diff --git a/SRC/dla_gerfsx_extended.f b/SRC/dla_gerfsx_extended.f
index cff131b8bf..92b0d76d44 100644
--- a/SRC/dla_gerfsx_extended.f
+++ b/SRC/dla_gerfsx_extended.f
@@ -53,7 +53,7 @@
 *> In addition to normwise error bound, the code provides maximum
 *> componentwise error bound if possible. See comments for ERRS_N
 *> and ERRS_C for details of the error bounds. Note that this
-*> subroutine is only resonsible for setting the second fields of
+*> subroutine is only responsible for setting the second fields of
 *> ERRS_N and ERRS_C.
 *> \endverbatim
 *
diff --git a/SRC/dla_porfsx_extended.f b/SRC/dla_porfsx_extended.f
index cabf631fd3..5c8850fef2 100644
--- a/SRC/dla_porfsx_extended.f
+++ b/SRC/dla_porfsx_extended.f
@@ -54,7 +54,7 @@
 *> In addition to normwise error bound, the code provides maximum
 *> componentwise error bound if possible. See comments for ERR_BNDS_NORM
 *> and ERR_BNDS_COMP for details of the error bounds. Note that this
-*> subroutine is only resonsible for setting the second fields of
+*> subroutine is only responsible for setting the second fields of
 *> ERR_BNDS_NORM and ERR_BNDS_COMP.
 *> \endverbatim
 *
diff --git a/SRC/dla_syrfsx_extended.f b/SRC/dla_syrfsx_extended.f
index 4ccac897aa..e1cde6fc2b 100644
--- a/SRC/dla_syrfsx_extended.f
+++ b/SRC/dla_syrfsx_extended.f
@@ -56,7 +56,7 @@
 *> In addition to normwise error bound, the code provides maximum
 *> componentwise error bound if possible. See comments for ERR_BNDS_NORM
 *> and ERR_BNDS_COMP for details of the error bounds. Note that this
-*> subroutine is only resonsible for setting the second fields of
+*> subroutine is only responsible for setting the second fields of
 *> ERR_BNDS_NORM and ERR_BNDS_COMP.
 *> \endverbatim
 *
diff --git a/SRC/sla_gbrfsx_extended.f b/SRC/sla_gbrfsx_extended.f
index 036f76b3b2..499d6bf900 100644
--- a/SRC/sla_gbrfsx_extended.f
+++ b/SRC/sla_gbrfsx_extended.f
@@ -54,7 +54,7 @@
 *> In addition to normwise error bound, the code provides maximum
 *> componentwise error bound if possible. See comments for ERR_BNDS_NORM
 *> and ERR_BNDS_COMP for details of the error bounds. Note that this
-*> subroutine is only resonsible for setting the second fields of
+*> subroutine is only responsible for setting the second fields of
 *> ERR_BNDS_NORM and ERR_BNDS_COMP.
 *> \endverbatim
 *
diff --git a/SRC/sla_gerfsx_extended.f b/SRC/sla_gerfsx_extended.f
index 5149b33937..de05d8eb41 100644
--- a/SRC/sla_gerfsx_extended.f
+++ b/SRC/sla_gerfsx_extended.f
@@ -54,7 +54,7 @@
 *> In addition to normwise error bound, the code provides maximum
 *> componentwise error bound if possible. See comments for ERRS_N
 *> and ERRS_C for details of the error bounds. Note that this
-*> subroutine is only resonsible for setting the second fields of
+*> subroutine is only responsible for setting the second fields of
 *> ERRS_N and ERRS_C.
 *> \endverbatim
 *
diff --git a/SRC/sla_porfsx_extended.f b/SRC/sla_porfsx_extended.f
index bb67bf4531..ada4cad21e 100644
--- a/SRC/sla_porfsx_extended.f
+++ b/SRC/sla_porfsx_extended.f
@@ -54,7 +54,7 @@
 *> In addition to normwise error bound, the code provides maximum
 *> componentwise error bound if possible. See comments for ERR_BNDS_NORM
 *> and ERR_BNDS_COMP for details of the error bounds. Note that this
-*> subroutine is only resonsible for setting the second fields of
+*> subroutine is only responsible for setting the second fields of
 *> ERR_BNDS_NORM and ERR_BNDS_COMP.
 *> \endverbatim
 *
diff --git a/SRC/sla_syrfsx_extended.f b/SRC/sla_syrfsx_extended.f
index 9bd878c255..d5096be024 100644
--- a/SRC/sla_syrfsx_extended.f
+++ b/SRC/sla_syrfsx_extended.f
@@ -56,7 +56,7 @@
 *> In addition to normwise error bound, the code provides maximum
 *> componentwise error bound if possible. See comments for ERR_BNDS_NORM
 *> and ERR_BNDS_COMP for details of the error bounds. Note that this
-*> subroutine is only resonsible for setting the second fields of
+*> subroutine is only responsible for setting the second fields of
 *> ERR_BNDS_NORM and ERR_BNDS_COMP.
 *> \endverbatim
 *
diff --git a/SRC/zla_gbrfsx_extended.f b/SRC/zla_gbrfsx_extended.f
index d8185ec073..fe4d635b12 100644
--- a/SRC/zla_gbrfsx_extended.f
+++ b/SRC/zla_gbrfsx_extended.f
@@ -54,7 +54,7 @@
 *> In addition to normwise error bound, the code provides maximum
 *> componentwise error bound if possible. See comments for ERR_BNDS_NORM
 *> and ERR_BNDS_COMP for details of the error bounds. Note that this
-*> subroutine is only resonsible for setting the second fields of
+*> subroutine is only responsible for setting the second fields of
 *> ERR_BNDS_NORM and ERR_BNDS_COMP.
 *> \endverbatim
 *
diff --git a/SRC/zla_gerfsx_extended.f b/SRC/zla_gerfsx_extended.f
index 355def45cf..9d618f2942 100644
--- a/SRC/zla_gerfsx_extended.f
+++ b/SRC/zla_gerfsx_extended.f
@@ -53,7 +53,7 @@
 *> In addition to normwise error bound, the code provides maximum
 *> componentwise error bound if possible. See comments for ERRS_N
 *> and ERRS_C for details of the error bounds. Note that this
-*> subroutine is only resonsible for setting the second fields of
+*> subroutine is only responsible for setting the second fields of
 *> ERRS_N and ERRS_C.
 *> \endverbatim
 *
diff --git a/SRC/zla_herfsx_extended.f b/SRC/zla_herfsx_extended.f
index 2c8c778af7..a55dd94318 100644
--- a/SRC/zla_herfsx_extended.f
+++ b/SRC/zla_herfsx_extended.f
@@ -55,7 +55,7 @@
 *> In addition to normwise error bound, the code provides maximum
 *> componentwise error bound if possible. See comments for ERR_BNDS_NORM
 *> and ERR_BNDS_COMP for details of the error bounds. Note that this
-*> subroutine is only resonsible for setting the second fields of
+*> subroutine is only responsible for setting the second fields of
 *> ERR_BNDS_NORM and ERR_BNDS_COMP.
 *> \endverbatim
 *
diff --git a/SRC/zla_porfsx_extended.f b/SRC/zla_porfsx_extended.f
index e2265c56f5..12e05e049e 100644
--- a/SRC/zla_porfsx_extended.f
+++ b/SRC/zla_porfsx_extended.f
@@ -54,7 +54,7 @@
 *> In addition to normwise error bound, the code provides maximum
 *> componentwise error bound if possible. See comments for ERR_BNDS_NORM
 *> and ERR_BNDS_COMP for details of the error bounds. Note that this
-*> subroutine is only resonsible for setting the second fields of
+*> subroutine is only responsible for setting the second fields of
 *> ERR_BNDS_NORM and ERR_BNDS_COMP.
 *> \endverbatim
 *
diff --git a/SRC/zla_syrfsx_extended.f b/SRC/zla_syrfsx_extended.f
index bdb7c05c3b..d6c241499a 100644
--- a/SRC/zla_syrfsx_extended.f
+++ b/SRC/zla_syrfsx_extended.f
@@ -55,7 +55,7 @@
 *> In addition to normwise error bound, the code provides maximum
 *> componentwise error bound if possible. See comments for ERR_BNDS_NORM
 *> and ERR_BNDS_COMP for details of the error bounds. Note that this
-*> subroutine is only resonsible for setting the second fields of
+*> subroutine is only responsible for setting the second fields of
 *> ERR_BNDS_NORM and ERR_BNDS_COMP.
 *> \endverbatim
 *
diff --git a/TESTING/EIG/cchkst2stg.f b/TESTING/EIG/cchkst2stg.f
index d2eaea4360..e4deb8ac85 100644
--- a/TESTING/EIG/cchkst2stg.f
+++ b/TESTING/EIG/cchkst2stg.f
@@ -1072,7 +1072,6 @@ SUBROUTINE CCHKST2STG( NSIZES, NN, NTYPES, DOTYPE, ISEED, THRESH,
                END IF
             END IF
 *
-*
 *           Do Tests 3 and 4 which are similar to 11 and 12 but with the
 *           D1 computed using the standard 1-stage reduction as reference
 *
@@ -1735,7 +1734,6 @@ SUBROUTINE CCHKST2STG( NSIZES, NN, NTYPES, DOTYPE, ISEED, THRESH,
                         END IF
                      END IF
 *
-*
 *                 Do test 28
 *
                      TEMP2 = TWO*( TWO*N-ONE )*ULP*
@@ -1792,7 +1790,6 @@ SUBROUTINE CCHKST2STG( NSIZES, NN, NTYPES, DOTYPE, ISEED, THRESH,
 *
 *           Do Tests 29 and 30
 *
-*
 *           Call CSTEMR to compute D2, do tests.
 *
 *           Compute D2
@@ -1832,7 +1829,6 @@ SUBROUTINE CCHKST2STG( NSIZES, NN, NTYPES, DOTYPE, ISEED, THRESH,
                   RESULT( 31 ) = TEMP2 / MAX( UNFL,
      $                           ULP*MAX( TEMP1, TEMP2 ) )
 *
-*
 *           Call CSTEMR(V,V) to compute D1 and Z, do tests.
 *
 *           Compute D1 and Z
@@ -1934,7 +1930,6 @@ SUBROUTINE CCHKST2STG( NSIZES, NN, NTYPES, DOTYPE, ISEED, THRESH,
                   RESULT( 34 ) = ZERO
                END IF
 *
-*
 *           Call CSTEMR(V,A) to compute D1 and Z, do tests.
 *
 *           Compute D1 and Z
@@ -1991,7 +1986,7 @@ SUBROUTINE CCHKST2STG( NSIZES, NN, NTYPES, DOTYPE, ISEED, THRESH,
                   END IF
                END IF
 *
-*           Do Test 34
+*           Do Test 37
 *
                TEMP1 = ZERO
                TEMP2 = ZERO
@@ -2010,7 +2005,6 @@ SUBROUTINE CCHKST2STG( NSIZES, NN, NTYPES, DOTYPE, ISEED, THRESH,
 *
 *           End of Loop -- Check for RESULT(j) > THRESH
 *
-*
 *           Print out tests which fail.
 *
             DO 290 JR = 1, NTEST
diff --git a/TESTING/EIG/dchkst2stg.f b/TESTING/EIG/dchkst2stg.f
index 90eeec33dc..2c98b802d3 100644
--- a/TESTING/EIG/dchkst2stg.f
+++ b/TESTING/EIG/dchkst2stg.f
@@ -1055,7 +1055,6 @@ SUBROUTINE DCHKST2STG( NSIZES, NN, NTYPES, DOTYPE, ISEED, THRESH,
                END IF
             END IF
 *
-*
 *           Do Tests 3 and 4 which are similar to 11 and 12 but with the
 *           D1 computed using the standard 1-stage reduction as reference
 *
@@ -1715,7 +1714,6 @@ SUBROUTINE DCHKST2STG( NSIZES, NN, NTYPES, DOTYPE, ISEED, THRESH,
                         END IF
                      END IF
 *
-*
 *                 Do test 28
 *
                      TEMP2 = TWO*( TWO*N-ONE )*ULP*
@@ -1814,7 +1812,6 @@ SUBROUTINE DCHKST2STG( NSIZES, NN, NTYPES, DOTYPE, ISEED, THRESH,
                   RESULT( 31 ) = TEMP2 / MAX( UNFL,
      $                           ULP*MAX( TEMP1, TEMP2 ) )
 *
-*
 *           Call DSTEMR(V,V) to compute D1 and Z, do tests.
 *
 *           Compute D1 and Z
@@ -1916,7 +1913,6 @@ SUBROUTINE DCHKST2STG( NSIZES, NN, NTYPES, DOTYPE, ISEED, THRESH,
                   RESULT( 34 ) = ZERO
                END IF
 *
-*
 *           Call DSTEMR(V,A) to compute D1 and Z, do tests.
 *
 *           Compute D1 and Z
@@ -1973,7 +1969,7 @@ SUBROUTINE DCHKST2STG( NSIZES, NN, NTYPES, DOTYPE, ISEED, THRESH,
                   END IF
                END IF
 *
-*           Do Test 34
+*           Do Test 37
 *
                TEMP1 = ZERO
                TEMP2 = ZERO
@@ -1992,7 +1988,6 @@ SUBROUTINE DCHKST2STG( NSIZES, NN, NTYPES, DOTYPE, ISEED, THRESH,
 *
 *           End of Loop -- Check for RESULT(j) > THRESH
 *
-*
 *           Print out tests which fail.
 *
             DO 290 JR = 1, NTEST
diff --git a/TESTING/EIG/schkst2stg.f b/TESTING/EIG/schkst2stg.f
index 62b93c08f8..ac5a3fc393 100644
--- a/TESTING/EIG/schkst2stg.f
+++ b/TESTING/EIG/schkst2stg.f
@@ -1055,7 +1055,6 @@ SUBROUTINE SCHKST2STG( NSIZES, NN, NTYPES, DOTYPE, ISEED, THRESH,
                END IF
             END IF
 *
-*
 *           Do Tests 3 and 4 which are similar to 11 and 12 but with the
 *           D1 computed using the standard 1-stage reduction as reference
 *
@@ -1715,7 +1714,6 @@ SUBROUTINE SCHKST2STG( NSIZES, NN, NTYPES, DOTYPE, ISEED, THRESH,
                         END IF
                      END IF
 *
-*
 *                 Do test 28
 *
                      TEMP2 = TWO*( TWO*N-ONE )*ULP*
@@ -1814,7 +1812,6 @@ SUBROUTINE SCHKST2STG( NSIZES, NN, NTYPES, DOTYPE, ISEED, THRESH,
                   RESULT( 31 ) = TEMP2 / MAX( UNFL,
      $                           ULP*MAX( TEMP1, TEMP2 ) )
 *
-*
 *           Call SSTEMR(V,V) to compute D1 and Z, do tests.
 *
 *           Compute D1 and Z
@@ -1916,7 +1913,6 @@ SUBROUTINE SCHKST2STG( NSIZES, NN, NTYPES, DOTYPE, ISEED, THRESH,
                   RESULT( 34 ) = ZERO
                END IF
 *
-*
 *           Call SSTEMR(V,A) to compute D1 and Z, do tests.
 *
 *           Compute D1 and Z
@@ -1973,7 +1969,7 @@ SUBROUTINE SCHKST2STG( NSIZES, NN, NTYPES, DOTYPE, ISEED, THRESH,
                   END IF
                END IF
 *
-*           Do Test 34
+*           Do Test 37
 *
                TEMP1 = ZERO
                TEMP2 = ZERO
@@ -1992,7 +1988,6 @@ SUBROUTINE SCHKST2STG( NSIZES, NN, NTYPES, DOTYPE, ISEED, THRESH,
 *
 *           End of Loop -- Check for RESULT(j) > THRESH
 *
-*
 *           Print out tests which fail.
 *
             DO 290 JR = 1, NTEST
diff --git a/TESTING/EIG/zchkst2stg.f b/TESTING/EIG/zchkst2stg.f
index 1e2dd624ca..b1ef808166 100644
--- a/TESTING/EIG/zchkst2stg.f
+++ b/TESTING/EIG/zchkst2stg.f
@@ -1072,7 +1072,6 @@ SUBROUTINE ZCHKST2STG( NSIZES, NN, NTYPES, DOTYPE, ISEED, THRESH,
                END IF
             END IF
 *
-*
 *           Do Tests 3 and 4 which are similar to 11 and 12 but with the
 *           D1 computed using the standard 1-stage reduction as reference
 *
@@ -1735,7 +1734,6 @@ SUBROUTINE ZCHKST2STG( NSIZES, NN, NTYPES, DOTYPE, ISEED, THRESH,
                         END IF
                      END IF
 *
-*
 *                 Do test 28
 *
                      TEMP2 = TWO*( TWO*N-ONE )*ULP*
@@ -1792,7 +1790,6 @@ SUBROUTINE ZCHKST2STG( NSIZES, NN, NTYPES, DOTYPE, ISEED, THRESH,
 *
 *           Do Tests 29 and 30
 *
-*
 *           Call ZSTEMR to compute D2, do tests.
 *
 *           Compute D2
@@ -1832,7 +1829,6 @@ SUBROUTINE ZCHKST2STG( NSIZES, NN, NTYPES, DOTYPE, ISEED, THRESH,
                   RESULT( 31 ) = TEMP2 / MAX( UNFL,
      $                           ULP*MAX( TEMP1, TEMP2 ) )
 *
-*
 *           Call ZSTEMR(V,V) to compute D1 and Z, do tests.
 *
 *           Compute D1 and Z
@@ -1934,7 +1930,6 @@ SUBROUTINE ZCHKST2STG( NSIZES, NN, NTYPES, DOTYPE, ISEED, THRESH,
                   RESULT( 34 ) = ZERO
                END IF
 *
-*
 *           Call ZSTEMR(V,A) to compute D1 and Z, do tests.
 *
 *           Compute D1 and Z
@@ -1991,7 +1986,7 @@ SUBROUTINE ZCHKST2STG( NSIZES, NN, NTYPES, DOTYPE, ISEED, THRESH,
                   END IF
                END IF
 *
-*           Do Test 34
+*           Do Test 37
 *
                TEMP1 = ZERO
                TEMP2 = ZERO
@@ -2010,7 +2005,6 @@ SUBROUTINE ZCHKST2STG( NSIZES, NN, NTYPES, DOTYPE, ISEED, THRESH,
 *
 *           End of Loop -- Check for RESULT(j) > THRESH
 *
-*
 *           Print out tests which fail.
 *
             DO 290 JR = 1, NTEST
diff --git a/TESTING/LIN/ctbt02.f b/TESTING/LIN/ctbt02.f
index 8c6a24b51e..bf58aa1da1 100644
--- a/TESTING/LIN/ctbt02.f
+++ b/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:
diff --git a/TESTING/LIN/ctpt02.f b/TESTING/LIN/ctpt02.f
index 12d8870c9d..9c32d4206a 100644
--- a/TESTING/LIN/ctpt02.f
+++ b/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:
diff --git a/TESTING/LIN/ctrt02.f b/TESTING/LIN/ctrt02.f
index 1addc5881c..72f44af4c7 100644
--- a/TESTING/LIN/ctrt02.f
+++ b/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:
diff --git a/TESTING/LIN/dtbt02.f b/TESTING/LIN/dtbt02.f
index 5f5192a003..347fbcb9c8 100644
--- a/TESTING/LIN/dtbt02.f
+++ b/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
 *
diff --git a/TESTING/LIN/dtpt02.f b/TESTING/LIN/dtpt02.f
index 9506a1d195..8b0b4a73d4 100644
--- a/TESTING/LIN/dtpt02.f
+++ b/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
 *
diff --git a/TESTING/LIN/dtrt02.f b/TESTING/LIN/dtrt02.f
index a7c95e136d..6d22f688ba 100644
--- a/TESTING/LIN/dtrt02.f
+++ b/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
 *
diff --git a/TESTING/LIN/stbt02.f b/TESTING/LIN/stbt02.f
index 31fbf9262a..f95b5d5f10 100644
--- a/TESTING/LIN/stbt02.f
+++ b/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
 *
diff --git a/TESTING/LIN/stpt02.f b/TESTING/LIN/stpt02.f
index 72936dd044..f81fa29b99 100644
--- a/TESTING/LIN/stpt02.f
+++ b/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
 *
diff --git a/TESTING/LIN/strt02.f b/TESTING/LIN/strt02.f
index df91501ded..26a95510a9 100644
--- a/TESTING/LIN/strt02.f
+++ b/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
 *
diff --git a/TESTING/LIN/ztbt02.f b/TESTING/LIN/ztbt02.f
index 3c9f119569..20917d9f7e 100644
--- a/TESTING/LIN/ztbt02.f
+++ b/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:
diff --git a/TESTING/LIN/ztpt02.f b/TESTING/LIN/ztpt02.f
index 8555053d14..d6b9dabafa 100644
--- a/TESTING/LIN/ztpt02.f
+++ b/TESTING/LIN/ztpt02.f
@@ -30,9 +30,10 @@
 *> ZTPT02 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:
diff --git a/TESTING/LIN/ztrt02.f b/TESTING/LIN/ztrt02.f
index 47e9f3463c..69f04c9dbc 100644
--- a/TESTING/LIN/ztrt02.f
+++ b/TESTING/LIN/ztrt02.f
@@ -30,10 +30,10 @@
 *>
 *> ZTRT02 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: