Physics – High Energy Physics – High Energy Physics - Lattice
Scientific paper
2009-10-06
Physics
High Energy Physics
High Energy Physics - Lattice
14 pages, 4 figures; as published in Comput. Phys. Commun., modified data in Figs. 2,3 and 4 for improved implementation of FO
Scientific paper
10.1016/j.cpc.2010.04.004
The overlap operator in lattice QCD requires the computation of the sign function of a matrix, which is non-Hermitian in the presence of a quark chemical potential. In previous work we introduced an Arnoldi-based Krylov subspace approximation, which uses long recurrences. Even after the deflation of critical eigenvalues, the low efficiency of the method restricts its application to small lattices. Here we propose new short-recurrence methods which strongly enhance the efficiency of the computational method. Using rational approximations to the sign function we introduce two variants, based on the restarted Arnoldi process and on the two-sided Lanczos method, respectively, which become very efficient when combined with multishift solvers. Alternatively, in the variant based on the two-sided Lanczos method the sign function can be evaluated directly. We present numerical results which compare the efficiencies of a restarted Arnoldi-based method and the direct two-sided Lanczos approximation for various lattice sizes. We also show that our new methods gain substantially when combined with deflation.
Bloch Jacques C. R.
Breu Tobias
Frommer Andreas
Heybrock Simon
Schäfer Katrin
No associations
LandOfFree
Short-recurrence Krylov subspace methods for the overlap Dirac operator at nonzero chemical potential 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 Short-recurrence Krylov subspace methods for the overlap Dirac operator at nonzero chemical potential, we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and Short-recurrence Krylov subspace methods for the overlap Dirac operator at nonzero chemical potential will most certainly appreciate the feedback.
Profile ID: LFWR-SCP-O-358078