2010
DOI: 10.48550/arxiv.1011.2775
|View full text |Cite
Preprint
|
Sign up to set email alerts
|

Multigrid solver for clover fermions

Abstract: We present an adaptive multigrid Dirac solver developed for Wilson clover fermions which offers order-of-magnitude reductions in solution time compared to conventional Krylov solvers. The solver incorporates even-odd preconditioning and mixed precision to solve the Dirac equation to double precision accuracy and shows only a mild increase in time to solution for decreasing quark mass. We show actual time to solution on production lattices in comparison to conventional Krylov solvers and will also discuss the s… Show more

Help me understand this report

Search citation statements

Order By: Relevance

Paper Sections

Select...
2
1
1
1

Citation Types

2
20
0

Year Published

2011
2011
2016
2016

Publication Types

Select...
4
2

Relationship

0
6

Authors

Journals

citations
Cited by 14 publications
(22 citation statements)
references
References 4 publications
2
20
0
Order By: Relevance
“…Second, we are leveraging a powerful new adaptive multigrid (MG) algorithm for inverting the Wilson-clover Dirac operator that is allowing us to compute the disconnected diagrams for both strange and light quarks at very little additional cost [64,65]. We are also taking advantage of clusters accelerated by graphics processing units (GPUs) using the QUDA library [66,67], which provides another substantial speedup.…”
Section: Discussionmentioning
confidence: 99%
“…Second, we are leveraging a powerful new adaptive multigrid (MG) algorithm for inverting the Wilson-clover Dirac operator that is allowing us to compute the disconnected diagrams for both strange and light quarks at very little additional cost [64,65]. We are also taking advantage of clusters accelerated by graphics processing units (GPUs) using the QUDA library [66,67], which provides another substantial speedup.…”
Section: Discussionmentioning
confidence: 99%
“…It is, therefore, efficient to precondition the matrix by deflating the low-eigenmodes. We implement such improvement using the multigrid solver [22,23] which has deflation built in. To obtain the LP estimate of the two-point function, we truncate the multigrid solver using a low-accuracy stopping criterion: the ratio (r LP ≡ |residue| LP /|source|) is chosen to be 10 −3 for all the ensembles.…”
Section: A Two-point Functionmentioning
confidence: 99%
“…These savings have been demonstrated in the original DD-HMC setup [14] as well as in mass preconditioned HMC [12]. A very similar idea is the adaptive multi-grid [15].…”
Section: Solvermentioning
confidence: 99%