guest ::
| Home > Publications database > Numerical Methods for the QCD Overlap Operator. III. Nested Iterations |
| Home > Publications database > Numerical Methods for the QCD Overlap Operator. III. Nested Iterations |
| Journal Article | PreJuSER-50021 |
; ; ; ; ;
2005
North Holland Publ. Co.
Amsterdam
This record in other databases: 
Please use a persistent id in citations: doi:10.1016/j.cpc.2004.10.005
Abstract: The numerical and computational aspects of chiral fermions in lattice quantum chromodynamics are extremely demanding. In the overlap framework, the computation of the fermion propagator leads to a nested iteration where the matrix vector multiplications in each step of an outer iteration have to be accomplished by an inner iteration; the latter approximates the product of the sign function of the hermitian Wilson fermion matrix with a vector.In this paper we investigate aspects of this nested paradigm. We examine several Krylov subspace methods to be used as an outer iteration for both propagator computations and the Hybrid Monte-Carlo scheme. We establish criteria on the accuracy of the inner iteration which allow to preserve an a priori given precision for the overall computation. It will turn out that the accuracy of the sign function can be relaxed as the outer iteration proceeds. Furthermore, we consider preconditioning strategies, where the preconditioner is built upon an inaccurate approximation to the sign function. Relaxation combined with preconditioning allows for considerable savings in computational efforts up to a factor of 4 as our numerical experiments illustrate. We also discuss the possibility of projecting the squared overlap operator into one chiral sector. (C) 2004 Elsevier B.V. All rights reserved.
Keyword(s): J ; lattice quantum chromodynamics (auto) ; overlap fermions (auto) ; matrix sign function (auto) ; inner-outer iterations (auto) ; relaxation (auto) ; flexible Krylov (auto) ; subspace methods (auto)
|
The record appears in these collections: |