.TH "TESTING/EIG/zunt03.f" 3 "Version 3.12.0" "LAPACK" \" -*- nroff -*-
.ad l
.nh
.SH NAME
TESTING/EIG/zunt03.f
.SH SYNOPSIS
.br
.PP
.SS "Functions/Subroutines"

.in +1c
.ti -1c
.RI "subroutine \fBzunt03\fP (rc, mu, mv, n, k, u, ldu, v, ldv, work, lwork, rwork, result, info)"
.br
.RI "\fBZUNT03\fP "
.in -1c
.SH "Function/Subroutine Documentation"
.PP 
.SS "subroutine zunt03 (character*( * ) rc, integer mu, integer mv, integer n, integer k, complex*16, dimension( ldu, * ) u, integer ldu, complex*16, dimension( ldv, * ) v, integer ldv, complex*16, dimension( * ) work, integer lwork, double precision, dimension( * ) rwork, double precision result, integer info)"

.PP
\fBZUNT03\fP 
.PP
\fBPurpose:\fP
.RS 4

.PP
.nf
 ZUNT03 compares two unitary matrices U and V to see if their
 corresponding rows or columns span the same spaces\&.  The rows are
 checked if RC = 'R', and the columns are checked if RC = 'C'\&.

 RESULT is the maximum of

    | V*V' - I | / ( MV ulp ), if RC = 'R', or

    | V'*V - I | / ( MV ulp ), if RC = 'C',

 and the maximum over rows (or columns) 1 to K of

    | U(i) - S*V(i) |/ ( N ulp )

 where abs(S) = 1 (chosen to minimize the expression), U(i) is the
 i-th row (column) of U, and V(i) is the i-th row (column) of V\&.
.fi
.PP
 
.RE
.PP
\fBParameters\fP
.RS 4
\fIRC\fP 
.PP
.nf
          RC is CHARACTER*1
          If RC = 'R' the rows of U and V are to be compared\&.
          If RC = 'C' the columns of U and V are to be compared\&.
.fi
.PP
.br
\fIMU\fP 
.PP
.nf
          MU is INTEGER
          The number of rows of U if RC = 'R', and the number of
          columns if RC = 'C'\&.  If MU = 0 ZUNT03 does nothing\&.
          MU must be at least zero\&.
.fi
.PP
.br
\fIMV\fP 
.PP
.nf
          MV is INTEGER
          The number of rows of V if RC = 'R', and the number of
          columns if RC = 'C'\&.  If MV = 0 ZUNT03 does nothing\&.
          MV must be at least zero\&.
.fi
.PP
.br
\fIN\fP 
.PP
.nf
          N is INTEGER
          If RC = 'R', the number of columns in the matrices U and V,
          and if RC = 'C', the number of rows in U and V\&.  If N = 0
          ZUNT03 does nothing\&.  N must be at least zero\&.
.fi
.PP
.br
\fIK\fP 
.PP
.nf
          K is INTEGER
          The number of rows or columns of U and V to compare\&.
          0 <= K <= max(MU,MV)\&.
.fi
.PP
.br
\fIU\fP 
.PP
.nf
          U is COMPLEX*16 array, dimension (LDU,N)
          The first matrix to compare\&.  If RC = 'R', U is MU by N, and
          if RC = 'C', U is N by MU\&.
.fi
.PP
.br
\fILDU\fP 
.PP
.nf
          LDU is INTEGER
          The leading dimension of U\&.  If RC = 'R', LDU >= max(1,MU),
          and if RC = 'C', LDU >= max(1,N)\&.
.fi
.PP
.br
\fIV\fP 
.PP
.nf
          V is COMPLEX*16 array, dimension (LDV,N)
          The second matrix to compare\&.  If RC = 'R', V is MV by N, and
          if RC = 'C', V is N by MV\&.
.fi
.PP
.br
\fILDV\fP 
.PP
.nf
          LDV is INTEGER
          The leading dimension of V\&.  If RC = 'R', LDV >= max(1,MV),
          and if RC = 'C', LDV >= max(1,N)\&.
.fi
.PP
.br
\fIWORK\fP 
.PP
.nf
          WORK is COMPLEX*16 array, dimension (LWORK)
.fi
.PP
.br
\fILWORK\fP 
.PP
.nf
          LWORK is INTEGER
          The length of the array WORK\&.  For best performance, LWORK
          should be at least N*N if RC = 'C' or M*M if RC = 'R', but
          the tests will be done even if LWORK is 0\&.
.fi
.PP
.br
\fIRWORK\fP 
.PP
.nf
          RWORK is DOUBLE PRECISION array, dimension (max(MV,N))
.fi
.PP
.br
\fIRESULT\fP 
.PP
.nf
          RESULT is DOUBLE PRECISION
          The value computed by the test described above\&.  RESULT is
          limited to 1/ulp to avoid overflow\&.
.fi
.PP
.br
\fIINFO\fP 
.PP
.nf
          INFO is INTEGER
          0  indicates a successful exit
          -k indicates the k-th parameter had an illegal value
.fi
.PP
 
.RE
.PP
\fBAuthor\fP
.RS 4
Univ\&. of Tennessee 
.PP
Univ\&. of California Berkeley 
.PP
Univ\&. of Colorado Denver 
.PP
NAG Ltd\&. 
.RE
.PP

.PP
Definition at line \fB160\fP of file \fBzunt03\&.f\fP\&.
.SH "Author"
.PP 
Generated automatically by Doxygen for LAPACK from the source code\&.