f07caf computes the solution to a real system of linear equations
where
is an
tridiagonal matrix and
and
are
matrices.
Note that equations
may be solved by interchanging the order of the arguments
du and
dl.
Anderson E, Bai Z, Bischof C, Blackford S, Demmel J, Dongarra J J, Du Croz J J, Greenbaum A, Hammarling S, McKenney A and Sorensen D (1999)
LAPACK Users' Guide (3rd Edition) SIAM, Philadelphia
https://www.netlib.org/lapack/lug
-
1:
– Integer
Input
-
On entry: , the number of linear equations, i.e., the order of the matrix .
Constraint:
.
-
2:
– Integer
Input
-
On entry: , the number of right-hand sides, i.e., the number of columns of the matrix .
Constraint:
.
-
3:
– Real (Kind=nag_wp) array
Input/Output
-
Note: the dimension of the array
dl
must be at least
.
On entry: must contain the subdiagonal elements of the matrix .
On exit: if no constraints are violated,
dl is overwritten by the (
) elements of the second superdiagonal of the upper triangular matrix
from the
factorization of
, in
.
-
4:
– Real (Kind=nag_wp) array
Input/Output
-
Note: the dimension of the array
d
must be at least
.
On entry: must contain the diagonal elements of the matrix .
On exit: if no constraints are violated,
d is overwritten by the
diagonal elements of the upper triangular matrix
from the
factorization of
.
-
5:
– Real (Kind=nag_wp) array
Input/Output
-
Note: the dimension of the array
du
must be at least
.
On entry: must contain the superdiagonal elements of the matrix .
On exit: if no constraints are violated,
du is overwritten by the
elements of the first superdiagonal of
.
-
6:
– Real (Kind=nag_wp) array
Input/Output
-
Note: the second dimension of the array
b
must be at least
.
to solve the equations
, where
is a single right-hand side,
b may be supplied as a one-dimensional array with length
.
On entry: the right-hand side matrix .
On exit: if , the solution matrix .
-
7:
– Integer
Input
-
On entry: the first dimension of the array
b as declared in the (sub)program from which
f07caf is called.
Constraint:
.
-
8:
– Integer
Output
-
On exit:
unless the routine detects an error (see
Section 6).
The computed solution for a single right-hand side,
, satisfies an equation of the form
where
and
is the
machine precision. An approximate error bound for the computed solution is given by
where
, the condition number of
with respect to the solution of the linear equations. See Section 4.4 of
Anderson et al. (1999) for further details.
Alternatives to
f07caf, which return condition and error estimates are
f04bcf and
f07cbf.
The complex analogue of this routine is
f07cnf.
This example solves the equations
where
is the tridiagonal matrix