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

NAME

6       ZTPCON  - the reciprocal of the condition number of a packed triangular
7       matrix A, in either the 1-norm or the infinity-norm
8

SYNOPSIS

10       SUBROUTINE ZTPCON( NORM, UPLO, DIAG, N, AP, RCOND, WORK, RWORK, INFO )
11
12           CHARACTER      DIAG, NORM, UPLO
13
14           INTEGER        INFO, N
15
16           DOUBLE         PRECISION RCOND
17
18           DOUBLE         PRECISION RWORK( * )
19
20           COMPLEX*16     AP( * ), WORK( * )
21

PURPOSE

23       ZTPCON estimates the reciprocal of the condition  number  of  a  packed
24       triangular matrix A, in either the 1-norm or the infinity-norm.
25
26       The norm of A is computed and an estimate is obtained for norm(inv(A)),
27       then the reciprocal of the condition number is computed as
28          RCOND = 1 / ( norm(A) * norm(inv(A)) ).
29
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       UPLO    (input) CHARACTER*1
39               = 'U':  A is upper triangular;
40               = 'L':  A is lower triangular.
41
42       DIAG    (input) CHARACTER*1
43               = 'N':  A is non-unit triangular;
44               = 'U':  A is unit triangular.
45
46       N       (input) INTEGER
47               The order of the matrix A.  N >= 0.
48
49       AP      (input) COMPLEX*16 array, dimension (N*(N+1)/2)
50               The  upper or lower triangular matrix A, packed columnwise in a
51               linear array.  The j-th column of A is stored in the  array  AP
52               as  follows:  if  UPLO  =  'U',  AP(i + (j-1)*j/2) = A(i,j) for
53               1<=i<=j; if UPLO = 'L', AP(i +  (j-1)*(2n-j)/2)  =  A(i,j)  for
54               j<=i<=n.   If  DIAG  =  'U', the diagonal elements of A are not
55               referenced and are assumed to be 1.
56
57       RCOND   (output) DOUBLE PRECISION
58               The reciprocal of the condition number of the  matrix  A,  com‐
59               puted as RCOND = 1/(norm(A) * norm(inv(A))).
60
61       WORK    (workspace) COMPLEX*16 array, dimension (2*N)
62
63       RWORK   (workspace) DOUBLE PRECISION array, dimension (N)
64
65       INFO    (output) INTEGER
66               = 0:  successful exit
67               < 0:  if INFO = -i, the i-th argument had an illegal value
68
69
70
71 LAPACK routine (version 3.1)    November 2006                       ZTPCON(1)
Impressum