The Tenth International Workshop on Lattice QFT and Numerical Analysis (QCDNA X)
from
Tuesday, June 27, 2017 (9:00 AM)
to
Thursday, June 29, 2017 (6:00 PM)
Monday, June 26, 2017
Tuesday, June 27, 2017
9:00 AM
Registration
Registration
9:00 AM  9:30 AM
9:30 AM
HMC: Geometry and Applications

Alessandro Barp
HMC: Geometry and Applications
Alessandro Barp
9:30 AM  10:15 AM
10:15 AM
Deflation for MonteCarlo estimation of the trace of a matrix inverse

Eloy Romero Alcalde
(
College of William and Mary
)
Deflation for MonteCarlo estimation of the trace of a matrix inverse
Eloy Romero Alcalde
(
College of William and Mary
)
10:15 AM  10:45 AM
In the context of computing disconnected diagrams, we investigate the efficient estimation of the trace of largescale matrix inverses. Our approach is based on the Hutchinson method (MonteCarlo averaging over matrix quadratures). Previous work showed that combining deflation against the lowest part of the spectrum with Hierarchical Probing can accelerate the convergence significantly. As the size of the matrix grows, however, the computation of enough singular vectors to achieve a reduction in the variance can be a bottleneck. In this work we take advantage of the fact that the singular vectors corresponding to the lowest modes can be represented well in a sparse basis. This allows us to compute and store efficiently the first 500~1000 lowest modes of the matrix spectrum, which is enough to obtain a good reduction in the variance. Moreover we discuss different projectors for deflating and the performance impact when the matrix is close to singular.
10:45 AM
Coffee break
Coffee break
10:45 AM  11:15 AM
11:15 AM
Progress and Challenge of Lattice Quantum Finite Elements (QFE) on Spheres

Richard Brower
(
Boston University
)
Progress and Challenge of Lattice Quantum Finite Elements (QFE) on Spheres
Richard Brower
(
Boston University
)
11:15 AM  12:15 PM
Extending lattice field to ultraviolet complete quantum field theory on any smooth Riemann manifolds is a challenging problem. By adapting element methods (FEM) and Regge geometry one recovers classical (IR) solution in the continuum. However to correctly handle UV divergences requires new counter terms to construct a what we call a "Quantum Finite Elements" (QFE) discrete Lagrangian on the simplicial complex. These UV counters for 2d phi 4th theory and free fermions on the two sphere (S2) have been tested numerically to high precession against the exact Ising solution. Methods to generalize the QFE construction to radial quantized 3d super renormalizable theories on R x S2 and challenges for asymptotical free 4d gauge theories on R x S3 will be presented.
12:30 PM
Lunch
Lunch
12:30 PM  2:00 PM
2:00 PM
Lattice Quantum Field Theory: Numerical Integration in an Infinite Number of Dimensions

Anthony Kennedy
(
Higgs Centre for Theoretical Physics, The University of Edinburgh
)
Lattice Quantum Field Theory: Numerical Integration in an Infinite Number of Dimensions
(Public Talks)
Anthony Kennedy
(
Higgs Centre for Theoretical Physics, The University of Edinburgh
)
2:00 PM  3:00 PM
Relativistic Quantum Field Theory (QFT) is the formalism upon which the Standard Model of particle physics is built, and it is remarkably successful and accurate. When the interactions are weak then the techniques of renormalized perturbation theory using Feynman diagrams works beautifully, but when the interactions are strong we have to turn to numerical evaluation of the functional integrals that define the quantum theory. The only way we know of to evaluate such infinitedimensional integrals is to use Markov Chain Monte Carlo (MCMC) techniques. I shall attempt to give an overview of the functional integral formalism of QFT, and how MCMC integration works. If time permits I will discuss how fermions (halfinteger spin) are dealt with, and maybe introduce the Hybrid (or Hamiltonian) Monte Carlo algorithm.
3:00 PM
Lattice QCD with mixed action  BoriciCreutz valence quarks on staggered sea

