PySCo (Python Simulations for Cosmology) runs cosmological N-body simulations across various cosmological models, including Newtonian gravity, MOND gravity (quasi-linear formulation), and parametrized gravity (scale independent). It is a multi-threaded Particle-Mesh code and uses multigrid or fast Fourier transform (FFT) methods. PySCo includes an initial condition generator and power spectrum estimation, and computes the background and growth of density perturbations.

Magrathea-Pathfinder propagates photons within cosmological simulations to construct observables. This high-performance framework uses a 3D Adaptive-Mesh Refinement and is built on top of the MAGRATHEA metalibrary (ascl:2203.023).