.TH "TESTING/EIG/chst01.f" 3 "Version 3.12.0" "LAPACK" \" -*- nroff -*- .ad l .nh .SH NAME TESTING/EIG/chst01.f .SH SYNOPSIS .br .PP .SS "Functions/Subroutines" .in +1c .ti -1c .RI "subroutine \fBchst01\fP (n, ilo, ihi, a, lda, h, ldh, q, ldq, work, lwork, rwork, result)" .br .RI "\fBCHST01\fP " .in -1c .SH "Function/Subroutine Documentation" .PP .SS "subroutine chst01 (integer n, integer ilo, integer ihi, complex, dimension( lda, * ) a, integer lda, complex, dimension( ldh, * ) h, integer ldh, complex, dimension( ldq, * ) q, integer ldq, complex, dimension( lwork ) work, integer lwork, real, dimension( * ) rwork, real, dimension( 2 ) result)" .PP \fBCHST01\fP .PP \fBPurpose:\fP .RS 4 .PP .nf CHST01 tests the reduction of a general matrix A to upper Hessenberg form: A = Q*H*Q'\&. Two test ratios are computed; RESULT(1) = norm( A - Q*H*Q' ) / ( norm(A) * N * EPS ) RESULT(2) = norm( I - Q'*Q ) / ( N * EPS ) The matrix Q is assumed to be given explicitly as it would be following CGEHRD + CUNGHR\&. In this version, ILO and IHI are not used, but they could be used to save some work if this is desired\&. .fi .PP .RE .PP \fBParameters\fP .RS 4 \fIN\fP .PP .nf N is INTEGER The order of the matrix A\&. N >= 0\&. .fi .PP .br \fIILO\fP .PP .nf ILO is INTEGER .fi .PP .br \fIIHI\fP .PP .nf IHI is INTEGER A is assumed to be upper triangular in rows and columns 1:ILO-1 and IHI+1:N, so Q differs from the identity only in rows and columns ILO+1:IHI\&. .fi .PP .br \fIA\fP .PP .nf A is COMPLEX array, dimension (LDA,N) The original n by n matrix A\&. .fi .PP .br \fILDA\fP .PP .nf LDA is INTEGER The leading dimension of the array A\&. LDA >= max(1,N)\&. .fi .PP .br \fIH\fP .PP .nf H is COMPLEX array, dimension (LDH,N) The upper Hessenberg matrix H from the reduction A = Q*H*Q' as computed by CGEHRD\&. H is assumed to be zero below the first subdiagonal\&. .fi .PP .br \fILDH\fP .PP .nf LDH is INTEGER The leading dimension of the array H\&. LDH >= max(1,N)\&. .fi .PP .br \fIQ\fP .PP .nf Q is COMPLEX array, dimension (LDQ,N) The orthogonal matrix Q from the reduction A = Q*H*Q' as computed by CGEHRD + CUNGHR\&. .fi .PP .br \fILDQ\fP .PP .nf LDQ is INTEGER The leading dimension of the array Q\&. LDQ >= max(1,N)\&. .fi .PP .br \fIWORK\fP .PP .nf WORK is COMPLEX array, dimension (LWORK) .fi .PP .br \fILWORK\fP .PP .nf LWORK is INTEGER The length of the array WORK\&. LWORK >= 2*N*N\&. .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) RESULT(1) = norm( A - Q*H*Q' ) / ( norm(A) * N * EPS ) RESULT(2) = norm( I - Q'*Q ) / ( N * EPS ) .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 \fB138\fP of file \fBchst01\&.f\fP\&. .SH "Author" .PP Generated automatically by Doxygen for LAPACK from the source code\&.