Skip to Main content Skip to Navigation

Automated Code Generation for Lattice QCD Simulation

Denis Barthou 1, 2 Gilbert Grosdidier 3 Konstantin Petrov 3, 4, 5 Michael Kruse 6, 7 Christine Eisenbeis 6, 7 Olivier Pène 4 Olivier Brand-Foissac 4 Claude Tadonki 8 Romain Dolbeau 9 
2 RUNTIME - Efficient runtime systems for parallel architectures
Inria Bordeaux - Sud-Ouest, UB - Université de Bordeaux, CNRS - Centre National de la Recherche Scientifique : UMR5800
6 GRAND-LARGE - Global parallel and distributed computing
LRI - Laboratoire de Recherche en Informatique, LIFL - Laboratoire d'Informatique Fondamentale de Lille, UP11 - Université Paris-Sud - Paris 11, Inria Saclay - Ile de France, CNRS - Centre National de la Recherche Scientifique : UMR8623
7 ParSys - LRI - Systèmes parallèles (LRI)
LRI - Laboratoire de Recherche en Informatique
Abstract : Quantum Chromodynamics (QCD) is the theory of strong nuclear force, responsible of the interactions between sub-nuclear particles. QCD simulations are typically performed through the lattice gauge theory approach, which provides a discrete analytical formalism called LQCD (Lattice Quantum Chromodynamics). LQCD simulations usually involve generating and then processing data on petabyte scale which demands multiple teraflop-years on supercomputers. Large parts of both, generation and analysis, can be reduced to the inversion of an extremely large matrix, the so-called Wilson-Dirac operator. For this purpose, and because this matrix is always sparse and structured, iterative methods are definitely considered. Therefore, the procedure of the application of this operator, resulting in a vector-matrix product, appears as a critical computation kernel that should be optimized as much as possible. Evaluating the Wilson-Dirac operator involves symmetric stencil calculations where each node has 8 neighbors. Such configuration is really hindering when it comes to memory accesses and data exchanges among processors. For current and future generation of supercomputers the hierarchical memory structure make it next to impossible for a physicist to write an efficient code. Addressing these issues in other to harvest an acceptable amount of computing cycles for the real need, which means reaching a good level of efficiency, is the main concern of this paper. We present here a Domain Specific Language and corresponding toolkit, called QIRAL, which is a complete solution from symbolic notation to simulation code.
Complete list of metadata
Contributor : Claude Tadonki Connect in order to contact the contributor
Submitted on : Sunday, December 15, 2013 - 9:45:35 AM
Last modification on : Sunday, June 26, 2022 - 12:00:19 PM
Long-term archiving on: : Tuesday, March 18, 2014 - 2:46:39 PM


Files produced by the author(s)


  • HAL Id : hal-00918812, version 1


Denis Barthou, Gilbert Grosdidier, Konstantin Petrov, Michael Kruse, Christine Eisenbeis, et al.. Automated Code Generation for Lattice QCD Simulation. [Research Report] RR-8417, INRIA. 2013, pp.13. ⟨hal-00918812⟩



Record views


Files downloads