相关论文: Algorithm Shootout: R versus RHMC

We study discretization effects of the Wilson and staggered Dirac operator with $N_{\rm c}>2$ using chiral random matrix theory (chRMT). We obtain analytical results for the joint probability density of Wilson-chRMT in terms of a…

We present a completed random matrix theory for staggered fermions which incorporates all taste symmetry breaking terms at their leading order from the staggered chiral Lagrangian. This is an extension of previous work which only included…

In the work we have considered the reduction techniques matrix element of the reaction interaction of elementary particles for massive fermions to spinor products . The our procedure is not more complex than CALCUL spinor techniques for…

In this paper, we evaluate the performance of four randomized optimization algorithms: Randomized Hill Climbing (RHC), Simulated Annealing (SA), Genetic Algorithms (GA), and MIMIC (Mutual Information Maximizing Input Clustering), across…

The randomized Kaczmarz (RK) method is an iterative method for approximating the least-squares solution of large linear systems of equations. The standard RK method uses sequential updates, making parallel computation difficult. Here, we…

We test the finite density algorithm in the canonical ensemble which combines the HMC update with the accept/reject step according to the ratio of the fermion number projected determinant to the unprojected one as a way of avoiding the…

We present a quantum algorithm for the simulation of molecular systems that is asymptotically more efficient than all previous algorithms in the literature in terms of the main problem parameters. As in previous work [Babbush et al., New…

We extend a technique for the chiral extrapolation of hadron masses calculated with dynamical fermions to those generated by quenched simulations. The method ensures the correct leading and next-to-leading non-analytic behaviour for either…

As a prerequisite to dynamical fermion simulations a detailed study of optimal parameters and scaling behavior is conducted for the quenched Schr\"odinger functional at fixed renormalized coupling. We compare standard hybrid overrelaxation…

We present a detailed design of a (P)HMC simulation algorithm for N_f = 2 + 1 + 1 maximally twisted Wilson quark flavours. The algorithm retains even/odd and mass-shift preconditionings combined with multiple Molecular Dynamics time scales…

Lattice chiral fermions are synonymous to the Ginsparg-Wilson relation. Indeed, this relation is satisfied by the overlap, domain wall and perfect action fermion kernel. In a recent work we have shown that it is possible to take a direct RG…

A new formulation of chiral fermions on the lattice is presented. It is a version of overlap fermions, but built from the computationally efficient staggered fermions rather than the previously used Wilson fermions. The construction reduces…

Randomized algorithms for very large matrix problems have received a great deal of attention in recent years. Much of this work was motivated by problems in large-scale data analysis, and this work was performed by individuals from many…

We propose a modification of the Hybrid-Monte-Carlo algorithm that allows for a larger step-size of the integration scheme at constant acceptance rate. The key ingredient is that the pseudo-fermion action is split into two parts. We test…

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)…

Let R be a discrete valuation ring (DVR) and K be its fraction field. If M is a matrix over R admitting a LU decomposition, it could happen that the entries of the factors L and U do not lie in R, but just in K. Having a good control on the…

I explore computer simulations of the dynamics of small multi-fermion lattice systems. The method is more general, but I concentrate on Hubbard type models where the fermions hop between a small number of connected sites. I use the natural…

In this work, we consider a rational approximation of the exponential function to design an algorithm for computing matrix exponential in the Hermitian case. Using partial fraction decomposition, we obtain a parallelizable method, where the…

We discuss the construction of a chiral random matrix model for staggered fermions. This model includes $O(a^2)$ corrections to the continuum limit of staggered fermions and is related to the zero momentum limit of the Lee-Sharpe Lagrangian…

Rank-revealing matrix decompositions provide an essential tool in spectral analysis of matrices, including the Singular Value Decomposition (SVD) and related low-rank approximation techniques. QR with Column Pivoting (QRCP) is usually…