Related papers: Comparison of Domain Wall Fermion Multigrid Method…
We construct the Schr\"odinger Functional (SF) setup for the M\"obius domain wall fermions (MDWF). The method is an extension of the method proposed by Takeda for the standard domain wall fermion. In order to fulfill the requirement that…
Stencil computations are widely used to simulate the change of state of physical systems across a multidimensional grid over multiple timesteps. The state-of-the-art techniques in this area fall into three groups: cache-aware tiled looping…
We show how the standard domain wall action can be simply modified to allow arbitrarily exact chiral symmetry at finite fifth dimensional extent $L_s$. We note that the method can be used for both quenched and dynamical calculations. We…
We describe our implementation of a multigrid solver for Wilson-clover fermions, which increases parallelism by solving for multiple right-hand sides (MRHS) simultaneously. The solver is based on Grid and thus runs on all computing…
We present the design and implementation details of a geometric multigrid method on adaptively refined meshes for massively parallel computations. The method uses local smoothing on the refined part of the mesh. Partitioning is achieved by…
We present a new multigrid scheme for solving the Poisson equation with Dirichlet boundary conditions on a Cartesian grid with irregular domain boundaries. This scheme was developed in the context of the Adaptive Mesh Refinement (AMR)…
Domain-wall fermions preserve chiral symmetry up to terms that decrease exponentially when the lattice size in the fifth dimension is taken to infinity. The associated rates of convergence are given by the low-lying eigenvalues of a simple…
Multigrid solvers are the standard in modern scientific computing simulations. Domain Decomposition Aggregation-Based Algebraic Multigrid, also known as the DD-$\alpha$AMG solver, is a successful realization of an algebraic multigrid solver…
We present a set of related Hybrid Monte Carlo methods to simulate an arbitrary number of dynamical overlap fermions. Each fermion is represented by a chiral pseudo-fermion field. The new algorithm reduces critical slowing down in the…
Methods for solving hyperbolic systems typically depend on unknown ordering (e.g., Gauss-Seidel, or sweep/wavefront/marching methods) to achieve good convergence. For many discretisations, mesh types or decompositions these methods do not…
Multigrid solvers for hierarchical hybrid grids (HHG) have been proposed to promote the efficient utilization of high performance computer architectures. These HHG meshes are constructed by uniformly refining a relatively coarse fully…
This paper is concerned with developing an efficient numerical algorithm for fast implementation of the sparse grid method for computing the $d$-dimensional integral of a given function. The new algorithm, called the MDI-SG ({\em multilevel…
We study a parameter optimization of domain-wall fermions to improve chiral symmetry based on machine learning. Domain-wall fermions involve coefficients along the fifth dimension, which can be treated as trainable parameters to reduce the…
We investigate the algorithms for dynamical overlap fermions aiming at improving the performance for large-scale simulations. We look for the best combination of Hybrid Monte Carlo options and iterative quark solvers with respect to the…
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…
This paper presents four novel domain decomposition algorithms integrated with nonlinear mapping techniques to address collocation-based solutions of eigenvalue problems involving sharp interfaces or steep gradients. The proposed methods…
We derive the exactly conserved vector, and almost conserved axial currents for rational approximations to the overlap operator with a general Mobius kernel. The approach maintains manifest Hermiticity, and allows matrix elements of the…
There has been much recent progress in the understanding and reduction of the computational cost of the Hybrid Monte Carlo algorithm for Lattice QCD as the quark mass parameter is reduced. In this letter we present a new solution to this…
We present a polynomial multigrid method for the nodal interior penalty formulation of the Poisson equation on three-dimensional Cartesian grids. Its key ingredient is a weighted overlapping Schwarz smoother operating on element-centered…
The inverse of the fermion matrix squared is used to define a transfer matrix for domain-wall fermions. When the domain-wall height $M$ is bigger than one, the transfer matrix is complex. Slowly suppressed chiral symmetry violations may…