mirror of
https://github.com/cosmo-sims/monofonIC.git
synced 2024-09-19 17:03:45 +02:00
49 lines
No EOL
1.9 KiB
C++
49 lines
No EOL
1.9 KiB
C++
#pragma once
|
|
|
|
#include <config_file.hh>
|
|
|
|
//! structure for cosmological parameters
|
|
struct CosmologyParameters
|
|
{
|
|
double
|
|
Omega_m, //!< baryon+dark matter density
|
|
Omega_b, //!< baryon matter density
|
|
Omega_DE, //!< dark energy density (cosmological constant or parameterised)
|
|
Omega_r, //!< photon + relativistic particle density
|
|
Omega_k, //!< curvature density
|
|
H0, //!< Hubble constant in km/s/Mpc
|
|
nspect, //!< long-wave spectral index (scale free is nspect=1)
|
|
sigma8, //!< power spectrum normalization
|
|
w_0, //!< dark energy equation of state parameter 1: w = w0 + a * wa
|
|
w_a, //!< dark energy equation of state parameter 2: w = w0 + a * wa
|
|
|
|
// below are helpers to store additional information
|
|
dplus, //!< linear perturbation growth factor
|
|
pnorm, //!< actual power spectrum normalisation factor
|
|
sqrtpnorm, //!< sqrt of power spectrum normalisation factor
|
|
vfact; //!< velocity<->displacement conversion factor in Zel'dovich approx.
|
|
|
|
explicit CosmologyParameters(ConfigFile cf)
|
|
{
|
|
Omega_b = cf.GetValue<double>("cosmology", "Omega_b");
|
|
Omega_m = cf.GetValue<double>("cosmology", "Omega_m");
|
|
Omega_DE = cf.GetValue<double>("cosmology", "Omega_L");
|
|
w_0 = cf.GetValueSafe<double>("cosmology", "w0", -1.0);
|
|
w_a = cf.GetValueSafe<double>("cosmology", "wa", 0.0);
|
|
|
|
Omega_r = cf.GetValueSafe<double>("cosmology", "Omega_r", 0.0); // no longer default to nonzero (8.3e-5)
|
|
Omega_k = 1.0 - Omega_m - Omega_DE - Omega_r;
|
|
|
|
H0 = cf.GetValue<double>("cosmology", "H0");
|
|
sigma8 = cf.GetValue<double>("cosmology", "sigma_8");
|
|
nspect = cf.GetValue<double>("cosmology", "nspec");
|
|
|
|
dplus = 0.0;
|
|
pnorm = 0.0;
|
|
vfact = 0.0;
|
|
}
|
|
|
|
CosmologyParameters(void)
|
|
{
|
|
}
|
|
}; |