Multi-level adaptive particle mesh (MLAPM): a c code for cosmological simulation. We present a computer code written in c that is designed to simulate structure formation from collisionless matter. The code is purely grid-based and uses a recursively refined Cartesian grid to solve Poisson’s equation for the potential, rather than obtaining the potential from a Green’s function. Refinements can have arbitrary shapes and in practice closely follow the complex morphology of the density field that evolves. The time-step shortens by a factor of 2 with each successive refinement. Competing approaches to N-body simulation are discussed from the point of view of the basic theory of N-body simulation. It is argued that an appropriate choice of softening length ε is of great importance and that ε should be at all points an appropriate multiple of the local interparticle separation. Unlike tree and P3M codes, multigrid codes automatically satisfy this requirement. We show that at early times and low densities in cosmological simulations, ε needs to be significantly smaller relative to the interparticle separation than in virialized regions. Tests of the ability of the code’s Poisson solver to recover the gravitational fields of both virialized haloes and Zel’dovich waves are presented, as are tests of the code’s ability to reproduce analytic solutions for plane-wave evolution. The times required to conduct a ΛCDM cosmological simulation for various configurations are compared with the times required to complete the same simulation with the ART, AP3M and GADGET codes. The power spectra, halo mass functions and halo–halo correlation functions of simulations conducted with different codes are compared.
Keywords for this software
References in zbMATH (referenced in 5 articles )
Showing results 1 to 5 of 5.
- Bramas, Bérenger; Hassan, Muhammad; Stamm, Benjamin: An integral equation formulation of the (N)-body dielectric spheres problem. II: Complexity analysis (2021)
- Hassan, Muhammad; Stamm, Benjamin: A linear scaling in accuracy numerical method for computing the electrostatic forces in the (N)-body dielectric spheres problem (2021)
- Sousbie, Thierry; Colombi, Stéphane: \textttColDICE: A parallel Vlasov-Poisson solver using moving adaptive simplicial tessellation (2016)
- Miniati, Francesco; Colella, Phillip: Block structured adaptive mesh and time refinement for hybrid, hyperbolic (+ N)-body systems (2007)
- Fujimoto, Keizo; Machida, Shinobu: Electromagnetic full particle code with adaptive mesh refinement technique: application to the current sheet evolution (2006)