F07PGF (DSPCON) estimates the condition number of a real symmetric indefinite matrix
$A$, where
$A$ has been factorized by
F07PDF (DSPTRF), using packed storage.
F07PGF (DSPCON) estimates the condition number (in the
$1$-norm) of a real symmetric indefinite matrix
$A$:
Since
$A$ is symmetric,
${\kappa}_{1}\left(A\right)={\kappa}_{\infty}\left(A\right)={\Vert A\Vert}_{\infty}{\Vert {A}^{-1}\Vert}_{\infty}$.
The routine should be preceded by a call to
F06RDF to compute
${\Vert A\Vert}_{1}$ and a call to
F07PDF (DSPTRF) to compute the Bunch–Kaufman factorization of
$A$. The routine then uses Higham's implementation of Hager's method (see
Higham (1988)) to estimate
${\Vert {A}^{-1}\Vert}_{1}$.
Higham N J (1988) FORTRAN codes for estimating the one-norm of a real or complex matrix, with applications to condition estimation ACM Trans. Math. Software 14 381–396
The computed estimate
RCOND is never less than the true value
$\rho $, and in practice is nearly always less than
$10\rho $, although examples can be constructed where
RCOND is much larger.
A call to F07PGF (DSPCON) involves solving a number of systems of linear equations of the form
$Ax=b$; the number is usually
$4$ or
$5$ and never more than
$11$. Each solution involves approximately
$2{n}^{2}$ floating point operations but takes considerably longer than a call to
F07PEF (DSPTRS) with one right-hand side, because extra care is taken to avoid overflow when
$A$ is approximately singular.
The complex analogues of this routine are
F07PUF (ZHPCON) for Hermitian matrices and
F07QUF (ZSPCON) for symmetric matrices.
This example estimates the condition number in the
$1$-norm (or
$\infty $-norm) of the matrix
$A$, where
Here
$A$ is symmetric indefinite, stored in packed form, and must first be factorized by
F07PDF (DSPTRF). The true condition number in the
$1$-norm is
$75.68$.