1ZGBCON(1)                LAPACK routine (version 3.2)                ZGBCON(1)
2
3
4

NAME

6       ZGBCON  - estimates the reciprocal of the condition number of a complex
7       general band matrix A, in either the 1-norm or the infinity-norm,
8

SYNOPSIS

10       SUBROUTINE ZGBCON( NORM, N, KL, KU, AB, LDAB, IPIV, ANORM, RCOND, WORK,
11                          RWORK, INFO )
12
13           CHARACTER      NORM
14
15           INTEGER        INFO, KL, KU, LDAB, N
16
17           DOUBLE         PRECISION ANORM, RCOND
18
19           INTEGER        IPIV( * )
20
21           DOUBLE         PRECISION RWORK( * )
22
23           COMPLEX*16     AB( LDAB, * ), WORK( * )
24

PURPOSE

26       ZGBCON  estimates  the  reciprocal of the condition number of a complex
27       general band matrix A, in either the 1-norm or the infinity-norm, using
28       the LU factorization computed by ZGBTRF.
29       An  estimate  is  obtained  for norm(inv(A)), and the reciprocal of the
30       condition number is computed as
31          RCOND = 1 / ( norm(A) * norm(inv(A)) ).
32

ARGUMENTS

34       NORM    (input) CHARACTER*1
35               Specifies whether the 1-norm condition number or the  infinity-
36               norm condition number is required:
37               = '1' or 'O':  1-norm;
38               = 'I':         Infinity-norm.
39
40       N       (input) INTEGER
41               The order of the matrix A.  N >= 0.
42
43       KL      (input) INTEGER
44               The number of subdiagonals within the band of A.  KL >= 0.
45
46       KU      (input) INTEGER
47               The number of superdiagonals within the band of A.  KU >= 0.
48
49       AB      (input) COMPLEX*16 array, dimension (LDAB,N)
50               Details  of  the LU factorization of the band matrix A, as com‐
51               puted by ZGBTRF.  U is  stored  as  an  upper  triangular  band
52               matrix  with KL+KU superdiagonals in rows 1 to KL+KU+1, and the
53               multipliers used during the factorization are  stored  in  rows
54               KL+KU+2 to 2*KL+KU+1.
55
56       LDAB    (input) INTEGER
57               The leading dimension of the array AB.  LDAB >= 2*KL+KU+1.
58
59       IPIV    (input) INTEGER array, dimension (N)
60               The  pivot  indices;  for  1 <= i <= N, row i of the matrix was
61               interchanged with row IPIV(i).
62
63       ANORM   (input) DOUBLE PRECISION
64               If NORM = '1' or 'O', the 1-norm of the original matrix A.   If
65               NORM = 'I', the infinity-norm of the original matrix A.
66
67       RCOND   (output) DOUBLE PRECISION
68               The  reciprocal  of  the condition number of the matrix A, com‐
69               puted as RCOND = 1/(norm(A) * norm(inv(A))).
70
71       WORK    (workspace) COMPLEX*16 array, dimension (2*N)
72
73       RWORK   (workspace) DOUBLE PRECISION array, dimension (N)
74
75       INFO    (output) INTEGER
76               = 0:  successful exit
77               < 0: if INFO = -i, the i-th argument had an illegal value
78
79
80
81 LAPACK routine (version 3.2)    November 2008                       ZGBCON(1)
Impressum