Seed methods for linear equations in lattice qcd problems with multiple right-hand sides
High Energy Physics - Lattice
2009-03-19 v2
Abstract
We consider three improvements to seed methods for Hermitian linear systems with multiple right-hand sides: only the Krylov subspace for the first system is used for seeding subsequent right-hand sides, the first right-hand side is solved past convergence, and periodic re-orthogonalization is used in order to control roundoff errors associated with the Conjugate Gradient algorithm. The method is tested for the case of Wilson fermions near kappa critical and a considerable speed up in the convergence is observed.
Cite
@article{arxiv.0901.3512,
title = {Seed methods for linear equations in lattice qcd problems with multiple right-hand sides},
author = {Abdou Abdel-Rehim and Ronald B. Morgan and Walter Wilcox},
journal= {arXiv preprint arXiv:0901.3512},
year = {2009}
}
Comments
7 pages, 2 figures. Presented at the 26th International Symposium on Lattice Field Theory, Williamsburg, VA, USA, 14-19 Jul 2008. Fixed a hyperlink to one of the references