#include <precision.h>

Data Fields
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

std::string	sBBN_file = "/bbn/sBBN_2017.dat"

double	recfast_z_initial = 1.0e4

int	recfast_Nz0 = 20000

double	thermo_z_initial_idm_dr = 1.0e9

double	tol_thermo_integration = 1.0e-7

int	recfast_Heswitch

double	recfast_fudge_He

int	recfast_Hswitch

double	recfast_fudge_H

double	recfast_delta_fudge_H

double	recfast_AGauss1 = -0.14

double	recfast_AGauss2

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

double	recfast_z_He_2 = 5000.0

double	recfast_delta_z_He_2

double	recfast_z_He_3 = 3500.0

double	recfast_delta_z_He_3

double	recfast_x_He0_trigger

double	recfast_x_He0_trigger_delta

double	recfast_x_H0_trigger

double	recfast_x_H0_trigger2

double	recfast_x_H0_trigger_delta

double	recfast_H_frac

double	reionization_z_start_max = 50.0

double	reionization_sampling = 5.0e-2

double	reionization_optical_depth_tol

double	reionization_start_factor

int	thermo_rate_smoothing_radius

std::string	hyrec_Alpha_inf_file

std::string	hyrec_R_inf_file

std::string	hyrec_two_photon_tables_file

double	k_min_tau0

double	k_max_tau0_over_l_max

double	k_step_sub

double	k_step_super

double	k_step_transition

double	k_step_super_reduction

double	k_per_decade_for_pk

double	idmdr_boost_k_per_decade_for_pk

double	k_per_decade_for_bao

double	k_bao_center

double	k_bao_width

double	start_small_k_at_tau_c_over_tau_h

double	start_large_k_at_tau_h_over_tau_k

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

int	tight_coupling_approximation

double	idm_dr_tight_coupling_trigger_tau_c_over_tau_k

double	idm_dr_tight_coupling_trigger_tau_c_over_tau_h

double	idm_drmd_tight_coupling_trigger_G_over_aH

int	l_max_g

int	l_max_pol_g

int	l_max_dr

int	l_max_dr_col

int	l_max_ur

int	l_max_idr

int	l_max_ncdm

int	l_max_g_ten

int	l_max_pol_g_ten

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

double	idr_streaming_trigger_tau_over_tau_k

double	idr_streaming_trigger_tau_c_over_tau

int	ur_fluid_approximation

double	ur_fluid_trigger_tau_over_tau_k = 30.0

int	ncdm_fluid_approximation

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

double	k_per_decade_primordial

double	primordial_inflation_ratio_min

double	primordial_inflation_ratio_max

int	primordial_inflation_phi_ini_maxit

double	primordial_inflation_pt_stepsize

double	primordial_inflation_bg_stepsize

double	primordial_inflation_tol_integration

double	primordial_inflation_attractor_precision_pivot

double	primordial_inflation_attractor_precision_initial

int	primordial_inflation_attractor_maxit

double	primordial_inflation_tol_curvature

double	primordial_inflation_aH_ini_target

double	primordial_inflation_end_dphi

double	primordial_inflation_end_logstep

double	primordial_inflation_small_epsilon

double	primordial_inflation_small_epsilon_tol = 0.01

double	primordial_inflation_extra_efolds

int	l_linstep

double	l_logstep

double	hyper_x_min

double	hyper_sampling_flat

double	hyper_sampling_curved_low_nu

double	hyper_sampling_curved_high_nu

double	hyper_nu_sampling_step

double	hyper_phi_min_abs

double	hyper_x_tol = 1.0e-4

double	hyper_flat_approximation_nu

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

double	k_max_limber_over_l_max_scalars

double	transfer_neglect_delta_k_S_t0

double	transfer_neglect_delta_k_S_t1

double	transfer_neglect_delta_k_S_t2

double	transfer_neglect_delta_k_S_e

double	transfer_neglect_delta_k_V_t1

double	transfer_neglect_delta_k_V_t2

double	transfer_neglect_delta_k_V_e

double	transfer_neglect_delta_k_V_b

double	transfer_neglect_delta_k_T_t2

double	transfer_neglect_delta_k_T_e

double	transfer_neglect_delta_k_T_b

double	transfer_neglect_late_source

double	l_switch_limber

double	l_switch_limber_for_nc_local_over_z

double	l_switch_limber_for_nc_los_over_z

double	selection_cut_at_sigma

double	selection_sampling

double	selection_sampling_bessel

double	selection_sampling_bessel_los

double	selection_tophat_edge

double	sigma_k_per_decade

double	nonlinear_min_k_max = 5.0

double	halofit_min_k_nonlinear

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

int	nsteps_for_p1h_integral = 256

int	accurate_lensing

int	num_mu_minus_lmax

int	delta_l_max = 500

double	tol_gauss_legendre

- general precision parameters
double	smallest_allowed_variation

void	ResolveDataPaths ()

void	parse (const FileContent &fc)

Detailed Description

All precision parameters.

Includes integrations steps, flags telling how the computation is to be performed, etc.

Member Function Documentation

◆ ResolveDataPaths()

void precision::ResolveDataPaths ( )

Prepend the runtime class_dir to the relative data-file path defaults (sBBN_file, hyrec_*). Call once after class_dir is set and before parse(), so a user-supplied absolute override in parse() still replaces verbatim.

◆ parse()

void precision::parse ( const FileContent & fc )

Parse precision parameters from a configuration file.

Field Documentation

◆ a_ini_over_a_today_default

double precision::a_ini_over_a_today_default = 1.e-14

Default initial value of scale factor used in the integration of background quantities. For models like ncdm, the code may decide to start the integration earlier.

◆ back_integration_stepsize

double precision::back_integration_stepsize = 7.e-3

