1ZPBCON(1)                LAPACK routine (version 3.1)                ZPBCON(1)
2
3
4

NAME

6       ZPBCON  -  the  reciprocal of the condition number (in the 1-norm) of a
7       complex Hermitian positive definite band matrix using the Cholesky fac‐
8       torization A = U**H*U or A = L*L**H computed by ZPBTRF
9

SYNOPSIS

11       SUBROUTINE ZPBCON( UPLO,  N,  KD,  AB, LDAB, ANORM, RCOND, WORK, RWORK,
12                          INFO )
13
14           CHARACTER      UPLO
15
16           INTEGER        INFO, KD, LDAB, N
17
18           DOUBLE         PRECISION ANORM, RCOND
19
20           DOUBLE         PRECISION RWORK( * )
21
22           COMPLEX*16     AB( LDAB, * ), WORK( * )
23

PURPOSE

25       ZPBCON estimates the reciprocal of the condition number (in the 1-norm)
26       of a complex Hermitian positive definite band matrix using the Cholesky
27       factorization A = U**H*U or A = L*L**H computed by ZPBTRF.
28
29       An estimate is obtained for norm(inv(A)), and  the  reciprocal  of  the
30       condition number is computed as RCOND = 1 / (ANORM * norm(inv(A))).
31
32

ARGUMENTS

34       UPLO    (input) CHARACTER*1
35               = 'U':  Upper triangular factor stored in AB;
36               = 'L':  Lower triangular factor stored in AB.
37
38       N       (input) INTEGER
39               The order of the matrix A.  N >= 0.
40
41       KD      (input) INTEGER
42               The  number of superdiagonals of the matrix A if UPLO = 'U', or
43               the number of sub-diagonals if UPLO = 'L'.  KD >= 0.
44
45       AB      (input) COMPLEX*16 array, dimension (LDAB,N)
46               The triangular factor U or L from the Cholesky factorization  A
47               =  U**H*U  or  A  =  L*L**H of the band matrix A, stored in the
48               first KD+1 rows of the array.  The j-th column of  U  or  L  is
49               stored  in  the j-th column of the array AB as follows: if UPLO
50               ='U', AB(kd+1+i-j,j) = U(i,j) for  max(1,j-kd)<=i<=j;  if  UPLO
51               ='L', AB(1+i-j,j)    = L(i,j) for j<=i<=min(n,j+kd).
52
53       LDAB    (input) INTEGER
54               The leading dimension of the array AB.  LDAB >= KD+1.
55
56       ANORM   (input) DOUBLE PRECISION
57               The 1-norm (or infinity-norm) of the Hermitian band matrix A.
58
59       RCOND   (output) DOUBLE PRECISION
60               The  reciprocal  of  the condition number of the matrix A, com‐
61               puted as RCOND = 1/(ANORM * AINVNM), where AINVNM is  an  esti‐
62               mate of the 1-norm of inv(A) computed in this routine.
63
64       WORK    (workspace) COMPLEX*16 array, dimension (2*N)
65
66       RWORK   (workspace) DOUBLE PRECISION array, dimension (N)
67
68       INFO    (output) INTEGER
69               = 0:  successful exit
70               < 0:  if INFO = -i, the i-th argument had an illegal value
71
72
73
74 LAPACK routine (version 3.1)    November 2006                       ZPBCON(1)
Impressum