f07jgf computes the reciprocal condition number of a real
symmetric positive definite tridiagonal matrix
, using the
factorization returned by
f07jdf.
f07jgf should be preceded by a call to
f07jdf, which computes a modified Cholesky factorization of the matrix
as
where
is a unit lower bidiagonal matrix and
is a diagonal matrix, with positive diagonal elements.
f07jgf then utilizes the factorization to compute
by a direct method, from which the reciprocal of the condition number of
,
is computed as
is returned, rather than
, since when
is singular
is infinite.
The computed condition number will be the exact condition number for a closely neighbouring matrix.
Background information to multithreading can be found in the
Multithreading documentation.
Please consult the
X06 Chapter Introduction for information on how to control and interrogate the OpenMP environment used within this routine. Please also consult the
Users' Note for your implementation for any additional implementation-specific information.
See Section 15.6 of
Higham (2002) for further details on computing the condition number of tridiagonal matrices.
The complex analogue of this routine is
f07juf.
This example computes the condition number of the symmetric positive definite tridiagonal matrix
given by