neutrinos_lra.h

Go to the documentation of this file.

#ifndef NEUTRINOS_LRA_H
#define NEUTRINOS_LRA_H
 
#include <bigfile-mpi.h>
#include "powerspectrum.h"
#include "cosmology.h"
 
struct _delta_tot_table {
    int nk;
    int nk_allocated;
    int namax;
    int ia;
    double delta_nu_prefac;
    int delta_tot_init_done;
    double **delta_tot;
    double * scalefact;
    double * delta_nu_init;
    double * delta_nu_last;
    double * wavenum;
    const _omega_nu * omnu;
    double Omeganonu;
    double light;
    double TimeTransfer;
};
typedef struct _delta_tot_table _delta_tot_table;
 
 
void init_neutrinos_lra(const int nk_in, const double TimeTransfer, const double TimeMax, const double Omega0, const _omega_nu * const omnu, const double UnitTime_in_s, const double UnitLength_in_cm);
 
/*Computes delta_nu from a CDM power spectrum.*/
void delta_nu_from_power(struct _powerspectrum * PowerSpectrum, Cosmology * CP, const double Time, const double TimeIC);
 
/*These functions save and load neutrino related data from the snapshots*/
void petaio_save_neutrinos(BigFile * bf, int ThisTask);
void petaio_read_neutrinos(BigFile * bf, int ThisTask);
/*Loads from the ICs*/
void petaio_read_icnutransfer(BigFile * bf, int ThisTask);
 
/*Save the neutrino power spectrum to a file*/
void powerspectrum_nu_save(struct _powerspectrum * PowerSpectrum, const char * OutputDir, const char * filename, const double Time);
#endif