Implementing an Algorithm for Solving Block Hessenberg Systems

This paper describes the implementation of a recursive descent methodfor solving block Hessenberg systems. Although the algorithm is
conceptually simple, its implementation in C (a natural choice of
language given the recursive nature of the algorithm and its data) is
nontrivial. Particularly important is the balance between ease of
use, computational efficiency, and flexibility.
(Also cross-referenced as UMIACS-TR-94-70)