Adaptive control in rollforward recovery for extreme scale multigrid

With the increasing number of compute components, failures in future exa-scale computer systems are expected to become more frequent. This motivates the study of novel resilience techniques. Here, we extend a recently proposed algorithm-based recovery method for multigrid iterations by introducing an adaptive control. After a fault, the healthy part of the system continues the iterative solution process, while the solution in the faulty domain is re-constructed by an asynchronous on-line recovery. The computations in both the faulty and healthy subdomains must be coordinated in a sensitive way, in particular, both under and over-solving must be avoided. Both of these waste computational resources and will therefore increase the overall time-to-solution. To control the local recovery and guarantee an optimal re-coupling, we introduce a stopping criterion based on a mathematical error estimator. It involves hierarchical weighted sums of residuals within the context of uniformly refined meshes and is well-suited in the context of parallel high-performance computing. The re-coupling process is steered by local contributions of the error estimator. We propose and compare two criteria which differ in their weights. Failure scenarios when solving up to $6.9\cdot10^{11}$ unknowns on more than 245\,766 parallel processes will be reported on a state-of-the-art peta-scale supercomputer demonstrating the robustness of the method

Spurious divergences in Dyson-Schwinger equations

We revisit the treatment of spurious ultraviolet divergences in the equation of motion of the gluon propagator caused by a momentum cutoff and the resulting violation of gauge invariance. With present continuum studies of the gluon propagator from its Dyson-Schwinger equation reaching the level of quantitatively accurate descriptions, it becomes increasingly important to understand how to subtract these spurious divergences in an unambiguous way. Here we propose such a method. It is based entirely on the asymptotic perturbative behavior of the QCD Green's functions without affecting non-perturbative aspects such as mass terms or the asymptotic infrared behavior. As a particular example, this allows us to assess the possible influence of the tadpole diagram beyond perturbation theory. Finally, we test this method numerically by solving the system of Dyson-Schwinger equations of the gluon and ghost propagators.Comment: 19 pages, 9 figs; agrees with published versio

Going beyond the propagators of Landau gauge Yang-Mills theory

We present results for the propagators and the ghost-gluon vertex of Landau gauge Yang-Mills theory obtained from Dyson-Schwinger equations. Solving these three quantities simultaneously constitutes a new step in truncating these equations. We also introduce a new model for the three-gluon vertex that is motivated by lattice results. It features a zero crossing which is confirmed a posteriori by a Dyson-Schwinger calculation. Within our setup we can reproduce lattice data very well. We establish that also for the ghost-gluon vertex a difference between decoupling and scaling solutions is present. For the scaling solution we discuss the possibility of modifying the infrared exponents via an angle dependence of the ghost-gluon vertex. However, no such dependence is found in our calculations. Finally, we calculate the Schwinger function for the gluon propagator.Comment: 8 pages, Confinement X proceeding

On the influence of three-point functions on the propagators of Landau gauge Yang-Mills theory

We solve the Dyson-Schwinger equations of the ghost and gluon propagators of Landau gauge Yang-Mills theory together with that of the ghost-gluon vertex. The latter plays a central role in many truncation schemes for functional equations. By including it dynamically we can determine its influence on the propagators. We also suggest a new model for the three-gluon vertex motivated by lattice data which plays a crucial role to obtain stable solutions when the ghost-gluon vertex is included. We find that both vertices have a sizable quantitative impact on the mid-momentum regime and contribute to the reduction of the gap between lattice and Dyson-Schwinger equation results. Furthermore, we establish that the three-gluon vertex dressing turns negative at low momenta as suggested by lattice results in three dimensions.Comment: 28 pages, 12 figures, matches published versio