primordial.h

Go to the documentation of this file.

 
#ifndef __PRIMORDIAL__
#define __PRIMORDIAL__
 
#include "perturbations.h"
 
enum primordial_spectrum_type {
                               analytic_Pk,
                               two_scales,
                               inflation_V,
                               inflation_H,
                               inflation_V_end,
                               external_Pk
};
 
enum linear_or_logarithmic {
                            linear,
                            logarithmic
};
 
enum potential_shape {
                      polynomial,
                      natural,
                      higgs_inflation
};
 
enum target_quantity {
                      _aH_,
                      _phi_,
                      _end_inflation_,
                      _a_
};
 
enum integration_direction {
                            backward,
                            forward
};
 
enum time_definition {
                      conformal,
                      proper
};
 
enum phi_pivot_methods {
                        N_star,
                        ln_aH_ratio,
                        ln_aH_ratio_auto
};
 
enum inflation_module_behavior {
                                numerical,
                                analytical
};
 
struct primordial {
 
 
  double k_pivot; 
  int has_k_max_for_primordial_pk;
  double k_max_for_primordial_pk; 
  enum primordial_spectrum_type primordial_spec_type; 
  /* - parameters describing the case primordial_spec_type = analytic_Pk : amplitudes, tilts, runnings, cross-correlations, ... */
 
  double A_s;  
  double n_s;  
  double alpha_s; 
  double beta_s;  
  double r;    
  double n_t;  
  double alpha_t; 
  double f_bi;  
  double n_bi;  
  double alpha_bi; 
  double f_cdi;  
  double n_cdi;  
  double alpha_cdi; 
  double f_nid;  
  double n_nid;  
  double alpha_nid; 
  double f_niv;  
  double n_niv;  
  double alpha_niv; 
  double c_ad_bi; 
  double n_ad_bi; 
  double alpha_ad_bi; 
  double c_ad_cdi; 
  double n_ad_cdi; 
  double alpha_ad_cdi; 
  double c_ad_nid; 
  double n_ad_nid; 
  double alpha_ad_nid; 
  double c_ad_niv; 
  double n_ad_niv; 
  double alpha_ad_niv; 
  double c_bi_cdi; 
  double n_bi_cdi; 
  double alpha_bi_cdi; 
  double c_bi_nid; 
  double n_bi_nid; 
  double alpha_bi_nid; 
  double c_bi_niv; 
  double n_bi_niv; 
  double alpha_bi_niv; 
  double c_cdi_nid; 
  double n_cdi_nid; 
  double alpha_cdi_nid; 
  double c_cdi_niv; 
  double n_cdi_niv; 
  double alpha_cdi_niv; 
  double c_nid_niv; 
  double n_nid_niv; 
  double alpha_nid_niv; 
  enum potential_shape potential;
 
  double V0;    
  double V1;    
  double V2;    
  double V3;    
  double V4;    
  /* parameters describing the case primordial_spec_type = inflation_H */
 
  double H0;    
  double H1;    
  double H2;    
  double H3;    
  double H4;    
  /* parameters describing inflation_V_end */
 
  double phi_end;   
  enum phi_pivot_methods phi_pivot_method; 
  double phi_pivot_target; 
  /* behavior of the inflation module */
  enum inflation_module_behavior behavior; 
  char*  command;  
  double custom1;  
  double custom2;  
  double custom3;  
  double custom4;  
  double custom5;  
  double custom6;  
  double custom7;  
  double custom8;  
  double custom9;  
  double custom10; 
 
 
  int md_size;      
  int * ic_size;    
  int * ic_ic_size; 
  int lnk_size;    
  double * lnk;    
  double ** lnpk;  
  double ** ddlnpk; 
  short ** is_non_zero; 
 
 
  double ** amplitude; 
  double ** tilt;      
  double ** running;   
 
 
  int index_in_a;       
  int index_in_phi;     
  int index_in_dphi;    
  int index_in_ksi_re;  
  int index_in_ksi_im;  
  int index_in_dksi_re; 
  int index_in_dksi_im; 
  int index_in_ah_re;   
  int index_in_ah_im;   
  int index_in_dah_re;  
  int index_in_dah_im;  
  int in_bg_size;       
  int in_size;          
 
 
  double phi_pivot;      
  double phi_min;        
  double phi_max;        
  double phi_stop;       
 
 
  short primordial_verbose; 
 
  ErrorMsg error_message; 
};
 
struct primordial_inflation_parameters_and_workspace {
 
  struct primordial * ppm;
  double N;
  double a2;
 
  double V;
  double dV;
  double ddV;
  double aH;
 
  double H;
  double dH;
  double ddH;
  double dddH;
 
  double zpp_over_z;
  double app_over_a;
 
  double k;
 
