Related papers: Performance of several Lanczos eigensolvers with H…
Quadratic minimization problems with orthogonality constraints (QMPO) play an important role in many applications of science and engineering. However, some existing methods may suffer from low accuracy or heavy workload for large-scale…
We use perturbative Symanzik improvement to create a new staggered-quark action (HISQ) that has greatly reduced one-loop taste-exchange errors, no tree-level order a^2 errors, and no tree-level order (am)^4 errors to leading order in the…
We study the lattice spacing dependence, or scaling, of physical quantities using the highly improved staggered quark (HISQ) action introduced by the HPQCD/UKQCD collaboration, comparing our results to similar simulations with the asqtad…
A modern Fortran implementation of three Dirac operators (Wilson, Brillouin, Susskind) in lattice QCD is presented, based on OpenMP shared-memory parallelization and SIMD pragmas. The main idea is to apply a Dirac operator to $N_v$ vectors…
The growing availability and usage of low precision foating point formats has attracts many interests of developing lower or mixed precision algorithms for scientific computing problems. In this paper we investigate the possibility of…
The fate of heavy quarkonia states in quark-gluon plasma is encoded in the temperature dependence of their spectral functions. Reconstruction of spectral functions from Euclidean lattice correlators is an ill-posed problem. Despite a…
We present the results of an effort to accelerate a Rational Hybrid Monte Carlo (RHMC) program for lattice quantum chromodynamics (QCD) simulation for 2 flavours of staggered fermions on multiple Kepler K20X GPUs distributed on different…
We have constructed the lowest few eigenvectors of the staggered Dirac operator on SU(3) gauge configurations, both quenched and dynamical. We use these modes to study the topological charge and to construct approximate hadronic…
We report on numerical experiments using deflation to compute quark propagators for the highly improved staggered quark (HISQ) action. The method is tested on HISQ gauge configurations, generated by the MILC collaboration, with lattice…
OpenQ$^\star$D code has been used by the RC$^\star$ collaboration for the generation of fully dynamical QCD+QED gauge configurations with C$^\star$ boundary conditions. In this talk, optimization of solvers provided with the openQ$^\star$D…
We report preliminary results on the chiral and deconfinement aspects of the QCD transition at finite temperature using the Highly Improved Staggered Quark (HISQ) action on lattices with temporal extent of N_{\tau}=6 and 8. The chiral…
We study spectra, localization properties and local chirality of eigenvectors of the lattice Dirac operator. We analyze ensembles of quenched SU(3) configurations on both sides of the QCD phase transition. Our Dirac operator is a systematic…
Complete spectra of the staggered Dirac operator $\Dirac$ are determined in four-dimensional $SU(2)$ gauge fields with and without dynamical fermions. An attempt is made to relate the performance of multigrid and conjugate gradient…
We calculate the one loop renormalisation parameters for the heavy-light axial-vector and vector currents using lattice perturbation theory. We use NonRelativistic QCD (NRQCD) heavy quarks and the Highly Improved Staggered Quark (HISQ)…
In this paper we present an improved dqds algorithm for computing all the singular values of a bidiagonal matrix to high relative accuracy. There are two key contributions: a novel deflation strategy that improves the convergence for badly…
The runtime of a Lattice QCD simulation is dominated by a small kernel, which calculates the product of a vector by a sparse matrix known as the "Dslash" operator. Therefore, this kernel is frequently optimized for various HPC…
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…
Computing eigenvalues is a computationally intensive task central to numerous applications in the natural sciences. Toward this end, we investigate the quantum block Krylov subspace projector (QBKSP) algorithm - a multireference quantum…
This paper proposes hybrid high-order eigensolvers for the computation of guaranteed lower eigenvalue bounds. These bounds display higher order convergence rates and are accessible to adaptive mesh-refining algorithms. The involved…
The overlap operator in lattice QCD requires the computation of the sign function of a matrix, which is non-Hermitian in the presence of a quark chemical potential. In previous work we introduced an Arnoldi-based Krylov subspace…