ZGTCON(3) LAPACK routine of NEC Numeric Library Collection ZGTCON(3)
NAME
ZGTCON
SYNOPSIS
SUBROUTINE ZGTCON (NORM, N, DL, D, DU, DU2, IPIV, ANORM, RCOND, WORK,
INFO)
PURPOSE
ZGTCON estimates the reciprocal of the condition number of a complex
tridiagonal matrix A using the LU factorization as computed by
ZGTTRF.
An estimate is obtained for norm(inv(A)), and the reciprocal of the
condition number is computed as RCOND = 1 / (ANORM * norm(inv(A))).
ARGUMENTS
NORM (input)
NORM is CHARACTER*1
Specifies whether the 1-norm condition number or the
infinity-norm condition number is required:
= '1' or 'O': 1-norm;
= 'I': Infinity-norm.
N (input)
N is INTEGER
The order of the matrix A. N >= 0.
DL (input)
DL is COMPLEX*16 array, dimension (N-1)
The (n-1) multipliers that define the matrix L from the
LU factorization of A as computed by ZGTTRF.
D (input)
D is COMPLEX*16 array, dimension (N)
The n diagonal elements of the upper triangular matrix U from
the LU factorization of A.
DU (input)
DU is COMPLEX*16 array, dimension (N-1)
The (n-1) elements of the first superdiagonal of U.
DU2 (input)
DU2 is COMPLEX*16 array, dimension (N-2)
The (n-2) elements of the second superdiagonal of U.
IPIV (input)
IPIV is INTEGER array, dimension (N)
The pivot indices; for 1 <= i <= n, row i of the matrix was
interchanged with row IPIV(i). IPIV(i) will always be either
i or i+1; IPIV(i) = i indicates a row interchange was not
required.
ANORM (input)
ANORM is DOUBLE PRECISION
If NORM = '1' or 'O', the 1-norm of the original matrix A.
If NORM = 'I', the infinity-norm of the original matrix A.
RCOND (output)
RCOND is DOUBLE PRECISION
The reciprocal of the condition number of the matrix A,
computed as RCOND = 1/(ANORM * AINVNM), where AINVNM is an
estimate of the 1-norm of inv(A) computed in this routine.
WORK (output)
WORK is COMPLEX*16 array, dimension (2*N)
INFO (output)
INFO is INTEGER
= 0: successful exit
< 0: if INFO = -i, the i-th argument had an illegal value
LAPACK routine 31 October 2017 ZGTCON(3)