Related papers: PLQCD library for Lattice QCD on multi-core machin…
Applying domain decomposition to the lattice Dirac operator and the associated quark propagator, we arrive at expressions which, with the proper insertion of random sources therein, can provide improvement to the estimation of the…
Development of scientific software involves tradeoffs between ease of use, generality, and performance. We describe the design of a general hyperbolic PDE solver that can be operated with the convenience of MATLAB yet achieves efficiency…
A numerical simulation algorithm for lattice QCD is described, in which the short- and long-distance effects of the sea quarks are treated separately. The algorithm can be regarded, to some extent, as an implementation at the quantum level…
We discuss a substantial update to the Grid software library for Lattice QCD, enabling it to port to multiple GPU architectures while retaining CPU vectorisation and SIMD execution within OpenMP threads. The GPU environments supported…
Since the development of the Beowulf project to build a parallel computer from commodity PC components, there have been many such clusters built. The MILC QCD code has been run on a variety of clusters and supercomputers. Key design…
This paper consists of a tutorial on the Augmented Lagrangian based Alternating Direction Inexact Newton method (ALADIN) and its application to distributed model predictive control (MPC). The focus is - for simplicity of presentation - on…
The lattice technique of studying the strong interaction of matter is used to obtain predictions of the hadronic spectrum. These simulations were performed by the UKQCD collaboration using full (unquenched) QCD. Details of the results, a…
Simulation of Lattice QCD is a challenging computational problem. Currently, technological trends in computation show multiple divergent models of computation. We are witnessing homogeneous multi-core architectures, the use of accelerator…
This paper reviews the most popular methods which are used in lattice QCD to compute the determinant of the lattice Dirac operator: Gaussian integral representation and noisy methods. Both of them lead naturally to matrix function problems.…
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…
We simulate Quantum Chromodynamics (QCD) in four Euclidean dimensions with two (degenerate mass) flavors of dynamical quarks. The Dirac operator we use is the so-called chirally improved Dirac operator. We discuss the algorithm used for the…
I define lattice fermions in five Euclidean dimensions and the corresponding effective theory in four dimensions. The main properties of these theories include the suppression of high momentum modes of the lattice Dirac operator and their…
Applying domain decomposition to the lattice Dirac operator and the associated quark propagator, we arrive at expressions which, with the proper insertion of random sources therein, can provide improvement to the estimation of the…
A simulation of quenched QCD with the overlap Dirac operator has been completed using 100 Wilson gauge configurations at beta = 6 on an 18^3 x 64 lattice and at beta = 5.85 on a 14^3 x 48 lattice, both in Landau gauge. We present results…
We describe the construction of a high performance parallel computer composed of PC components, as well as the performance test in lattice QCD.
We suggest that the lattice Dirac spectra in QCD at finite temperature may be understood using a gaussian unitary ensemble for Wilson fermions, and a chiral gaussian unitary ensemble for Kogut-Susskind fermions. For Kogut-Susskind fermions,…
We present an update of BQCD, our Hybrid Monte Carlo program for simulating lattice QCD. BQCD is one of the main production codes of the QCDSF collaboration and is used by CSSM and in some Japanese finite temperature and finite density…
This paper introduces QCDLAB, a design and research tool for lattice QCD algorithms. The tool, a collection of MATLAB functions, is based on a ``small-code'' and a ``minutes-run-time'' algorithmic design philosophy. The present version uses…
We investigate the state-of-the-art Lanczos eigensolvers available in the Grid and QUDA libraries. They include Implicitly Restarted Lanczos, Thick-Restart Lanczos, and Block Lanczos. We measure and analyze their performance for the Highly…
The increasing imbalance between the computing capabilities of individual nodes and the internode bandwidth makes it highly desirable for any Lattice QCD algorithm to minimize the amount of internode communication. One of the relatively new…