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,pk_numerical_nowiggle,pk_analytic_nowiggle};
 
enum source_extrapolation {extrap_zero,extrap_only_max,extrap_only_max_units,extrap_max_scaled,extrap_hmcode,extrap_user_defined};
 
enum hmcode_baryonic_feedback_model {hmcode_emu_dmonly, hmcode_owls_dmonly, hmcode_owls_ref, hmcode_owls_agn, hmcode_owls_dblim, hmcode_user_defined};
enum hmcode_version {hmcode_version_2015, hmcode_version_2020, hmcode_version_2020_unfitted, hmcode_version_2020_baryonic};
 
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; 
  enum hmcode_version hm_version; 
  double c_min;      
  double eta_0;      
  double z_infinity; 
  short has_pk_eq;   
  int nk_wiggle;     
  double log10T_heat_hmcode; 
  short has_pk_analytic_nowiggle; 
  short has_pk_numerical_nowiggle; 
 
 
  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 redundant 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 ln_tau_size_nl;  
  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;
  double ** ln_pk_l_extra;   
  double ** ddln_pk_l_extra; 
  double * ln_pk_l_an_extra; 
  double * ddln_pk_l_an_extra; 
  int * pk_l_nw_index;  
  double * ln_pk_l_nw_extra; 
  double * ddln_pk_l_nw_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;   
  short is_allocated; 
};
 
/********************************************************************************/
 
/* @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 primordial *ppm,
                         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_pk_analytic_nowiggle(
                                   struct precision *ppr,
                                   struct background *pba,
                                   struct primordial * ppm,
                                   struct fourier *pfo
                                   );
 
  int fourier_wnw_split(
                        struct precision *ppr,
                        struct background *pba,
                        struct primordial * ppm,
                        struct fourier *pfo
                        );
 
  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
                         );
 
#ifdef __cplusplus
}
#endif
 
#endif
/* @endcond */