cptcon.f(3) LAPACK cptcon.f(3)NAMEcptcon.f-
subroutine cptcon (N, D, E, ANORM, RCOND, RWORK, INFO)
subroutine cptcon (integerN, real, dimension( * )D, complex, dimension( *
)E, realANORM, realRCOND, real, dimension( * )RWORK, integerINFO)
CPTCON computes the reciprocal of the condition number (in the
1-norm) of a complex Hermitian positive definite tridiagonal matrix
using the factorization A = L*D*L**H or A = U**H*D*U computed by
Norm(inv(A)) is computed by a direct method, and the reciprocal of
the condition number is computed as
RCOND = 1 / (ANORM * norm(inv(A))).
N is INTEGER
The order of the matrix A. N >= 0.
D is REAL array, dimension (N)
The n diagonal elements of the diagonal matrix D from the
factorization of A, as computed by CPTTRF.
E is COMPLEX array, dimension (N-1)
The (n-1) off-diagonal elements of the unit bidiagonal factor
U or L from the factorization of A, as computed by CPTTRF.
ANORM is REAL
The 1-norm of the original matrix A.
RCOND is REAL
The reciprocal of the condition number of the matrix A,
computed as RCOND = 1/(ANORM * AINVNM), where AINVNM is the
1-norm of inv(A) computed in this routine.
RWORK is REAL array, dimension (N)
INFO is INTEGER
= 0: successful exit
< 0: if INFO = -i, the i-th argument had an illegal value
Univ. of Tennessee
Univ. of California Berkeley
Univ. of Colorado Denver
The method used is described in Nicholas J. Higham, "Efficient
Algorithms for Computing the Condition Number of a Tridiagonal
Matrix", SIAM J. Sci. Stat. Comput., Vol. 7, No. 1, January 1986.
Definition at line 120 of file cptcon.f.
Generated automatically by Doxygen for LAPACK from the source code.
Version 3.4.2 Tue Sep 25 2012 cptcon.f(3)