primordial.h

Go to the documentation of this file.

 
#ifndef __PRIMORDIAL__
#define __PRIMORDIAL__
 
#include <string>
 
#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 = 0.05; 
  enum primordial_spectrum_type primordial_spec_type =
      analytic_Pk; 
  /* - parameters describing the case primordial_spec_type = analytic_Pk : amplitudes, tilts, runnings, cross-correlations, ... */
 
  double A_s    = 2.215e-9; 
  double sigma8 = 0.;       
  double n_s = 0.9619; 
  double alpha_s = 0.; 
  double r = 1.; 
  double n_t =
      -r / 8. *
      (2. - r / 8. -
       n_s); 
  double alpha_t = r / 8. *
                   (r / 8. + n_s -
                    1.); 
  double f_bi     = 1.; 
  double n_bi     = 1.; 
  double alpha_bi = 0.; 
  double f_cdi     = 1.; 
  double n_cdi     = 1.; 
  double alpha_cdi = 0.; 
  double f_nid =
      1.; 
  double n_nid     = 1.; 
  double alpha_nid = 0.; 
  double f_niv =
      1.; 
  double n_niv     = 1.; 
  double alpha_niv = 0.; 
  double c_ad_bi     = 0.; 
  double n_ad_bi     = 0.; 
  double alpha_ad_bi = 0.; 
  double c_ad_cdi     = 0.; 
  double n_ad_cdi     = 0.; 
  double alpha_ad_cdi = 0.; 
  double c_ad_nid     = 0.; 
  double n_ad_nid     = 0.; 
  double alpha_ad_nid = 0.; 
  double c_ad_niv     = 0.; 
  double n_ad_niv     = 0.; 
  double alpha_ad_niv = 0.; 
  double c_bi_cdi     = 0.; 
  double n_bi_cdi     = 0.; 
  double alpha_bi_cdi = 0.; 
  double c_bi_nid     = 0.; 
  double n_bi_nid     = 0.; 
  double alpha_bi_nid = 0.; 
  double c_bi_niv     = 0.; 
  double n_bi_niv     = 0.; 
  double alpha_bi_niv = 0.; 
  double c_cdi_nid     = 0.; 
  double n_cdi_nid     = 0.; 
  double alpha_cdi_nid = 0.; 
  double c_cdi_niv     = 0.; 
  double n_cdi_niv     = 0.; 
  double alpha_cdi_niv = 0.; 
  double c_nid_niv     = 0.; 
  double n_nid_niv     = 0.; 
  double alpha_nid_niv = 0.; 
  enum potential_shape potential = polynomial;
 
  double V0 = 1.25e-13;  
  double V1 = -1.12e-14; 
  double V2 = -6.95e-14; 
  double V3 = 0.;        
  double V4 = 0.;        
  /* parameters describing the case primordial_spec_type = inflation_H */
 
  double H0 = 3.69e-6;  
  double H1 = -5.84e-7; 
  double H2 = 0.;       
  double H3 = 0.;       
  double H4 = 0.;       
  /* parameters describing inflation_V_end */
 
  double phi_end = 0.; 
  enum phi_pivot_methods phi_pivot_method =
      N_star; 
  double phi_pivot_target =
      60; 
  /* behavior of the inflation module */
  enum inflation_module_behavior behavior =
      numerical; 
  std::string command;  
  double custom1  = 0.; 
  double custom2  = 0.; 
  double custom3  = 0.; 
  double custom4  = 0.; 
  double custom5  = 0.; 
  double custom6  = 0.; 
  double custom7  = 0.; 
  double custom8  = 0.; 
  double custom9  = 0.; 
  double custom10 = 0.; 
 
 
  short primordial_verbose =
      0; 
};
 
class PrimordialModule;
struct primordial_inflation_parameters_and_workspace {
  primordial_inflation_parameters_and_workspace(PrimordialModule* p_m) : primordial_module(p_m) {}
  const PrimordialModule* primordial_module;
  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;
};
 
 
#define _K_PER_DECADE_PRIMORDIAL_MIN_ 1.
 
 
#endif