Default stepsize in conformal time for the background integration, in units for the conformal Hubble time. dtau = back_integration_stepsize/aH

◆ tol_background_integration

double precision::tol_background_integration = 1.e-6

Tolerance of the background integration, giving the allowed relative integration error. Passed to the evolver in background_solve_evolver(). The default is 1e-6 (not the historical 1e-2): the removed legacy RK solver had no dense output, so the requested output points capped the step size and the step count acted as the effective precision knob; the evolver instead relies on this tolerance directly, so it must be tight.

◆ tol_initial_Omega_r

double precision::tol_initial_Omega_r = 1.e-4

Tolerance of the deviation of $\Omega_r$ from 1 for which to start integration: The starting point of integration will be chosen, such that the Omega of radiation at that point is close to 1 within tolerance. (Class starts background integration during complete radiation domination)

◆ tol_M_ncdm

double precision::tol_M_ncdm = 1.e-7

Tolerance of relative deviation of the used non-cold dark matter mass compared to that which would give the correct density. The dark matter mass is estimated from the dark matter density using a Newton-Method. In the nonrelativistic limit, this could be estimated using M=density/number density

◆ tol_ncdm_synchronous

double precision::tol_ncdm_synchronous = 3.4e-3

Tolerance on the relative precision of the integration over non-cold dark matter phase-space distributions in the synchronous gauge. The convenience input key "tol_ncdm" sets both this and tol_ncdm_newtonian.

◆ tol_ncdm_newtonian

double precision::tol_ncdm_newtonian = 9.05e-5

Tolerance on the relative precision of the integration over non-cold dark matter phase-space distributions in the newtonian gauge. Newtonian gauge needs a finer ncdm momentum grid than synchronous for the same P(k) accuracy, hence the tighter default.

◆ tol_ncdm_bg

double precision::tol_ncdm_bg = 1.e-5

Tolerance on the relative precision of the integration over non-cold dark matter phase-space distributions during the background evolution.

◆ tol_ncdm_initial_w

double precision::tol_ncdm_initial_w = 1.e-3

Tolerance on the initial deviation of non-cold dark matter from being fully relativistic. Using w = pressure/density, this quantifies the maximum deviation from 1/3. (for relativistic species)

◆ tol_tau_eq

double precision::tol_tau_eq = 1.e-6

Tolerance on the deviation of the conformal time of equality from the true value in 1/Mpc.

◆ Omega0_cdm_min_synchronous

double precision::Omega0_cdm_min_synchronous = 1.e-10

Minimum amount of cdm to allow calculations in synchronous gauge comoving with cdm.

◆ sBBN_file

std::string precision::sBBN_file = "/bbn/sBBN_2017.dat"

Big Bang Nucleosynthesis file path. The file specifies the predictions for $Y_\mathrm{He}$ for given $\omega_b$ and $N_\mathrm{eff}$ .

◆ recfast_z_initial

double precision::recfast_z_initial = 1.0e4

The initial z for the recfast calculation of the recombination history, e.g. 10^4

◆ recfast_Nz0

int precision::recfast_Nz0 = 20000

Number of recfast integration steps, e.g. if this is 1.10^4 and the previous one is 10^4, the step will be Delta z = 0.5

◆ thermo_z_initial_idm_dr

double precision::thermo_z_initial_idm_dr = 1.0e9

If there is interacting DM, we want the thermodynamics table to start at a much larger z, in order to capture the possible non-trivial behavior of the dark matter interaction rate at early times:

The new initial redshift will be thermo_z_initial_idm_dr
the highest redhsift will be sampled with thermo_Nz1_idm_dr values, and the step will be Delta z = (thermo_z_initial_idm_dr-recfast_z_initial)/thermo_Nz1_idm_dr For instance, if the previous value is 10^9 and this value is 10^4, then Delta z simeq 10^5
But the first interval after recfast_z_initial will be better sampled with thermo_Nz2_idm_dr values, in order to ensure a smoother transition from a small step to a large step. The intermediate stepsize will then be Delta z = (thermo_z_initial_idm_dr-recfast_z_initial)/thermo_Nz1_idm_dr/thermo_Nz1_idm_dr. For instance, if the three values are (10^9, 10^4, 10^2), then the intermediate timestep is Delta z simeq 10^3

◆ tol_thermo_integration

double precision::tol_thermo_integration = 1.0e-7

Tolerance for thermodynamical integration. RECFAST now uses the shared evolver interface with dense output, so this is a true ODE tolerance rather than a legacy per-redshift-bin RK correction tolerance.

◆ recfast_Heswitch

int precision::recfast_Heswitch

Initial value:

=

6

from recfast 1.4, specifies how accurate the Helium recombination should be handled

◆ recfast_fudge_He

double precision::recfast_fudge_He

Initial value:

=

0.86

from recfast 1.4, fugde factor for Peeble's equation coefficient of Helium

◆ recfast_Hswitch

int precision::recfast_Hswitch

Initial value:

=

1

from recfast 1.5, specifies how accurate the Hydrogen recombination should be handled

◆ recfast_fudge_H

double precision::recfast_fudge_H

Initial value:

=

1.14

from recfast 1.4, fudge factor for Peeble's equation coeffient of Hydrogen

◆ recfast_delta_fudge_H

double precision::recfast_delta_fudge_H

Initial value:

=

-0.015

from recfast 1.5.2, increasing Hydrogen fudge factor if Hswitch is enabled

◆ recfast_AGauss1

double precision::recfast_AGauss1 = -0.14

from recfast 1.5, Gaussian Peeble prefactor fit, amplitude

◆ recfast_AGauss2

double precision::recfast_AGauss2

Initial value:

=

0.079

from recfast 1.5.2, Gaussian Peeble prefactor fit, amplitude

◆ recfast_zGauss1

