1DGTTRF(1)                LAPACK routine (version 3.2)                DGTTRF(1)
2
3
4

NAME

6       DGTTRF  -  computes  an LU factorization of a real tridiagonal matrix A
7       using elimination with partial pivoting and row interchanges
8

SYNOPSIS

10       SUBROUTINE DGTTRF( N, DL, D, DU, DU2, IPIV, INFO )
11
12           INTEGER        INFO, N
13
14           INTEGER        IPIV( * )
15
16           DOUBLE         PRECISION D( * ), DL( * ), DU( * ), DU2( * )
17

PURPOSE

19       DGTTRF computes an LU factorization of  a  real  tridiagonal  matrix  A
20       using elimination with partial pivoting and row interchanges.  The fac‐
21       torization has the form
22          A = L * U
23       where L is a product of permutation and unit lower bidiagonal  matrices
24       and  U  is upper triangular with nonzeros in only the main diagonal and
25       first two superdiagonals.
26

ARGUMENTS

28       N       (input) INTEGER
29               The order of the matrix A.
30
31       DL      (input/output) DOUBLE PRECISION array, dimension (N-1)
32               On entry, DL must contain the (n-1) sub-diagonal elements of A.
33               On exit, DL is overwritten by the (n-1) multipliers that define
34               the matrix L from the LU factorization of A.
35
36       D       (input/output) DOUBLE PRECISION array, dimension (N)
37               On entry, D must contain the diagonal elements of A.  On  exit,
38               D is overwritten by the n diagonal elements of the upper trian‐
39               gular matrix U from the LU factorization of A.
40
41       DU      (input/output) DOUBLE PRECISION array, dimension (N-1)
42               On entry, DU must contain the (n-1) super-diagonal elements  of
43               A.   On  exit,  DU  is overwritten by the (n-1) elements of the
44               first super-diagonal of U.
45
46       DU2     (output) DOUBLE PRECISION array, dimension (N-2)
47               On exit, DU2 is overwritten by the (n-2) elements of the second
48               super-diagonal of U.
49
50       IPIV    (output) INTEGER array, dimension (N)
51               The  pivot  indices;  for  1 <= i <= n, row i of the matrix was
52               interchanged with row IPIV(i).  IPIV(i) will always be either i
53               or  i+1;  IPIV(i)  =  i  indicates  a  row  interchange was not
54               required.
55
56       INFO    (output) INTEGER
57               = 0:  successful exit
58               < 0:  if INFO = -k, the k-th argument had an illegal value
59               > 0:  if INFO = k, U(k,k) is exactly  zero.  The  factorization
60               has  been  completed, but the factor U is exactly singular, and
61               division by zero will occur if it is used to solve a system  of
62               equations.
63
64
65
66 LAPACK routine (version 3.2)    November 2008                       DGTTRF(1)
Impressum