TESTING/LIN/csyt03.f(3)

NAME
TESTING/LIN/csyt03.f

SYNOPSIS

  Functions/Subroutines
subroutine csyt03 (uplo, n, a, lda, ainv, ldainv, work, ldwork, rwork, rcond, resid)
CSYT03

Function/Subroutine Documentation

  subroutine csyt03 (character uplo, integer n, complex, dimension( lda, * ) a, integer lda, complex, dimension( ldainv, * ) ainv, integer ldainv, complex, dimension( ldwork, * ) work, integer ldwork, real, dimension( * ) rwork, real rcond, real resid)

CSYT03

Purpose:

    !>
    !> CSYT03 computes the residual for a complex symmetric matrix times
    !> its inverse:
    !>    norm( I - A*AINV ) / ( N * norm(A) * norm(AINV) * EPS )
    !> where EPS is the machine epsilon.
    !>

Parameters
    UPLO

    !>          UPLO is CHARACTER*1
    !>          Specifies whether the upper or lower triangular part of the
    !>          complex symmetric matrix A is stored:
    !>          = 'U':  Upper triangular
    !>          = 'L':  Lower triangular
    !>

    N

    !>          N is INTEGER
    !>          The number of rows and columns of the matrix A.  N >= 0.
    !>

    A

    !>          A is COMPLEX array, dimension (LDA,N)
    !>          The original complex symmetric matrix A.
    !>

    LDA

    !>          LDA is INTEGER
    !>          The leading dimension of the array A.  LDA >= max(1,N)
    !>

    AINV

    !>          AINV is COMPLEX array, dimension (LDAINV,N)
    !>          On entry, the inverse of the matrix A, stored as a symmetric
    !>          matrix in the same format as A.
    !>          In this version, AINV is expanded into a full matrix and
    !>          multiplied by A, so the opposing triangle of AINV will be
    !>          changed; i.e., if the upper triangular part of AINV is
    !>          stored, the lower triangular part will be used as work space.
    !>

    LDAINV

    !>          LDAINV is INTEGER
    !>          The leading dimension of the array AINV.  LDAINV >= max(1,N).
    !>

    WORK

    !>          WORK is COMPLEX array, dimension (LDWORK,N)
    !>

    LDWORK

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

    RWORK

    !>          RWORK is REAL array, dimension (N)
    !>

    RCOND

    !>          RCOND is REAL
    !>          The reciprocal of the condition number of A, computed as
    !>          RCOND = 1/ (norm(A) * norm(AINV)).
    !>

    RESID

    !>          RESID is REAL
    !>          norm(I - A*AINV) / ( N * norm(A) * norm(AINV) * EPS )
    !>

Author
    Univ. of Tennessee

    Univ. of California Berkeley

    Univ. of Colorado Denver

    NAG Ltd.

Definition at line 124 of file csyt03.f.

Author

Generated automatically by Doxygen for LAPACK from the source code.