spinosaurus computes two‑point correlation functions of galaxy shapes and densities using Lagrangian perturbation theory (LPT). Given cosmological input (e.g., matter power spectrum) and bias parameters, it returns shape–shape and density–shape power spectra (and related statistics) including the effects of long-wavelength displacement resummation. The code relies on numpy, scipy, and pyFFTW (ascl:2109.009) for efficient Fourier transforms and supports both a library interface and example notebooks for common use cases.

The anzu package offers two independent codes for hybrid Lagrangian bias models in large-scale structure. The first code measures the hybrid "basis functions"; the second takes measurements of these basis functions and constructs an emulator to obtain predictions from them at any cosmology (within the bounds of the training set). anzu is self-contained; given a set of N-body simulations used to build emulators, it measures the basis functions. Alternatively, given measurements of the basis functions, anzu should in principle be useful for constructing a custom emulator.

SkyLine generates mock line-intensity maps (both in 3D and 2D) in a lightcone from a halo catalog, accounting for the evolution of clustering and astrophysical properties, and observational effects such as spectral and angular resolutions, line-interlopers, and galactic foregrounds. Using a given astrophysical model for the luminosity of each line, the code paints the signal for each emitter and generates the map, adding coherently all contributions of interest. In addition, SkyLine can generate maps with the distribution of Luminous Red Galaxies and Emitting Line Galaxies.