PolyBin3D is a Python code that estimates the binned power spectrum and bispectrum for 3D fields such as the distributions of matter and galaxies. It is a sister code to PolyBin, which computes the polyspectra of data on the two-sphere. For each statistic, two estimators are available: the standard (ideal) estimators, which do not take into account the mask, and window-deconvolved estimators. In the second case, we require computation of a Fisher matrix; this depends on binning and the mask, but does not need to be recomputed for each new simulation. The code also supports GPU acceleration using JAX.
Link to source: https://github.com/oliverphilcox/PolyBin3D
Related papers: