Related papers: Adaptive Multigrid Algorithm for Lattice QCD
We have developed an adaptive multigrid code for solving the Poisson equation in gravitational simulations. Finer rectangular subgrids are adaptively created in locations where the density exceeds a local level-dependent threshold. We…
We perform a comparison between different lattice regularizations of the Dirac operator for massless fermions in the framework of the single and two flavor Schwinger model. We consider a) the Wilson-Dirac operator at the critical value of…
The computational complexity of naive, sampling-based uncertainty quantification for 3D partial differential equations is extremely high. Multilevel approaches, such as multilevel Monte Carlo (MLMC), can reduce the complexity significantly,…
At physical light quark masses, efficient linear solvers are crucial for carrying out the millions of inversions of the Dirac matrix required for obtaining high statistics in quark correlation functions. Adaptive algebraic multi-grid…
A new class of lattice Dirac operators which satisfy the index theorem have been recently proposed on the basis of the algebraic relation $\gamma_{5}(\gamma_{5}D) + (\gamma_{5}D)\gamma_{5} = 2a^{2k+1}(\gamma_{5}D)^{2k+2}$. Here $k$ stands…
The gauge covariant lattice Dirac operator D which has recently been proposed by Neuberger satisfies the Ginsparg-Wilson relation and thus preserves chiral symmetry. The operator also avoids a doubling of fermion species, but its locality…
Overlap fermions are particularly well suited to study the finite temperature dynamics of the chiral symmetry restoration transition of QCD, which might be just an analytic crossover. Using gauge field configurations on a 24^3x10 lattice…
The first order condition of the constrained minimization problem leads to a saddle point problem. A multigrid method using a multiplicative Schwarz smoother for saddle point problems can thus be interpreted as a successive subspace…
We establish global existence and derive sharp pointwise decay estimates of solutions to cubic Dirac and Dirac-Klein-Gordon systems on a curved background, close to the Minkowski spacetime. By squaring the Dirac operator, we reduce the…
Computing disconnected diagrams in Lattice QCD (operator insertion in a quark loop) entails the computationally demanding problem of taking the trace of the all to all quark propagator. We first outline the basic algorithm used to compute a…
We investigate exact symmetries of a staggered fermion in D dimensions. The Dirac operator is reformulated by SO(2D) Clifford algebra. The chiral symmetry, rotational invariance and parity symmetries are clarified in any dimension. Local…
I review recent theoretical developments and numerical results of twisted mass QCD. I argue that, combined with an efficient algorithm, twisted mass QCD can be an attractive QCD lattice action, to perform large scale simulations at small…
In master-field simulations of lattice QCD, the expectation values of interest are obtained from a single or at most a few representative gauge-field configurations on very large lattices. If the light quarks are included, the generation of…
This work proposes a scheme for significantly reducing the computational complexity of discretized problems involving the non-smooth forward propagation of uncertainty by combining the adaptive hierarchical sparse grid stochastic…
An idealized multigrid algorithm for the computation of propagators of staggered fermions is investigated. Exemplified in four-dimensional $SU(2)$ gauge fields, it is shown that the idealized algorithm preserves criticality under…
Adaptive gradient methods, such as AdaGrad, are among the most successful optimization algorithms for neural network training. While these methods are known to achieve better dimensional dependence than stochastic gradient descent (SGD) for…
We propose the first optimal geometric multigrid solver for hybrid high-order discretizations that can handle arbitrary polytopal agglomeration hierarchies in both two and three dimensions. The key ingredient is the use of modified skeleton…
It is shown that the 21264 Alpha processor can reach about 20% sustained efficiency for the inversion of the Wilson-Dirac operator. Since fast ethernet is not sufficient to get balancing between computation and communication on reasonable…
It is well known that the block Krylov subspace solvers work efficiently for some cases of the solution of differential equations with multiple right-hand sides. In lattice QCD calculation of physical quantities on a given configuration…
I review the lattice formulations of vector-like gauge theories (e.g. QCD) with domain-wall/overlap fermions, and discuss how to optimize the chiral symmetry for any finite $ N_s $ (sites in the fifth dimension). In this formulation, quark…