The goal of large-scale structure (LSS) analysis is to extract cosmological information from the nonlinear matter density field.
Nearly all modern cosmological analyses are built upon two-point statistics that probe this field [1]}, [2]}.
Increasingly precise measurements of two-point correlators in modern galaxy surveys demand percent-level accuracy of theoretical models of these correlators.
This is especially true on smaller scales where the effect of survey sample variance is minimal and non-perturbative effects dominate.
However, models of two-point statistics require a tradeoff between the scales they access and the amount of theoretical control they have.
