fourier.h

Go to the documentation of this file.

 
#include "primordial.h"
#include "trigonometric_integrals.h"
 
#ifndef __FOURIER__
#define __FOURIER__
 
#define _M_EV_TOO_BIG_FOR_HALOFIT_ 10. 
#define _M_SUN_ 1.98847e30 
#define _MAX_NUM_EXTRAPOLATION_ 100000
 
enum non_linear_method {nl_none,nl_halofit,nl_HMcode};
enum pk_outputs {pk_linear,pk_nonlinear};
 
enum source_extrapolation {extrap_zero,extrap_only_max,extrap_only_max_units,extrap_max_scaled,extrap_hmcode,extrap_user_defined};
 
enum halofit_integral_type {halofit_integral_one, halofit_integral_two, halofit_integral_three};
 
enum hmcode_baryonic_feedback_model {nl_emu_dmonly, nl_owls_dmonly, nl_owls_ref, nl_owls_agn, nl_owls_dblim, nl_user_defined};
enum out_sigmas {out_sigma,out_sigma_prime,out_sigma_disp};
 
struct fourier {
 
 
  enum non_linear_method method; 
  enum source_extrapolation extrapolation_method; 
  enum hmcode_baryonic_feedback_model feedback; 
  double c_min;      
  double eta_0;      
  double z_infinity; 
  short has_pk_eq;  
 
 
  int index_md_scalars; 
  int ic_size;         
  int ic_ic_size;      
  short * is_non_zero; 
 
 
  short has_pk_m;  
  short has_pk_cb; 
  int index_pk_m;  
  int index_pk_cb; 
  /* and two redundent but useful indices: */
 
  int index_pk_total;      
  int index_pk_cluster;    
  int pk_size;     
 
 
  short has_pk_matter; 
  int k_size;      
  int k_size_pk;   
  double * k;      
  double * ln_k;   
  double * ln_tau;     
  int ln_tau_size;     
  int index_ln_tau_pk; 
  double ** ln_pk_ic_l;   
  double ** ddln_pk_ic_l; 
  double ** ln_pk_l;   
  double ** ddln_pk_l; 
  double ** ln_pk_nl;   
  double ** ddln_pk_nl; 
  double * sigma8;   
 
 
  int k_size_extra;
  int tau_size;    
  double * tau;    
  double ** nl_corr_density;   
  double ** k_nl;              
  int index_tau_min_nl;        
 
 
  int index_pk_eq_w;                
  int index_pk_eq_Omega_m;          
  int pk_eq_size;                   
  int pk_eq_tau_size;               
  double * pk_eq_tau;               
  double * pk_eq_w_and_Omega;       
  double * pk_eq_ddw_and_ddOmega;   
 
 
  short fourier_verbose;    
  ErrorMsg error_message;   
};
 
struct fourier_workspace {
 
 
  double * rtab; 
  double * stab; 
  double * ddstab; 
  double * growtable;
  double * ztable;
  double * tautable;
 
  double ** sigma_8;
  double ** sigma_disp;
  double ** sigma_disp_100;
  double ** sigma_prime;
 
  double dark_energy_correction; 
 
};
 
/********************************************************************************/
 
/* @cond INCLUDE_WITH_DOXYGEN */
/*
 * Boilerplate for C++
 */
