1997
DOI: 10.1016/s0550-3213(96)00579-2
|View full text |Cite
|
Sign up to set email alerts
|

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… Show more

Help me understand this report
View preprint versions

Search citation statements

Order By: Relevance

Paper Sections

Select...
2
2
1

Citation Types

2
94
0

Year Published

2000
2000
2014
2014

Publication Types

Select...
5
3
1

Relationship

1
8

Authors

Journals

citations
Cited by 82 publications
(96 citation statements)
references
References 16 publications
2
94
0
Order By: Relevance
“…While the intermediate mass term introduces additional pseudofermion fields, it helps to reduce the overall fermion force in the MD steps, making it possible to use larger step sizes while maintaining a good acceptance rate [39]. Another reduction in computational cost comes from the use of the chronological inverter [40]. To avoid the loss of reversibility, we set stringent stopping conditions for the conjugate gradient matrix inversion: 10 −9 for the MD evolution, and 10 −10 for the Metropolis step.…”
Section: Simulation Details a Parameters And Algorithmsmentioning
confidence: 99%
“…While the intermediate mass term introduces additional pseudofermion fields, it helps to reduce the overall fermion force in the MD steps, making it possible to use larger step sizes while maintaining a good acceptance rate [39]. Another reduction in computational cost comes from the use of the chronological inverter [40]. To avoid the loss of reversibility, we set stringent stopping conditions for the conjugate gradient matrix inversion: 10 −9 for the MD evolution, and 10 −10 for the Metropolis step.…”
Section: Simulation Details a Parameters And Algorithmsmentioning
confidence: 99%
“…We also give the number, N csg i , of solutions of the Dirac equation we save for the chronological solver guess [35] with the purpose of evaluating the two force terms (i = 1, 2) associated to pseudofermion integration (i = 0 refers to the pure gauge force). The notation N csg 1,2 = 0 means that no chronological solver guess was used there.…”
Section: Simulation Algorithm and Error Analysismentioning
confidence: 99%
“…[35]. In our GRAL simulations the depth of the ex- The parameters of our simulations within the GRAL project together with those of previous SESAM and TχL runs which we have compiled here for reference (see also Table II).…”
Section: Simulationmentioning
confidence: 99%