double precision::recfast_zGauss1 = 7.28

from recfast 1.5, Gaussian Peeble prefactor fit, center

◆ recfast_zGauss2

double precision::recfast_zGauss2 = 6.73

from recfast 1.5.2, Gaussian Peeble prefactor fit, center

◆ recfast_wGauss1

double precision::recfast_wGauss1 = 0.18

from recfast 1.5, Gaussian Peeble prefactor fit, width

◆ recfast_wGauss2

double precision::recfast_wGauss2 = 0.33

from recfast 1.5, Gaussian Peeble prefactor fit, width

◆ recfast_z_He_1

double precision::recfast_z_He_1 = 8000.0

from recfast 1.4, Starting value of Helium recombination 1

◆ recfast_delta_z_He_1

double precision::recfast_delta_z_He_1

Initial value:

=

50.0

Smoothing factor for recombination approximation switching, found to be OK on 3.09.10

◆ recfast_z_He_2

double precision::recfast_z_He_2 = 5000.0

from recfast 1.4, Ending value of Helium recombination 1

◆ recfast_delta_z_He_2

double precision::recfast_delta_z_He_2

Initial value:

=

100.0

Smoothing factor for recombination approximation switching, found to be OK on 3.09.10

◆ recfast_z_He_3

double precision::recfast_z_He_3 = 3500.0

from recfast 1.4, Starting value of Helium recombination 2

◆ recfast_delta_z_He_3

double precision::recfast_delta_z_He_3

Initial value:

=

50.0

Smoothing factor for recombination approximation switching, found to be OK on 3.09.10

◆ recfast_x_He0_trigger

double precision::recfast_x_He0_trigger

Initial value:

=

0.995

Switch for Helium full calculation during reco, raised from 0.99 to 0.995 for smoother Helium

◆ recfast_x_He0_trigger_delta

double precision::recfast_x_He0_trigger_delta

Initial value:

=

0.05

Smoothing factor for recombination approximation switching, found to be OK on 3.09.10

◆ recfast_x_H0_trigger

double precision::recfast_x_H0_trigger

Initial value:

=

0.995

Switch for Hydrogen full calculation during reco, raised from 0.99 to 0.995 for smoother Hydrogen

◆ recfast_x_H0_trigger2

double precision::recfast_x_H0_trigger2

Initial value:

=

0.995

Switch for Hydrogen full calculation during reco, for changing Hydrogen flag, raised from 0.98 to same as previous one for smoother Hydrogen

◆ recfast_x_H0_trigger_delta

double precision::recfast_x_H0_trigger_delta

Initial value:

=

0.05

Smoothing factor for recombination approximation switching, found to be OK on 3.09.10

◆ recfast_H_frac

double precision::recfast_H_frac

Initial value:

=

1.0e-3

from recfast 1.4, specifies the time at which the temperature evolution is calculated by the more precise equation

◆ reionization_z_start_max

double precision::reionization_z_start_max = 50.0

Maximum starting value in z for reionization

◆ reionization_sampling

double precision::reionization_sampling = 5.0e-2

Sampling density in z during reionization

◆ reionization_optical_depth_tol

double precision::reionization_optical_depth_tol

Initial value:

=

1.0e-4

Relative tolerance on finding the user-given optical depth of reionization given a certain redshift of reionization

◆ reionization_start_factor

double precision::reionization_start_factor

Initial value:

=

8.0

Searching optical depth corresponding to the redshift is started from an initial offset beyond z_reionization_start, multiplied by reionization_width

◆ thermo_rate_smoothing_radius

int precision::thermo_rate_smoothing_radius

Initial value:

=

50

Smoothing in redshift of the variation rate of $\exp(-\kappa)$ , g, and $\frac{dg}{d\tau}$ that is used as a timescale afterwards

◆ hyrec_Alpha_inf_file

std::string precision::hyrec_Alpha_inf_file

Initial value:

=

"/hyrec/Alpha_inf.dat"

File containing the alpha parameter of hyrec

◆ hyrec_R_inf_file

std::string precision::hyrec_R_inf_file

Initial value:

=

"/hyrec/R_inf.dat"

File containing the R_inf parameter of hyrec

◆ hyrec_two_photon_tables_file

std::string precision::hyrec_two_photon_tables_file

Initial value:

=

"/hyrec/two_photon_tables.dat"

File containing the two-photon interaction parameter of hyrec

◆ k_min_tau0

double precision::k_min_tau0

Initial value:

=

0.1

number defining k_min for the computation of Cl's and P(k)'s (dimensionless): (k_min tau_0), usually chosen much smaller than one

◆ k_max_tau0_over_l_max

double precision::k_max_tau0_over_l_max

Initial value:

=

1.8

number defining k_max for the computation of Cl's (dimensionless): (k_max tau_0)/l_max, usually chosen around two. Since v3.2.2, the separate full-Limber grid keeps CMB lensing accurate at high l, allowing the standard transfer grid to stop at 1.8 instead of 2.4.

◆ k_step_sub

double precision::k_step_sub

Initial value:

=

0.05

step in k space, in units of one period of acoustic oscillation at decoupling, for scales inside sound horizon at decoupling

◆ k_step_super

double precision::k_step_super

Initial value:

=

0.002

step in k space, in units of one period of acoustic oscillation at decoupling, for scales above sound horizon at decoupling

◆ k_step_transition

double precision::k_step_transition

Initial value:

=

0.2

dimensionless number regulating the transition from 'sub' steps to 'super' steps. Decrease for more precision.

◆ k_step_super_reduction

double precision::k_step_super_reduction

Initial value:

=

0.1

the step k_step_super is reduced by this amount in the k-->0 limit (below scale of Hubble and/or curvature radius)

◆ k_per_decade_for_pk

double precision::k_per_decade_for_pk

Initial value:

=

