--- rpl/lapack/lapack/zla_lin_berr.f 2010/08/07 13:22:37 1.2
+++ rpl/lapack/lapack/zla_lin_berr.f 2017/06/17 10:54:17 1.13
@@ -1,15 +1,111 @@
- SUBROUTINE ZLA_LIN_BERR ( N, NZ, NRHS, RES, AYB, BERR )
+*> \brief \b ZLA_LIN_BERR computes a component-wise relative backward error.
*
-* -- LAPACK routine (version 3.2.2) --
-* -- Contributed by James Demmel, Deaglan Halligan, Yozo Hida and --
-* -- Jason Riedy of Univ. of California Berkeley. --
-* -- June 2010 --
+* =========== DOCUMENTATION ===========
*
-* -- LAPACK is a software package provided by Univ. of Tennessee, --
-* -- Univ. of California Berkeley and NAG Ltd. --
+* Online html documentation available at
+* http://www.netlib.org/lapack/explore-html/
+*
+*> \htmlonly
+*> Download ZLA_LIN_BERR + dependencies
+*>
+*> [TGZ]
+*>
+*> [ZIP]
+*>
+*> [TXT]
+*> \endhtmlonly
+*
+* Definition:
+* ===========
+*
+* SUBROUTINE ZLA_LIN_BERR ( N, NZ, NRHS, RES, AYB, BERR )
+*
+* .. Scalar Arguments ..
+* INTEGER N, NZ, NRHS
+* ..
+* .. Array Arguments ..
+* DOUBLE PRECISION AYB( N, NRHS ), BERR( NRHS )
+* COMPLEX*16 RES( N, NRHS )
+* ..
+*
+*
+*> \par Purpose:
+* =============
+*>
+*> \verbatim
+*>
+*> ZLA_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 componentwise absolute value of the matrix
+*> or vector Z.
+*> \endverbatim
+*
+* Arguments:
+* ==========
+*
+*> \param[in] N
+*> \verbatim
+*> N is INTEGER
+*> The number of linear equations, i.e., the order of the
+*> matrix A. N >= 0.
+*> \endverbatim
+*>
+*> \param[in] NZ
+*> \verbatim
+*> NZ is INTEGER
+*> We add (NZ+1)*SLAMCH( 'Safe minimum' ) to R(i) in the numerator to
+*> guard against spuriously zero residuals. Default value is N.
+*> \endverbatim
+*>
+*> \param[in] NRHS
+*> \verbatim
+*> NRHS is INTEGER
+*> The number of right hand sides, i.e., the number of columns
+*> of the matrices AYB, RES, and BERR. NRHS >= 0.
+*> \endverbatim
+*>
+*> \param[in] RES
+*> \verbatim
+*> RES is COMPLEX*16 array, dimension (N,NRHS)
+*> The residual matrix, i.e., the matrix R in the relative backward
+*> error formula above.
+*> \endverbatim
+*>
+*> \param[in] AYB
+*> \verbatim
+*> AYB is DOUBLE PRECISION 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 zla_gerfsx_extended.f).
+*> \endverbatim
+*>
+*> \param[out] BERR
+*> \verbatim
+*> BERR is DOUBLE PRECISION array, dimension (NRHS)
+*> The componentwise relative backward error from the formula above.
+*> \endverbatim
+*
+* Authors:
+* ========
+*
+*> \author Univ. of Tennessee
+*> \author Univ. of California Berkeley
+*> \author Univ. of Colorado Denver
+*> \author NAG Ltd.
+*
+*> \date June 2016
+*
+*> \ingroup complex16OTHERcomputational
+*
+* =====================================================================
+ SUBROUTINE ZLA_LIN_BERR ( N, NZ, NRHS, RES, AYB, BERR )
+*
+* -- LAPACK computational routine (version 3.7.0) --
+* -- LAPACK is a software package provided by Univ. of Tennessee, --
+* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
+* June 2016
*
- IMPLICIT NONE
-* ..
* .. Scalar Arguments ..
INTEGER N, NZ, NRHS
* ..
@@ -18,39 +114,6 @@
COMPLEX*16 RES( N, NRHS )
* ..
*
-* Purpose
-* =======
-*
-* ZLA_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 componentwise absolute value of the matrix
-* or vector Z.
-*
-* 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) DOUBLE PRECISION array, dimension (N,NRHS)
-* The residual matrix, i.e., the matrix R in the relative backward
-* error formula above.
-*
-* AYB (input) DOUBLE PRECISION 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 zla_gerfsx_extended.f).
-*
-* BERR (output) COMPLEX*16 array, dimension (NRHS)
-* The componentwise relative backward error from the formula above.
-*
* =====================================================================
*
* .. Local Scalars ..