background_neutrino_horizon.h

 
#ifndef __BACKGROUND__
#define __BACKGROUND__
 
#include "common.h"
#include "quadrature.h"
#include "growTable.h"
#include "arrays.h"
#include "dei_rkck.h"
#include "parser.h"
 
enum equation_of_state {CLP,EDE};
 
 
enum varconst_dependence {varconst_none,varconst_instant};
 
enum vecback_format {short_info, normal_info, long_info};
 
enum interpolation_method {inter_normal, inter_closeby};
 
struct background
{
 
  double H0; 
  double h;  
  double Omega0_g; 
  double T_cmb;    
  double Omega0_b; 
  double Omega0_ur; 
  double Omega0_cdm;      
  double Omega0_idm; 
  double Omega0_idr; 
  double T_idr;      
  double Omega0_dcdmdr;   
  double Omega_ini_dcdm;  
  double Gamma_dcdm;      
  double tau_dcdm;
 
  int N_ncdm;                            
  /* the following parameters help to define tabulated ncdm p-s-d passed in file */
  char * ncdm_psd_files;                 
  int * got_files;                       
  /* the following parameters help to define the analytical ncdm phase space distributions (p-s-d) */
  double * ncdm_psd_parameters;          
  double * M_ncdm;                       
  double * m_ncdm_in_eV;                 
  double * Omega0_ncdm, Omega0_ncdm_tot; 
  double * T_ncdm,T_ncdm_default;        
  double * ksi_ncdm, ksi_ncdm_default;   
  double * deg_ncdm, deg_ncdm_default;    
  int * ncdm_input_q_size; 
  double * ncdm_qmax;      
  double Omega0_k;         
  double Omega0_lambda;    
  double Omega0_fld;       
  double Omega0_scf;       
  short use_ppf; 
  double c_gamma_over_c_fld; 
  enum equation_of_state fluid_equation_of_state; 
  double w0_fld;   
  double wa_fld;   
  double cs2_fld;  
  double Omega_EDE;        
  double * scf_parameters; 
  short attractor_ic_scf;  
  int scf_tuning_index;    
  double phi_ini_scf;      
  double phi_prime_ini_scf;
  int scf_parameters_size; 
  double varconst_alpha; 
  double varconst_me; 
  enum varconst_dependence varconst_dep; 
  double varconst_transition_redshift; 
 
 
 
  double age; 
  double conformal_age; 
  double K; 
  int sgnK; 
  double Neff; 
  double Omega0_dcdm; 
  double Omega0_dr; 
  double Omega0_m;  
  double Omega0_r;  
  double Omega0_de; 
  double Omega0_nfsm; 
  double a_eq;      
  double H_eq;      
  double z_eq;      
  double tau_eq;    
 
 
 
  int index_bg_a;             
  int index_bg_H;             
  int index_bg_H_prime;       
  /* end of vector in short format, now quantities in normal format */
 
  int index_bg_rho_g;         
  int index_bg_rho_b;         
  int index_bg_rho_cdm;       
  int index_bg_rho_idm;       
  int index_bg_rho_lambda;    
  int index_bg_rho_fld;       
  int index_bg_w_fld;         
  int index_bg_rho_idr;       
  int index_bg_rho_ur;        
  int index_bg_rho_dcdm;      
  int index_bg_rho_dr;        
  int index_bg_phi_scf;       
  int index_bg_phi_prime_scf; 
  int index_bg_V_scf;         
  int index_bg_dV_scf;        
  int index_bg_ddV_scf;       
  int index_bg_rho_scf;       
  int index_bg_p_scf;         
  int index_bg_p_prime_scf;         
  int index_bg_rho_ncdm1;     
  int index_bg_p_ncdm1;       
  int index_bg_pseudo_p_ncdm1;
  int index_bg_rho_tot;       
  int index_bg_p_tot;         
  int index_bg_p_tot_prime;   
  int index_bg_Omega_r;       
  /* end of vector in normal format, now quantities in long format */
 
