Numerical Methods for the QCD Overlap Operator:III. Nested Iterations

Details Numerical Methods for the QCD Overlap Operator:III. Nested Iterations Numerical Methods for the QCD Overlap Operator:III. Nested Iterations

2004-05-05

arxiv.org/abs/hep-lat/0405003v3

Comput.Phys.Commun. 165 (2005) 221-242

Physics

High Energy Physics

High Energy Physics - Lattice

33 Pages; citations added

Scientific paper

10.1016/j.cpc.2004.10.005

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.

Affiliated with

Also associated with

No associations

Rating

Numerical Methods for the QCD Overlap Operator:III. Nested Iterations does not yet have a rating. At this time, there are no reviews or comments for this scientific paper.

If you have personal experience with Numerical Methods for the QCD Overlap Operator:III. Nested Iterations, we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and Numerical Methods for the QCD Overlap Operator:III. Nested Iterations will most certainly appreciate the feedback.

Rate now

Comments { 0 }

Profile ID: LFWR-SCP-O-467176