omega_nu_single.h

Go to the documentation of this file.

#ifndef OMEGA_NU_SINGLE_H
#define OMEGA_NU_SINGLE_H
#include <gsl/gsl_interp.h>
 
#define TNUCMB     (pow(4/11.,1/3.)*1.00328)
#define NUSPECIES 3
 
struct _rho_nu_single {
    double * loga;
    double * rhonu;
    gsl_interp * interp;
    gsl_interp_accel * acc;
    /*Neutrino mass for this structure*/
    double mnu;
};
typedef struct _rho_nu_single _rho_nu_single;
 
void rho_nu_init(_rho_nu_single * rho_nu_tab, double a0, const double mnu, const double kBtnu);
 
double rho_nu(const _rho_nu_single * rho_nu_tab, const double a, const double kT);
 
struct _hybrid_nu {
    /*True if hybrid neutrinos are enabled*/
    int enabled;
    /* This is the fraction of neutrino mass not followed by the analytic integrator.
    The analytic method is cutoff at q < qcrit (specified using vcrit, below) and use
    particles for the slower neutrinos.*/
    double nufrac_low[NUSPECIES];
    /* Time at which to turn on the particle neutrinos.
     * Ultimately we want something better than this.*/
    double nu_crit_time;
    /*Critical velocity as a fraction of lightspeed*/
    double vcrit;
};
typedef struct _hybrid_nu _hybrid_nu;
 
void init_hybrid_nu(_hybrid_nu * const hybnu, const double mnu[], const double vcrit, const double light, const double nu_crit_time, const double kBtnu);
 
double particle_nu_fraction(const _hybrid_nu * const hybnu, const double a, int i);
 
double nufrac_low(const double qc);
 
struct _omega_nu {
    /*Pointers to the array of structures we use to store rho_nu*/
    _rho_nu_single RhoNuTab[NUSPECIES];
    /* Which species have the same mass and can thus be counted together.*/
    int nu_degeneracies[NUSPECIES];
    /* Prefactor to turn density into matter density omega*/
    double rhocrit;
    /*neutrino temperature times Boltzmann constant*/
    double kBtnu;
    /*CMB temperature*/
    double tcmb0;
    /* Pointer to structure for hybrid neutrinos. */
    _hybrid_nu hybnu;
};
typedef struct _omega_nu _omega_nu;
 
void init_omega_nu(_omega_nu * const omnu, const double MNu[], const double a0, const double HubbleParam, const double tcmb0);
 
double get_omega_nu(const _omega_nu * const omnu, const double a);
 
double get_omega_nu_nopart(const _omega_nu * const omnu, const double a);
 
double get_omegag(const _omega_nu * const omnu, const double a);
 
double omega_nu_single(const _omega_nu * const rho_nu_tab, const double a, const int i);
 
#endif