core_power.h

Go to the documentation of this file.

#pragma once
#include "precision.hpp"
#include <gsl/gsl_spline.h>

/*******************/
class Cosmo_Param; 
class Sim_Param; 
template <typename T, size_t N> class Data_Vec; 

void norm_pwr(Cosmo_Param& cosmo);

FTYPE_t norm_growth_factor(const Cosmo_Param& cosmo);

FTYPE_t growth_factor(FTYPE_t a, const Cosmo_Param& cosmo);

FTYPE_t growth_rate(FTYPE_t a, const Cosmo_Param& cosmo);

FTYPE_t growth_change(FTYPE_t a, const Cosmo_Param& cosmo);

FTYPE_t Omega_lambda(FTYPE_t a, const Cosmo_Param& cosmo);

FTYPE_t lin_pow_spec(FTYPE_t a, FTYPE_t k, const Cosmo_Param& cosmo);

FTYPE_t non_lin_pow_spec(FTYPE_t a, FTYPE_t k, const Cosmo_Param& cosmo);

class Interp_obj
{
public:
    Interp_obj(): is_init(false) {}
    ~Interp_obj();
    double operator()(double x) const;
    template <typename T, size_t N>
    void init(const Data_Vec<T, N>& data);
    double x_min, x_max;

private:
    bool is_init;
    gsl_spline* spline;
    gsl_interp_accel* acc;
};

template <typename T, size_t N>
class Extrap_Pk : public Interp_obj
{
public:
    Extrap_Pk(const Data_Vec<T, N>& data, const Sim_Param& sim);
    Extrap_Pk(const Data_Vec<T, N>& data, const Sim_Param& sim, const size_t m_l, const size_t n_u);
    Extrap_Pk(const Data_Vec<T, N>& data, const Sim_Param& sim, const size_t m_l, const size_t n_l,
              const size_t m_u, const size_t n_u);
    double operator()(double k) const;

    void fit_lin(const Data_Vec<T, N>& data, const size_t m, const size_t n, double& A);
    void fit_power_law(const Data_Vec<T, N>& data, const size_t m, const size_t n, double& A, double& n_s);

    double A_low; 
    const Cosmo_Param& cosmo;
    double A_up, n_s; 
    T k_min, k_max; 
};

template <typename T, size_t N>
class Extrap_Pk_Nl : public Extrap_Pk<T, N>
{
public:
    Extrap_Pk_Nl(const Data_Vec<T, N>& data, const Sim_Param &sim, T A_nl, T a_eff);
    const T A_nl, a_eff, k_split;
    double operator()(double k) const;
};

template<class P>
void gen_corr_func_binned_gsl_qawf(const Sim_Param &sim, const P& P_k, Data_Vec<FTYPE_t, 2>& corr_func_binned);

void gen_corr_func_binned_gsl_qawf_lin(const Sim_Param &sim, FTYPE_t a, Data_Vec<FTYPE_t, 2>& corr_func_binned);

void gen_corr_func_binned_gsl_qawf_nl(const Sim_Param &sim, FTYPE_t a, Data_Vec<FTYPE_t, 2>& corr_func_binned);

template<class P>
void gen_sigma_binned_gsl_qawf(const Sim_Param &sim, const P& P_k, Data_Vec<FTYPE_t, 2>& sigma_binned);

void gen_sigma_func_binned_gsl_qawf_lin(const Sim_Param &sim, FTYPE_t a, Data_Vec<FTYPE_t, 2>& sigma_binned);

void gen_sigma_func_binned_gsl_qawf_nl(const Sim_Param &sim, FTYPE_t a, Data_Vec<FTYPE_t, 2>& sigma_binned);