Performance Characteristics of TreePM codes

Abstract

We present a detailed analysis of the error budget for the TreePM method for doing cosmological N-Body simulations. It is shown that the choice of filter for splitting the inverse square force into short and long range components suggested in Bagla (2002) is close to optimum. We show that the error in the long range component of the force contributes very little to the total error in force. Errors introduced by the tree approximation for the short range force are different from those for the inverse square force, and these errors dominate the total error in force. We calculate the distribution function for error in force for clustered and unclustered particle distributions. This gives an idea of the error in realistic situations for different choices of parameters of the TreePM algorithm. We test the code by simulating a few power law models and checking for scale invariance.

19 pages, 10 postscript figures

Cited by in corpus (30)

• The cosmological simulation code GADGET-2
• The FLAMINGO project: cosmological hydrodynamical simulations for large-scale structure and galaxy cluster surveys
• Comparing AMR and SPH Cosmological Simulations: I. Dark Matter & Adiabatic Simulations
• Large-scale dark matter simulations
• Grid Based Linear Neutrino Perturbations in Cosmological N-body Simulations
• SWIFT: A modern highly-parallel gravity and smoothed particle hydrodynamics solver for astrophysical and cosmological applications
• GRAPE-6A: A single-card GRAPE-6 for parallel PC-GRAPE cluster system
• HI as a Probe of the Large Scale Structure in the Post-Reionization Universe
• Fractal Dimension as a measure of the scale of Homogeneity
• Simulation techniques for cosmological simulations
• Numerical Simulations of the Lyman-alpha forest - A comparison of Gadget-2 and Enzo
• A fast recursive coordinate bisection tree for neighbour search and gravity
• Comments on the size of the simulation box in cosmological N-Body simulations
• Effects of the size of cosmological N-Body simulations on physical quantities -- I: Mass Function
• Lyman-α Forest and Cosmic Weak Lensing in a Warm Dark Matter Universe
• Cosmological simulation in tides: power spectra and halo shape responses, and shape assembly bias
• The FLAMINGO project: Baryon effects on the matter power spectrum
• Effects of the size of cosmological N-Body simulations on physical quantities - III: Skewness
• A Hybrid Fast Multipole Method for Cosmological N-body Simulations
• The Adaptive TreePM: An Adaptive Resolution Code for Cosmological N-body Simulations
• Power transfer in nonlinear gravitational clustering and asymptotic universality
• Gravitational collapse in an expanding background and the role of substructure II: Excess power at small scales and its effect of collapse of structures at larger scales
• A Modified TreePM Code
• Gravitational collapse in an expanding background and the role of substructure I: Planar collapse
• Gravitational Collapse in an Expanding Universe: Scaling Relations for Two-Dimensional Collapse Revisited
• Dynamic Zoom Simulations: a fast, adaptive algorithm for simulating lightcones
• The FLAMINGO simulations data release
• TreePM Method for Two-Dimensional Cosmological Simulations
• Challenges in Non-Linear Gravitational Clustering
• Summary of ICGC04 Cosmology Workshop