Chronological inversion method for the Dirac matrix in hybrid Monte Carlo
Year: 1997
Type: article
Abstract: In hybrid Monte Carlo simulations for full QCD, the gauge fields evolve smoothly as a function of molecular dynamics time. Here we investigate improved methods of estimating the trial or starting solutions for the Dirac matrix inversion as superpositions of a chronological sequence of solutions in the recent past. By taking as the trial solution the vector which minimizes the residual in the linear space spanned by the past solutions, the number of conjugate gradient iterations per unit MD time is decreased by at least a factor of 2. Extensions of this basic approach to precondition the conjugate gradient iterations are also discussed.
Source: Nuclear Physics B
Authors: R.C. Brower, T. Ivanenko, A.R. Levi, K.N. Orginos
Language: English
Subfield: Nuclear and High Energy Physics
Field: Physics and Astronomy
Domain: Physical Sciences
Open Access status: hybrid
Estimated APC paid:
OpenAlex ID: W2003143357
