TESTING/EIG/zget22.f(3)

NAME
TESTING/EIG/zget22.f

SYNOPSIS

  Functions/Subroutines
subroutine zget22 (transa, transe, transw, n, a, lda, e, lde, w, work, rwork, result)
ZGET22

Function/Subroutine Documentation

  subroutine zget22 (character transa, character transe, character transw, integer n, complex*16, dimension( lda, * ) a, integer lda, complex*16, dimension( lde, * ) e, integer lde, complex*16, dimension( * ) w, complex*16, dimension( * ) work, double precision, dimension( * ) rwork, double precision, dimension( 2 ) result)

ZGET22

Purpose:

    !>
    !> ZGET22 does an eigenvector check.
    !>
    !> The basic test is:
    !>
    !>    RESULT(1) = | A E  -  E W | / ( |A| |E| ulp )
    !>
    !> using the 1-norm.  It also tests the normalization of E:
    !>
    !>    RESULT(2) = max | m-norm(E(j)) - 1 | / ( n ulp )
    !>                 j
    !>
    !> where E(j) is the j-th eigenvector, and m-norm is the max-norm of a
    !> vector.  The max-norm of a complex n-vector x in this case is the
    !> maximum of |re(x(i)| + |im(x(i)| over i = 1, ..., n.
    !>

Parameters
    TRANSA

    !>          TRANSA is CHARACTER*1
    !>          Specifies whether or not A is transposed.
    !>          = 'N':  No transpose
    !>          = 'T':  Transpose
    !>          = 'C':  Conjugate transpose
    !>

    TRANSE

    !>          TRANSE is CHARACTER*1
    !>          Specifies whether or not E is transposed.
    !>          = 'N':  No transpose, eigenvectors are in columns of E
    !>          = 'T':  Transpose, eigenvectors are in rows of E
    !>          = 'C':  Conjugate transpose, eigenvectors are in rows of E
    !>

    TRANSW

    !>          TRANSW is CHARACTER*1
    !>          Specifies whether or not W is transposed.
    !>          = 'N':  No transpose
    !>          = 'T':  Transpose, same as TRANSW = 'N'
    !>          = 'C':  Conjugate transpose, use -WI(j) instead of WI(j)
    !>

    N

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

    A

    !>          A is COMPLEX*16 array, dimension (LDA,N)
    !>          The matrix whose eigenvectors are in E.
    !>

    LDA

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

    E

    !>          E is COMPLEX*16 array, dimension (LDE,N)
    !>          The matrix of eigenvectors. If TRANSE = 'N', the eigenvectors
    !>          are stored in the columns of E, if TRANSE = 'T' or 'C', the
    !>          eigenvectors are stored in the rows of E.
    !>

    LDE

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

    W

    !>          W is COMPLEX*16 array, dimension (N)
    !>          The eigenvalues of A.
    !>

    WORK

    !>          WORK is COMPLEX*16 array, dimension (N*N)
    !>

    RWORK

    !>          RWORK is DOUBLE PRECISION array, dimension (N)
    !>

    RESULT

    !>          RESULT is DOUBLE PRECISION array, dimension (2)
    !>          RESULT(1) = | A E  -  E W | / ( |A| |E| ulp )
    !>          RESULT(2) = max | m-norm(E(j)) - 1 | / ( n ulp )
    !>                       j
    !>

Author
    Univ. of Tennessee

    Univ. of California Berkeley

    Univ. of Colorado Denver

    NAG Ltd.

Definition at line 142 of file zget22.f.

Author

Generated automatically by Doxygen for LAPACK from the source code.