Software Open Access

Shift Block Recycled GMRES: A Sylvester equation based minimum residual Krylov subspace method for solving multiple shifted systems

Soodhalter, Kirk M.

We propose a block Krylov subspace method for solving a family of shifted linear systems based on the GCRO-DR method proposed in [Parks et al. SISC 2005], which is an iterative method allowing for the efficient minimization of the the residual over an augmented block Krylov subspace. It has been previously shown that using augmented methods to solve multiple shifted systems simultaneously can been problematic for non-symmetric problems.  Here, this difficulty is overcome by using the Sylvester equation formulation of a family of shifted systems and applying a block Krylov subspace method to solve this Sylvester equation (a strategy proposed by Simoncini in 1996).  This strategy is fully compatible with augmentation and, specifically, with a GCRO-type minimization method.

Associated to the paper: Block Krylov Subspace Recycling for Shifted Systems with Unrelated Right-Hand Sides (Kirk M. Soodhalter), In SIAM Journal on Scientific Computing, volume 38, 2016.
Files (769.2 kB)
Name Size
conf5_0-4x4-10.mat
md5:a9e3659abe093ed34bebd79e7b899fba
717.6 kB Download
rgmres_bl_sh.m
md5:bae0cb7eadd78ed2ee4b772ef5dfe78d
42.1 kB Download
shifted_block_gcrodr_examples.m
md5:74d8137a230fbaf4fc746fc6fe92fbe9
9.5 kB Download
71
5
views
downloads
All versions This version
Views 7171
Downloads 55
Data volume 853.4 kB853.4 kB
Unique views 7070
Unique downloads 22

Share

Cite as