1CGBCON(1) LAPACK routine (version 3.1) CGBCON(1)
2
3
4
6 CGBCON - the reciprocal of the condition number of a complex general
7 band matrix A, in either the 1-norm or the infinity-norm,
8
10 SUBROUTINE CGBCON( 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 REAL ANORM, RCOND
18
19 INTEGER IPIV( * )
20
21 REAL RWORK( * )
22
23 COMPLEX AB( LDAB, * ), WORK( * )
24
26 CGBCON 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 CGBTRF.
29
30 An estimate is obtained for norm(inv(A)), and the reciprocal of the
31 condition number is computed as
32 RCOND = 1 / ( norm(A) * norm(inv(A)) ).
33
34
36 NORM (input) CHARACTER*1
37 Specifies whether the 1-norm condition number or the infinity-
38 norm condition number is required:
39 = '1' or 'O': 1-norm;
40 = 'I': Infinity-norm.
41
42 N (input) INTEGER
43 The order of the matrix A. N >= 0.
44
45 KL (input) INTEGER
46 The number of subdiagonals within the band of A. KL >= 0.
47
48 KU (input) INTEGER
49 The number of superdiagonals within the band of A. KU >= 0.
50
51 AB (input) COMPLEX array, dimension (LDAB,N)
52 Details of the LU factorization of the band matrix A, as com‐
53 puted by CGBTRF. U is stored as an upper triangular band
54 matrix with KL+KU superdiagonals in rows 1 to KL+KU+1, and the
55 multipliers used during the factorization are stored in rows
56 KL+KU+2 to 2*KL+KU+1.
57
58 LDAB (input) INTEGER
59 The leading dimension of the array AB. LDAB >= 2*KL+KU+1.
60
61 IPIV (input) INTEGER array, dimension (N)
62 The pivot indices; for 1 <= i <= N, row i of the matrix was
63 interchanged with row IPIV(i).
64
65 ANORM (input) REAL
66 If NORM = '1' or 'O', the 1-norm of the original matrix A. If
67 NORM = 'I', the infinity-norm of the original matrix A.
68
69 RCOND (output) REAL
70 The reciprocal of the condition number of the matrix A, com‐
71 puted as RCOND = 1/(norm(A) * norm(inv(A))).
72
73 WORK (workspace) COMPLEX array, dimension (2*N)
74
75 RWORK (workspace) REAL array, dimension (N)
76
77 INFO (output) INTEGER
78 = 0: successful exit
79 < 0: if INFO = -i, the i-th argument had an illegal value
80
81
82
83 LAPACK routine (version 3.1) November 2006 CGBCON(1)