10.0

if values needed between kmax inferred from k_oscillations and k_kmax_for_pk, this gives the number of k per decade outside the BAO region

◆ idmdr_boost_k_per_decade_for_pk

double precision::idmdr_boost_k_per_decade_for_pk

Initial value:

=

1.0

boost factor for the case of DAO in idm-idr models

◆ k_per_decade_for_bao

double precision::k_per_decade_for_bao

Initial value:

=

70.0

if values needed between kmax inferred from k_oscillations and k_kmax_for_pk, this gives the number of k per decade inside the BAO region (for finer sampling)

◆ k_bao_center

double precision::k_bao_center

Initial value:

=

3.0

in ln(k) space, the central value of the BAO region where sampling is finer is defined as k_rec times this number (recommended: 3, i.e. finest sampling near 3rd BAO peak)

◆ k_bao_width

double precision::k_bao_width

Initial value:

=

4.0

in ln(k) space, width of the BAO region where sampling is finer: this number gives roughly the number of BAO oscillations well resolved on both sides of the central value (recommended: 4, i.e. finest sampling from before first up to 3+4=7th peak)

◆ start_small_k_at_tau_c_over_tau_h

double precision::start_small_k_at_tau_c_over_tau_h

Initial value:

=

0.0015

largest wavelengths start being sampled when universe is sufficiently opaque. This is quantified in terms of the ratio of thermo to hubble time scales, $\tau_c/\tau_H$ . Start when start_largek_at_tau_c_over_tau_h equals this ratio. Decrease this value to start integrating the wavenumbers earlier in time.

◆ start_large_k_at_tau_h_over_tau_k

double precision::start_large_k_at_tau_h_over_tau_k

Initial value:

=

0.07

largest wavelengths start being sampled when mode is sufficiently outside Hubble scale. This is quantified in terms of the ratio of hubble time scale to wavenumber time scale, $\tau_h/\tau_k$ which is roughly equal to (k*tau). Start when this ratio equals start_large_k_at_tau_k_over_tau_h. Decrease this value to start integrating the wavenumbers earlier in time.

◆ tight_coupling_trigger_tau_c_over_tau_h

double precision::tight_coupling_trigger_tau_c_over_tau_h = 0.015

when to switch off tight-coupling approximation: first condition: $\tau_c/\tau_H$ > tight_coupling_trigger_tau_c_over_tau_h. Decrease this value to switch off earlier in time. If this number is larger than start_sources_at_tau_c_over_tau_h, the code returns an error, because the source computation requires tight-coupling to be switched off.

◆ tight_coupling_trigger_tau_c_over_tau_k

double precision::tight_coupling_trigger_tau_c_over_tau_k = 0.01

when to switch off tight-coupling approximation: second condition: $\tau_c/\tau_k \equiv k \tau_c$ < tight_coupling_trigger_tau_c_over_tau_k. Decrease this value to switch off earlier in time.

◆ start_sources_at_tau_c_over_tau_h

double precision::start_sources_at_tau_c_over_tau_h

Initial value:

=

0.008

sources start being sampled when universe is sufficiently opaque. This is quantified in terms of the ratio of thermo to hubble time scales, $\tau_c/\tau_H$ . Start when start_sources_at_tau_c_over_tau_h equals this ratio. Decrease this value to start sampling the sources earlier in time.

◆ tight_coupling_approximation

int precision::tight_coupling_approximation

Initial value:

= static_cast<int>(

tca_method::compromise_CLASS)

method for tight coupling approximation

◆ idm_dr_tight_coupling_trigger_tau_c_over_tau_k

double precision::idm_dr_tight_coupling_trigger_tau_c_over_tau_k

Initial value:

=

0.01

when to switch off the dark-tight-coupling approximation, first condition (see normal tca for full definition)

◆ idm_dr_tight_coupling_trigger_tau_c_over_tau_h

double precision::idm_dr_tight_coupling_trigger_tau_c_over_tau_h

Initial value:

=

0.015

when to switch off the dark-tight-coupling approximation, second condition (see normal tca for full definition)

◆ idm_drmd_tight_coupling_trigger_G_over_aH

double precision::idm_drmd_tight_coupling_trigger_G_over_aH

Initial value:

=

100000

when to switch off the dark-tight-coupling approximation in DRMD, should be larger than at least 100 (currently set to a very high number as the code runs perfectly fine without the approximation.)

◆ l_max_g

int precision::l_max_g

Initial value:

=

12

number of momenta in Boltzmann hierarchy for photon temperature (scalar), at least 4

◆ l_max_pol_g

int precision::l_max_pol_g

Initial value:

=

10

number of momenta in Boltzmann hierarchy for photon polarization (scalar), at least 4

◆ l_max_dr

int precision::l_max_dr

Initial value:

=

17

number of momenta in Boltzmann hierarchy for decay radiation, at least 4

◆ l_max_dr_col

int precision::l_max_dr_col

Initial value:

=

17

number of collision terms in Boltzmann hierarchy for decay radiation, at least 2

◆ l_max_ur

int precision::l_max_ur

Initial value:

=

17

number of momenta in Boltzmann hierarchy for relativistic neutrino/relics (scalar), at least 4

◆ l_max_idr

int precision::l_max_idr

Initial value:

=

17

number of momenta in Boltzmann hierarchy for interacting dark radiation

◆ l_max_ncdm

int precision::l_max_ncdm

Initial value:

=

17

number of momenta in Boltzmann hierarchy for relativistic neutrino/relics (scalar), at least 4

◆ l_max_g_ten

int precision::l_max_g_ten

Initial value:

=

5

number of momenta in Boltzmann hierarchy for photon temperature (tensor), at least 4

◆ l_max_pol_g_ten

int precision::l_max_pol_g_ten

Initial value:

=

5

