Chronological inversion method for the Dirac matrix in hybrid Monte Carlo
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.
 Publication:

Nuclear Physics B
 Pub Date:
 February 1997
 DOI:
 10.1016/S05503213(96)005792
 arXiv:
 arXiv:heplat/9509012
 Bibcode:
 1997NuPhB.484..353B
 Keywords:

 High Energy Physics  Lattice
 EPrint:
 35 pages, 18 EPS figures A new "preconditioning" method, derived from the Chronological Inversion, is described. Some new figures are appended. Some reorganization of the material has taken place