Related papers: Multiscale Gaussian Random Fields for Cosmological…
This technical paper describes a software package that was designed to produce initial conditions for large cosmological simulations in the context of the Horizon collaboration. These tools generalize E. Bertschinger's Grafic1 software to…
We present a generic algorithm for generating Gaussian random initial conditions for cosmological simulations on periodic rectangular lattices. We show that imposing periodic boundary conditions on the real-space correlator and choosing…
We provide theoretical procedures and practical recipes to simulate non-Gaussian correlated, homogeneous random fields with prescribed marginal distributions and cross-correlation structure, either in a N-dimensional Cartesian space or on…
We discuss a new algorithm to generate multi-scale initial conditions with multiple levels of refinements for cosmological "zoom-in" simulations. The method uses an adaptive convolution of Gaussian white noise with a real space transfer…
We describe how to define an extremely large discrete realisation of a Gaussian white noise field that has a hierarchical structure and the property that the value of any part of the field can be computed quickly. Tiny subregions of such a…
The large-scale structure in cosmology is highly non-Gaussian at late times and small length scales, making it difficult to describe analytically. Parameter inference, data reconstruction, and data generation tasks in cosmology are greatly…
A computer code is described for the simulation of gravitational lensing data. The code incorporates adaptive mesh refinement in choosing which rays to shoot based on the requirements of the source size, location and surface brightness…
Reconstructing cosmological initial conditions (ICs) from late-time observations is a difficult task, which relies on the use of computationally expensive simulators alongside sophisticated statistical methods to navigate multi-million…
In this paper we will consider the problem of the numerical simulation of non-Gaussian, scalar random fields with a prescribed correlation structure provided either by a theoretical model or computed on a set of observational data.…
Gaussian processes are popular and flexible models for spatial, temporal, and functional data, but they are computationally infeasible for large datasets. We discuss Gaussian-process approximations that use basis functions at multiple…
In this contribution a broad overview of the methodologies of cosmological N-body simulations and a short introduction explaining the general idea behind such simulations is presented. After explaining how to set up the initial conditions…
We introduce statistical techniques required to handle complex computer models with potential applications to astronomy. Computer experiments play a critical role in almost all fields of scientific research and engineering. These computer…
Reconstructing astrophysical and cosmological fields from observations is challenging. It requires accounting for non-linear transformations, mixing of spatial structure, and noise. In contrast, forward simulators that map fields to…
Numerical simulations are becoming a more effective tool for conducting detailed investigations into the evolution of our universe. In this article, we show how the framework of numerical relativity can be used for studying cosmological…
We have developed an adaptive multigrid code for solving the Poisson equation in gravitational simulations. Finer rectangular subgrids are adaptively created in locations where the density exceeds a local level-dependent threshold. We…
The standard paradigm of cosmology assumes General Relativity (GR) is a valid theory for gravity at scales in which it has not been properly tested. Developing novel tests of GR and its alternatives is crucial if we want to give strength to…
To numerically evolve the full Einstein equations (or modifications thereof), simulations of cosmological spacetimes must rely on a particular formulation of the field equations combined with a specific gauge/frame choice. Yet truly…
There has been a lot of research interest in modified gravity theories which utilise the Vainshtein mechanism to recover standard general relativity in regions with high matter density, such as the Dvali-Gabadadze-Porrati and Galileon…
We consider problems related to initial meshing and adaptive mesh refinement for the electromagnetic simulation of various structures. The quality of the initial mesh and the performance of the adaptive refinement are of great importance…
We make a very large realisation of a Gaussian white noise field, called PANPHASIA, public by releasing software that computes this field. Panphasia is designed specifically for setting up Gaussian initial conditions for cosmological…