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

NAME

6       ZTPCON  -  estimates the reciprocal of the condition number of a packed
7       triangular 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.  The
25       norm of A is computed and an estimate  is  obtained  for  norm(inv(A)),
26       then the reciprocal of the condition number is computed as
27          RCOND = 1 / ( norm(A) * norm(inv(A)) ).
28

ARGUMENTS

30       NORM    (input) CHARACTER*1
31               Specifies  whether the 1-norm condition number or the infinity-
32               norm condition number is required:
33               = '1' or 'O':  1-norm;
34               = 'I':         Infinity-norm.
35
36       UPLO    (input) CHARACTER*1
37               = 'U':  A is upper triangular;
38               = 'L':  A is lower triangular.
39
40       DIAG    (input) CHARACTER*1
41               = 'N':  A is non-unit triangular;
42               = 'U':  A is unit triangular.
43
44       N       (input) INTEGER
45               The order of the matrix A.  N >= 0.
46
47       AP      (input) COMPLEX*16 array, dimension (N*(N+1)/2)
48               The upper or lower triangular matrix A, packed columnwise in  a
49               linear  array.   The j-th column of A is stored in the array AP
50               as follows: if UPLO = 'U',  AP(i  +  (j-1)*j/2)  =  A(i,j)  for
51               1<=i<=j;  if  UPLO  =  'L', AP(i + (j-1)*(2n-j)/2) = A(i,j) for
52               j<=i<=n.  If DIAG = 'U', the diagonal elements  of  A  are  not
53               referenced and are assumed to be 1.
54
55       RCOND   (output) DOUBLE PRECISION
56               The  reciprocal  of  the condition number of the matrix A, com‐
57               puted as RCOND = 1/(norm(A) * norm(inv(A))).
58
59       WORK    (workspace) COMPLEX*16 array, dimension (2*N)
60
61       RWORK   (workspace) DOUBLE PRECISION array, dimension (N)
62
63       INFO    (output) INTEGER
64               = 0:  successful exit
65               < 0:  if INFO = -i, the i-th argument had an illegal value
66
67
68
69 LAPACK routine (version 3.2)    November 2008                       ZTPCON(1)
Impressum