DHST01 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 DGEHRD + DORGHR.
 In this version, ILO and IHI are not used and are assumed to be 1 and
 N, respectively.

N

          N is INTEGER
          The order of the matrix A.  N >= 0.

ILO

          ILO is INTEGER

IHI

          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.

A

          A is DOUBLE PRECISION array, dimension (LDA,N)
          The original n by n matrix A.

LDA

          LDA is INTEGER
          The leading dimension of the array A.  LDA >= max(1,N).

H

          H is DOUBLE PRECISION array, dimension (LDH,N)
          The upper Hessenberg matrix H from the reduction A = Q*H*Q'
          as computed by DGEHRD.  H is assumed to be zero below the
          first subdiagonal.

LDH

          LDH is INTEGER
          The leading dimension of the array H.  LDH >= max(1,N).

Q

          Q is DOUBLE PRECISION array, dimension (LDQ,N)
          The orthogonal matrix Q from the reduction A = Q*H*Q' as
          computed by DGEHRD + DORGHR.

LDQ

          LDQ is INTEGER
          The leading dimension of the array Q.  LDQ >= max(1,N).

WORK

          WORK is DOUBLE PRECISION array, dimension (LWORK)

LWORK

          LWORK is INTEGER
          The length of the array WORK.  LWORK >= 2*N*N.

RESULT

          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 )

Univ. of Tennessee

Univ. of California Berkeley

Univ. of Colorado Denver

NAG Ltd.