halomod.tools.corr_to_power_ogata

halomod.tools.corr_to_power_ogata(corr, r, k, h=0.005, power_pos=(True, False), atol=1e-15, rtol=0.001)[source]

Convert an isotropic 3D correlation function to a power spectrum.

Uses Ogata’s method (Ogata 2005) for Hankel Transforms in 3D.

Parameters
  • corr (np.ndarray) – The correlation function to convert as a function of r.

  • r (np.ndarray) – Array of same length as corr, giving the real-space co-ordinates.

  • k (np.ndarray) – The fourier-space co-ordinates to which to transform.

  • n (int, optional) – The number of subdivisions in the integral.

  • h (int, optional) – Controls the spacing of the intervals (note the intervals are not equispaced). Smaller numbers give smaller intervals.

Notes

See the hankel documentation for details on the implementation here. This particular function is restricted to a spherical transform.