.TH "TESTING/EIG/zhst01.f" 3 "Version 3.12.0" "LAPACK" \" -*- nroff -*- .ad l .nh .SH NAME TESTING/EIG/zhst01.f .SH SYNOPSIS .br .PP .SS "Functions/Subroutines" .in +1c .ti -1c .RI "subroutine \fBzhst01\fP (n, ilo, ihi, a, lda, h, ldh, q, ldq, work, lwork, rwork, result)" .br .RI "\fBZHST01\fP " .in -1c .SH "Function/Subroutine Documentation" .PP .SS "subroutine zhst01 (integer n, integer ilo, integer ihi, complex*16, dimension( lda, * ) a, integer lda, complex*16, dimension( ldh, * ) h, integer ldh, complex*16, dimension( ldq, * ) q, integer ldq, complex*16, dimension( lwork ) work, integer lwork, double precision, dimension( * ) rwork, double precision, dimension( 2 ) result)" .PP \fBZHST01\fP .PP \fBPurpose:\fP .RS 4 .PP .nf !> !> ZHST01 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 ZGEHRD + ZUNGHR\&. !> !> 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*16 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*16 array, dimension (LDH,N) !> The upper Hessenberg matrix H from the reduction A = Q*H*Q' !> as computed by ZGEHRD\&. 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*16 array, dimension (LDQ,N) !> The orthogonal matrix Q from the reduction A = Q*H*Q' as !> computed by ZGEHRD + ZUNGHR\&. !> .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*16 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 DOUBLE PRECISION array, dimension (N) !> .fi .PP .br \fIRESULT\fP .PP .nf !> RESULT is DOUBLE PRECISION 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 \fBzhst01\&.f\fP\&. .SH "Author" .PP Generated automatically by Doxygen for LAPACK from the source code\&.