SARABANDE: 3/4 Point Correlation Functions with Fast Fourier Transforms

We present a new python package sarabande for measuring 3 & 4 Point Correlation Functions (3/4 PCFs) in $\mathcal {O} (N_{\mathrm{g}}\log N_{\mathrm{g}})$ time using Fast Fourier Transforms (FFTs), with Ng the number of grid points used for the FFT. sarabande can measure both projected and full 3 and 4 PCFs on gridded 2D and 3D datasets. The general technique is to generate suitable angular basis functions on an underlying grid, radially bin these to create kernels, and convolve these kernels with the original gridded data to obtain expansion coefficients about every point simultaneously. These coefficients are then combined to give us the 3/4 PCF as expanded in our basis. We apply sarabande to simulations of the Interstellar Medium (ISM) to show the results and scaling of calculating both the full and projected 3/4 PCFs.