CCL
ccl_neutrinos.h
Go to the documentation of this file.
1 
2 #ifdef __cplusplus
3 extern "C" {
4 #endif
5 
6 //#pragma once
7 #include "ccl_core.h"
8 #include "gsl/gsl_spline.h"
9 #include "gsl/gsl_const_mksa.h"
10 
11 // maximum number of species
12 #define CCL_MAX_NU_SPECIES 3
13 // limits for the precomputed spline of phase
14 // space diagram in MNU/T
15 #define CCL_NU_MNUT_MIN 1e-4
16 #define CCL_NU_MNUT_MAX 500
17 // and number of points
18 #define CCL_NU_MNUT_N 1000
19 
20 // The combination of constants required in Omeganuh2
21 #define NU_CONST (8. * pow(M_PI,5) *pow((KBOLTZ/ HPLANCK),3)* KBOLTZ/(15. *pow( CLIGHT,3))* (8. * M_PI * GNEWT) / (3. * 100.*100.*1000.*1000. /MPC_TO_METER /MPC_TO_METER * CLIGHT * CLIGHT))
22 
30 gsl_spline* calculate_nu_phasespace_spline(int *status);
31 
44 double Omeganuh2 (double a, double Neff, double mnu, double TCMB, gsl_interp_accel* accel, int * status);
45 
46 #ifdef __cplusplus
47 }
48 #endif
gsl_spline * calculate_nu_phasespace_spline(int *status)
Definition: ccl_neutrinos.c:28
double Omeganuh2(double a, double Neff, double mnu, double TCMB, gsl_interp_accel *accel, int *status)
Definition: ccl_neutrinos.c:96