1ZPOCON(1) LAPACK routine (version 3.2) ZPOCON(1)
2
3
4
6 ZPOCON - estimates the reciprocal of the condition number (in the
7 1-norm) of a complex Hermitian positive definite matrix using the
8 Cholesky factorization A = U**H*U or A = L*L**H computed by ZPOTRF
9
11 SUBROUTINE ZPOCON( UPLO, N, A, LDA, ANORM, RCOND, WORK, RWORK, INFO )
12
13 CHARACTER UPLO
14
15 INTEGER INFO, LDA, N
16
17 DOUBLE PRECISION ANORM, RCOND
18
19 DOUBLE PRECISION RWORK( * )
20
21 COMPLEX*16 A( LDA, * ), WORK( * )
22
24 ZPOCON estimates the reciprocal of the condition number (in the 1-norm)
25 of a complex Hermitian positive definite matrix using the Cholesky fac‐
26 torization A = U**H*U or A = L*L**H computed by ZPOTRF. An estimate is
27 obtained for norm(inv(A)), and the reciprocal of the condition number
28 is computed as RCOND = 1 / (ANORM * norm(inv(A))).
29
31 UPLO (input) CHARACTER*1
32 = 'U': Upper triangle of A is stored;
33 = 'L': Lower triangle of A is stored.
34
35 N (input) INTEGER
36 The order of the matrix A. N >= 0.
37
38 A (input) COMPLEX*16 array, dimension (LDA,N)
39 The triangular factor U or L from the Cholesky factorization A
40 = U**H*U or A = L*L**H, as computed by ZPOTRF.
41
42 LDA (input) INTEGER
43 The leading dimension of the array A. LDA >= max(1,N).
44
45 ANORM (input) DOUBLE PRECISION
46 The 1-norm (or infinity-norm) of the Hermitian matrix A.
47
48 RCOND (output) DOUBLE PRECISION
49 The reciprocal of the condition number of the matrix A, com‐
50 puted as RCOND = 1/(ANORM * AINVNM), where AINVNM is an esti‐
51 mate of the 1-norm of inv(A) computed in this routine.
52
53 WORK (workspace) COMPLEX*16 array, dimension (2*N)
54
55 RWORK (workspace) DOUBLE PRECISION array, dimension (N)
56
57 INFO (output) INTEGER
58 = 0: successful exit
59 < 0: if INFO = -i, the i-th argument had an illegal value
60
61
62
63 LAPACK routine (version 3.2) November 2008 ZPOCON(1)