number of momenta in Boltzmann hierarchy for photon polarization (tensor), at least 4

◆ curvature_ini

double precision::curvature_ini = 1.0

initial condition for curvature for adiabatic

◆ entropy_ini

double precision::entropy_ini = 1.0

initial condition for entropy perturbation for isocurvature

◆ gw_ini

double precision::gw_ini = 1.0

initial condition for tensor metric perturbation h

◆ perturb_integration_stepsize

double precision::perturb_integration_stepsize = 0.5

default step $d \tau$ in perturbation integration, in units of the timescale involved in the equations (usually, the min of , , $1/\dot{\kappa}$ )

◆ perturb_sampling_stepsize

double precision::perturb_sampling_stepsize = 0.1

default step $d \tau$ for sampling the source function, in units of the timescale involved in the sources: $(\dot{\kappa}- \ddot{\kappa}/\dot{\kappa})^{-1}$

◆ perturbations_sampling_boost_above_age_fraction

double precision::perturbations_sampling_boost_above_age_fraction = 0.9

Age fraction above which source sampling is twice as fine. This improves the low-l CMB lensing line-of-sight integral; 1.0 disables the boost.

◆ tol_perturb_integration

double precision::tol_perturb_integration = 1.0e-5

control parameter for the precision of the perturbation integration, IMPORTANT FOR SETTING THE STEPSIZE OF NDF15

◆ c_gamma_k_H_square_max

double precision::c_gamma_k_H_square_max = 1.0e3

cutoff relevant for controlling stiffness in the PPF scheme. It is neccessary for the Runge-Kutta evolver, but not for ndf15. However, the approximation is excellent for a cutoff value of 1000, so we leave it on for both evolvers. (CAMB uses a cutoff value of 30.)

◆ tol_tau_approx

double precision::tol_tau_approx = 1.0e-10

precision with which the code should determine (by bisection) the times at which sources start being sampled, and at which approximations must be switched on/off (units of Mpc)

◆ radiation_streaming_approximation

int precision::radiation_streaming_approximation = static_cast<int>(rsa_method::rsa_MD_with_reio)

method for switching off photon perturbations

◆ radiation_streaming_trigger_tau_over_tau_k

double precision::radiation_streaming_trigger_tau_over_tau_k = 45.0

when to switch off photon perturbations, ie when to switch on photon free-streaming approximation (keep density and thtau, set shear and higher momenta to zero): first condition: $k \tau$ > radiation_streaming_trigger_tau_h_over_tau_k

◆ radiation_streaming_trigger_tau_c_over_tau

double precision::radiation_streaming_trigger_tau_c_over_tau = 5.0

when to switch off photon perturbations, ie when to switch on photon free-streaming approximation (keep density and theta, set shear and higher momenta to zero): second condition:

◆ idr_streaming_approximation

int precision::idr_streaming_approximation

Initial value:

= static_cast<int>(

rsa_idr_method::rsa_idr_none)

method for dark radiation free-streaming approximation

◆ idr_streaming_trigger_tau_over_tau_k

double precision::idr_streaming_trigger_tau_over_tau_k

Initial value:

=

50.0

when to switch on dark radiation (idr) free-streaming approximation, first condition

◆ idr_streaming_trigger_tau_c_over_tau

double precision::idr_streaming_trigger_tau_c_over_tau

Initial value:

=

10.0

when to switch on dark radiation (idr) free-streaming approximation, second condition

◆ ur_fluid_approximation

int precision::ur_fluid_approximation

Initial value:

= static_cast<int>(

ufa_method::ufa_CLASS)

method for ultra relativistic fluid approximation

◆ ur_fluid_trigger_tau_over_tau_k

double precision::ur_fluid_trigger_tau_over_tau_k = 30.0

when to switch off ur (massless neutrinos / ultra-relativistic relics) fluid approximation

◆ ncdm_fluid_approximation

int precision::ncdm_fluid_approximation

Initial value:

= static_cast<int>(

ncdmfa_method::ncdmfa_CLASS)

method for non-cold dark matter fluid approximation

◆ ncdm_fluid_trigger_tau_over_tau_k

double precision::ncdm_fluid_trigger_tau_over_tau_k = 31.0

when to switch off ncdm (massive neutrinos / non-cold relics) fluid approximation

◆ neglect_CMB_sources_below_visibility

double precision::neglect_CMB_sources_below_visibility = 1.0e-3

whether CMB source functions can be approximated as zero when visibility function g(tau) is tiny

◆ evolver

evolver_type precision::evolver = evolver_type::ndf15

The type of evolver to use: options are ndf15 or rk

◆ k_per_decade_primordial

double precision::k_per_decade_primordial

Initial value:

=

10.0

logarithmic sampling for primordial spectra (number of points per decade in k space)

◆ primordial_inflation_ratio_min

double precision::primordial_inflation_ratio_min

Initial value:

=

100.0

for each k, start following wavenumber when aH = k/primordial_inflation_ratio_min

◆ primordial_inflation_ratio_max

double precision::primordial_inflation_ratio_max

Initial value:

=
      1.0 /
      50.0

for each k, stop following wavenumber, at the latest, when aH = k/primordial_inflation_ratio_max

◆ primordial_inflation_phi_ini_maxit

int precision::primordial_inflation_phi_ini_maxit

Initial value:

=

10000

maximum number of iteration when searching a suitable initial field value phi_ini (value reached when no long-enough slow-roll period before the pivot scale)

◆ primordial_inflation_pt_stepsize

double precision::primordial_inflation_pt_stepsize

Initial value:

=

0.01

controls the integration timestep for inflaton perturbations

◆ primordial_inflation_bg_stepsize

double precision::primordial_inflation_bg_stepsize

Initial value:

=

0.005

controls the integration timestep for inflaton background

