Fast Methods for Cosmological Simulations
FastSim serves as a tool for quick N-body simulations in modified gravity.
ccl_constants.h
Go to the documentation of this file.
1 #ifndef __CCL_CONSTANTS_H_INCLUDED__
2 #define __CCL_CONSTANTS_H_INCLUDED__
3 
4 #include <gsl/gsl_const_mksa.h>
5 
6 //Spline types
7 #define A_SPLINE_TYPE gsl_interp_akima
8 #define K_SPLINE_TYPE gsl_interp_akima
9 #define L_SPLINE_TYPE gsl_interp_akima
10 #define M_SPLINE_TYPE gsl_interp_akima
11 #define D_SPLINE_TYPE gsl_interp_akima
12 #define PNL_SPLINE_TYPE gsl_interp2d_bicubic
13 #define PLIN_SPLINE_TYPE gsl_interp2d_bicubic
14 #define CORR_SPLINE_TYPE gsl_interp_akima
15 
18 #ifndef M_PI
19 
22 #define M_PI 3.14159265358979323846
23 #endif
24 
28 #define K_PIVOT 0.05
29 
33 #define CLIGHT_HMPC 2997.92458 //H0^-1 in Mpc/h
34 
38 //#define GNEWT 6.6738e-11 /(from PDG 2013) in m^3/Kg/s^2
39 //#define GNEWT 6.67428e-11 // CLASS VALUE
40 #define GNEWT 6.67408e-11 // from CODATA 2014
41 
45 //#define SOLAR_MASS GSL_CONST_MKSA_SOLAR_MASS
46 //#define SOLAR_MASS 1.9885e30 //(from PDG 2015) in Kg
47 #define SOLAR_MASS 1.9884754153381438e+30 //from IAU 2015
48 
51 #define MPC_TO_METER 3.08567758149e22
52 
56 #define PC_TO_METER 3.08567758149e16
57 
61 #define RHO_CRITICAL ((3*100*100)/(8*M_PI*GNEWT)) * (1000*1000*MPC_TO_METER/SOLAR_MASS)
62 
66 //#define KBOLTZ GSL_CONST_MKSA_BOLTZMANN
67 #define KBOLTZ 1.38064852e-23 //from CODATA 2014
68 
72 //#define STBOLTZ GSL_CONST_MKSA_STEFAN_BOLTZMANN_CONSTANT
73 #define STBOLTZ 5.670367e-8 //from CODATA 2014
74 
77 //#define HPLANCK GSL_CONST_MKSA_PLANCKS_CONSTANT_H
78 #define HPLANCK 6.626070040e-34 //from CODATA 2014
79 
83 //#define CLIGHT GSL_CONST_MKSA_SPEED_OF_LIGHT
84 #define CLIGHT 299792458.0 //from CODATA 2014
85 
89 //#define EV_IN_J GSL_CONST_MKSA_ELECTRON_VOLT
90 #define EV_IN_J 1.6021766208e-19 //from CODATA 2014
91 
95 #define TCMB 2.725
96 //#define TCMB 2.7255 // CLASS value
97 
101 #define TNCDM 0.71611
102 
109 #define DELTAM12_sq 7.62E-5
110 #define DELTAM13_sq_pos 2.55E-3
111 #define DELTAM13_sq_neg -2.43E-3
112 
113 
114 //Precision parameters
118 #define GSL_EPSREL 1E-4
119 
124 #define GSL_N_ITERATION 1000
125 
130 #define GSL_INTEGRATION_GAUSS_KRONROD_POINTS GSL_INTEG_GAUSS41
131 
135 #define GSL_EPSABS_NU 1E-7
136 
140 #define GSL_EPSREL_NU 1E-7
141 
145 #define GSL_N_ITERATION_NU 1000
146 
150 #define GSL_EPSREL_SIGMAR 1E-5
151 
155 #define GSL_EPSREL_DIST 1E-6
156 
160 #define GSL_EPSREL_GROWTH 1E-6
161 
165 #define GSL_EPSREL_DNDZ 1E-6
166 
170 #define EPS_SCALEFAC_GROWTH 1E-6
171 
172 
173 #endif