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.