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

NAME

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

SYNOPSIS

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

PURPOSE

24       DGBCON  estimates the reciprocal of the condition number of a real gen‐
25       eral band matrix A, in either the 1-norm or  the  infinity-norm,  using
26       the LU factorization computed by DGBTRF.
27       An  estimate  is  obtained  for norm(inv(A)), and the reciprocal of the
28       condition number is computed as
29          RCOND = 1 / ( norm(A) * norm(inv(A)) ).
30

ARGUMENTS

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