Subhasish Basak
(
NISER
)
Lattice QCD with mixed action  BoriciCreutz valence quarks on staggered sea
Subhasish Basak
(
NISER
)
3:00 PM  3:30 PM
We used BoriciCreutz fermion to study discrete chiral symmetry breaking at strong coupling in 2dim GrossNeveu model and mass spectra in 2dim field theories. Mixed action lattice QCD study with BoriciCreutz valence quarks on staggered sea quarks is carried out. The counter terms are fixed nonperturbatively to restore the broken symmetries. The effect of partial quenching and unitarity violation is investigated.
3:30 PM
Quarkantiquark excited flux tube

Nuno Cardoso
(
IST
)
Quarkantiquark excited flux tube
Nuno Cardoso
(
IST
)
3:30 PM  4:00 PM
We present color field profiles for some of the first SU(3) gluonic excitations of the flux tube in the presence of a static quarkantiquark pair.
4:00 PM
Coffee break
Coffee break
4:00 PM  4:30 PM
4:30 PM
DDαAMG on QPACE 3

Daniel Ritchmann
(
University of Regensburg
)
Peter Georg
(
University of Regensburg
)
DDαAMG on QPACE 3
Daniel Ritchmann
(
University of Regensburg
)
Peter Georg
(
University of Regensburg
)
4:30 PM  5:30 PM
We describe our experiences porting the Regensburg implementation of the DDαAMG solver from QPACE 2 to QPACE 3. We first review porting from the first generation Intel Xeon Phi processor (Knights Corner) to its successor (Knights Landing). Secondly, we present an overview of OmniPath, comparing it to the wellknown competitor InfiniBand. Finally, we present the performance of the code on a single processor as well as the scaling on many nodes.
5:30 PM
MILC code performance on high end CPU and GPU supercomputer clusters

Ruizi Li
MILC code performance on high end CPU and GPU supercomputer clusters
Ruizi Li
5:30 PM  6:00 PM
With recent developments in parallel supercomputing architecture, many core, multicore, and GPU processors are now commonplace resulting in more levels of parallelism, memory hierarchy, and programming complexity. It has been necessary to adapt the MILC code to these new processors starting with NVIDIA GPUs and more recently the Intel Xeon Phi processors. We report on our efforts to port and optimize our code for the Intel Knights Landing architecture. We consider performance of the MILC code with MPI and OpenMP, and optimizations with QOPQDP and QPhiX. For the latter approach we concentrate on the staggered conjugate gradient and gauge force. We also consider performance on recent NVIDIA GPUs using the QUDA library.
7:00 PM
Fado ao Centro
Fado ao Centro
7:00 PM  7:50 PM
8:15 PM
Workshop Dinner at "Solar do Bacalhau"
Workshop Dinner at "Solar do Bacalhau"
8:15 PM  10:30 PM
Wednesday, June 28, 2017
9:30 AM
Algorithmic advances in NSPT

Mattia Dalla Brida
(
Universita' di MilanoBicocca & INFN
)
Algorithmic advances in NSPT
Mattia Dalla Brida
(
Universita' di MilanoBicocca & INFN
)
9:30 AM  10:15 AM
Numerical stochastic perturbation theory (NSPT) is a powerful tool that allows perturbation expansions in QCD and other interesting theories to be estimated to high order in the interactions. The standard algorithms on which NSPT is based on, however, suffer from several limitations which in practice restrict the potential of these techniques. In this talk I will review the recent algorithmic advances in this field and show how these significantly reduce the computational effort for precise and accurate determinations. This opens up the way to tackle challenging and interesting new problems, as will be illustrated by a highly nontrivial computation.
10:15 AM
Numerical Stochastic Perturbation Theory in $\varphi^4$ Theory

Marco Garofalo
(
Higgs Centre for Theoretical Physics, The University of Edinburgh
)
Numerical Stochastic Perturbation Theory in $\varphi^4$ Theory
Marco Garofalo
(
Higgs Centre for Theoretical Physics, The University of Edinburgh
)
10:15 AM  10:45 AM
Numerical stochastic perturbation theory (NSPT) is a powerful tool for estimating highorder perturbative expansions in lattice field theory. The standard NSPT is based on the Langevin equation. In this contribution, we investigate in $\varphi^4$ theory some alternative methods. In particular, we present a study of the recently proposed Instantaneous Stochastic Perturbation Theory, as well as a formulation of numerical stochastic perturbation theory based on Generalised Hybrid Molecular Dynamics algorithms.
10:45 AM
Coffee break
Coffee break
10:45 AM  11:15 AM
11:15 AM
Nonperturbative anisotropy calibration in lattice QCD at strong coupling

