Augmented unprojected Krylov subspace methods from an alternative view of an existing framework
share
Augmented Krylov subspace methods aid in accelerating the convergence of a standard Krylov subspace method by including additional vectors in the search space. These methods are commonly used in High-Performance Computing (HPC) applications to considerably reduce the number of matrix vector products required when building a basis for the Krylov subspace. In a recent survey on subspace recycling iterative methods [Soodhalter et al, GAMM-Mitt. 2020], a framework was presented which describes a wide class of such augmented methods. The framework describes these methods broadly in a two step process. Step one involves solving a projected problem via a standard Krylov subspace, or other projection method, and step two then performs an additional projection into the augmentation subspace. In this work we show that the projected problem one must solve in step one has an equivalent unprojected formulation. We then show how this observation allows the framework to be adapted to describe unprojected augmented Krylov subspace methods. We demonstrate this with two examples. We first show one can recover the R^3GMRES algorithm from this view of the framework, and then we use the framework to derive the first unprojected augmented Full Orthogonalization Method (FOM). This method allows one to recycle information from previous system solves, and we thus denote it as unproj rFOM (unprojected recycled FOM).
READ FULL TEXT