TESTING/LIN/zsyt03.f(3)    Library Functions Manual    TESTING/LIN/zsyt03.f(3)

NAME
       TESTING/LIN/zsyt03.f

SYNOPSIS
   Functions/Subroutines
       subroutine zsyt03 (uplo, n, a, lda, ainv, ldainv, work, ldwork, rwork,
           rcond, resid)
           ZSYT03

Function/Subroutine Documentation
   subroutine zsyt03 (character uplo, integer n, complex*16, dimension( lda, *
       ) a, integer lda, complex*16, dimension( ldainv, * ) ainv, integer
       ldainv, complex*16, dimension( ldwork, * ) work, integer ldwork, double
       precision, dimension( * ) rwork, double precision rcond, double
       precision resid)
       ZSYT03

       Purpose:


           !>
           !> ZSYT03 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*16 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*16 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*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)
           !>

           RCOND

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

           RESID

           !>          RESID is DOUBLE PRECISION
           !>          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 zsyt03.f.

Author
       Generated automatically by Doxygen for LAPACK from the source code.

LAPACK                          Version 3.12.0         TESTING/LIN/zsyt03.f(3)