Related papers: Warm starting the projected Gauss-Seidel algorithm…
The recent accelerated growth in the computing power has generated popularization of experimentation with dynamic computer models in various physical and engineering applications. Despite the extensive statistical research in computer…
A Gaussian process has been one of the important approaches for emulating computer simulations. However, the stationarity assumption for a Gaussian process and the intractability for large-scale dataset limit its availability in practice.…
We investigate the adiabatic index $\Gamma_{\mathrm{th}}$ of hot and dense nuclear matter from chiral effective field theory and find that the results are systematically larger than from typical mean field models. We start by constructing…
We present several numerical results on granular mixtures. In particular, we examine the efficiency of diffusion as a mixing mechanism in these systems. The collisions are inelastic and to compensate the energy loss, we thermalize the…
The behavior of the average velocity, its deviation and average squared velocity are characterized using three techniques for a 1-D dissipative impact system. The system -- a particle, or an ensemble of non interacting particles, moving in…
We consider the Gibbs measure of a general interacting particle system for a certain class of ``weakly interacting" kernels. In particular, we show that the local point process converges to a Poisson point process as long as the inverse…
We present methods for emulating the matter power spectrum by combining information from cosmological $N$-body simulations at different resolutions. An emulator allows estimation of simulation output by interpolating across the parameter…
This paper studies the effective convergence of iterative methods for solving convex minimization problems using block Gauss--Seidel algorithms. It investigates whether it is always possible to algorithmically terminate the iteration in…
One of important questions concerning particles which compose the Dark Matter (DM) is their average speed. We consider the model of relativistic weakly interacting massive particles and try to impose an upper bound on their actual and past…
Coarse-grained models that preserve hydrodynamics provide a natural approach to study collective properties of soft-matter systems. Here, we demonstrate that commonly used integration schemes in dissipative particle dynamics give rise to…
A method is presented for the computation of the one-loop effective action at finite temperature and density. The method is based on an expansion in the number of spatial covariant derivatives. It applies to general background field…
We present a simplified and fast method for simulating minor mergers between galaxy clusters. Instead of following the evolution of the dark matter halos directly by the N-body method, we employ a rigid potential approximation for both…
To address the time-consuming and computationally intensive issues of traditional ART algorithms for flame combustion diagnosis, inspired by flame simulation technology, we propose a novel representation method for flames. By modeling the…
The method of Alternating Projections (AP) is a fundamental iterative technique with applications to problems in machine learning, optimization and signal processing. Examples include the Gauss-Seidel algorithm which is used to solve…
We introduce a broad class of models called semiparametric spatial point process for making inference between spatial point patterns and spatial covariates. These models feature an intensity function with both parametric and nonparametric…
Global optimization heuristics are popular to optimize hard non-convex problems. Despite their irrefutably large cost-to-solution, in the lack of other working greedy or convex approaches, global optimization algorithms remain the…
For many survey-based spatial modelling problems, responses are observed as spatially aggregated over survey regions due to limited resources. Covariates, from weather models and satellite imageries, can be observed at many different…
A method for computing the thermopower in interacting systems is proposed. This approach, which relies on Monte Carlo simulations, is illustrated first for a diatomic chain of hard-point elastically colliding particles and then in the case…
We have presented some practical consequences on the molecular-dynamics simulations arising from the numerical algorithm published recently in paper Int. J. Mod. Phys. C 16, 413 (2005). The algorithm is not a finite-difference method and…
Many processes can correspond to reactive impregnation in porous solids. These processes are usually numerically computed by classical methods like finite element method, finite volume method, etc. The disadvantage of these methods remains…