  int index_bg_rho_crit;      
  int index_bg_Omega_m;       
  int index_bg_conf_distance; 
  int index_bg_ang_distance;  
  int index_bg_lum_distance;  
  int index_bg_time;          
  int index_bg_rs;            
  int index_bg_nu;            
  int index_bg_D;             
  int index_bg_f;             
  int index_bg_varc_alpha;    
  int index_bg_varc_me;      
  int bg_size_short;  
  int bg_size_normal; 
  int bg_size;        
 
 
 
  int bt_size;               
  double * loga_table;       
  double * tau_table;        
  double * z_table;          
  double * background_table; 
 
 
 
  double * d2tau_dz2_table; 
  double * d2z_dtau2_table; 
  double * d2background_dloga2_table; 
 
 
 
  int index_bi_rho_dcdm;
  int index_bi_rho_dr;  
  int index_bi_rho_fld; 
  int index_bi_phi_scf;       
  int index_bi_phi_prime_scf; 
  int index_bi_time;    
  int index_bi_rs;      
  int index_bi_nu;      
  int index_bi_tau;     
  int index_bi_D;       
  int index_bi_D_prime; 
  int bi_B_size;        
  int bi_size;          
 
 
  short has_cdm;       
  short has_idm;       
  short has_dcdm;      
  short has_dr;        
  short has_scf;       
  short has_ncdm;      
  short has_lambda;    
  short has_fld;       
  short has_ur;        
  short has_idr;       
  short has_curvature; 
  short has_varconst;  
 
 
 
  int * ncdm_quadrature_strategy; 
  double ** q_ncdm_bg;  
  double ** w_ncdm_bg;  
  double ** q_ncdm;     
  double ** w_ncdm;     
  double ** dlnf0_dlnq_ncdm; 
  int * q_size_ncdm_bg; 
  int * q_size_ncdm;    
  double * factor_ncdm; 
 
 
  short shooting_failed;  
  ErrorMsg shooting_error; 
  short background_verbose; 
  ErrorMsg error_message; 
  short is_allocated; 
};
 
 
struct background_parameters_and_workspace {
 
  /* structures containing fixed input parameters (indices, ...) */
  struct background * pba;
 
  /* workspace */
  double * pvecback;
 
};
 
struct background_parameters_for_distributions {
 
  /* structures containing fixed input parameters (indices, ...) */
  struct background * pba;
 
  /* Additional parameters */
 
  /* Index of current distribution function */
  int n_ncdm;
 
  /* Used for interpolating in file of tabulated p-s-d: */
  int tablesize;
  double *q;
  double *f0;
  double *d2f0;
  int last_index;
 
};
 
/**************************************************************/
/* @cond INCLUDE_WITH_DOXYGEN */
/*
 * Boilerplate for C++
 */
