function reduces a real, nonsymmetric
to upper Hessenberg form. The result is an upper Hessenberg array with
eigenvalues that are identical to those of the original array
Hessenberg array is stored in elements (
+ 1. Elements with
+ 1 are to be thought of as zero, but are returned with random values.
ELMHES is based on the routine
described in section 11.5 of
Numerical Recipes in C: The Art of Scientific Computing
(Second Edition), published by Cambridge University Press, and is used by permission.
Result = ELMHES(
real, nonsymmetric array.
Set this keyword to force the computation to be done in double-precision arithmetic.
Set this keyword if the input array
is in column-major format (composed of column vectors) rather than in row-major format (composed of row vectors).
Set this keyword to disable balancing. By default, a balancing algorithm is applied to
. Balancing a nonsymmetric array is recommended to reduce the sensitivity of eigenvalues to rounding errors.
See the description of HQR for an example using this function.