comments unified: add an EPS description in xQRT17

Author: jip

TESTING/LIN/cqrt17.f

@@ -28,9 +28,10 @@
 *>
 *> CQRT17 computes the ratio
 *>
-*>    norm(R**H * op(A))/(norm(A)*alpha*max(M,N,NRHS)*eps),
+*>    norm(R**H * op(A)) / ( norm(A) * alpha * max(M,N,NRHS) * EPS ),
 *>
-*> where R = B - op(A)*X, op(A) is A or A**H, depending on TRANS, and
+*> where R = B - op(A)*X, op(A) is A or A**H, depending on TRANS, EPS
+*> is the machine epsilon, and
 *>
 *>    alpha = norm(B) if IRESID = 1 (zero-residual problem)
 *>    alpha = norm(R) if IRESID = 2 (otherwise).

TESTING/LIN/dqrt17.f

@@ -28,9 +28,10 @@
 *>
 *> DQRT17 computes the ratio
 *>
-*>    norm(R**T * op(A))/(norm(A)*alpha*max(M,N,NRHS)*eps),
+*>    norm(R**T * op(A)) / ( norm(A) * alpha * max(M,N,NRHS) * EPS ),
 *>
-*> where R = B - op(A)*X, op(A) is A or A**T, depending on TRANS, and
+*> where R = B - op(A)*X, op(A) is A or A**T, depending on TRANS, EPS
+*> is the machine epsilon, and
 *>
 *>    alpha = norm(B) if IRESID = 1 (zero-residual problem)
 *>    alpha = norm(R) if IRESID = 2 (otherwise).

TESTING/LIN/sqrt17.f

@@ -28,9 +28,10 @@
 *>
 *> SQRT17 computes the ratio
 *>
-*>    norm(R**T * op(A))/(norm(A)*alpha*max(M,N,NRHS)*eps),
+*>    norm(R**T * op(A)) / ( norm(A) * alpha * max(M,N,NRHS) * EPS ),
 *>
-*> where R = B - op(A)*X, op(A) is A or A**T, depending on TRANS, and
+*> where R = B - op(A)*X, op(A) is A or A**T, depending on TRANS, EPS
+*> is the machine epsilon, and
 *>
 *>    alpha = norm(B) if IRESID = 1 (zero-residual problem)
 *>    alpha = norm(R) if IRESID = 2 (otherwise).

TESTING/LIN/zqrt17.f

@@ -28,9 +28,10 @@
 *>
 *> ZQRT17 computes the ratio
 *>
-*>    norm(R**H * op(A))/(norm(A)*alpha*max(M,N,NRHS)*eps),
+*>    norm(R**H * op(A)) / ( norm(A) * alpha * max(M,N,NRHS) * EPS ),
 *>
-*> where R = B - op(A)*X, op(A) is A or A**H, depending on TRANS, and
+*> where R = B - op(A)*X, op(A) is A or A**H, depending on TRANS, EPS
+*> is the machine epsilon, and
 *>
 *>    alpha = norm(B) if IRESID = 1 (zero-residual problem)
 *>    alpha = norm(R) if IRESID = 2 (otherwise).