◆ primordial_inflation_tol_integration

double precision::primordial_inflation_tol_integration

Initial value:

=

1.0e-3

controls the precision of the ODE integration during inflation

◆ primordial_inflation_attractor_precision_pivot

double precision::primordial_inflation_attractor_precision_pivot

Initial value:

=

0.001

targeted precision when searching attractor solution near phi_pivot

◆ primordial_inflation_attractor_precision_initial

double precision::primordial_inflation_attractor_precision_initial

Initial value:

=

0.1

targeted precision when searching attractor solution near phi_ini

◆ primordial_inflation_attractor_maxit

int precision::primordial_inflation_attractor_maxit

Initial value:

=

10

maximum number of iteration when searching attractor solution

◆ primordial_inflation_tol_curvature

double precision::primordial_inflation_tol_curvature

Initial value:

=

1.0e-3

for each k, stop following wavenumber, at the latest, when curvature perturbation R is stable up to to this tolerance

◆ primordial_inflation_aH_ini_target

double precision::primordial_inflation_aH_ini_target

Initial value:

=

0.9

control the step size in the search for a suitable initial field value

◆ primordial_inflation_end_dphi

double precision::primordial_inflation_end_dphi

Initial value:

=

1.0e-10

first bracketing width, when trying to bracket the value phi_end at which inflation ends naturally

◆ primordial_inflation_end_logstep

double precision::primordial_inflation_end_logstep

Initial value:

=

10.0

logarithmic step for updating the bracketing width, when trying to bracket the value phi_end at which inflation ends naturally

◆ primordial_inflation_small_epsilon

double precision::primordial_inflation_small_epsilon

Initial value:

=

0.1

value of slow-roll parameter epsilon used to define a field value phi_end close to the end of inflation (doesn't need to be exactly at the end): epsilon(phi_end)=small_epsilon (should be smaller than one)

◆ primordial_inflation_small_epsilon_tol

double precision::primordial_inflation_small_epsilon_tol = 0.01

tolerance in the search for phi_end

◆ primordial_inflation_extra_efolds

double precision::primordial_inflation_extra_efolds

Initial value:

=

2.0

a small number of efolds, irrelevant at the end, used in the search for the pivot scale (backward from the end of inflation)

◆ l_linstep

int precision::l_linstep

Initial value:

=

40

factor for logarithmic spacing of values of l over which bessel and transfer functions are sampled

◆ l_logstep

double precision::l_logstep

Initial value:

=

1.12

maximum spacing of values of l over which Bessel and transfer functions are sampled (so, spacing becomes linear instead of logarithmic at some point)

◆ hyper_x_min

double precision::hyper_x_min

Initial value:

=

1.0e-5

flat case: lower bound on the smallest value of x at which we sample $\Phi_l^{\nu}(x)$ or

◆ hyper_sampling_flat

double precision::hyper_sampling_flat

Initial value:

=

8.0

flat case: number of sampled points x per approximate wavelength $2\pi$ , should remain >7.5

◆ hyper_sampling_curved_low_nu

double precision::hyper_sampling_curved_low_nu

Initial value:

=

7.0

open/closed cases: number of sampled points x per approximate wavelength $2\pi/\nu$ , when $\nu$ smaller than hyper_nu_sampling_step

◆ hyper_sampling_curved_high_nu

double precision::hyper_sampling_curved_high_nu

Initial value:

=

3.0

open/closed cases: number of sampled points x per approximate wavelength $2\pi/\nu$ , when $\nu$ greater than hyper_nu_sampling_step

◆ hyper_nu_sampling_step

double precision::hyper_nu_sampling_step

Initial value:

=

1000.0

open/closed cases: value of nu at which sampling changes

◆ hyper_phi_min_abs

double precision::hyper_phi_min_abs

Initial value:

=

1.0e-10

small value of Bessel function used in calculation of first point x ( $\Phi_l^{\nu}(x)$ equals hyper_phi_min_abs)

◆ hyper_x_tol

double precision::hyper_x_tol = 1.0e-4

tolerance parameter used to determine first value of x

◆ hyper_flat_approximation_nu

double precision::hyper_flat_approximation_nu

Initial value:

=

4000.0

value of nu below which the flat approximation is used to compute Bessel function

◆ q_linstep

double precision::q_linstep = 0.45

asymptotic linear sampling step in q UNHANDLED: space, in units of $2\pi/r_a(\tau_rec)$ UNHANDLED: (comoving angular diameter distance to UNHANDLED: recombination), very important for CMB

◆ q_logstep_spline

double precision::q_logstep_spline = 170.0

initial logarithmic sampling step in q UNHANDLED: space, in units of $2\pi/r_a(\tau_{rec})$ UNHANDLED: (comoving angular diameter distance to UNHANDLED: recombination), very important for CMB and LSS

◆ q_logstep_open

double precision::q_logstep_open = 6.0

in open models, the value of UNHANDLED: q_logstep_spline must be decreased UNHANDLED: according to curvature. Increasing UNHANDLED: this number will make the calculation UNHANDLED: more accurate for large positive UNHANDLED: Omega_k

◆ q_logstep_trapzd

double precision::q_logstep_trapzd = 20.0

initial logarithmic sampling step in q UNHANDLED: space, in units of $2\pi/r_a(\tau_{rec})$ UNHANDLED: (comoving angular diameter distance to UNHANDLED: recombination), in the case of small UNHANDLED: q's in the closed case, for which one UNHANDLED: must used trapezoidal integration UNHANDLED: instead of spline (the number of q's UNHANDLED: for which this is the case decreases UNHANDLED: with curvature and vanishes in the UNHANDLED: flat limit)

◆ q_numstep_transition

double precision::q_numstep_transition = 250.0

