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)