A Python code to compute the non-Gaussian (regular trispectrum and its shot noise) part of the analytic covariance matrix of the redshift-space galaxy power spectrum multipoles, using an FFTLog-based method proposed in Kobayashi 2023. The galaxy trispectrum is based on tree-level standard perturbation theory as described in Wadekar & Scoccimarro 2020, but with a slightly different galaxy bias expansion (see Appendix A of Kobayashi 2023).
The code computes the non-Gaussian covariance of the power spectrum monopole, quadrupole, hexadecapole, and their cross-covariance up to kmax ~ 0.4 h/Mpc.
The following packages need to be installed.
- numpy
- scipy
- sympy
You can find an example Jupyter notebook here.
- Yosuke Kobayashi (yosukekobayashi@arizona.edu)
- Yosuke Kobayashi, Fast computation of non-Gaussian covariance of redshift-space galaxy power spectrum multipoles (arXiv)