number of steps for the transition UNHANDLED: from q_logstep_trapzd steps to UNHANDLED: q_logstep_spline steps (transition UNHANDLED: must be smooth for spline)

◆ q_logstep_limber

double precision::q_logstep_limber

Initial value:

=

1.025

logarithmic q-step ratio for the separate full-Limber CMB lensing grid

◆ k_max_limber_over_l_max_scalars

double precision::k_max_limber_over_l_max_scalars

Initial value:

=

0.001

full-Limber perturbation source cutoff k_max/l_max_scalars in 1/Mpc

◆ transfer_neglect_delta_k_S_t0

double precision::transfer_neglect_delta_k_S_t0

Initial value:

=

0.15

for temperature source function T0 of scalar mode, range of k values (in 1/Mpc) taken into account in transfer function: for l < (k-delta_k)*tau0, ie for k > (l/tau0 + delta_k), the transfer function is set to zero

◆ transfer_neglect_delta_k_S_t1

double precision::transfer_neglect_delta_k_S_t1

Initial value:

=

0.04

same for temperature source function T1 of scalar mode

◆ transfer_neglect_delta_k_S_t2

double precision::transfer_neglect_delta_k_S_t2

Initial value:

=

0.15

same for temperature source function T2 of scalar mode

◆ transfer_neglect_delta_k_S_e

double precision::transfer_neglect_delta_k_S_e

Initial value:

=

0.11

same for polarization source function E of scalar mode

◆ transfer_neglect_delta_k_V_t1

double precision::transfer_neglect_delta_k_V_t1

Initial value:

=

1.0

same for temperature source function T1 of vector mode

◆ transfer_neglect_delta_k_V_t2

double precision::transfer_neglect_delta_k_V_t2

Initial value:

=

1.0

same for temperature source function T2 of vector mode

◆ transfer_neglect_delta_k_V_e

double precision::transfer_neglect_delta_k_V_e

Initial value:

=

1.0

same for polarization source function E of vector mode

◆ transfer_neglect_delta_k_V_b

double precision::transfer_neglect_delta_k_V_b

Initial value:

=

1.0

same for polarization source function B of vector mode

◆ transfer_neglect_delta_k_T_t2

double precision::transfer_neglect_delta_k_T_t2

Initial value:

=

0.2

same for temperature source function T2 of tensor mode

◆ transfer_neglect_delta_k_T_e

double precision::transfer_neglect_delta_k_T_e

Initial value:

=

0.25

same for polarization source function E of tensor mode

◆ transfer_neglect_delta_k_T_b

double precision::transfer_neglect_delta_k_T_b

Initial value:

=

0.1

same for polarization source function B of tensor mode

◆ transfer_neglect_late_source

double precision::transfer_neglect_late_source

Initial value:

=

400.0

value of l below which the CMB source functions can be neglected at late time, excepted when there is a Late ISW contribution

◆ l_switch_limber

double precision::l_switch_limber

Initial value:

=

10.

when to use the Limber approximation for project gravitational potential cl's

◆ l_switch_limber_for_nc_local_over_z

double precision::l_switch_limber_for_nc_local_over_z

Initial value:

=

100.0

when to use the Limber approximation for local number count contributions to cl's (relative to central redshift of each bin)

◆ l_switch_limber_for_nc_los_over_z

double precision::l_switch_limber_for_nc_los_over_z

Initial value:

=

30.0

when to use the Limber approximation for number count contributions to cl's integrated along the line-of-sight (relative to central redshift of each bin)

◆ selection_cut_at_sigma

double precision::selection_cut_at_sigma

Initial value:

=

5.0

in sigma units, where to cut gaussian selection functions

◆ selection_sampling

double precision::selection_sampling

Initial value:

=

50.0

controls sampling of integral over time when selection functions vary quicker than Bessel functions. Increase for better sampling.

◆ selection_sampling_bessel

double precision::selection_sampling_bessel

Initial value:

=

20.0

controls sampling of integral over time when selection functions vary slower than Bessel functions. Increase for better sampling. IMPORTANT for lensed contributions.

◆ selection_sampling_bessel_los

double precision::selection_sampling_bessel_los

Initial value:

=

20.0

controls sampling of integral over time when selection functions vary slower than Bessel functions. This parameter is specific to number counts contributions to Cl integrated along the line of sight. Increase for better sampling

◆ selection_tophat_edge

double precision::selection_tophat_edge

Initial value:

=

0.1

controls how smooth are the edge of top-hat window function (<<1 for very sharp, 0.1 for sharp)

◆ sigma_k_per_decade

double precision::sigma_k_per_decade

Initial value:

=

80.

logarithmic stepsize controlling the precision of integrals for sigma(R,k) and similar quantitites

◆ nonlinear_min_k_max

double precision::nonlinear_min_k_max = 5.0

when UNHANDLED: using an algorithm to compute nonlinear UNHANDLED: corrections, like halofit or hmcode, UNHANDLED: k_max must be at least equal to this UNHANDLED: value. Calculations are done internally UNHANDLED: until this k_max, but the P(k,z) output UNHANDLED: is still controlled by P_k_max_1/Mpc or UNHANDLED: P_k_max_h/Mpc even if they are UNHANDLED: smaller

◆ halofit_min_k_nonlinear

double precision::halofit_min_k_nonlinear

Initial value:

=

1.0e-4

parameters relevant for HALOFIT computation value of k in 1/Mpc below which non-linear corrections will be neglected

◆ halofit_min_k_max

double precision::halofit_min_k_max = 5.0

DEPRECATED: should use instead nonlinear_min_k_max

◆ halofit_k_per_decade

double precision::halofit_k_per_decade = 80.0

halofit needs to evalute integrals UNHANDLED: (linear power spectrum times some UNHANDLED: kernels). They are sampled using UNHANDLED: this logarithmic step size.

