Diff for /rpl/lapack/lapack/zla_lin_berr.f between versions 1.4 and 1.5

version 1.4, 2010/12/21 13:53:48 version 1.5, 2011/11/21 20:43:14
Line 1 Line 1
       SUBROUTINE ZLA_LIN_BERR ( N, NZ, NRHS, RES, AYB, BERR )  *> \brief \b ZLA_LIN_BERR
 *  *
 *     -- LAPACK routine (version 3.2.2)                                 --  *  =========== DOCUMENTATION ===========
 *     -- Contributed by James Demmel, Deaglan Halligan, Yozo Hida and --  
 *     -- Jason Riedy of Univ. of California Berkeley.                 --  
 *     -- June 2010                                                    --  
 *  *
 *     -- LAPACK is a software package provided by Univ. of Tennessee, --  * Online html documentation available at 
 *     -- Univ. of California Berkeley and NAG Ltd.                    --  *            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 ..  *     .. Scalar Arguments ..
       INTEGER            N, NZ, NRHS        INTEGER            N, NZ, NRHS
 *     ..  *     ..
Line 18 Line 114
       COMPLEX*16         RES( N, NRHS )        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 ..  *     .. Local Scalars ..

Removed from v.1.4  
changed lines
  Added in v.1.5


CVSweb interface <joel.bertrand@systella.fr>