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

.in +1c
.ti -1c
.RI "subroutine \fBchbt21\fP (uplo, n, ka, ks, a, lda, d, e, u, ldu, work, rwork, result)"
.br
.RI "\fBCHBT21\fP "
.in -1c
.SH "Function/Subroutine Documentation"
.PP 
.SS "subroutine chbt21 (character uplo, integer n, integer ka, integer ks, complex, dimension( lda, * ) a, integer lda, real, dimension( * ) d, real, dimension( * ) e, complex, dimension( ldu, * ) u, integer ldu, complex, dimension( * ) work, real, dimension( * ) rwork, real, dimension( 2 ) result)"

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

.PP
.nf
 CHBT21  generally checks a decomposition of the form

         A = U S U**H

 where **H means conjugate transpose, A is hermitian banded, U is
 unitary, and S is diagonal (if KS=0) or symmetric
 tridiagonal (if KS=1)\&.

 Specifically:

         RESULT(1) = | A - U S U**H | / ( |A| n ulp ) and
         RESULT(2) = | I - U U**H | / ( n ulp )
.fi
.PP
 
.RE
.PP
\fBParameters\fP
.RS 4
\fIUPLO\fP 
.PP
.nf
          UPLO is CHARACTER
          If UPLO='U', the upper triangle of A and V will be used and
          the (strictly) lower triangle will not be referenced\&.
          If UPLO='L', the lower triangle of A and V will be used and
          the (strictly) upper triangle will not be referenced\&.
.fi
.PP
.br
\fIN\fP 
.PP
.nf
          N is INTEGER
          The size of the matrix\&.  If it is zero, CHBT21 does nothing\&.
          It must be at least zero\&.
.fi
.PP
.br
\fIKA\fP 
.PP
.nf
          KA is INTEGER
          The bandwidth of the matrix A\&.  It must be at least zero\&.  If
          it is larger than N-1, then max( 0, N-1 ) will be used\&.
.fi
.PP
.br
\fIKS\fP 
.PP
.nf
          KS is INTEGER
          The bandwidth of the matrix S\&.  It may only be zero or one\&.
          If zero, then S is diagonal, and E is not referenced\&.  If
          one, then S is symmetric tri-diagonal\&.
.fi
.PP
.br
\fIA\fP 
.PP
.nf
          A is COMPLEX array, dimension (LDA, N)
          The original (unfactored) matrix\&.  It is assumed to be
          hermitian, and only the upper (UPLO='U') or only the lower
          (UPLO='L') will be referenced\&.
.fi
.PP
.br
\fILDA\fP 
.PP
.nf
          LDA is INTEGER
          The leading dimension of A\&.  It must be at least 1
          and at least min( KA, N-1 )\&.
.fi
.PP
.br
\fID\fP 
.PP
.nf
          D is REAL array, dimension (N)
          The diagonal of the (symmetric tri-) diagonal matrix S\&.
.fi
.PP
.br
\fIE\fP 
.PP
.nf
          E is REAL array, dimension (N-1)
          The off-diagonal of the (symmetric tri-) diagonal matrix S\&.
          E(1) is the (1,2) and (2,1) element, E(2) is the (2,3) and
          (3,2) element, etc\&.
          Not referenced if KS=0\&.
.fi
.PP
.br
\fIU\fP 
.PP
.nf
          U is COMPLEX array, dimension (LDU, N)
          The unitary matrix in the decomposition, expressed as a
          dense matrix (i\&.e\&., not as a product of Householder
          transformations, Givens transformations, etc\&.)
.fi
.PP
.br
\fILDU\fP 
.PP
.nf
          LDU is INTEGER
          The leading dimension of U\&.  LDU must be at least N and
          at least 1\&.
.fi
.PP
.br
\fIWORK\fP 
.PP
.nf
          WORK is COMPLEX array, dimension (N**2)
.fi
.PP
.br
\fIRWORK\fP 
.PP
.nf
          RWORK is REAL array, dimension (N)
.fi
.PP
.br
\fIRESULT\fP 
.PP
.nf
          RESULT is REAL array, dimension (2)
          The values computed by the two tests described above\&.  The
          values are currently limited to 1/ulp, to avoid overflow\&.
.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 \fB150\fP of file \fBchbt21\&.f\fP\&.
.SH "Author"
.PP 
Generated automatically by Doxygen for LAPACK from the source code\&.