f04lef solves a system of tridiagonal equations following the factorization by
f01lef. This routine is intended for applications such as inverse iteration as well as straightforward linear equation applications.
Following the factorization of the
$n\times n$ tridiagonal matrix
$(T-\lambda I)$ as
by
f01lef,
f04lef may be used to solve any of the equations
for
$x$, the choice of equation being controlled by the argument
job. In each case there is an option to perturb zero or very small diagonal elements of
$U$, this option being intended for use in applications such as inverse iteration.
If on entry
${\mathbf{ifail}}=0$ or
$\mathrm{-1}$, explanatory error messages are output on the current error message unit (as defined by
x04aaf).
The computed solution of the equations
$(T-\lambda I)x=y$, say
$\overline{x}$, will satisfy an equation of the form
where
$E$ can be expected to satisfy a bound of the form
$\alpha $ being a modest constant and
$\epsilon $ being the
machine precision. The computed solution of the equations
${(T-\lambda I)}^{\mathrm{T}}x=y$ and
$Ux=y$ will satisfy similar results. The above result implies that the relative error in
$\overline{x}$ satisfies
where
$c(T-\lambda I)$ is the condition number of
$(T-\lambda I)$ with respect to inversion. Thus if
$(T-\lambda I)$ is nearly singular,
$\overline{x}$ can be expected to have a large relative error. Note that
f01lef incorporates a test for near singularity.
If you have single systems of tridiagonal equations to solve you are advised that
f07caf requires less storage and will normally be faster than the combination of
f01lef and
f04lef, but
f07caf does not incorporate a test for near singularity.
This example solves the two sets of tridiagonal equations
where
The example program first calls
f01lef to factorize
$T$ and then two calls are made to
f04lef, one to solve
$Tx=y$ and the second to solve
${T}^{\mathrm{T}}x=y$.