Massively Parallel Cosmological SPH Simulation Software - MP-Gadget.

Source code browser (may be slightly out-sync from current master branch)

This version of Gadget is derived from main P-Gadget / Gadget-2. It is the source code used to run the BlueTides simulation (http://bluetides-project.org).

The infrastructure is heavily reworked. As a summary:

• A better PM solver for long range force with Pencil FFT.
• A better Tree solver with faster threading and less redundant code.
• A better Domain decomposition that scales to half a million cores.
• A easier to use IO module with a Python binding.
• A more intuitive parameter file parser with schema and docstrings.
• A cleaner code base with less conditional compilation flags.

Physics models:

• Pressure Entropy SPH and Density Entropy SPH
• Radiation background in the expansion
• Massive neutrinos
• Dark energy
• ICs have species dependent density and velocity transfer functions
• Generic halo tracer particle seeding
• Various wind feedback and blackhole feedback models
• Various star formation criteria
• Primordial and metal cooling using updated recombination rates from the Sherwood simulation.
• Helium reionization
• Fluctuating UV background

First time users:

git clone https://github.com/MP-Gadget/MP-Gadget.git
cd MP-Gadget
make -j

We will need gsl. On HPC systems with the modules command, usually it can be loaded with

module load gsl

env | grep GSL  # check if GSL path is reasonable

On a common PC/Linux system, refer to your package vendor how to install gsl and gsl-devel.

If you wish to perform compile-time customisation (to, eg, change optimizations or use different compilers), you need an Options.mk file. The initial defaults are stored in Options.mk.example.

For other systems you should use the customised Options.mk file in the platform-options directory. For example, for Stampede 2 you should do:

cp platform-options/Options.mk.stampede2 Options.mk

For generic intel compiler based clusters, start with platform-options/Options.mk.icc

Compile-time options may be set in Options.mk. The remaining compile time options are generally only useful for development or debugging. All science options are set using a parameter file at runtime.

• DEBUG which enables various internal code consistency checks for debugging.
• VALGRIND which if set disables the internal memory allocator and allocates memory from the system. This is required for debugging memory allocation errors with valgrind of the address sanitizer.
• NO_ISEND_IRECV_IN_DOMAIN disables the use of asynchronous send and receive in our custom MPI_Alltoallv implementation, for buggy MPI libraries.
• NO_OPENMP_SPINLOCK uses the OpenMP default locking routines. These are often much slower than the default pthread spinlocks. However, they are necessary for Mac, which does not provide pthreads.

If compilation fails with errors related to the GSL, you may also need to set the GSL_INC or GSL_LIB variables in Options.mk to the filesystem path containing the GSL headers and libraries.

To run a N-Body sim, use IC files with no gas particles.

Now we are ready to build

make -j

In the end, we will have 2 binaries:

ls gadget/MP-Gadget genic/MP-GenIC

1. MP-Gadget is the main simulation program.
2. MP-GenIC is the initial condition generator.

Most options are configured at runtime with options in the config files. The meaning of these options are documented in the params.c files in the gadget/ and genic/ subdirectories.

Find examples in examples/.

• dm-only : Dark Matter only
• lya : Lyman Alpha only
• hydro : hydro
• small : hydro with low resolution

Control number of threads with OMP_NUM_THREADS. A good value is 10-20 threads.

A longer user guide in LaTeX can be found here: https://www.overleaf.com/read/kzksrgnzhtnh

The snapshot is in bigfile format. For data analysis in Python, use

pip install bigfile

Refer to https://github.com/rainwoodman/bigfile for usage.

Otherwise directly open the blocks with Fortran or C, noting the data-type information and attributes in header and attrs files (in plain text)

Cray updated their GLIBC to 2.22+ recently. A good move but it happens to be a buggy version of GLIBC: https://sourceware.org/bugzilla/show_bug.cgi?id=19590 causing non-existing symbols like _ZGVcN4v___log_finite. Adding -lmvec -lmvec_nonshared to GSL_LIBS works around the issue.

Bigfile is incorporated using git-subtree, in the depends/bigfile prefix. The command to update it (squash is currently mandatory) is:

git subtree pull --prefix depends/bigfile "https://github.com/rainwoodman/bigfile.git" master --squash

Gadget-2 was authored by Volker Springel. The original P-GADGET3 was maintained by Volker Springel

MP-Gadget is maintained by Yu Feng and Simeon Bird.

Contributors to MP-Gadget include:

Nicholas Battaglia, James Davies, Nishikanta Khandai, Yueying Ni, Karime Maamari, Chris Pederson and Lauren Anderson.

Pull requests should ideally be reviewed. Here are some links on how to conduct review:

https://smartbear.com/learn/code-review/best-practices-for-peer-code-review/ http://web.mit.edu/6.005/www/fa15/classes/04-code-review/

We never get around to write a proper code paper on MP-Gadget.

For usage of the code, here is a DOI for this repository that you can cite

It helps us to keep track of uses.

Most files are licensed under GPLv2+.

Please refer to the source files for details.

master branch status: