Actions of noncompact groups and algorithm design: a case study
Numerical matrix computations involving actions of noncompact transformation groups are known to produce numerical problems since the elements of the pertaining matrix representations are inherently unbounded. In this case study we analyse numerical problems occurring in a class of algorithms that is based on actions of the pseudo-orthogonal group O/sub n,m/-a group that is noncompact (hyperbolic geometry) and well established in signal processing (Schur methods). As a major result, it is shown how to exploit the additional degrees of freedom in defining coordinate frames in a Grassmannian setting in order to impose an a priori bound on the norm of the transformation matrices. This way, numerically disastrous situations can be circumvented systematically. Hence, it becomes possible to develop modified algorithms which exhibit superior numerical performance for a large class of problems based on, for example, hyperbolic transformations.