[BACK]Return to zla_gbrcond_x.f CVS log [TXT] [DIR] Up to [local] / rpl / lapack / lapack

Diff for /rpl/lapack/lapack/zla_gbrcond_x.f between versions 1.5 and 1.6

version 1.5, 2011/07/22 07:38:16 version 1.6, 2011/11/21 20:43:13
Line 1 Line 1
   *> \brief \b ZLA_GBRCOND_X
   *
   *  =========== DOCUMENTATION ===========
   *
   * Online html documentation available at 
   *            http://www.netlib.org/lapack/explore-html/ 
   *
   *> \htmlonly
   *> Download ZLA_GBRCOND_X + dependencies 
   *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.tgz?format=tgz&filename=/lapack/lapack_routine/zla_gbrcond_x.f"> 
   *> [TGZ]</a> 
   *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.zip?format=zip&filename=/lapack/lapack_routine/zla_gbrcond_x.f"> 
   *> [ZIP]</a> 
   *> <a href="http://www.netlib.org/cgi-bin/netlibfiles.txt?format=txt&filename=/lapack/lapack_routine/zla_gbrcond_x.f"> 
   *> [TXT]</a>
   *> \endhtmlonly 
   *
   *  Definition:
   *  ===========
   *
   *       DOUBLE PRECISION FUNCTION ZLA_GBRCOND_X( TRANS, N, KL, KU, AB,
   *                                                LDAB, AFB, LDAFB, IPIV,
   *                                                X, INFO, WORK, RWORK )
   * 
   *       .. Scalar Arguments ..
   *       CHARACTER          TRANS
   *       INTEGER            N, KL, KU, KD, KE, LDAB, LDAFB, INFO
   *       ..
   *       .. Array Arguments ..
   *       INTEGER            IPIV( * )
   *       COMPLEX*16         AB( LDAB, * ), AFB( LDAFB, * ), WORK( * ),
   *      $                   X( * )
   *       DOUBLE PRECISION   RWORK( * )
   *  
   *  
   *
   *> \par Purpose:
   *  =============
   *>
   *> \verbatim
   *>
   *>    ZLA_GBRCOND_X Computes the infinity norm condition number of
   *>    op(A) * diag(X) where X is a COMPLEX*16 vector.
   *> \endverbatim
   *
   *  Arguments:
   *  ==========
   *
   *> \param[in] TRANS
   *> \verbatim
   *>          TRANS is CHARACTER*1
   *>     Specifies the form of the system of equations:
   *>       = 'N':  A * X = B     (No transpose)
   *>       = 'T':  A**T * X = B  (Transpose)
   *>       = 'C':  A**H * X = B  (Conjugate Transpose = Transpose)
   *> \endverbatim
   *>
   *> \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] KL
   *> \verbatim
   *>          KL is INTEGER
   *>     The number of subdiagonals within the band of A.  KL >= 0.
   *> \endverbatim
   *>
   *> \param[in] KU
   *> \verbatim
   *>          KU is INTEGER
   *>     The number of superdiagonals within the band of A.  KU >= 0.
   *> \endverbatim
   *>
   *> \param[in] AB
   *> \verbatim
   *>          AB is COMPLEX*16 array, dimension (LDAB,N)
   *>     On entry, the matrix A in band storage, in rows 1 to KL+KU+1.
   *>     The j-th column of A is stored in the j-th column of the
   *>     array AB as follows:
   *>     AB(KU+1+i-j,j) = A(i,j) for max(1,j-KU)<=i<=min(N,j+kl)
   *> \endverbatim
   *>
   *> \param[in] LDAB
   *> \verbatim
   *>          LDAB is INTEGER
   *>     The leading dimension of the array AB.  LDAB >= KL+KU+1.
   *> \endverbatim
   *>
   *> \param[in] AFB
   *> \verbatim
   *>          AFB is COMPLEX*16 array, dimension (LDAFB,N)
   *>     Details of the LU factorization of the band matrix A, as
   *>     computed by ZGBTRF.  U is stored as an upper triangular
   *>     band matrix with KL+KU superdiagonals in rows 1 to KL+KU+1,
   *>     and the multipliers used during the factorization are stored
   *>     in rows KL+KU+2 to 2*KL+KU+1.
   *> \endverbatim
   *>
   *> \param[in] LDAFB
   *> \verbatim
   *>          LDAFB is INTEGER
   *>     The leading dimension of the array AFB.  LDAFB >= 2*KL+KU+1.
   *> \endverbatim
   *>
   *> \param[in] IPIV
   *> \verbatim
   *>          IPIV is INTEGER array, dimension (N)
   *>     The pivot indices from the factorization A = P*L*U
   *>     as computed by ZGBTRF; row i of the matrix was interchanged
   *>     with row IPIV(i).
   *> \endverbatim
   *>
   *> \param[in] X
   *> \verbatim
   *>          X is COMPLEX*16 array, dimension (N)
   *>     The vector X in the formula op(A) * diag(X).
   *> \endverbatim
   *>
   *> \param[out] INFO
   *> \verbatim
   *>          INFO is INTEGER
   *>       = 0:  Successful exit.
   *>     i > 0:  The ith argument is invalid.
   *> \endverbatim
   *>
   *> \param[in] WORK
   *> \verbatim
   *>          WORK is COMPLEX*16 array, dimension (2*N).
   *>     Workspace.
   *> \endverbatim
   *>
   *> \param[in] RWORK
   *> \verbatim
   *>          RWORK is DOUBLE PRECISION array, dimension (N).
   *>     Workspace.
   *> \endverbatim
   *
   *  Authors:
   *  ========
   *
   *> \author Univ. of Tennessee 
   *> \author Univ. of California Berkeley 
   *> \author Univ. of Colorado Denver 
   *> \author NAG Ltd. 
   *
   *> \date November 2011
   *
   *> \ingroup complex16GBcomputational
   *
   *  =====================================================================
       DOUBLE PRECISION FUNCTION ZLA_GBRCOND_X( TRANS, N, KL, KU, AB,        DOUBLE PRECISION FUNCTION ZLA_GBRCOND_X( TRANS, N, KL, KU, AB,
      $                                         LDAB, AFB, LDAFB, IPIV,       $                                         LDAB, AFB, LDAFB, IPIV,
      $                                         X, INFO, WORK, RWORK )       $                                         X, INFO, WORK, RWORK )
 *  *
 *     -- LAPACK routine (version 3.2.1)                               --  *  -- LAPACK computational routine (version 3.4.0) --
 *     -- Contributed by James Demmel, Deaglan Halligan, Yozo Hida and --  *  -- LAPACK is a software package provided by Univ. of Tennessee,    --
 *     -- Jason Riedy of Univ. of California Berkeley.                 --  *  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
 *     -- April 2009                                                   --  *     November 2011
 *  
 *     -- LAPACK is a software package provided by Univ. of Tennessee, --  
 *     -- Univ. of California Berkeley and NAG Ltd.                    --  
 *  *
       IMPLICIT NONE  
 *     ..  
 *     .. Scalar Arguments ..  *     .. Scalar Arguments ..
       CHARACTER          TRANS        CHARACTER          TRANS
       INTEGER            N, KL, KU, KD, KE, LDAB, LDAFB, INFO        INTEGER            N, KL, KU, KD, KE, LDAB, LDAFB, INFO
