Householder's method to reduce a real, symmetric array to
TRIRED is based on the routine
described in section 11.2 of
Numerical Recipes in C: The Art of Scientific Computing
(Second Edition), published by Cambridge University Press, and is used by permission.
A, D, E
real, symmetric array that is replaced, on exit, by the orthogonal array
effecting the transformation. The routine TRIQL can use this result to find the eigenvectors of the array
-element output vector containing the diagonal elements of the tridiagonal array.
-element output vector containing the off-diagonal elements.
Set this keyword to force the computation to be done in double-precision arithmetic.
See the description of
for an example using this function.