precision.h

Go to the documentation of this file.

 
#ifndef CLASS_PRECISION_H
#define CLASS_PRECISION_H
 
#include <float.h>  // DBL_EPSILON
#include <string>
 
#include "constants.h"  // physics constants used in precision defaults
#include "errors.h"     // status codes
 
class FileContent;  // forward decl for precision::parse
 
enum class evolver_type {
  rk,    /* Runge-Kutta integrator */
  ndf15, /* stiff integrator */
  rkdp45 /* Dormand-Prince 4(5) explicit adaptive integrator */
};
 
enum class pk_def {
  delta_m_squared, 
  delta_tot_squared, 
  delta_bc_squared, 
  delta_tot_from_poisson_squared 
};
enum class file_format { class_format, camb_format };
 
enum class tca_method {
  first_order_MB,
  first_order_CAMB,
  first_order_CLASS,
  second_order_CRS,
  second_order_CLASS,
  compromise_CLASS
};
enum class rsa_method { rsa_null, rsa_MD, rsa_MD_with_reio, rsa_none };
enum class rsa_idr_method { rsa_idr_none, rsa_idr_MD };
enum class ufa_method { ufa_mb, ufa_hu, ufa_CLASS, ufa_none };
enum class ncdmfa_method { ncdmfa_mb, ncdmfa_hu, ncdmfa_CLASS, ncdmfa_none };
 
struct precision {
  std::string class_dir;
 
  /*
   * Background Quantities
   * */
 
  double a_ini_over_a_today_default = 1.e-14;
  double back_integration_stepsize = 7.e-3;
  double tol_background_integration = 1.e-6;
  double tol_initial_Omega_r = 1.e-4;
  double tol_M_ncdm = 1.e-7;
  double tol_ncdm_synchronous = 3.4e-3;
  double tol_ncdm_newtonian = 9.05e-5;
  double tol_ncdm_bg = 1.e-5;
  double tol_ncdm_initial_w = 1.e-3;
  double tol_tau_eq = 1.e-6;
  double Omega0_cdm_min_synchronous = 1.e-10;
  /*
   * Currently unused parameter.
   */
  //class_precision_parameter(safe_phi_scf,double,0.0)
  std::string sBBN_file = "/bbn/sBBN_2017.dat";
 
  /*
   *  Thermodynamical quantities
   * */
 
  double recfast_z_initial = 1.0e4;
  int recfast_Nz0 = 20000;
  double thermo_z_initial_idm_dr = 1.0e9;
  int thermo_Nz1_idm_dr          = 10000;
  int thermo_Nz2_idm_dr          = 100;
  double tol_thermo_integration = 1.0e-7;
  /*
   * Recfast 1.4 switch parameters
   */
  int recfast_Heswitch =
      6; 
  double recfast_fudge_He =
      0.86; 
  /*
   * Recfast 1.5 parameters
   */
  int recfast_Hswitch =
      1; 
  double recfast_fudge_H =
      1.14; 
  double recfast_delta_fudge_H =
      -0.015; 
  double recfast_AGauss1 = -0.14; 
  double recfast_AGauss2 =
      0.079; 
  double recfast_zGauss1 = 7.28; 
  double recfast_zGauss2 = 6.73; 
  double recfast_wGauss1 = 0.18; 
  double recfast_wGauss2 = 0.33; 
  double recfast_z_He_1 = 8000.0; 
  double recfast_delta_z_He_1 =
      50.0; 
  double recfast_z_He_2 = 5000.0; 
  double recfast_delta_z_He_2 =
      100.0; 
  double recfast_z_He_3 = 3500.0; 
  double recfast_delta_z_He_3 =
      50.0; 
  double recfast_x_He0_trigger =
      0.995; 
  double recfast_x_He0_trigger_delta =
      0.05; 
  double recfast_x_H0_trigger =
      0.995; 
  double recfast_x_H0_trigger2 =
      0.995; 
  double recfast_x_H0_trigger_delta =
      0.05; 
  double recfast_H_frac =
      1.0e-3; 
  double reionization_z_start_max = 50.0;   
  double reionization_sampling    = 5.0e-2; 
  double reionization_optical_depth_tol =
      1.0e-4; 
  double reionization_start_factor =
      8.0; 
  int thermo_rate_smoothing_radius =
      50; 
  std::string hyrec_Alpha_inf_file =
      "/hyrec/Alpha_inf.dat"; 
  std::string hyrec_R_inf_file =
      "/hyrec/R_inf.dat"; 
  std::string hyrec_two_photon_tables_file =
      "/hyrec/two_photon_tables.dat"; 
  double k_min_tau0 =
      0.1; 
  double k_max_tau0_over_l_max =
      1.8; 
  double k_step_sub =
      0.05; 
  double k_step_super =
      0.002; 
  double k_step_transition =
      0.2; 
  double k_step_super_reduction =
      0.1; 
  double k_per_decade_for_pk =
      10.0; 
  double idmdr_boost_k_per_decade_for_pk =
      1.0; 
  double k_per_decade_for_bao =
      70.0; 
  double k_bao_center =
      3.0; 
  double k_bao_width =
      4.0; 
  double start_small_k_at_tau_c_over_tau_h =
      0.0015; 
  double start_large_k_at_tau_h_over_tau_k =
      0.07; 
  double tight_coupling_trigger_tau_c_over_tau_h = 0.015;
 
