Related papers: COSMICA: a novel parallel GPU code for Cosmic Rays…
Numerical simulations of models and theories that describe complex systems such as spin glasses are becoming increasingly important. Beyond fundamental research, these computational methods also find practical applications in fields like…
Context. The diffusive shock acceleration mechanism has been widely accepted as the acceleration mechanism for galactic cosmic rays. While self-consistent hybrid simulations have shown how power-law spectra are produced, detailed…
For the understanding of the origin and propagation of ultra-high energy cosmic rays (UHECR) we developed a new approach to simulating UHECRs from an arbitrary number of sources based on Monte Carlo technique. The method consists of a…
We introduce a GPU-accelerated Monte Carlo framework for nonconvex, free-final-time trajectory optimization problems. This framework makes use of the prox-linear method, which belongs to the larger family of sequential convex programming…
Multidimensional nucleosynthesis studies with hundreds of nuclei linked through thousands of nuclear processes are still computationally prohibitive. To date, most nucleosynthesis studies rely either on hydrostatic/hydrodynamic simulations…
The Morse-Smale complex is a well studied topological structure that represents the gradient flow behavior between critical points of a scalar function. It supports multi-scale topological analysis and visualization of feature-rich…
Matrix multiplication is a foundational operation in scientific computing and machine learning, yet its computational complexity makes it a significant bottleneck for large-scale applications. The shift to parallel architectures, primarily…
Massively parallel architectures offer the potential to significantly accelerate an application relative to their serial counterparts. However, not all applications exhibit an adequate level of data and/or task parallelism to exploit such…
Some experiments have found multi-core events in cosmic rays extensive air showers which should be interpreted by hadronic interaction theory. In this paper, the multi-core events are reproduced by Monte Carlo simulation with CORSIKA. The…
$N$-body simulations are computationally expensive, so machine-learning (ML)-based emulation techniques have emerged as a way to increase their speed. Although fast, surrogate models have limited trustworthiness due to potentially…
Solving the question of the origin of ultra-high energy cosmic rays (UHECRs) requires the development of detailed simulation tools in order to interpret the experimental data and draw conclusions on the UHECR universe. CRPropa is a public…
FLUKA is a general purpose Monte Carlo transport and interaction code used for fundamental physics and for a wide range of applications. These include Cosmic Ray Physics (muons, neutrinos, EAS, underground physics), both for basic research…
We describe a hybrid Direct Simulation Monte Carlo (DSMC) code for simultaneously solving the collisional Boltzmann equation for gas and the collisionless Boltzmann equation for stars and dark matter for problems important to galaxy…
This article describes a new, fully adaptive Particle-Multiple-Mesh numerical simulation code developed primarily for simulations of small regions (such as a group of galaxies) in a cosmological context. It integrates the equations of…
We present a general method for accelerating by more than an order of magnitude the convolution of pixelated function on the sphere with a radially-symmetric kernel. Our method splits the kernel into a compact real-space, and a compact…
Starting from a high-level problem description in terms of partial differential equations using abstract tensor notation, the Chemora framework discretizes, optimizes, and generates complete high performance codes for a wide range of…
We describe a new parallel N-body code for cosmological simulations. The code is based on a work- and data sharing scheme, and is implemented within the Cray Research Corporation's CRAFT programming environment. Different data distribution…
The FLUKA Monte Carlo code, presently used in cosmic ray physics, contains packages to sample soft hadronic processes which are built according to the Dual Parton Model. This is a phenomenological model capable of reproducing many of the…
Understanding the universe is hampered by the elusiveness of its most common constituent, cold dark matter. Almost impossible to observe, dark matter can be studied effectively by means of simulation and there is probably no other research…
We describe a novel simulation method that eliminates the slowing-down problem in the Monte Carlo simulations of imaginary-time path integrals near the continuum limit. This method combines a stochastic blocking procedure with the multigrid…