Updated documentation for EPIR routines. Changed ERRS_{N,C} variable names. Other cosmetic changes.

1 files changed, 29 insertions, 2 deletions

diff --git a/SRC/sla_lin_berr.f b/SRC/sla_lin_berr.f
index 74cb2c9f..7a92d17a 100644
--- a/SRC/sla_lin_berr.f
+++ b/SRC/sla_lin_berr.f
@@ -21,12 +21,39 @@
 *  Purpose
 *  =======
 *
-*     SLA_LIN_BERR computes component-wise relative backward error from
+*     SLA_LIN_BERR computes componentwise relative backward error from
 *     the formula
 *         max(i) ( abs(R(i)) / ( abs(op(A_s))*abs(Y) + abs(B_s) )(i) )
-*     where abs(Z) is the component-wise absolute value of the matrix
+*     where abs(Z) is the componentwise absolute value of the matrix
 *     or vector Z.
 *
+*  Arguments
+*  ==========
+*
+*     N       (input) INTEGER
+*     The number of linear equations, i.e., the order of the
+*     matrix A.  N >= 0.
+*
+*     NZ      (input) INTEGER
+*     We add (NZ+1)*SLAMCH( 'Safe minimum' ) to R(i) in the numerator to
+*     guard against spuriously zero residuals. Default value is N.
+*
+*     NRHS    (input) INTEGER
+*     The number of right hand sides, i.e., the number of columns
+*     of the matrices AYB, RES, and BERR.  NRHS >= 0.
+*
+*     RES    (input) REAL array, dimension (N,NRHS)
+*     The residual matrix, i.e., the matrix R in the relative backward
+*     error formula above.
+*
+*     AYB    (input) REAL array, dimension (N, NRHS)
+*     The denominator in the relative backward error formula above, i.e.,
+*     the matrix abs(op(A_s))*abs(Y) + abs(B_s). The matrices A, Y, and B
+*     are from iterative refinement (see sla_gerfsx_extended.f).
+*     
+*     RES    (output) REAL array, dimension (NRHS)
+*     The componentwise relative backward error from the formula above.
+*
 *  =====================================================================
 *
 *     .. Local Scalars ..