Related papers: A general-purpose timestep criterion for simulatio…
We propose a faster algorithm for individual based simulations for adaptive dynamics based on a simple modification to the standard Gillespie Algorithm for simulating stochastic birth-death processes. We provide an analytical explanation…
Gradient algorithms are classical in adaptive control and parameter estimation. For instantaneous quadratic cost functions they lead to a linear time-varying dynamic system that converges exponentially under persistence of excitation…
We study the variational principle over an Hilbert-Einstein like action for an extended geometry taking into account torsion and non-metricity. By extending the semi-Riemannian geometry, we obtain an effective energy-momentum tensor which…
Implementing particle-interaction configuration and time integration are performance intensive essentials of particle-based methods. In this paper, a dual-criteria time-stepping method is proposed to improve the computational efficiency of…
In simulations of fluid motion time accuracy has proven to be elusive. We seek highly accurate methods with strong enough stability properties to deal with the richness of scales of many flows. These methods must also be easy to implement…
The growing field of large-scale time domain astronomy requires methods for probabilistic data analysis that are computationally tractable, even with large datasets. Gaussian Processes are a popular class of models used for this purpose…
Recently, tensor algebra have witnessed significant applications across various domains. Each operator in tensor algebra features different computational workload and precision. However, current general accelerators, such as VPU, GPGPU, and…
It is shown in this article that if the Einstein Equivalence Principle is valid on a particular metric theory of gravitation in a spherically symmetric space-time, then the time metric component is not equal to the negative of the inverse…
We present a new implicit asymptotic preserving time integration scheme for charged-particle orbit computation in arbitrary electromagnetic fields. The scheme is built on the Crank-Nicolson integrator and continues to recover full-orbit…
Computational efficiency demands discretised, hierarchically organised, and individually adaptive time-step sizes (known as the block-step scheme) for the time integration of N-body models. However, most existing N-body codes adapt…
This paper describes a fourth-order integration algorithm for the gravitational N-body problem based on discrete Lagrangian mechanics. When used with shared timesteps, the algorithm is momentum conserving and symplectic. We generalize the…
With the improvement in sensitivity of gravitational wave (GW) detectors and the increasing diversity of GW sources, there is a strong need for accurate GW waveform models for data analysis. While the current model accuracy assessments…
We introduce a scalable approach to Gaussian process inference that combines spatio-temporal filtering with natural gradient variational inference, resulting in a non-conjugate GP method for multivariate data that scales linearly with…
One of the main targets for ground-based gravitational wave (GW) detectors such as Advanced LIGO (Laser Interferometer Gravitational wave Observatory) and Virgo is coalescences of neutron star (NS) binaries. Even though a NS's macroscopic…
Error entropy is a important nonlinear similarity measure, and it has received increasing attention in many practical applications. The default kernel function of error entropy criterion is Gaussian kernel function, however, which is not…
The strong equivalence principle is extended in application to averaged dynamical fields in cosmology to include the role of the average density in the determination of inertial frames. The resulting cosmological equivalence principle is…
We consider in this paper a class of composite optimization problems whose objective function is given by the summation of a general smooth and nonsmooth component, together with a relatively simple nonsmooth term. We present a new class of…
We introduce the public version of the cosmological magnetohydrodynamical moving-mesh simulation code Arepo. This version contains a finite-volume magnetohydrodynamics algorithm on an unstructured, dynamic Voronoi tessellation coupled to a…
We review various approaches to approximating general relativistic effects in hydrodynamic simulations of stellar core collapse and post-bounce evolution. Different formulations of a modified Newtonian gravitational potential are presented.…
Galaxy clusters are the largest gravitationally bound systems, and they continue their growth through mergers in a hierarchical {\Lambda}CDM Universe. Therefore, we can describe the merger stage of a cluster as the dynamical state of…