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

.in +1c
.ti -1c
.RI "subroutine \fBdstt21\fP (n, kband, ad, ae, sd, se, u, ldu, work, result)"
.br
.RI "\fBDSTT21\fP "
.in -1c
.SH "Function/Subroutine Documentation"
.PP 
.SS "subroutine dstt21 (integer n, integer kband, double precision, dimension( * ) ad, double precision, dimension( * ) ae, double precision, dimension( * ) sd, double precision, dimension( * ) se, double precision, dimension( ldu, * ) u, integer ldu, double precision, dimension( * ) work, double precision, dimension( 2 ) result)"

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

.PP
.nf
!>
!> DSTT21 checks a decomposition of the form
!>
!>    A = U S U'
!>
!> where ' means transpose, A is symmetric tridiagonal, U is orthogonal,
!> and S is diagonal (if KBAND=0) or symmetric tridiagonal (if KBAND=1)\&.
!> Two tests are performed:
!>
!>    RESULT(1) = | A - U S U' | / ( |A| n ulp )
!>
!>    RESULT(2) = | I - UU' | / ( n ulp )
!> 
.fi
.PP
 
.RE
.PP
\fBParameters\fP
.RS 4
\fIN\fP 
.PP
.nf
!>          N is INTEGER
!>          The size of the matrix\&.  If it is zero, DSTT21 does nothing\&.
!>          It must be at least zero\&.
!> 
.fi
.PP
.br
\fIKBAND\fP 
.PP
.nf
!>          KBAND is INTEGER
!>          The bandwidth of the matrix S\&.  It may only be zero or one\&.
!>          If zero, then S is diagonal, and SE is not referenced\&.  If
!>          one, then S is symmetric tri-diagonal\&.
!> 
.fi
.PP
.br
\fIAD\fP 
.PP
.nf
!>          AD is DOUBLE PRECISION array, dimension (N)
!>          The diagonal of the original (unfactored) matrix A\&.  A is
!>          assumed to be symmetric tridiagonal\&.
!> 
.fi
.PP
.br
\fIAE\fP 
.PP
.nf
!>          AE is DOUBLE PRECISION array, dimension (N-1)
!>          The off-diagonal of the original (unfactored) matrix A\&.  A
!>          is assumed to be symmetric tridiagonal\&.  AE(1) is the (1,2)
!>          and (2,1) element, AE(2) is the (2,3) and (3,2) element, etc\&.
!> 
.fi
.PP
.br
\fISD\fP 
.PP
.nf
!>          SD is DOUBLE PRECISION array, dimension (N)
!>          The diagonal of the (symmetric tri-) diagonal matrix S\&.
!> 
.fi
.PP
.br
\fISE\fP 
.PP
.nf
!>          SE is DOUBLE PRECISION array, dimension (N-1)
!>          The off-diagonal of the (symmetric tri-) diagonal matrix S\&.
!>          Not referenced if KBSND=0\&.  If KBAND=1, then AE(1) is the
!>          (1,2) and (2,1) element, SE(2) is the (2,3) and (3,2)
!>          element, etc\&.
!> 
.fi
.PP
.br
\fIU\fP 
.PP
.nf
!>          U is DOUBLE PRECISION array, dimension (LDU, N)
!>          The orthogonal matrix in the decomposition\&.
!> 
.fi
.PP
.br
\fILDU\fP 
.PP
.nf
!>          LDU is INTEGER
!>          The leading dimension of U\&.  LDU must be at least N\&.
!> 
.fi
.PP
.br
\fIWORK\fP 
.PP
.nf
!>          WORK is DOUBLE PRECISION array, dimension (N*(N+1))
!> 
.fi
.PP
.br
\fIRESULT\fP 
.PP
.nf
!>          RESULT is DOUBLE PRECISION array, dimension (2)
!>          The values computed by the two tests described above\&.  The
!>          values are currently limited to 1/ulp, to avoid overflow\&.
!>          RESULT(1) is always modified\&.
!> 
.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 \fB125\fP of file \fBdstt21\&.f\fP\&.
.SH "Author"
.PP 
Generated automatically by Doxygen for LAPACK from the source code\&.