Related papers: Linear systems solvers - recent developments and i…
Master-field simulations offer an approach to lattice QCD in which calculations are performed on a small number of large-volume gauge-field configurations. The latter is advantageous for simulations in which the global topological charge is…
In the context of isogeometric analysis, we consider two discretization approaches that make the resulting stiffness matrix nonsymmetric even if the differential operator is self-adjoint. These are the collocation method and the…
Bivariate matrix functions provide a unified framework for various tasks in numerical linear algebra, including the solution of linear matrix equations and the application of the Fr\'echet derivative. In this work, we propose a novel…
We propose a method to improve lattice operators composed of Wilson fermions which allows the removal of all corrections of $O(a)$, including those proportional to the quark mass, leaving only errors of $O(a^2)$. The method exploits the…
Chiral perturbation theory for eigenvalue distributions, and equivalently random matrix theory, has recently been extended to include lattice effects for Wilson fermions. We test the predictions by comparison to eigenvalue distributions of…
Recently developed analytic methods in the framework of the Field Correlator Method are reviewed in this series of four lectures and results of calculations are compared to lattice data and experiment. Recent lattice data demonstrating the…
The novel contribution of this paper relies in the proposal of a fully implicit numerical method designed for nonlinear degenerate parabolic equations, in its convergence/stability analysis, and in the study of the related computational…
It is interesting to superimpose the Pauli-Villars regularization on the lattice regularization. We illustrate how this scheme works by evaluating the axial anomaly in a simple lattice fermion model, the Pauli-Villars Lagrangian with a…
We study the expectation value of the phase of the fermion determinant for Wilson lattice fermions with chemical potential. We use quenched SU(3) ensembles and implement a recently proposed exact dimensional reduction of the fermion…
We present and analyze a class of nonsymmetric preconditioners within a normal (weighted least-squares) matrix form for use in GMRES to solve nonsymmetric matrix problems that typically arise in finite element discretizations. An example of…
We discuss the successes and limitations of statistical sampling for a sequence of models studied in the context of lattice QCD and emphasize the need for new methods to deal with finite-density and real-time evolution. We show that these…
Most efficient linear solvers use composable algorithmic components, with the most common model being the combination of a Krylov accelerator and one or more preconditioners. A similar set of concepts may be used for nonlinear algebraic…
Wilson's formulation of lattice QCD is attractive for many reasons, but perhaps mainly because of its simplicity and conceptual clarity. Numerical simulations of the Wilson theory (and of its improved versions) tend to be extremely…
We consider the iterative solution of large linear systems of equations in which the coefficient matrix is the sum of two terms, a sparse matrix $A$ and a possibly dense, rank deficient matrix of the form $\gamma UU^T$, where $\gamma > 0$…
We propose an augmented Lagrangian-based preconditioner to accelerate the convergence of Krylov subspace methods applied to linear systems of equations with a block three-by-three structure such as those arising from mixed finite element…
The extreme computational costs of calculating the sign of the Wilson matrix within the overlap operator have so far prevented four dimensional dynamical overlap simulations on realistic lattice sizes, because the computational power…
New solution method for the systems of linear equations in commutative integral domains is proposed. Its complexity is the same that the complexity of the matrix multiplication.
We study lattice formulations of the two-dimensional N=2 Wess-Zumino model with a cubic superpotential. Discretizations with and without lattice supersymmetries are compared. We observe that the "Nicolai improvement" introduces new problems…
The worldline casting of a gauge field system with spin-1/2 matter fields has provided a, particle-based, first quantization formalism in the framework of which the Bern-Kosower algorithms for efficient computations in QCD acquire a simple…
We discuss a number of lattice fermion actions solving the Ginsparg-Wilson relation. We also consider short ranged approximate solutions. In particular, we are interested in reducing the lattice artifacts, while avoiding (or suppressing)…