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

NAME

6       SGTCON  -  estimates  the  reciprocal of the condition number of a real
7       tridiagonal matrix A using the LU factorization as computed by SGTTRF
8

SYNOPSIS

10       SUBROUTINE SGTCON( NORM, N, DL, D, DU, DU2, IPIV, ANORM,  RCOND,  WORK,
11                          IWORK, INFO )
12
13           CHARACTER      NORM
14
15           INTEGER        INFO, N
16
17           REAL           ANORM, RCOND
18
19           INTEGER        IPIV( * ), IWORK( * )
20
21           REAL           D( * ), DL( * ), DU( * ), DU2( * ), WORK( * )
22

PURPOSE

24       SGTCON  estimates  the  reciprocal  of  the  condition number of a real
25       tridiagonal matrix A using the LU factorization as computed by  SGTTRF.
26       An  estimate  is  obtained  for norm(inv(A)), and the reciprocal of the
27       condition number is computed as RCOND = 1 / (ANORM * 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       N       (input) INTEGER
37               The order of the matrix A.  N >= 0.
38
39       DL      (input) REAL array, dimension (N-1)
40               The (n-1) multipliers that define the matrix L from the LU fac‐
41               torization of A as computed by SGTTRF.
42
43       D       (input) REAL array, dimension (N)
44               The n diagonal elements of the upper triangular matrix  U  from
45               the LU factorization of A.
46
47       DU      (input) REAL array, dimension (N-1)
48               The (n-1) elements of the first superdiagonal of U.
49
50       DU2     (input) REAL array, dimension (N-2)
51               The (n-2) elements of the second superdiagonal of U.
52
53       IPIV    (input) INTEGER array, dimension (N)
54               The  pivot  indices;  for  1 <= i <= n, row i of the matrix was
55               interchanged with row IPIV(i).  IPIV(i) will always be either i
56               or  i+1;  IPIV(i)  =  i  indicates  a  row  interchange was not
57               required.
58
59       ANORM   (input) REAL
60               If NORM = '1' or 'O', the 1-norm of the original matrix A.   If
61               NORM = 'I', the infinity-norm of the original matrix A.
62
63       RCOND   (output) REAL
64               The  reciprocal  of  the condition number of the matrix A, com‐
65               puted as RCOND = 1/(ANORM * AINVNM), where AINVNM is  an  esti‐
66               mate of the 1-norm of inv(A) computed in this routine.
67
68       WORK    (workspace) REAL array, dimension (2*N)
69
70       IWORK   (workspace) INTEGER array, dimension (N)
71
72       INFO    (output) INTEGER
73               = 0:  successful exit
74               < 0:  if INFO = -i, the i-th argument had an illegal value
75
76
77
78 LAPACK routine (version 3.2)    November 2008                       SGTCON(1)
Impressum