ZGBCON(3)      LAPACK routine of NEC Numeric Library Collection      ZGBCON(3)



NAME
       ZGBCON

SYNOPSIS
       SUBROUTINE ZGBCON (NORM, N, KL, KU, AB, LDAB, IPIV, ANORM, RCOND, WORK,
           RWORK, INFO)



PURPOSE
            ZGBCON estimates the reciprocal of the condition number of a complex
            general band matrix A, in either the 1-norm or the infinity-norm,
            using the LU factorization computed by ZGBTRF.

            An estimate is obtained for norm(inv(A)), and the reciprocal of the
            condition number is computed as
               RCOND = 1 / ( norm(A) * norm(inv(A)) ).




ARGUMENTS
           NORM      (input)
                     NORM is CHARACTER*1
                     Specifies whether the 1-norm condition number or the
                     infinity-norm condition number is required:
                     = '1' or 'O':  1-norm;
                     = 'I':         Infinity-norm.

           N         (input)
                     N is INTEGER
                     The order of the matrix A.  N >= 0.

           KL        (input)
                     KL is INTEGER
                     The number of subdiagonals within the band of A.  KL >= 0.

           KU        (input)
                     KU is INTEGER
                     The number of superdiagonals within the band of A.  KU >= 0.

           AB        (input)
                     AB is COMPLEX*16 array, dimension (LDAB,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.

           LDAB      (input)
                     LDAB is INTEGER
                     The leading dimension of the array AB.  LDAB >= 2*KL+KU+1.

           IPIV      (input)
                     IPIV is INTEGER array, dimension (N)
                     The pivot indices; for 1 <= i <= N, row i of the matrix was
                     interchanged with row IPIV(i).

           ANORM     (input)
                     ANORM is DOUBLE PRECISION
                     If NORM = '1' or 'O', the 1-norm of the original matrix A.
                     If NORM = 'I', the infinity-norm of the original matrix A.

           RCOND     (output)
                     RCOND is DOUBLE PRECISION
                     The reciprocal of the condition number of the matrix A,
                     computed as RCOND = 1/(norm(A) * norm(inv(A))).

           WORK      (output)
                     WORK is COMPLEX*16 array, dimension (2*N)

           RWORK     (output)
                     RWORK is DOUBLE PRECISION array, dimension (N)

           INFO      (output)
                     INFO is INTEGER
                     = 0:  successful exit
                     < 0: if INFO = -i, the i-th argument had an illegal value



LAPACK routine                  31 October 2017                      ZGBCON(3)