background.h

Go to the documentation of this file.

 
#ifndef __BACKGROUND__
#define __BACKGROUND__
 
#include <memory>
#include <vector>
 
#include "arrays.h"
#include "common.h"
#include "dei_rkck.h"
#include "parser.h"
#include "quadrature.h"
 
class NonColdDarkMatter;
class DarkRadiation;
 
enum equation_of_state { CLP, EDE, PhenoAxion };
 
enum class ClosureSpecies { None, Lambda, Fluid, ScalarField };
 
struct background {
 
  double T_cmb = 2.7255; 
  double h = 0.67556; 
  double H0 =
      h * 1.e5 /
      _c_; 
  double Omega0_g =
      Omega0gFromTcmb(T_cmb,
                      h); 
  double Omega0_b = 0.022032 /
                    (h * h); 
  // Per-species Omega0_* fields removed: each species owns its density
  // through its constructor (queryable via species->GetOmega0()).  Only
  // photons (Omega0_g), baryons (Omega0_b), and curvature (Omega0_k) live
  // on pba.
 
  // Fluid physics params (fluid_equation_of_state, w0_fld, wa_fld, cs2_fld,
  // Omega_EDE) live on FluidSpecies; PPF-only params (use_ppf, c_gamma_over_c_fld)
  // live on PpfFluid — see FluidSpecies::CreateAll for the parser and
  // species/fluid.h / species/ppf_fluid.h for the accessors.
 
  // DCDM physics params (Gamma_dcdm, Omega_ini_dcdm) live on DCDMSpecies — see
  // DCDM_DR_Species::CreateAll for the parser and species/dcdm.h for the
  // accessors used by cross-module readers.
  // IDR physics params (T_idr, l_max_idr) live on IDRSpecies — see
  // IDM_DR_IDR_Species::CreateAll for the parser and species/idr.h for the
  // accessors used by cross-module readers.
 
  double Omega0_k = 0.; 
 
  double K = 0.; 
  int sgnK = 0;  
 
  ClosureSpecies closure_species = ClosureSpecies::None;
 
 
  short short_info;  
  short normal_info; 
  short long_info;   
  short
      inter_normal; 
  short
      inter_closeby; 
 
 
  unsigned int number_of_threads;
 
  short background_verbose =
      0; 
};
 
 
#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 
 
#endif