Related papers: RHMC with Block Solvers and Multiple Pseudofermion…
Robust and efficient solvers for coupled-adjoint linear systems are crucial to successful aerostructural optimization. Monolithic and partitioned strategies can be applied. The monolithic approach is expected to offer better robustness and…
We consider the problem of estimating the probability of a large loss from a financial portfolio, where the future loss is expressed as a conditional expectation. Since the conditional expectation is intractable in most cases, one may…
We investigate the application of Krylov space methods to the solution of shifted linear systems of the form (A+\sigma) x - b = 0 for several values of \sigma simultaneously, using only as many matrix-vector operations as the solution of a…
Machine learning and deep learning have revolutionized computational physics, particularly the simulation of complex systems. Equivariance is essential for simulating physical systems because it imposes a strong inductive bias on the…
Lattice QCD calculations including the effects of one or more non-degenerate sea quark flavors are conventionally performed using the Rational Hybrid Monte Carlo (RHMC) algorithm, which computes the square root of the determinant of…
I review recent machine trends and algorithmic developments for dynamical lattice QCD simulations with the HMC algorithm for Wilson-type fermions. The topics include the trend toward multi-core processors and general purpose GPU (GPGPU)…
We present preliminary results of a calculation of the QCD renormalization constants (RCs) for quark bilinear operators, evaluated non-perturbatively on the lattice in the RI'-MOM scheme. The calculation is performed by using dedicated…
This article gives an introduction to the multilevel blocking (MLB) approach to both the fermion and the dynamical sign problem in path-integral Monte Carlo simulations. MLB is able to substantially relieve the sign problem in many…
We report about on-going simulations of $N_f=2+1$ lattice QCD. We use a tadpole improved Symanzik gauge action and stout link smeared Wilson fermions with a clover term. We employ the Hasenbusch trick for the degenerate u- and d-quarks, and…
Frequent Directions, as a deterministic matrix sketching technique, has been proposed for tackling low-rank approximation problems. This method has a high degree of accuracy and practicality, but experiences a lot of computational cost for…
Neutronic calculations for reactors are a daunting task when using Monte Carlo (MC) methods. As high-performance computing has advanced, the simulation of a reactor is nowadays more readily done, but design and optimization with multiple…
We review a number of topics related to block variable renormalisation group transformations of quantum fields on the lattice, and to the emerging perfect lattice actions. We first illustrate this procedure by considering scalar fields.…
We investigate the performance of the hybrid Monte Carlo algorithm, the standard algorithm used for lattice QCD simulations involving fermions, in updating non-trivial global topological structures. We find that the hybrid Monte Carlo…
In this work we investigate replacing standard quadrature techniques used in deterministic linear solvers with a fixed-seed Quasi-Monte Carlo calculation to obtain more accurate and efficient solutions to the neutron transport equation…
Clover improvement is the standard choice for lattice simulations of QCD, when the lattice artefacts coming from Wilson fermions have to be reduced. However, the clover improvement is not limited to QCD, but can be applied to a wider range…
The quark chromoelectric dipole (qCEDM) operator is a CP-violating operator describing, at hadronic energies, beyond-the-standard-model contributions to the electric dipole moment of particles with nonzero spin. In this paper we define…
Complex Langevin simulations provide an alternative to sample path integrals with complex weights and therefore are suited to determine the phase diagram of QCD from first principles. We use our proposed method of Dynamic Stabilisation (DS)…
Considering Ginsparg-Wilson type fermions dynamically in lattice QCD simulations is a challenging task. The hope is to be able to approach smaller pion masses and to eventually reach physical situations. The price to pay is substantially…
We present simulation results for lattice QCD with chiral fermions in small volumes, where the epsilon-expansion of chiral perturbation theory applies. Our data for the low lying Dirac eigenvalues, as well as mesonic correlation functions,…
The fermion determinant is a highly non-local object and its logarithm is an extensive quantity. For these reasons it is widely believed that the determinant cannot be treated in acceptance steps of gauge link configurations that differ in…