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

NAME

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

SYNOPSIS

10       SUBROUTINE CTRCON( NORM, UPLO, DIAG, N, A,  LDA,  RCOND,  WORK,  RWORK,
11                          INFO )
12
13           CHARACTER      DIAG, NORM, UPLO
14
15           INTEGER        INFO, LDA, N
16
17           REAL           RCOND
18
19           REAL           RWORK( * )
20
21           COMPLEX        A( LDA, * ), WORK( * )
22

PURPOSE

24       CTRCON estimates the reciprocal of the condition number of a triangular
25       matrix A, in either the 1-norm or the infinity-norm.
26
27       The norm of A is computed and an estimate is obtained for norm(inv(A)),
28       then the reciprocal of the condition number is computed as
29          RCOND = 1 / ( norm(A) * norm(inv(A)) ).
30
31

ARGUMENTS

33       NORM    (input) CHARACTER*1
34               Specifies  whether the 1-norm condition number or the infinity-
35               norm condition number is required:
36               = '1' or 'O':  1-norm;
37               = 'I':         Infinity-norm.
38
39       UPLO    (input) CHARACTER*1
40               = 'U':  A is upper triangular;
41               = 'L':  A is lower triangular.
42
43       DIAG    (input) CHARACTER*1
44               = 'N':  A is non-unit triangular;
45               = 'U':  A is unit triangular.
46
47       N       (input) INTEGER
48               The order of the matrix A.  N >= 0.
49
50       A       (input) COMPLEX array, dimension (LDA,N)
51               The triangular matrix A.  If UPLO =  'U',  the  leading  N-by-N
52               upper  triangular part of the array A contains the upper trian‐
53               gular matrix, and the strictly lower triangular part  of  A  is
54               not referenced.  If UPLO = 'L', the leading N-by-N lower trian‐
55               gular part of the array A contains the lower triangular matrix,
56               and  the strictly upper triangular part of A is not referenced.
57               If DIAG = 'U', the diagonal elements of A are also  not  refer‐
58               enced and are assumed to be 1.
59
60       LDA     (input) INTEGER
61               The leading dimension of the array A.  LDA >= max(1,N).
62
63       RCOND   (output) REAL
64               The  reciprocal  of  the condition number of the matrix A, com‐
65               puted as RCOND = 1/(norm(A) * norm(inv(A))).
66
67       WORK    (workspace) COMPLEX array, dimension (2*N)
68
69       RWORK   (workspace) REAL array, dimension (N)
70
71       INFO    (output) INTEGER
72               = 0:  successful exit
73               < 0:  if INFO = -i, the i-th argument had an illegal value
74
75
76
77 LAPACK routine (version 3.1)    November 2006                       CTRCON(1)
Impressum