An adaptive mesh, massively-parallel, cosmological simulation code
Nyx code solves equations of compressible hydrodynamics on an adaptive grid hierarchy coupled with an N-body treatment of dark matter. The gas dynamics in Nyx uses a finite volume methodology on a set of 3-D Eulerian grids; dark matter is represented as discrete particles moving under the influence of gravity. Particles are evolved via a particle-mesh method, using Cloud-in-Cell deposition/interpolation scheme. Both baryonic and dark matter contribute to the gravitational field. In addition, Nyx includes physics needed to accurately model the intergalactic medium: in optically thin limit and assuming ionization equilibrium, the code calculates heating and cooling processes of the primordial-composition gas in an ionizing ultraviolet background radiation field. Additional physics capabilities are under development.
While Nyx can run on any Linux system in general, we particularly focus on supercomputer systems. Nyx is parallelized with MPI + X, where X can be OpenMP on multicore architectures and CUDA/HIP/DPC++ on hybrid CPU/GPU architectures. In the OpenMP regime, Nyx has been successfully run at parallel concurrency of up to 2,097,152 on NERSC's Cori-KNL. With Cuda implementation, it was run on up to 13,824 GPUs on OLCF's Summit.
More information on Nyx can be found at the main web page and the online documentation.
To compile the code we require C++11 compliant compilers that support MPI-2 or higher implementation. If threads or accelerators are used, we require OpenMP 4.5 or higher, Cuda 9 or higher, or HIP-Clang. To use Nyx, you also need AMReX.
For example, to compile the Lyman alpha (LyA) executable on Summit:
$ module load gcc/6.4.0 cuda/11.0.3
$ git clone https://github.com/AMReX-Codes/amrex.git
$ git clone https://github.com/AMReX-astro/Nyx.git
$ cd Nyx/Exec/LyA
$ make -j 12 USE_CUDA=TRUESee the Getting Started section for more information.
Please see CONTRIBUTING.md for details on how to contribute to AMReX.
Nyx outputs certain global diagnostics at each timestep and plot files at regular intervals, or at user-specified redshifts. Visualization packages VisIt, Paraview, yt, and Amrvis have built-in support for the AMReX file format used by Nyx.
In addition, Nyx interfaces with two post-processing suites, Reeber and Gimlet. Reeber uses topological methods to construct merge trees of scalar fields, which is in turn used to find halos. Gimlet computes a variety of quantities related to the Lyman-alpha forest science. These suites are fully MPI-parallel and can be run either in situ or in-transit, or with a combination of both (see Friesen et al. 2016).
Nyx Copyright (c) 2017, The Regents of the University of California, through Lawrence Berkeley National Laboratory (subject to receipt of any required approvals from the U.S. Dept. of Energy). All rights reserved.
Details of the license can be found in license.txt file.
If you have questions about your rights to use or distribute this software, please contact Berkeley Lab's Innovation & Partnerships Office at [email protected].
For questions, comments, suggestions, contact Jean Sexton ([email protected]) or Zarija Lukic ([email protected]).