VPAStab(J,LJ,L): An iterative method with look-ahead for the solution of large sparse linear systems

VPAStab(J,LJ,L) denotes an iterative method designed to give an accurate numerical solution of a large sparse system of linear equations. The iteration comprises LL-steps and possibly jump-steps. An LL-step involves stabilisation (minimisation) over an LL-dimensional subspace following LL levels of vector-Padé approximation. A jump-step (look-ahead stage) avoids near-breakdowns by jumping past them by JJ levels simultaneously. A numerical implementation of VPAStab(2,L2,L) is considered in the context of some well-known examples.