Deflation for inversion with multiple right-hand sides in QCD
Most calculations in lattice Quantum Chromodynamics (QCD) involve the solution of a series of linear systems of equations with exceedingly large matrices and a large number of right hand sides. Iterative methods for these problems can be sped up significantly if we deflate approximations of appropriate invariant spaces from the initial guesses. Recently we have developed eigCG, a modification of the Conjugate Gradient (CG) method, which while solving a linear system can reuse a window of the CG vectors to compute eigenvectors almost as accurately as the Lanczos method. The number of approximate eigenvectors can increase as more systems are solved. In this paper we review some of the characteristics of eigCG and show how it helps remove the critical slowdown in QCD calculations. Moreover, we study scaling with lattice volume and an extension of the technique to nonsymmetric problems.
- Research Organization:
- Thomas Jefferson National Accelerator Facility, Newport News, VA (United States)
- Sponsoring Organization:
- USDOE Office of Science (SC)
- DOE Contract Number:
- AC05-06OR23177
- OSTI ID:
- 1003867
- Report Number(s):
- JLAB-THY-09-977; DOE/OR/23177-1412
- Journal Information:
- J. Phys., Conf. Ser., Journal Name: J. Phys., Conf. Ser. Vol. 180
- Country of Publication:
- United States
- Language:
- English
Similar Records
Extending the eigCG algorithm to non-symmetric linear systems with multiple right-hand sides
Multigrid deflation for Lattice QCD