  enum integration_direction integrate;
  enum time_definition time;
 
};
 
 
/*************************************************************************************************************/
/* @cond INCLUDE_WITH_DOXYGEN */
/*
 * Boilerplate for C++
 */
#ifdef __cplusplus
extern "C" {
#endif
 
  int primordial_spectrum_at_k(
                               struct primordial * ppm,
                               int index_md,
                               enum linear_or_logarithmic mode,
                               double k,
                               double * pk
                               );
 
  int primordial_init(
                      struct precision  * ppr,
                      struct perturbations   * ppt,
                      struct primordial * ppm
                      );
 
  int primordial_free(
                      struct primordial * ppm
                      );
 
  int primordial_indices(
                         struct perturbations   * ppt,
                         struct primordial * ppm
                         );
 
  int primordial_get_lnk_list(
                              struct primordial * ppm,
                              double kmin,
                              double kmax,
                              double k_per_decade
                              );
 
  int primordial_analytic_spectrum_init(
                                        struct perturbations   * ppt,
                                        struct primordial * ppm
                                        );
 
  int primordial_analytic_spectrum(
                                   struct primordial * ppm,
                                   int index_md,
                                   int index_ic1_ic2,
                                   double k,
                                   double * pk
                                   );
 
  int primordial_inflation_potential(
                                     struct primordial * ppm,
                                     double phi,
                                     double * V,
                                     double * dV,
                                     double * ddV
                                     );
 
  int primordial_inflation_hubble(
                                  struct primordial * ppm,
                                  double phi,
                                  double * H,
                                  double * dH,
                                  double * ddH,
                                  double * dddH
                                  );
 
  int primordial_inflation_indices(
                                   struct primordial * ppm
                                   );
 
  int primordial_inflation_solve_inflation(
                                           struct perturbations * ppt,
                                           struct primordial * ppm,
                                           struct precision * ppr
                                           );
 
  int primordial_inflation_analytic_spectra(
                                            struct perturbations * ppt,
                                            struct primordial * ppm,
                                            struct precision * ppr,
                                            double * y_ini
                                            );
 
  int primordial_inflation_spectra(
                                   struct perturbations * ppt,
                                   struct primordial * ppm,
                                   struct precision * ppr,
                                   double * y_ini
                                   );
 
  int primordial_inflation_one_wavenumber(
                                          struct perturbations * ppt,
                                          struct primordial * ppm,
                                          struct precision * ppr,
                                          double * y_ini,
                                          int index_k
                                          );
 
  int primordial_inflation_one_k(
                                 struct primordial * ppm,
                                 struct precision * ppr,
                                 double k,
                                 double * y,
                                 double * dy,
                                 double * curvature,
                                 double * tensor
                                 );
 
  int primordial_inflation_find_attractor(
                                          struct primordial * ppm,
                                          struct precision * ppr,
                                          double phi_0,
                                          double precision,
                                          double * y,
                                          double * dy,
                                          double * H_0,
                                          double * dphidt_0
                                          );
 
  int primordial_inflation_evolve_background(
                                             struct primordial * ppm,
                                             struct precision * ppr,
                                             double * y,
                                             double * dy,
                                             enum target_quantity target,
                                             double stop,
                                             short check_epsilon,
                                             enum integration_direction direction,
                                             enum time_definition time
                                             );
 
  int primordial_inflation_check_potential(
                                           struct primordial * ppm,
                                           double phi,
                                           double * V,
                                           double * dV,
                                           double * ddV
                                           );
 
  int primordial_inflation_check_hubble(
                                        struct primordial * ppm,
                                        double phi,
                                        double *H,
                                        double * dH,
                                        double * ddH,
                                        double * dddH
                                        );
 
  int primordial_inflation_get_epsilon(
                                       struct primordial * ppm,
                                       double phi,
                                       double * epsilon
                                       );
 
  int primordial_inflation_find_phi_pivot(
                                          struct primordial * ppm,
                                          struct precision * ppr,
                                          double * y,
                                          double * dy
                                          );
 
  int primordial_inflation_derivs(
                                  double tau,
                                  double * y,
                                  double * dy,
                                  void * parameters_and_workspace,
                                  ErrorMsg error_message
                                  );
 
  int primordial_external_spectrum_init(
                                        struct perturbations * ppt,
                                        struct primordial * ppm
                                        );
 
  int primordial_output_titles(struct perturbations * ppt,
                               struct primordial * ppm,
                               char titles[_MAXTITLESTRINGLENGTH_]
                               );
 
  int primordial_output_data(struct perturbations * ppt,
                             struct primordial * ppm,
                             int number_of_titles,
                             double *data);
#ifdef __cplusplus
}
#endif
 
/**************************************************************/
 
 
#define _K_PER_DECADE_PRIMORDIAL_MIN_ 1.
 
 
#endif
/* @endcond */