DGTSL(3F)DGTSL(3F)NAMEDGTSL - DGTSL given a general tridiagonal matrix and a right hand side
will find the solution.
SYNOPSYS
SUBROUTINE DGTSL(N,C,D,E,B,INFO)
DESCRIPTION
On Entry
N INTEGER
is the order of the tridiagonal matrix.
C DOUBLE PRECISION(N)
is the subdiagonal of the tridiagonal matrix.
C(2) through C(N) should contain the subdiagonal.
On output C is destroyed.
D DOUBLE PRECISION(N)
is the diagonal of the tridiagonal matrix.
On output D is destroyed.
E DOUBLE PRECISION(N)
is the superdiagonal of the tridiagonal matrix.
E(1) through E(N-1) should contain the superdiagonal.
On output E is destroyed.
B DOUBLE PRECISION(N)
is the right hand side vector. On Return
B is the solution vector.
INFO INTEGER
= 0 normal value.
= K if the K-th element of the diagonal becomes
exactly zero. The subroutine returns when
this is detected. LINPACK. This version dated 08/14/78 . Jack
Dongarra, Argonne National Laboratory. No externals Fortran DABS
Page 1