DGTT02 computes the residual for the solution to a tridiagonal
 system of equations:
    RESID = norm(B - op(A)*X) / (norm(op(A)) * norm(X) * EPS),
 where EPS is the machine epsilon.
 The norm used is the 1-norm.

TRANS

          TRANS is CHARACTER
          Specifies the form of the residual.
          = 'N':  B - A    * X  (No transpose)
          = 'T':  B - A**T * X  (Transpose)
          = 'C':  B - A**H * X  (Conjugate transpose = Transpose)

N

          N is INTEGER
          The order of the matrix A.  N >= 0.

NRHS

          NRHS is INTEGER
          The number of right hand sides, i.e., the number of columns
          of the matrices B and X.  NRHS >= 0.

DL

          DL is DOUBLE PRECISION array, dimension (N-1)
          The (n-1) sub-diagonal elements of A.

D

          D is DOUBLE PRECISION array, dimension (N)
          The diagonal elements of A.

DU

          DU is DOUBLE PRECISION array, dimension (N-1)
          The (n-1) super-diagonal elements of A.

X

          X is DOUBLE PRECISION array, dimension (LDX,NRHS)
          The computed solution vectors X.

LDX

          LDX is INTEGER
          The leading dimension of the array X.  LDX >= max(1,N).

B

          B is DOUBLE PRECISION array, dimension (LDB,NRHS)
          On entry, the right hand side vectors for the system of
          linear equations.
          On exit, B is overwritten with the difference B - op(A)*X.

LDB

          LDB is INTEGER
          The leading dimension of the array B.  LDB >= max(1,N).

RESID

          RESID is DOUBLE PRECISION
          norm(B - op(A)*X) / (norm(op(A)) * norm(X) * EPS)

Univ. of Tennessee

Univ. of California Berkeley

Univ. of Colorado Denver

NAG Ltd.