#ifdef __cplusplus
extern "C" {
#endif
 
  int background_at_z(
                      struct background *pba,
                      double a_rel,
                      enum vecback_format return_format,
                      enum interpolation_method inter_mode,
                      int * last_index,
                      double * pvecback
                      );
 
  int background_at_tau(
                        struct background *pba,
                        double tau,
                        enum vecback_format return_format,
                        enum interpolation_method inter_mode,
                        int * last_index,
                        double * pvecback
                        );
 
  int background_tau_of_z(
                          struct background *pba,
                          double z,
                          double * tau
                          );
 
  int background_z_of_tau(
                          struct background *pba,
                          double tau,
                          double * z
                          );
 
  int background_functions(
                           struct background *pba,
                           double a_rel,
                           double * pvecback_B,
                           enum vecback_format return_format,
                           double * pvecback
                           );
 
  int background_w_fld(
                       struct background * pba,
                       double a,
                       double * w_fld,
                       double * dw_over_da_fld,
                       double * integral_fld);
 
  int background_varconst_of_z(
                               struct background* pba,
                               double z,
                               double* alpha,
                               double* me
                               );
 
  int background_init(
                      struct precision *ppr,
                      struct background *pba
                      );
 
  int background_free(
                      struct background *pba
                      );
 
  int background_free_noinput(
                              struct background *pba
                              );
 
  int background_free_input(
                            struct background *pba
                            );
 
  int background_indices(
                         struct background *pba
                         );
 
  int background_ncdm_distribution(
                                   void *pba,
                                   double q,
                                   double * f0
                                   );
 
  int background_ncdm_test_function(
                                    void *pba,
                                    double q,
                                    double * test
                                    );
 
  int background_ncdm_init(
                           struct precision *ppr,
                           struct background *pba
                           );
 
  int background_ncdm_momenta(
                              double * qvec,
                              double * wvec,
                              int qsize,
                              double M,
                              double factor,
                              double z,
                              double * n,
                              double * rho,
                              double * p,
                              double * drho_dM,
                              double * pseudo_p
                              );
 
  int background_ncdm_M_from_Omega(
                                   struct precision *ppr,
                                   struct background *pba,
                                   int species
                                   );
 
  int background_checks(
                        struct precision * ppr,
                        struct background *pba
                        );
 
  int background_solve(
                       struct precision *ppr,
                       struct background *pba
                       );
 
  int background_initial_conditions(
                                    struct precision *ppr,
                                    struct background *pba,
                                    double * pvecback,
                                    double * pvecback_integration,
                                    double * loga_ini
                                    );
 
  int background_find_equality(
                               struct precision *ppr,
                               struct background *pba
                               );
 
 
  int background_output_titles(struct background * pba,
                               char titles[_MAXTITLESTRINGLENGTH_]
                               );
 
  int background_output_data(
                             struct background *pba,
                             int number_of_titles,
                             double *data);
 
  int background_derivs(
                        double loga,
                        double * y,
                        double * dy,
                        void * parameters_and_workspace,
                        ErrorMsg error_message
                        );
 
  int background_sources(
                         double loga,
                         double * y,
                         double * dy,
                         int index_loga,
                         void * parameters_and_workspace,
                         ErrorMsg error_message
                         );
 
  int background_timescale(
                           double loga,
                           void * parameters_and_workspace,
                           double * timescale,
                           ErrorMsg error_message
                           );
 
  int background_output_budget(
                               struct background* pba
                               );
 
  double V_scf(
               struct background *pba,
               double phi
               );
 
  double dV_scf(
                struct background *pba,
                double phi
                );
 
  double ddV_scf(
                 struct background *pba,
                 double phi
                 );
 
  double Q_scf(
               struct background *pba,
               double phi,
               double phi_prime
               );
 
#ifdef __cplusplus
}
#endif
 
/**************************************************************/
 
 
#define _Mpc_over_m_ 3.085677581282e22  
/* remark: CAMB uses 3.085678e22: good to know if you want to compare  with high accuracy */
 
#define _Gyr_over_Mpc_ 3.06601394e2 
#define _c_ 2.99792458e8            
#define _G_ 6.67428e-11             
#define _eV_ 1.602176487e-19        
/* parameters entering in Stefan-Boltzmann constant sigma_B */
#define _k_B_ 1.3806504e-23
#define _h_P_ 6.62606896e-34
/* remark: sigma_B = 2 pi^5 k_B^4 / (15h^3c^2) = 5.670400e-8
   = Stefan-Boltzmann constant in W/m^2/K^4 = Kg/K^4/s^3 */
 
 
 
#define _h_BIG_ 1.5            
#define _h_SMALL_ 0.3         
#define _omegab_BIG_ 0.039    
#define _omegab_SMALL_ 0.005  
 
 
#define _SCALE_BACK_ 0.1  
#define _PSD_DERIVATIVE_EXP_MIN_ -30 
#define _PSD_DERIVATIVE_EXP_MAX_ 2  
#define _zeta3_ 1.2020569031595942853997381615114499907649862923404988817922 
#define _zeta5_ 1.0369277551433699263313654864570341680570809195019128119741 
 
 
#endif
/* @endcond */