--- rpl/lapack/lapack/zla_lin_berr.f	2010/12/21 13:53:48	1.4
+++ rpl/lapack/lapack/zla_lin_berr.f	2011/11/21 20:43:14	1.5
@@ -1,15 +1,111 @@
-      SUBROUTINE ZLA_LIN_BERR ( N, NZ, NRHS, RES, AYB, BERR )
+*> \brief \b ZLA_LIN_BERR
 *
-*     -- 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 
+*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/zla_lin_berr.f"> 
+*> [TGZ]</a> 
+*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/zla_lin_berr.f"> 
+*> [ZIP]</a> 
+*> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/zla_lin_berr.f"> 
+*> [TXT]</a>
+*> \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 DOUBLE PRECISION 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 COMPLEX*16 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 November 2011
+*
+*> \ingroup complex16OTHERcomputational
+*
+*  =====================================================================
+      SUBROUTINE ZLA_LIN_BERR ( N, NZ, NRHS, RES, AYB, BERR )
+*
+*  -- LAPACK computational routine (version 3.4.0) --
+*  -- LAPACK is a software package provided by Univ. of Tennessee,    --
+*  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
+*     November 2011
 *
-      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 ..