◆ halofit_sigma_precision

double precision::halofit_sigma_precision = 0.05

a smaller value will lead to a UNHANDLED: more precise halofit result at the highest UNHANDLED: redshift at which halofit can make computations, UNHANDLED: at the expense of requiring a larger k_max; but UNHANDLED: this parameter is not relevant for the UNHANDLED: precision on P_nl(k,z) at other redshifts, so UNHANDLED: there is normally no need to change it

◆ halofit_tol_sigma

double precision::halofit_tol_sigma = 1.0e-6

tolerance required on sigma(R) when UNHANDLED: matching the condition sigma(R_nl)=1, UNHANDLED: whcih defines the wavenumber of UNHANDLED: non-linearity, k_nl=1./R_nl

◆ pk_eq_z_max

double precision::pk_eq_z_max = 5.0

Maximum z for the pk_eq method

◆ pk_eq_tol

double precision::pk_eq_tol = 1.0e-7

Tolerance on the pk_eq method for finding the pk

◆ hmcode_max_k_extra

double precision::hmcode_max_k_extra = 1.e6

Parameters relevant for HMcode computation parameter specifying the maximum k value for UNHANDLED: the extrapolation of the linear power spectrum UNHANDLED: (needed for the sigma computation)

◆ hmcode_min_k_max

double precision::hmcode_min_k_max = 5.

DEPRECATED: should use instead nonlinear_min_k_max

◆ hmcode_tol_sigma

double precision::hmcode_tol_sigma = 1.e-6

tolerance required on sigma(R) when matching the UNHANDLED: condition sigma(R_nl)=1, which defines the wavenumber UNHANDLED: of non-linearity, k_nl=1./R_nl

◆ n_hmcode_tables

int precision::n_hmcode_tables = 64

parameters controlling stepsize and min/max r & a values for sigma(r) & grow table

◆ nsteps_for_p1h_integral

int precision::nsteps_for_p1h_integral = 256

parameters controlling stepsize and min/max halomass values for the 1-halo-power integral

◆ accurate_lensing

int precision::accurate_lensing

Initial value:

=

0

switch between Gauss-Legendre quadrature integration and simple quadrature on a subdomain of angles

◆ num_mu_minus_lmax

int precision::num_mu_minus_lmax

Initial value:

=

70

difference between num_mu and l_max, increase for more precision

◆ delta_l_max

int precision::delta_l_max = 500

difference between l_max in unlensed and lensed spectra

◆ tol_gauss_legendre

double precision::tol_gauss_legendre

Initial value:

=

DBL_EPSILON

tolerance with which quadrature points are found: must be very small for an accurate integration (if not entered manually, set automatically to match machine precision)

◆ smallest_allowed_variation

double precision::smallest_allowed_variation

Initial value:

=

DBL_EPSILON

machine-dependent, assigned automatically by the code

The documentation for this struct was generated from the following files:

Data Fields

- general precision parameters

Detailed Description

Member Function Documentation

◆ ResolveDataPaths()

◆ parse()

Field Documentation

◆ a_ini_over_a_today_default

◆ back_integration_stepsize

◆ tol_background_integration

◆ tol_initial_Omega_r

◆ tol_M_ncdm

◆ tol_ncdm_synchronous

◆ tol_ncdm_newtonian

◆ tol_ncdm_bg

◆ tol_ncdm_initial_w

◆ tol_tau_eq

◆ Omega0_cdm_min_synchronous

◆ sBBN_file

◆ recfast_z_initial

◆ recfast_Nz0

◆ thermo_z_initial_idm_dr

◆ tol_thermo_integration

◆ recfast_Heswitch

◆ recfast_fudge_He

◆ recfast_Hswitch

◆ recfast_fudge_H

◆ recfast_delta_fudge_H

◆ recfast_AGauss1

◆ recfast_AGauss2

◆ recfast_zGauss1

◆ recfast_zGauss2

◆ recfast_wGauss1

◆ recfast_wGauss2

◆ recfast_z_He_1

◆ recfast_delta_z_He_1

◆ recfast_z_He_2

◆ recfast_delta_z_He_2

◆ recfast_z_He_3

◆ recfast_delta_z_He_3

◆ recfast_x_He0_trigger

◆ recfast_x_He0_trigger_delta

◆ recfast_x_H0_trigger

◆ recfast_x_H0_trigger2

◆ recfast_x_H0_trigger_delta

◆ recfast_H_frac

◆ reionization_z_start_max

◆ reionization_sampling

◆ reionization_optical_depth_tol

◆ reionization_start_factor

◆ thermo_rate_smoothing_radius

◆ hyrec_Alpha_inf_file

◆ hyrec_R_inf_file

◆ hyrec_two_photon_tables_file

◆ k_min_tau0

◆ k_max_tau0_over_l_max

◆ k_step_sub

◆ k_step_super

◆ k_step_transition

◆ k_step_super_reduction

◆ k_per_decade_for_pk

◆ idmdr_boost_k_per_decade_for_pk

◆ k_per_decade_for_bao

◆ k_bao_center

◆ k_bao_width

◆ start_small_k_at_tau_c_over_tau_h

◆ start_large_k_at_tau_h_over_tau_k

◆ tight_coupling_trigger_tau_c_over_tau_h

◆ tight_coupling_trigger_tau_c_over_tau_k

◆ start_sources_at_tau_c_over_tau_h

◆ tight_coupling_approximation

◆ idm_dr_tight_coupling_trigger_tau_c_over_tau_k

◆ idm_dr_tight_coupling_trigger_tau_c_over_tau_h

◆ idm_drmd_tight_coupling_trigger_G_over_aH

◆ l_max_g

◆ l_max_pol_g

◆ l_max_dr

◆ l_max_dr_col

◆ l_max_ur

◆ l_max_idr