TESTING/LIN/zgtt01.f(3)

NAME
TESTING/LIN/zgtt01.f

SYNOPSIS

  Functions/Subroutines
subroutine zgtt01 (n, dl, d, du, dlf, df, duf, du2, ipiv, work, ldwork, rwork, resid)
ZGTT01

Function/Subroutine Documentation

  subroutine zgtt01 (integer n, complex*16, dimension( * ) dl, complex*16, dimension( * ) d, complex*16, dimension( * ) du, complex*16, dimension( * ) dlf, complex*16, dimension( * ) df, complex*16, dimension( * ) duf, complex*16, dimension( * ) du2, integer, dimension( * ) ipiv, complex*16, dimension( ldwork, * ) work, integer ldwork, double precision, dimension( * ) rwork, double precision resid)

ZGTT01

Purpose:

     ZGTT01 reconstructs a tridiagonal matrix A from its LU factorization
     and computes the residual
        norm(L*U - A) / ( norm(A) * EPS ),
     where EPS is the machine epsilon.

Parameters
    N

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

    DL

              DL is COMPLEX*16 array, dimension (N-1)
              The (n-1) sub-diagonal elements of A.

    D

              D is COMPLEX*16 array, dimension (N)
              The diagonal elements of A.

    DU

              DU is COMPLEX*16 array, dimension (N-1)
              The (n-1) super-diagonal elements of A.

    DLF

              DLF is COMPLEX*16 array, dimension (N-1)
              The (n-1) multipliers that define the matrix L from the
              LU factorization of A.

    DF

              DF is COMPLEX*16 array, dimension (N)
              The n diagonal elements of the upper triangular matrix U from
              the LU factorization of A.

    DUF

              DUF is COMPLEX*16 array, dimension (N-1)
              The (n-1) elements of the first super-diagonal of U.

    DU2

              DU2 is COMPLEX*16 array, dimension (N-2)
              The (n-2) elements of the second super-diagonal of U.

    IPIV

              IPIV is INTEGER array, dimension (N)
              The pivot indices; for 1 <= i <= n, row i of the matrix was
              interchanged with row IPIV(i).  IPIV(i) will always be either
              i or i+1; IPIV(i) = i indicates a row interchange was not
              required.

    WORK

              WORK is COMPLEX*16 array, dimension (LDWORK,N)

    LDWORK

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

    RWORK

              RWORK is DOUBLE PRECISION array, dimension (N)

    RESID

              RESID is DOUBLE PRECISION
              The scaled residual:  norm(L*U - A) / (norm(A) * EPS)

Author
    Univ. of Tennessee

    Univ. of California Berkeley

    Univ. of Colorado Denver

    NAG Ltd.

Definition at line 132 of file zgtt01.f.

Author

Generated automatically by Doxygen for LAPACK from the source code.