Related papers: Deflation and Flexible SAP-Preconditioning of GMRE…
We discuss the usage and applicability of deflation methods for the overlap lattice Dirac operator, focussing on calculating the eigenvalues using a method similar to the eigCG algorithm used for other Dirac operators. The overlap operator,…
In lattice QCD computations a substantial amount of work is spent in solving discretized versions of the Dirac equation. Conventional Krylov solvers show critical slowing down for large system sizes and physically interesting parameter…
Numerical simulations of quantum chromodynamics (QCD) on a lattice require the frequent solution of linear systems of equations with large, sparse and typically ill-conditioned matrices. Algebraic multigrid methods are meanwhile the…
Computing the trace of the inverse of large matrices is typically addressed through statistical methods. Deflating out the lowest eigenvectors or singular vectors of the matrix reduces the variance of the trace estimator. This work…
Lattice QCD simulations at small lattice spacings and quark masses close to their physical values are technically challenging. In particular, the simulations can get trapped in the topological charge sectors of field space or may run into…
Recent developments have shown that a lot can be gained for QCD simulations from GPU hardware. This can be exploited especially in the case of Ginsparg-Wilson fermions when the com putational costs are particularly high. In this work, we…
We propose improved estimators to compute the reweighting factors which are needed for lattice QCD calculations that rely on twisted-mass reweighting for the light quark contribution and the Rational Hybrid Monte Carlo (RHMC) algorithm for…
In lattice QCD computations a substantial amount of work is spent in solving linear systems arising in Wilson's discretization of the Dirac equations. We show first numerical results of the extension of the two-level DD-\alpha AMG method to…
Versions of GMRES with deflation of eigenvalues are applied to lattice QCD problems. Approximate eigenvectors corresponding to the smallest eigenvalues are generated at the same time that linear equations are solved. The eigenvectors…
In order to develop fast inversion algorithms we have used overlap solvers in two dimensions. Lattice QED theory with U(1) group symmetry in two dimensional space-times dimensions has always been a testing ground for algorithms. By the…
The overlap Dirac operator, which satisfies the Ginsparg-Wilson relation, realizes exact chiral symmetry on the lattice without any unphysical doubler modes. To perform the path integrals, one should, however, note that the overlap fermion…
We compute the low lying eigenvalues of the Hermitian Dirac operator in lattice QCD with $N_{\rm f} = 2+1+1$ twisted mass fermions. We discuss whether these eigenvalues are in the $\epsilon$-regime or the $p$-regime of Wilson chiral…
The microscopic spectral density of the Wilson Dirac operator for two flavor lattice QCD is analyzed. The computation includes the leading order a^2 corrections of the chiral Lagrangian in the microscopic limit. The result is used to…
We present a new multigrid solver that is suitable for the Dirac operator in the presence of disordered gauge fields. The key behind the success of the algorithm is an adaptive projection onto the coarse grids that preserves the near null…
Generalized inverses play a fundamental role in numerical linear algebra, particularly when matrices are rectangular, singular, or rank deficient. Even when the input matrix is sparse, generalized inverses such as the M-P pseudoinverse are…
We present numerical results for the 2-flavour Schwinger model with dynamical chiral lattice fermions. We insert an approximately chiral hypercube Dirac operator into the overlap formula to construct the overlap hypercube operator. This is…
Recently various new concepts for the construction of Dirac operators in lattice Quantum Chromodynamics (QCD) have been introduced. These operators satisfy the so-called Ginsparg-Wilson condition (GWC), thus obeying the Atiyah-Singer index…
Individual agents in a multi-agent system (MAS) may have decoupled open-loop dynamics, but a cooperative control objective usually results in coupled closed-loop dynamics thereby making the control design computationally expensive. The…
Solving discretized versions of the Dirac equation represents a large share of execution time in lattice Quantum Chromodynamics (QCD) simulations. Many high-performance computing (HPC) clusters use graphics processing units (GPUs) to offer…
This talk presents results of a comparitive study of iterative algorithms like minimal residue ($MR$) and conjugate gradient ($CG$, $BiCG\gamma_5$, and \bicgstab) used for inverting the Dirac matrix $M$. The tests were done on the…