A Parallel SSOR Preconditioner for Lattice QCD
arXiv:hep-lat/9602019 · doi:10.1016/0010-4655(96)00089-6
Abstract
We present a parallelizable SSOR preconditioning scheme for Krylov subspace iterative solvers which proves to be efficient in lattice QCD applications involving Wilson fermions. Our preconditioner is based on a locally lexicographic ordering of the lattice points. In actual hybrid Monte Carlo applications with the bi-conjugate gradient stabilized method BiCGstab, we achieve a gain factor of about 2 in the number of iterations compared to conventional odd-even preconditioning. Whether this translates into similar reductions in run time will depend on the parallel computer in use. We discuss implementation issues using the `Eisenstat-trick' and machine specific advantages of the method for the APE100/Quadrics parallel computer. In a full QCD simulation with Wilson fermions on a 512-processor Quadrics QH4 we find a gain in cpu-time of a factor of 1.7 over odd-even preconditioning for a 24^3 x 40 lattice.
21 pages, Tex-file + Postscript figures
References in corpus (4)
- A New Approach to the Problem of Dynamical Quarks in Numerical Simulations of Lattice QCD
- Accelerating Wilson Fermion Matrix Inversions by Means of the Stabilized Biconjugate Gradient Algorithm
- Chronological Inversion Method for the Dirac Matrix in Hybrid Monte Carlo
- Fast Krylov Space Methods for Calculation of Quark Propogator
Cited by in corpus (36)
- 2+1 Flavor Lattice QCD toward the Physical Point
- Schwarz-preconditioned HMC algorithm for two-flavour lattice QCD
- HMC algorithm with multiple time scale integration and mass preconditioning
- Numerical Methods for the QCD Overlap Operator: I. Sign-Function and Error Bounds
- Computing and deflating eigenvalues while solving multiple right hand side linear systems in Quantum Chromodynamics
- Light and Strange Hadron Spectroscopy with Dynamical Wilson Fermions
- Hadron masses and matrix elements from the QCD Schr"odinger functional
- Unqueched Kogut-Susskind quark propagator in Lattice Landau Gauge QCD
- Non-perturbative renormalization of static-light four-fermion operators in quenched lattice QCD
- Finite-Size Effects in Lattice QCD with Dynamical Wilson Fermions
- Non-perturbatively renormalised light quark masses from a lattice simulation with N_f=2
- SESAM and TXL Results for Wilson Action--A Status Report
- Non-perturbative renormalization of the static axial current in quenched QCD
- Comparative Benchmarks of full QCD Algorithms
- Continuous external momenta in non-perturbative lattice simulations: a computation of renormalization factors
- The PHMC algorithm for simulations of dynamical fermions: II - Performance analysis
- Progress in lattice algorithms
- Non-perturbative running of the average momentum of non-singlet parton densities
- Preconditioning of Improved and ``Perfect'' Fermion Actions
- Modified Block BiCGSTAB for Lattice QCD
- Efficiency of the UV-filtered Multiboson algorithm
- The LHMC Algorithm for Free Field Theory: Reexamining Overrelaxation
- Computing the eta and eta' Mesons in Lattice QCD
- Critical Dynamics of the Hybrid Monte Carlo Algorithm
- Full QCD Algorithms towards the Chiral Limit
- Iterative methods for overlap and twisted mass fermions
- A Parallel SSOR Preconditioner for Lattice QCD
- Optimised Dirac Operators on the Lattice: Construction, Properties and Applications
- Neural-network preconditioners for solving the Dirac equation in lattice gauge theory
- SSOR Preconditioning of Improved Actions
- Cost of QCD simulations with n_f=2 dynamical Wilson fermions
- Recent development and perspectives of machines for lattice QCD
- A Preconditioner for Improved Fermion Actions
- The Hybrid Monte Carlo Algorithm for Quantum Chromodynamics
- Block Algorithms for Quark Propagator Calculation
- QCD on α-Clusters