  double tight_coupling_trigger_tau_c_over_tau_k = 0.01;
 
  double start_sources_at_tau_c_over_tau_h =
      0.008; 
  int tight_coupling_approximation = static_cast<int>(
      tca_method::compromise_CLASS); 
  double idm_dr_tight_coupling_trigger_tau_c_over_tau_k =
      0.01; 
  double idm_dr_tight_coupling_trigger_tau_c_over_tau_h =
      0.015; 
  double idm_drmd_tight_coupling_trigger_G_over_aH =
      100000; 
  int l_max_g =
      12; 
  int l_max_pol_g =
      10; 
  int l_max_dr =
      17; 
  int l_max_dr_col =
      17; 
  int l_max_ur =
      17; 
  int l_max_idr =
      17; 
  int l_max_ncdm =
      17; 
  int l_max_g_ten =
      5; 
  int l_max_pol_g_ten =
      5; 
  double curvature_ini = 1.0; 
  double entropy_ini   = 1.0; 
  double gw_ini        = 1.0; 
  double perturb_integration_stepsize = 0.5;
 
  double perturb_sampling_stepsize = 0.1;
  double perturbations_sampling_boost_above_age_fraction = 0.9;
 
  double tol_perturb_integration = 1.0e-5;
 
  double c_gamma_k_H_square_max = 1.0e3;
 
  double tol_tau_approx = 1.0e-10;
 
  int radiation_streaming_approximation = static_cast<int>(rsa_method::rsa_MD_with_reio);
 
  double radiation_streaming_trigger_tau_over_tau_k = 45.0;
 
  double radiation_streaming_trigger_tau_c_over_tau = 5.0;
 
  int idr_streaming_approximation = static_cast<int>(
      rsa_idr_method::rsa_idr_none); 
  double idr_streaming_trigger_tau_over_tau_k =
      50.0; 
  double idr_streaming_trigger_tau_c_over_tau =
      10.0; 
  int ur_fluid_approximation = static_cast<int>(
      ufa_method::ufa_CLASS); 
  double ur_fluid_trigger_tau_over_tau_k = 30.0;
 
  int ncdm_fluid_approximation = static_cast<int>(
      ncdmfa_method::ncdmfa_CLASS); 
  double ncdm_fluid_trigger_tau_over_tau_k = 31.0;
 
  double neglect_CMB_sources_below_visibility = 1.0e-3;
 
  evolver_type evolver = evolver_type::ndf15;
 
  /*
   * Primordial parameters
   * */
 
  double k_per_decade_primordial =
      10.0; 
  double primordial_inflation_ratio_min =
      100.0; 
  double primordial_inflation_ratio_max =
      1.0 /
      50.0; 
  int primordial_inflation_phi_ini_maxit =
      10000; 
  double primordial_inflation_pt_stepsize =
      0.01; 
  double primordial_inflation_bg_stepsize =
      0.005; 
  double primordial_inflation_tol_integration =
      1.0e-3; 
  double primordial_inflation_attractor_precision_pivot =
      0.001; 
  double primordial_inflation_attractor_precision_initial =
      0.1; 
  int primordial_inflation_attractor_maxit =
      10; 
  double primordial_inflation_tol_curvature =
      1.0e-3; 
  double primordial_inflation_aH_ini_target =
      0.9; 
  double primordial_inflation_end_dphi =
      1.0e-10; 
  double primordial_inflation_end_logstep =
      10.0; 
  double primordial_inflation_small_epsilon =
      0.1; 
  double primordial_inflation_small_epsilon_tol = 0.01; 
  double primordial_inflation_extra_efolds =
      2.0; 
  /*
   * Transfer function parameters
   * */
 