Helvio Jose Caleiro Vairinhos
(
Fisica Teorica e Met. Mat.Inst. Nacional de Invest. Cientif.Un
)
Nonperturbative anisotropy calibration in lattice QCD at strong coupling
Helvio Jose Caleiro Vairinhos
(
Fisica Teorica e Met. Mat.Inst. Nacional de Invest. Cientif.Un
)
11:15 AM  11:45 AM
We propose a simple criterion for the nonperturbative renormalization of the anisotropy coupling in lattice QCD with massless staggered fermions, in the strong coupling limit. We compute numerically and to high precision the renormalised anisotropy, and the analogue of Karsch’s coefficients, using diagrammatic Monte Carlo algorithms and multihistogram reweighting. We observe a large nonperturbative correction to the mean field anisotropy, and we analyse the implications of such a correction on the continuous time limits of the phase diagram of lattice QCD at strong coupling, and of the baryon mass.
11:45 AM
Algebraic Multigrid: Theory and Practice

James Brannick
Algebraic Multigrid: Theory and Practice
James Brannick
11:45 AM  12:30 PM
This talk gives an overview of recent progress made in the design and analysis of algebraic multigrid methods. The focus is on the setup algorithm that automatically constructs the multilevel hierarchy used in the solve phase. A sharp twogrid theory is introduced and then used to derive various quality measures of the coarse spaces constructed by the setup algorithm, based on the ideas of compatible relaxation, a related identity that assumes the use of the socalled ideal interpolation operator, and an optimal form of classical algebraic multigrid interpolation that gives the best possible twogrid convergence rate. Various numerical results are presented to illustrate these theoretical results. As a test problem, we focus on a finite volume discretization of a scalar diffusion problem with highly varying (discontinuous) diffusion coefficient.
12:30 PM
Lunch
Lunch
12:30 PM  2:00 PM
2:00 PM
GPU Computing to the Exascale and Beyond

Kate Clark
(
NVIDIA
)
GPU Computing to the Exascale and Beyond
(Public Talks)
Kate Clark
(
NVIDIA
)
2:00 PM  3:00 PM
With the demise of Denard scaling, it is well known that we cannot go faster for greater throughput, rather we have to go wider. However, with the imminent demise of Moore's Law, there lies continued challenges in reaching and exceeding the Exascale. We discuss how and why GPU computing provides a solution to take computational science to this next level. We consider some of the software and algorithm challenges that the Exascale will bring. Lastly, we take a dive into the newlylaunched Volta GPU architecture, the architecture that will drive the nextgeneration CORAL Supercomputers.
3:00 PM
HMC, high level structures and architecture support in Grid

Guido Cossu
(
The University of Edinburgh
)
HMC, high level structures and architecture support in Grid
Guido Cossu
(
The University of Edinburgh
)
3:00 PM  3:45 PM
I will present the strategies we developed in Grid for the HMC sector, in order to support a variety of behaviours without code replication. I will also discuss the current status of the architecture support in view of the upcoming machines.
3:45 PM
Analyzing AMA data on 48^3 × 96 lattices

Oliver Witsel
Analyzing AMA data on 48^3 × 96 lattices
Oliver Witsel
3:45 PM  4:10 PM
As an example for analyzing data on RBCUKQCD's 48^3 x 96 ensemble with physical pions, I investigate correlations between propagator sources placed on the same configuration and discuss options for analyzing this data set to extract the Bsmeson decay constant. The data are generated taking advantage of allmode averaging (AMA) to reduce the numerical costs.
4:10 PM
Coffee break
Coffee break
4:10 PM  4:40 PM
4:45 PM
Visit to the Old University
Visit to the Old University
4:45 PM  7:00 PM
Thursday, June 29, 2017
9:30 AM
Scaling Multigrid to the Exascale