Line 23 Line 171
       DOUBLE PRECISION   RWORK( * )        DOUBLE PRECISION   RWORK( * )
 *  *
 *  *
 *  Purpose  
 *  =======  
 *  
 *     ZLA_GBRCOND_X Computes the infinity norm condition number of  
 *     op(A) * diag(X) where X is a COMPLEX*16 vector.  
 *  
 *  Arguments  
 *  =========  
 *  
 *     TRANS   (input) CHARACTER*1  
 *     Specifies the form of the system of equations:  
 *       = 'N':  A * X = B     (No transpose)  
 *       = 'T':  A**T * X = B  (Transpose)  
 *       = 'C':  A**H * X = B  (Conjugate Transpose = Transpose)  
 *  
 *     N       (input) INTEGER  
 *     The number of linear equations, i.e., the order of the  
 *     matrix A.  N >= 0.  
 *  
 *     KL      (input) INTEGER  
 *     The number of subdiagonals within the band of A.  KL >= 0.  
 *  
 *     KU      (input) INTEGER  
 *     The number of superdiagonals within the band of A.  KU >= 0.  
 *  
 *     AB      (input) COMPLEX*16 array, dimension (LDAB,N)  
 *     On entry, the matrix A in band storage, in rows 1 to KL+KU+1.  
 *     The j-th column of A is stored in the j-th column of the  
 *     array AB as follows:  
 *     AB(KU+1+i-j,j) = A(i,j) for max(1,j-KU)<=i<=min(N,j+kl)  
 *  
 *     LDAB    (input) INTEGER  
 *     The leading dimension of the array AB.  LDAB >= KL+KU+1.  
 *  
 *     AFB     (input) COMPLEX*16 array, dimension (LDAFB,N)  
 *     Details of the LU factorization of the band matrix A, as  
 *     computed by ZGBTRF.  U is stored as an upper triangular  
 *     band matrix with KL+KU superdiagonals in rows 1 to KL+KU+1,  
 *     and the multipliers used during the factorization are stored  
 *     in rows KL+KU+2 to 2*KL+KU+1.  
 *  
 *     LDAFB   (input) INTEGER  
 *     The leading dimension of the array AFB.  LDAFB >= 2*KL+KU+1.  
 *  
 *     IPIV    (input) INTEGER array, dimension (N)  
 *     The pivot indices from the factorization A = P*L*U  
 *     as computed by ZGBTRF; row i of the matrix was interchanged  
 *     with row IPIV(i).  
 *  
 *     X       (input) COMPLEX*16 array, dimension (N)  
 *     The vector X in the formula op(A) * diag(X).  
 *  
 *     INFO    (output) INTEGER  
 *       = 0:  Successful exit.  
 *     i > 0:  The ith argument is invalid.  
 *  
 *     WORK    (input) COMPLEX*16 array, dimension (2*N).  
 *     Workspace.  
 *  
 *     RWORK   (input) DOUBLE PRECISION array, dimension (N).  
 *     Workspace.  
 *  
 *  =====================================================================  *  =====================================================================
 *  *
 *     .. Local Scalars ..  *     .. Local Scalars ..
Removed from v.1.5  
changed lines
  Added in v.1.6

CVSweb interface <joel.bertrand@systella.fr>