  int l_linstep =
      40; 
  double l_logstep =
      1.12; 
  double hyper_x_min =
      1.0e-5; 
  double hyper_sampling_flat =
      8.0; 
  double hyper_sampling_curved_low_nu =
      7.0; 
  double hyper_sampling_curved_high_nu =
      3.0; 
  double hyper_nu_sampling_step =
      1000.0; 
  double hyper_phi_min_abs =
      1.0e-10; 
  double hyper_x_tol = 1.0e-4; 
  double hyper_flat_approximation_nu =
      4000.0; 
  double q_linstep = 0.45; 
  double q_logstep_spline = 170.0; 
  double q_logstep_open = 6.0; 
  double q_logstep_trapzd = 20.0; 
  double q_numstep_transition = 250.0; 
  double q_logstep_limber =
      1.025; 
  double k_max_limber_over_l_max_scalars =
      0.001; 
  double transfer_neglect_delta_k_S_t0 =
      0.15; 
  double transfer_neglect_delta_k_S_t1 =
      0.04; 
  double transfer_neglect_delta_k_S_t2 =
      0.15; 
  double transfer_neglect_delta_k_S_e =
      0.11; 
  double transfer_neglect_delta_k_V_t1 =
      1.0; 
  double transfer_neglect_delta_k_V_t2 =
      1.0; 
  double transfer_neglect_delta_k_V_e =
      1.0; 
  double transfer_neglect_delta_k_V_b =
      1.0; 
  double transfer_neglect_delta_k_T_t2 =
      0.2; 
  double transfer_neglect_delta_k_T_e =
      0.25; 
  double transfer_neglect_delta_k_T_b =
      0.1; 
  double transfer_neglect_late_source =
      400.0; 
  double l_switch_limber =
      10.; 
  // For density Cl, we recommend not to use the Limber approximation
  // at all, and hence to put here a very large number (e.g. 10000); but
  // if you have wide and smooth selection functions you may wish to
  // use it; then 100 might be OK
  double l_switch_limber_for_nc_local_over_z =
      100.0; 
  // For terms integrated along the line-of-sight involving spherical
  // Bessel functions (but not their derivatives), Limber
  // approximation works well. High precision can be reached with 2000
  // only. But if you have wide and smooth selection functions you may
  // reduce to e.g. 30.
  double l_switch_limber_for_nc_los_over_z =
      30.0; 
  double selection_cut_at_sigma =
      5.0; 
  double selection_sampling =
      50.0; 
  double selection_sampling_bessel =
      20.0; 
  double selection_sampling_bessel_los =
      20.0; 
  double selection_tophat_edge =
      0.1; 
  /*
   * Nonlinear module precision parameters
   * */
 
  double sigma_k_per_decade =
      80.; 
  double nonlinear_min_k_max = 5.0; 
  double halofit_min_k_nonlinear =
      1.0e-4; 
  double halofit_min_k_max = 5.0; 
  double halofit_k_per_decade = 80.0; 
  double halofit_sigma_precision = 0.05; 
  double halofit_tol_sigma = 1.0e-6; 
  double pk_eq_z_max = 5.0;    
  double pk_eq_tol   = 1.0e-7; 
  double hmcode_max_k_extra = 1.e6; 
  double hmcode_min_k_max = 5.; 
  double hmcode_tol_sigma = 1.e-6; 
  int n_hmcode_tables      = 64;
  double rmin_for_sigtab   = 1.e-5;
  double rmax_for_sigtab   = 1.e3;
  double ainit_for_growtab = 1.e-3;
  double amax_for_growtab  = 1.;
 
  int nsteps_for_p1h_integral  = 256;
  double mmin_for_p1h_integral = 1.e3;
  double mmax_for_p1h_integral = 1.e18;
 
  /*
   * Lensing precision parameters
   * */
 
  int accurate_lensing =
      0; 
  int num_mu_minus_lmax =
      70;                
  int delta_l_max = 500; 
  double tol_gauss_legendre =
      DBL_EPSILON; 
  double smallest_allowed_variation =
      DBL_EPSILON; 
 
  void ResolveDataPaths();
 
  void parse(const FileContent& fc);
};
 
#endif  // CLASS_PRECISION_H