Kate Clark
(
NVIDIA
)
Scaling Multigrid to the Exascale
Kate Clark
(
NVIDIA
)
9:30 AM  10:15 AM
Owing to its success in removing the critical slowing down of Dirac linear systems, adaptive multigrid is now a standard solver in the arsenal of tools that the lattice field theorist expects. In this work we report on the latest progress in improving the strong scaling of adaptive multigrid algorithms when running on GPUaccelerated architectures using the QUDA library. Techniques include Schwarz preconditioning, pipelined solvers, precision truncation and RDMAenabled MPI. Furthermore, we report on progress on optimizing the adaptive setup process in order to increase its applicability to Hybrid Monte Carlo. Finally we discuss the challenges in scaling multigrid to the Exascalegeneration of supercomputers.
10:15 AM
Novel Approaches for Staggered Multigrid Algorithms

Evan Weinberg
Novel Approaches for Staggered Multigrid Algorithms
Evan Weinberg
10:15 AM  10:45 AM
Critical slowing down in the fermion sector is a leading obstacle facing the approach to the continuum in lattice gauge theory simulations. Adaptive multigrid ($\alpha$MG) methods offer a solution to the resulting superlinear growth in the cost of iterative Krylov solves. Ongoing research has suggested that previously developed $\alpha$MG methods, such as the successful formulation for Wilsonclover fermions, cannot be applied in a blackbox fashion to other discretizations. In this talk I will expand on subtle issues offered by the staggered discretization. I will also discuss novel approaches that may sidestep these issues.
10:45 AM
Coffee break
Coffee break
10:45 AM  11:15 AM
11:15 AM
An Eigensolver for the Hermitian Dirac Operator with Multigrid Acceleration

Artur Strebel
(
University of Wuppertal
)
An Eigensolver for the Hermitian Dirac Operator with Multigrid Acceleration
Artur Strebel
(
University of Wuppertal
)
11:15 AM  11:45 AM
In this talk we present a Davidson type eigensolver combined with the DD$\alpha$AMG multigrid solver library. The basic Davidson method is adjusted to our multigrid method and the structure of the hermitian Dirac operator in a way that both methods benefit from each other. We compare the resulting eigensolver with a Chebychev filtered Arnoldi method (PARPACK) and the multi purpose eigensolver library PRIMME based on a variety of scaling and performance studies.
11:45 AM
Local Adaptive Refinement on Lattice Gauge Fields

Edward White, Jr.
Local Adaptive Refinement on Lattice Gauge Fields
Edward White, Jr.
11:45 AM  12:15 PM
Adaptive Mesh Reﬁnement (AMR) has been widely used in computational ﬂuid dynamics, shock hydrodynamics, astrophysics, turbulence modeling and combustion to improve the performance of algorithms running large, complex problems. The results are often impressive. However, for various reasons which will be discussed, traditional AMR is not directly applicable to lattice QCD. This talk proposes a new numerical method, Local Adaptive Reﬁnement (LAR), for lattice QCD. LAR is strongly motivated by AMR, and is designed to work with the aggregationbased algebraic multigrid framework of lattice QCD. The premise behind all adaptive reﬁnement methods is that certain numerical problems have solutions that exhibit local variations, which can be examined more fully at higher resolution (more closely spaced grid or lattice points) than the base level of resolution used for the simulation. As a result, problems that would normally require very high resolution over a large domain (which is very computationally expensive) could instead be solved at lower resolution over the large domain, with automated local reﬁnement to higher resolution as needed. The scale and nature of the problems commonly found at the leading edge of today’s lattice QCD work suggest that the field could be a beneﬁt from adaptive reﬁnement. Among the topics that will be discussed are the need for an algebraic variablebased approach to reﬁnement, as opposed to a geometric gridbased approach, methods of determining the regions in which reﬁnement should occur, and ways to implement LAR in a parallel machine environment.
12:15 PM
End of the workshop
End of the workshop
12:15 PM  1:00 PM