Related papers: A comparison of updating algorithms for large N re…
Diagrammatic Monte Carlo methods provide robust routines for accurate computations of correlated electronic systems in the thermodynamical limit. Recently, its versatility was extended to SU(N) Hubbard model, where the core is a novel…
Hidden Markov models (HMMs) are popular models to identify a finite number of latent states from sequential data. However, fitting them to large data sets can be computationally demanding because most likelihood maximization techniques…
This paper introduces fast R updating algorithms specifically designed for statistical applications, including regression, filtering, and model selection, where data structures change frequently. Although traditional QR decomposition is…
Using a dual representation of lattice fermion models that is based on spin-charge transformation and fermionisation of the original description, I derive an algorithm for diagrammatic Monte Carlo simulation of strongly correlated systems.…
We present an exact Monte Carlo algorithm designed to sample theories where the energy is a sum of many couplings of decreasing strength. The algorithm avoids the computation of almost all non-leading terms. Its use is illustrated by…
This paper proposes the use of particle swarm optimization method (PSO) for finite element (FE) model updating. The PSO method is compared to the existing methods that use simulated annealing (SA) or genetic algorithms (GA) for FE model for…
We have continued our systematic investigations of the numerical simulations of lattice gauge theories in the dual formulation. These include: i) a more practical implementation of the quasi-local updating technique, ii) a thorough…
The accuracy of Green Function Monte Carlo (GFMC) simulations can be greatly improved by a clever choice of the approximate ground state wave function that controls configuration sampling. This trial wave function typically depends on many…
Scalable decentralized optimization in large-scale systems hinges on efficient communication. A common way to reduce communication overhead is to perform multiple local updates between two communication rounds, as in federated learning.…
We present a study on the performance of Wang-Landau algorithm in a lattice model of liquid crystals which is a continuous lattice spin model. We propose a novel method of the spin update scheme in a continuous lattice spin model. The…
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…
Nonnegative matrix factorization (NMF) is a data analysis technique used in a great variety of applications such as text mining, image processing, hyperspectral data analysis, computational biology, and clustering. In this paper, we…
We improve the algorithm to noninvasively update the response matrix using information from the orbit-feedback system, described in [1]. The new version is capable of adapting to slow changes of the lattice, albeit at the expense of…
At small lattice spacing, or when using overlap fermions, lattice QCD simulations tend to become stuck in a single topological sector. Physical observables, e.g.\ hadron masses, then differ from their full QCD counterparts by $1/V$…
We study and derive algorithms for nonlinear eigenvalue problems, where the system matrix depends on the eigenvector, or several eigenvectors (or their corresponding invariant subspace). The algorithms are derived from an implicit…
The results of numerical simulation using a modified Monte Carlo method with a heat bath algorithm for the pseudospin model of cuprates are presented. The temperature phase diagrams are constructed for various degrees of doping and for…
Monte-Carlo sampling of lattice model Hamiltonians is a well-established technique in statistical mechanics for studying the configurational entropy of crystalline materials. When species to be distributed on the lattice model carry charge,…
We present a simulation algorithm for Hamiltonian fermion lattice models. A guiding trial wave function is adaptively optimized during Monte Carlo evolution. We apply the method to the two dimensional Gross-Neveu model and analyze systematc…
A new algorithm for simulation of theories with dynamical fermions is presented. The algorithm is based on obtaining the new configuration U' from the old one U by solving the equation M(U')\eta= \omega M(U)\eta, where M is fermionic…
Gauge fixing is an essential step in lattice QCD calculations, particularly for studying gauge-dependent observables. Traditional iterative algorithms are computationally expensive and often suffer from critical slowing down and scaling…