#ifdef __cplusplus
extern "C" {
#endif
 
  /* external functions (meant to be called from other modules) */
 
  int fourier_pk_at_z(
                      struct background * pba,
                      struct fourier *pfo,
                      enum linear_or_logarithmic mode,
                      enum pk_outputs pk_output,
                      double z,
                      int index_pk,
                      double * out_pk,
                      double * out_pk_ic
                      );
 
  int fourier_pks_at_z(
                       struct background * pba,
                       struct fourier *pfo,
                       enum linear_or_logarithmic mode,
                       enum pk_outputs pk_output,
                       double z,
                       double * out_pk,
                       double * out_pk_ic,
                       double * out_pk_cb,
                       double * out_pk_cb_ic
                       );
 
  int fourier_pk_at_k_and_z(
                            struct background * pba,
                            struct primordial * ppm,
                            struct fourier *pfo,
                            enum pk_outputs pk_output,
                            double k,
                            double z,
                            int index_pk,
                            double * out_pk,
                            double * out_pk_ic
                            );
 
  int fourier_pks_at_k_and_z(
                             struct background * pba,
                             struct primordial * ppm,
                             struct fourier *pfo,
                             enum pk_outputs pk_output,
                             double k,
                             double z,
                             double * out_pk,
                             double * out_pk_ic,
                             double * out_pk_cb,
                             double * out_pk_cb_ic
                             );
 
  int fourier_pks_at_kvec_and_zvec(
                                   struct background * pba,
                                   struct fourier * pfo,
                                   enum pk_outputs pk_output,
                                   double * kvec,
                                   int kvec_size,
                                   double * zvec,
                                   int zvec_size,
                                   double * out_pk,
                                   double * out_pk_cb
                                   );
 
  int fourier_sigmas_at_z(
                          struct precision * ppr,
                          struct background * pba,
                          struct fourier * pfo,
                          double R,
                          double z,
                          int index_pk,
                          enum out_sigmas sigma_output,
                          double * result
                          );
 
  int fourier_pk_tilt_at_k_and_z(
                                 struct background * pba,
                                 struct primordial * ppm,
                                 struct fourier * pfo,
                                 enum pk_outputs pk_output,
                                 double k,
                                 double z,
                                 int index_pk,
                                 double * pk_tilt
                                 );
 
  int fourier_k_nl_at_z(
                        struct background *pba,
                        struct fourier * pfo,
                        double z,
                        double * k_nl,
                        double * k_nl_cb
                        );
 
  /* internal functions */
 
  int fourier_init(
                   struct precision *ppr,
                   struct background *pba,
                   struct thermodynamics *pth,
                   struct perturbations *ppt,
                   struct primordial *ppm,
                   struct fourier *pfo
                   );
 
  int fourier_free(
                   struct fourier *pfo
                   );
 
  int fourier_indices(
                      struct precision *ppr,
                      struct background *pba,
                      struct perturbations * ppt,
                      struct primordial * ppm,
                      struct fourier * pfo
                      );
 
  int fourier_get_k_list(
                         struct precision *ppr,
                         struct perturbations * ppt,
                         struct fourier * pfo
                         );
 
  int fourier_get_tau_list(
                           struct perturbations * ppt,
                           struct fourier * pfo
                           );
 
  int fourier_get_source(
                         struct background * pba,
                         struct perturbations * ppt,
                         struct fourier * pfo,
                         int index_k,
                         int index_ic,
                         int index_tp,
                         int index_tau,
                         double ** sources,
                         double * source);
 
  int fourier_pk_linear(
                        struct background *pba,
                        struct perturbations *ppt,
                        struct primordial *ppm,
                        struct fourier *pfo,
                        int index_pk,
                        int index_tau,
                        int k_size,
                        double * lnpk,
                        double * lnpk_ic
                        );
 
  int fourier_sigmas(
                     struct fourier * pfo,
                     double R,
                     double *lnpk_l,
                     double *ddlnpk_l,
                     int k_size,
                     double k_per_decade,
                     enum out_sigmas sigma_output,
                     double * result
                     );
 
  int fourier_sigma_at_z(
                         struct background * pba,
                         struct fourier * pfo,
                         double R,
                         double z,
                         int index_pk,
                         double k_per_decade,
                         double * result
                         );
 
  int fourier_halofit(
                      struct precision *ppr,
                      struct background *pba,
                      struct perturbations *ppt,
                      struct primordial *ppm,
                      struct fourier *pfo,
                      int index_pk,
                      double tau,
                      double *pk_nl,
                      double *lnpk_l,
                      double *ddlnpk_l,
                      double *k_nl,
                      short * halofit_found_k_max
                      );
 
  int fourier_halofit_integrate(
                                struct fourier *pfo,
                                double * integrand_array,
                                int integrand_size,
                                int ia_size,
                                int index_ia_k,
                                int index_ia_pk,
                                int index_ia_sum,
                                int index_ia_ddsum,
                                double R,
                                enum halofit_integral_type type,
                                double * sum
                                );
 
  int fourier_hmcode(
                     struct precision *ppr,
                     struct background *pba,
                     struct perturbations *ppt,
                     struct primordial *ppm,
                     struct fourier *pfo,
                     int index_pk,
                     int index_tau,
                     double tau,
                     double *pk_nl,
                     double **lnpk_l,
                     double **ddlnpk_l,
                     double *k_nl,
                     short * halofit_found_k_max,
                     struct fourier_workspace * pnw
                     );
 
  int fourier_hmcode_workspace_init(
                                    struct precision *ppr,
                                    struct background *pba,
                                    struct fourier *pfo,
                                    struct fourier_workspace * pnw
                                    );
 
  int fourier_hmcode_workspace_free(
                                    struct fourier *pfo,
                                    struct fourier_workspace * pnw
                                    );
 
  int fourier_hmcode_dark_energy_correction(
                                            struct precision *ppr,
                                            struct background *pba,
                                            struct fourier *pfo,
                                            struct fourier_workspace * pnw
                                            );
 
  int fourier_hmcode_baryonic_feedback(
                                       struct fourier *pfo
                                       );
 
 
  int fourier_hmcode_fill_sigtab(
                                 struct precision *ppr,
                                 struct background * pba,
                                 struct perturbations *ppt,
                                 struct primordial * ppm,
                                 struct fourier * pfo,
                                 int index_tau,
                                 double *lnpk_l,
                                 double *ddlnpk_l,
                                 struct fourier_workspace * pnw
                                 );
 
  int fourier_hmcode_fill_growtab(
                                  struct precision *ppr,
                                  struct background * pba,
                                  struct fourier * pfo,
                                  struct fourier_workspace * pnw
                                  );
 
  int fourier_hmcode_growint(
                             struct precision *ppr,
                             struct background * pba,
                             struct fourier * pfo,
                             double a,
                             double w,
                             double wa,
                             double * growth
                             );
 
  int fourier_hmcode_window_nfw(
                                struct fourier * pfo,
                                double k,
                                double rv,
                                double c,
                                double *window_nfw
                                );
 
  int fourier_hmcode_halomassfunction(
                                      double nu,
                                      double *hmf
                                      );
 
  int fourier_hmcode_sigma8_at_z(
                                 struct background *pba,
                                 struct fourier * pfo,
                                 double z,
                                 double * sigma_8,
                                 double * sigma_8_cb,
                                 struct fourier_workspace * pnw
                                 );
 
  int fourier_hmcode_sigmadisp_at_z(
                                    struct background *pba,
                                    struct fourier * pfo,
                                    double z,
                                    double * sigma_disp,
                                    double * sigma_disp_cb,
                                    struct fourier_workspace * pnw
                                    );
 
  int fourier_hmcode_sigmadisp100_at_z(
                                       struct background *pba,
                                       struct fourier * pfo,
                                       double z,
                                       double * sigma_disp_100,
                                       double * sigma_disp_100_cb,
                                       struct fourier_workspace * pnw
                                       );
 
  int fourier_hmcode_sigmaprime_at_z(
                                     struct background *pba,
                                     struct fourier * pfo,
                                     double z,
                                     double * sigma_prime,
                                     double * sigma_prime_cb,
                                     struct fourier_workspace * pnw
                                     );
 
#ifdef __cplusplus
}
#endif
 
#endif
/* @endcond */