Fast and stable solution of banded-plus-semiseparable linear systems

摘要

We describe a stable algorithm, having linear complexity, for the solution of banded-plus-semiseparable linear systems. The algorithm exploits the structural properties of the inverse of a semiseparable matrix. Stability is achieved by combining these properties with partial pivoting techniques. Several numerical experiments are shown to confirm the effectiveness of the proposed approach.