Doc/html/cosmo_8h_source.html

#include <iostream>

#include <vector>

#include <memory>

#include <cstdlib>

#include <cstddef>

#include <math.h>


#ifdef ENABLE_GSL

#include <gsl/gsl_integration.h>

#include <gsl/gsl_deriv.h>

#include <cmath>

#endif


#ifndef PI

#define PI  3.141593

#endif


#ifndef Grav

#define Grav 4.7788e-20  // Newton's constant over c^2 in Mpc / Msun

#endif


#ifndef lightspeed

#define lightspeed 2.99792458e5  // km / s

#endif


#ifndef error_message

#define error_message

#define ERROR_MESSAGE() std::cerr << "ERROR: file: " << __FILE__ << " line: " << __LINE__ << std::endl;

#endif


#define CRITD 2.49783e18    /* critical density / Ho^2 solar/Mpc */

#define CRITD2 2.7752543e11 /* critical density / h^2 M_sun/Mpc^3 */


#ifndef cosmo_declare


enum class CosmoParamSet {WMAP5yr,Millennium,Planck1yr,Planck15,Planck18,BigMultiDark,Uchuu,none};


std::string to_string(const CosmoParamSet &p);

std::ostream &operator<<(std::ostream &os,const CosmoParamSet &p);


class COSMOLOGY{

public:


  COSMOLOGY(CosmoParamSet cosmo_p = CosmoParamSet::WMAP5yr);

  COSMOLOGY(double omegam,double omegal,double h, double w = -1,double wa = 0

            ,bool justdistances=false

  );

  COSMOLOGY(const COSMOLOGY &cosmo);

  ~COSMOLOGY();


  COSMOLOGY& operator=(const COSMOLOGY &cosmo);


  // returns the parameter set if any

  CosmoParamSet ParamSet(){return cosmo_set;}

  void ParamSet(std::string &s);


  void PrintCosmology(short physical = 0) const;

  std::string print(short physical= 0) const;


  // Lengths

  double rcurve() const;

  //double emptyDist(double zo,double z);

  double coorDist(double zo,double z) const;

  double coorDist(double z) const {return coorDist(0,z);}

  double radDist(double zo,double z) const;

  double radDist(double z) const {return radDist(0,z);}

  double angDist(double zo,double z) const;

  double angDist(double z) const {return angDist(0,z);}

  double lumDist(double zo,double z) const;

  double lumDist(double z) const {return lumDist(0,z);}


  double d_coorDist_dOmo(double zo,double z) const;

  double d_coorDist_dw(double zo,double z) const;

  double d_coorDist_dw1(double zo,double z) const;


  double DRradius(double zo,double z,double pfrac);

  double DRradius2(double zo,double z);


  double invCoorDist(double d) const;

  double invRadDist(double d) const;

  double invComovingDist(double d) const;


  double scalefactor(double rad) const;

  double Omegam(double z) const;

  double rho_crit_comoving(double z) const;

  //double drdz_empty(double x);

  double drdz(double x) const;

  //double adrdz(double x) const;

  double drdz_dark(double x) const;

  inline double dark_factor(double x) const;

  //double adrdz_dark(double x) const;


  double delta_c() const;


  double DeltaVir(double z,int caseunit=0) const;

  double Deltao(double m) const;

  double time(double z) const;

  double nonlinMass(double z) const;


  // Stuff having to do with the power spectrum

  double power_normalize(double sigma8);

  double power_linear(double k,double z);

  double Dgrowth(double z) const;


  double powerCDMz(double k,double z);

  double psdfdm(double z,double m,int caseunit=0);

  double halo_bias (double m, double z, int t=0);

  double stdfdm(double z,double m,int caseunit=0);

  double powerlawdfdm(double z,double m,double alpha,int caseunit=0);

  double haloNumberDensity(double m,double z,double a, int t,double alpha = 0.0);


  double CorrelationFunction(double radius,double redshift

                      ,double k_max = 100,double k_min = 1.0e-3);


  struct CorrFunctorType{

    CorrFunctorType(COSMOLOGY *cosmo,double radius,double redshift)

    : cosmology(cosmo),z(redshift),r(radius)

    {

      norm = 0.5/PI/PI/(1+z)/(1+z);

    };


    COSMOLOGY *cosmology;

    double z;

    double r;


    double norm;


    double operator () (double k) {


      double rk = r*k;

      double jo = sin(rk)/rk;

      if(rk < 1.0e-3) jo = 1.0;


      return norm*jo*k*k*cosmology->powerCDMz(k,z);

    }

  };


  double haloNumberDensityOnSky (double m,double z1,double z2,int t,double alpha = 0.0);

  double haloNumberInBufferedCone (double mass ,double z1,double z2,double fov,double buffer ,int type ,double alpha=0.0);

  double haloMassInBufferedCone (double mass ,double z1,double z2,double fov,double buffer ,int type ,double alpha=0.0);


  double TopHatVariance(double m) const;

  double TopHatVarianceR(double R,double z);

  double TopHatVarianceM(double M,double z);

  //double gradius(double R,double rd);


  double getTimefromDeltaC(double dc);

  double getZfromDeltaC(double dc);


  void sethubble(double ht){ h = ht; TFmdm_set_cosm();

    power_normalize(sig8);  calc_interp_dist(); cosmo_set = CosmoParamSet::none;}


  double gethubble() const {return h;}

  double getHubble() const {return 100*h/lightspeed;}


  void setindex(double nn){ n = nn;  power_normalize(sig8); cosmo_set = CosmoParamSet::none;}

  double getindex() const { return n;}


  void setOmega_matter(double Omega_matter,bool FLAT = false){Omo = Omega_matter; if(FLAT) Oml = 1-Omo ; TFmdm_set_cosm(); power_normalize(sig8); calc_interp_dist(); cosmo_set = CosmoParamSet::none;}

  double getOmega_matter() const {return Omo;}


  void setOmega_lambda(double Omega_lambda,bool FLAT = false){Oml = Omega_lambda;  if(FLAT) Oml = 1-Omo ; TFmdm_set_cosm(); power_normalize(sig8); calc_interp_dist(); cosmo_set = CosmoParamSet::none;}

  double getOmega_lambda() const {return Oml;}


  void setOmega_baryon(double Omega_baryon){Omb = Omega_baryon; TFmdm_set_cosm(); power_normalize(sig8); calc_interp_dist(); cosmo_set = CosmoParamSet::none;

  }


  double getOmega_baryon() const {return Omb;}


  void setOmega_neutrino(double Omega_neutrino){Omnu = Omega_neutrino; TFmdm_set_cosm(); power_normalize(sig8);  calc_interp_dist(); cosmo_set = CosmoParamSet::none;

  }


  double getOmega_neutrino() const {return Omnu;}


  void setNneutrino(double Nneutrino){Nnu = Nneutrino; TFmdm_set_cosm(); power_normalize(sig8); cosmo_set = CosmoParamSet::none;}

  double getNneutrino() const {return Nnu;}


  void setW(double w){ww = w;  calc_interp_dist(); cosmo_set = CosmoParamSet::none;}

  double getW() const {return ww;}

  void setW1(double w){ww1 = w;  calc_interp_dist(); cosmo_set = CosmoParamSet::none;}

  double getW1() const {return ww1;}


  void setdndlnk(double w){dndlnk = w; power_normalize(sig8); cosmo_set = CosmoParamSet::none;}

  double getdndlnk() const {return dndlnk;}


  void setgamma(double gamm){gamma = gamm; cosmo_set = CosmoParamSet::none;}

  double getgamma() const {return gamma;}

  // If physical = 1 all Omega are Omega*h^2 otherwise they have the usual definitions.

  // 2 gamma parameterization is used for dark energy

  void setDEtype(short tt){darkenergy = tt; cosmo_set = CosmoParamSet::none;}

  short getDEtype() const {return darkenergy;}


  void setSigma8(double my_sig8){power_normalize(my_sig8); cosmo_set = CosmoParamSet::none;}

  double getSigma8() const {return sig8;}


  void dzdangDist(double D,double z[],double dzdD[]);


  double totalMassDensityinHalos(int t,double alpha,double m_min,double z,double z1,double z2);


  void setInterpolation(double z_interp, std::size_t n_interp);


  double SigmaCrit(double zlens,double zsource) const;


protected:

  void SetConcordenceCosmology(CosmoParamSet cosmo_p);


  CosmoParamSet cosmo_set;


  // structure rappers to make integration thread safe

  struct drdz_struct{

    drdz_struct(COSMOLOGY const &cosmo):cos(cosmo){};

    double operator()(double x) const {return cos.drdz(x);}

    COSMOLOGY const &cos;

  };

  struct drdzdark_struct{

    drdzdark_struct(COSMOLOGY const &cosmo):cos(cosmo){};

    double operator()(double x) const {return cos.drdz_dark(x);}

    COSMOLOGY const &cos;

  };

  struct adrdz_struct{

    adrdz_struct(COSMOLOGY const &cosmo):cos(cosmo){};

    double operator()(double x) const {return cos.drdz(x)/x;}

    COSMOLOGY const &cos;

  };

  struct adrdzdark_struct{

    adrdzdark_struct(COSMOLOGY const &cosmo):cos(cosmo){};

    double operator()(double x) const {return cos.drdz_dark(x)/x;}

    COSMOLOGY const &cos;

  };


  struct normL_struct{

    normL_struct(COSMOLOGY &cosmo):cos(cosmo){};

    double operator()(double x) const {return cos.normL(x);}

    COSMOLOGY &cos;

  };


  double ddrdzdOmo(double x) const;

  double ddrdzdw(double x) const;

  double ddrdzdw1(double x) const;


  struct ddrdzdOmo_struct{

    ddrdzdOmo_struct(COSMOLOGY const &cosmo):cos(cosmo){};

    double operator()(double x) const {return cos.ddrdzdOmo(x);}

    COSMOLOGY const &cos;

  };

  struct ddrdzdw_struct{

    ddrdzdw_struct(COSMOLOGY const &cosmo):cos(cosmo){};

    double operator()(double x) const {return cos.ddrdzdw(x);}

    COSMOLOGY const &cos;

  };

  struct ddrdzdw1_struct{

    ddrdzdw1_struct(COSMOLOGY const &cosmo):cos(cosmo){};

    double operator()(double x) const {return cos.ddrdzdw1(x);}

    COSMOLOGY const &cos;

  };


  bool init_structure_functions;


  double h;

  double n;

  double Omo;

  double Oml;

  double Omb;

  double Omnu;

  double Nnu;

  double ww;

  double ww1;

  double dndlnk;


  double gamma;

  // If physical = 1 all Omega are Omega*h^2 otherwise they have the usual definitions.

  //short physical;

  // 2 gamma parameterization is used for dark energy

  short darkenergy;


  /* table for growth parameter */

  //std::auto_ptr<double> *aa;

  //std::auto_ptr<double> *growth;

  //int Ntable;


  double A;

  double sig8;  /* do not access these normalization outside */


  double Rtophat;

  double ztmp;


  double powerEH(double k,double z);

  double powerEHv2(double k);

  double powerloc(double k,double z);

  double npow(double k);

  double normL(double lgk);

  double gradius(double R,double rd) const;


  int nbin;  // number of points when binning to interpolate

  std:: vector<double> vDeltaCz,vlz,vt;

  double DpropDz(double z);

  double dsigdM(double m);

  double timeEarly(double a) const;

  double dNdz(double z);


  typedef double (COSMOLOGY::*pt2MemFunc)(double) const;

  typedef double (COSMOLOGY::*pt2MemFuncNonConst)(double);


  double nintegrateDcos(pt2MemFunc func,double a,double b,double tols) const;

  double trapzdDcoslocal(pt2MemFunc func, double a, double b, int n, double *s2) const;

  void polintD(double xa[], double ya[], int n, double x, double *y, double *dy) const;

  double nintegrateDcos(pt2MemFuncNonConst func,double a,double b,double tols);

  double trapzdDcoslocal(pt2MemFuncNonConst func, double a, double b, int n, double *s2);

  double dfridrDcos(pt2MemFunc func, double x, double h, double *err);

  double f4(double u) const;


  int ni;

  std::vector<float> xf, wf;


  // The parameters used in Eisenstein & Hu power spectrum

  /* in powerEH.c */

  short TFmdm_set_cosm_change_z(double redshift);

  short TFmdm_set_cosm();

  double TFmdm_onek_mpc(double kk);

  double TFmdm_onek_hmpc(double kk);


  double f_baryon; // Baryon fraction

  double f_bnu;  // Baryon + Massive Neutrino fraction

  double f_cb;  // Baryon + CDM fraction

  double f_cdm;  // CDM fraction

  double f_hdm;  // Massive Neutrino fraction


  double alpha_gamma;  // sqrt(alpha_nu)

  double alpha_nu;      // The small-scale suppression

  double beta_c;   // The correction to the log in the small-scale

  double num_degen_hdm;  // Number of degenerate massive neutrino species

  double growth_k0;      // D_1(z) -- the growth function as k->0

  double growth_to_z0;  // D_1(z)/D_1(0) -- the growth relative to z=0

  double k_equality;  // The comoving wave number of the horizon at equality

  double obhh;       // Omega_baryon * hubble^2

  double omega_curv;  // = 1 - omega_matter - omega_lambda

  double omega_lambda_z; // Omega_lambda at the given redshift

  double omega_matter_z; // Omega_matter at the given redshift

  double omhh;       // Omega_matter * hubble^2

  double onhh;       // Omega_hdm * hubble^2

  double p_c;       // The correction to the exponent before drag epoch

  double p_cb;       // The correction to the exponent after drag epoch

  double sound_horizon_fit;  // The sound horizon at the drag epoch

  double theta_cmb;      // The temperature of the CMB, in units of 2.7 K

  double y_drag;   // Ratio of z_equality to z_drag

  double z_drag;   // Redshift of the drag epoch

  double z_equality;  // Redshift of matter-radiation equality


  void setinternals();

  // in powerEHv2.c

  void TFset_parameters(double omega0hh, double f_baryon, double Tcmb);

  double TFfit_onek(double k, double *tf_baryon, double *tf_cdm);


  double R_drag;  // Photon-baryon ratio at drag epoch

  double R_equality; // Photon-baryon ratio at equality epoch

  double sound_horizon; // Sound horizon at drag epoch, in Mpc

  double k_silk;  // Silk damping scale, in Mpc^-1

  double alpha_c;     // CDM suppression

  double alpha_b;     // Baryon suppression

  double beta_b;  // Baryon envelope shift

  double beta_node; // Sound horizon shift

  double k_peak;  // Fit to wavenumber of first peak, in Mpc^-1


  // temporary variables for doing interations

  int tmp_type;

  double tmp_alpha;

  double tmp_mass;

  double tmp_a;


  // interpolation of functions

  double z_interp;

  std::size_t n_interp;

  void calc_interp_dist();

  double interp(const std::vector<double>& table, double z) const;

  double invert(const std::vector<double>& table, double f_z) const;

  std::vector<double> redshift_interp;

  std::vector<double> coorDist_interp;

  std::vector<double> radDist_interp;

};


typedef COSMOLOGY *CosmoHndl;


class NFW_Utility {

public:

  NFW_Utility(){}

  ~NFW_Utility(){};


  // methods for NFW profile

  double NFW_V200(double M200,double R200);

  double NFW_Vmax(double cons,double M200,double R200);

  double NFW_Vr(double x,double cons,double M200,double R200);

  double NFW_M(double x,double cons,double M200);

  double NFW_deltac(double cons);

  double NFW_Concentration(double Vmax,double M200,double R200);

  double NFW_rho(double cons,double x);


  void match_nfw(float Vmax,float R_half,float mass,float *cons,float *Rsize);

  float Rsize(float cons,float Vmax,float mass);

  float g_func(float x);


private:

  float Vmax,R_half,mass;


  typedef float (NFW_Utility::*MemFunc)(float);

  float zbrentD(MemFunc func, float a,float b,float tols);

  float nfwfunc(float cons);

  float funcforconcentration(float cons);

};


double drdz_wrapper(double x, void *params);

double drdz_dark_wrapper(double x, void *params);

double Deltao_wrapper(double m, void *params);


#define cosmo_declare

#endif


/*** in cosmo.c ***/

int cosmo_compare(CosmoHndl cos1,CosmoHndl cos2);

void cosmo_copy(CosmoHndl cos1,CosmoHndl cos2);

void ders(double z,double Da[],double dDdz[]);

void dir(double r,double a[],double dadr[]);

double arctanh(double x);

double fmini(double a,double b);

double fmaxi(double a,double b);

double **dmatrixcos(long nrl, long nrh, long ncl, long nch);

void free_dmatrixcos(double **m, long nrl, long nrh, long ncl, long nch);


/* in nfw.c */


/* in powerEHv2.c */

void TFset_parameters(double omega0hh, double f_baryon, double Tcmb);

double TFfit_onek(double k, double *tf_baryon, double *tf_cdm);

COSMOLOGY
The cosmology and all the functions required to calculated quantities based on the cosmology.
Definition cosmo.h:52

COSMOLOGY::DeltaVir
double DeltaVir(double z, int caseunit=0) const
Virial overdensity.
Definition cosmo.cpp:1092

COSMOLOGY::setOmega_neutrino
void setOmega_neutrino(double Omega_neutrino)
Omega neutrino, renormalizes P(k) to keep sig8 fixed.
Definition cosmo.h:216

COSMOLOGY::stdfdm
double stdfdm(double z, double m, int caseunit=0)
Sheth-Tormen mass function.
Definition cosmo.cpp:803

COSMOLOGY::setInterpolation
void setInterpolation(double z_interp, std::size_t n_interp)
set interpolation range and number of points
Definition cosmo.cpp:1465

COSMOLOGY::haloNumberDensityOnSky
double haloNumberDensityOnSky(double m, double z1, double z2, int t, double alpha=0.0)
Number of halos with mass larger than m (in solar masses) between redshifts z1 and z2 per square degr...
Definition cosmo.cpp:953

COSMOLOGY::powerloc
double powerloc(double k, double z)
The linear power spectrum without growth factor growth factor should be normalized to 1 at z=0.
Definition powerCDM.cpp:144

COSMOLOGY::radDist
double radDist(double zo, double z) const
Non-comoving radial distance in units Mpc also known as the lookback time. This is coorDist only inte...
Definition cosmo.cpp:687

COSMOLOGY::npow
double npow(double k)
Logorithmic slope of the power spectrum.
Definition powerCDM.cpp:124

COSMOLOGY::d_coorDist_dOmo
double d_coorDist_dOmo(double zo, double z) const
derivative of coordinate distance with respect to Omo in flat cosmology. Useful for Fisher matrix cal...
Definition cosmo.cpp:672

COSMOLOGY::setNneutrino
void setNneutrino(double Nneutrino)
Number of neutrino species, renormalizes P(k) to keep sig8 fixed.
Definition cosmo.h:221

COSMOLOGY::rcurve
double rcurve() const
Curvature radius in Mpc.
Definition cosmo.cpp:461

COSMOLOGY::Nnu
double Nnu
Number of neutrino species.
Definition cosmo.h:324

COSMOLOGY::print
std::string print(short physical=0) const
Print cosmological parameters.
Definition cosmo.cpp:390

COSMOLOGY::setOmega_matter
void setOmega_matter(double Omega_matter, bool FLAT=false)
Omega matter, renormalizes P(k) to keep sig8 fixed.
Definition cosmo.h:203

COSMOLOGY::PrintCosmology
void PrintCosmology(short physical=0) const
Print cosmological parameters.
Definition cosmo.cpp:359

COSMOLOGY::SetConcordenceCosmology
void SetConcordenceCosmology(CosmoParamSet cosmo_p)
Sets cosmology to WMAP 2009 model. This is done automatically in the constructor.
Definition cosmo.cpp:198

COSMOLOGY::powerlawdfdm
double powerlawdfdm(double z, double m, double alpha, int caseunit=0)
Power-law mass function.
Definition cosmo.cpp:846

COSMOLOGY::gradius
double gradius(double R, double rd) const
Incorporates curvature for angular size distance.
Definition cosmo.cpp:755

COSMOLOGY::powerEH
double powerEH(double k, double z)
The power spectrum from Eisinstein & Hu with neutrinos but no BAO.
Definition powerCDM.cpp:218

COSMOLOGY::invComovingDist
double invComovingDist(double d) const
The inverse of the angular size distance in units Mpc, works within interpolation range.
Definition cosmo.cpp:733

COSMOLOGY::psdfdm
double psdfdm(double z, double m, int caseunit=0)
Press-Schechter halo mass function - default  i.e. fraction of mass or if caseunit==1 the .
Definition cosmo.cpp:770

COSMOLOGY::power_normalize
double power_normalize(double sigma8)
Set the linear normalization for the power spectrum.
Definition powerCDM.cpp:155

COSMOLOGY::scalefactor
double scalefactor(double rad) const
The scale factor, a = 1/(1+z), as a function of radius in Mpc.
Definition powerCDM.cpp:84

COSMOLOGY::delta_c
double delta_c() const
the critical over density
Definition cosmo.cpp:1242

COSMOLOGY::Omnu
double Omnu
Omega neutrino.
Definition cosmo.h:322

COSMOLOGY::getHubble
double getHubble() const
Hubble parameter in 1/Mpc units.
Definition cosmo.h:196

COSMOLOGY::SigmaCrit
double SigmaCrit(double zlens, double zsource) const
The lensing critical density in Msun / Mpc^2.
Definition cosmo.cpp:1540

COSMOLOGY::invCoorDist
double invCoorDist(double d) const
The inverse of the coordinate distance in units Mpc, returning redshift. It works within interpolatio...
Definition cosmo.cpp:726

COSMOLOGY::DRradius
double DRradius(double zo, double z, double pfrac)
Comoving Dyer-Roeder angular size distance for lambda=0.
Definition cosmo.cpp:485

COSMOLOGY::setW
void setW(double w)
Dark energy equation of state parameter p/rho = w + w_1 (1+z)
Definition cosmo.h:225

COSMOLOGY::setgamma
void setgamma(double gamm)
Alternative to w for dark energy/ alt. grav. structure evolution.
Definition cosmo.h:235

COSMOLOGY::rho_crit_comoving
double rho_crit_comoving(double z) const
comoving critical density in M_sun/Mpc^3
Definition cosmo.cpp:607

COSMOLOGY::gamma
double gamma
Alternative to w for dark energy/ alt. grav. structure evolution.
Definition cosmo.h:333

COSMOLOGY::d_coorDist_dw1
double d_coorDist_dw1(double zo, double z) const
derivative of coordinate distance with respect to w1. Useful for Fisher matrix calculations
Definition cosmo.cpp:680

COSMOLOGY::lumDist
double lumDist(double zo, double z) const
The bolometric luminosity distance in units Mpc.
Definition cosmo.cpp:719

COSMOLOGY::dndlnk
double dndlnk
Running of primordial spectral index, P(k)_primordial \propto pow(k/h,n+dndlnk*log(k))
Definition cosmo.h:330

COSMOLOGY::Omegam
double Omegam(double z) const
Matter density parameter at redshift z.
Definition cosmo.cpp:478

COSMOLOGY::setOmega_lambda
void setOmega_lambda(double Omega_lambda, bool FLAT=false)
Omega lambda, renormalizes P(k) to keep sig8 fixed.
Definition cosmo.h:207

COSMOLOGY::dsigdM
double dsigdM(double m)
Derivative of  in CDM model.
Definition cosmo.cpp:1073

COSMOLOGY::halo_bias
double halo_bias(double m, double z, int t=0)
Halo bias, uses formalism by Mo-White.
Definition cosmo.cpp:1271

COSMOLOGY::setindex
void setindex(double nn)
Primordial spectral index, renormalizes P(k) to keep sig8 fixed.
Definition cosmo.h:199

COSMOLOGY::ww
double ww
Dark energy equation of state parameter p/rho = ww + ww_1 (1+z)
Definition cosmo.h:326

COSMOLOGY::Dgrowth
double Dgrowth(double z) const
Linear growth factor normalized to 1 at z=0.
Definition cosmo.cpp:588

COSMOLOGY::Oml
double Oml
Omega lambda.
Definition cosmo.h:318

COSMOLOGY::coorDist
double coorDist(double zo, double z) const
The coordinate distance in units Mpc. This is the radial distance found by integrating 1/H(z)....
Definition cosmo.cpp:657

COSMOLOGY::time
double time(double z) const
Return the time from the Big Bang in Gyr at a given redshift z.
Definition cosmo.cpp:1182

COSMOLOGY::TopHatVariance
double TopHatVariance(double m) const
Mass variance . This uses a fitting formula for the CDM model which might not be perfectly accurate....
Definition cosmo.cpp:1062

COSMOLOGY::angDist
double angDist(double zo, double z) const
The angular size distance in units Mpc.
Definition cosmo.cpp:704

COSMOLOGY::h
double h
Hubble paremters in units of 100 km/s/Mpc.
Definition cosmo.h:312

COSMOLOGY::setdndlnk
void setdndlnk(double w)
Running of primordial spectral index, P(k)_primordial \propto pow(k/h,n+dndlnk*log(k)),...
Definition cosmo.h:231

COSMOLOGY::invRadDist
double invRadDist(double d) const
The inverse of the coordinate distance in units Mpc, works within interpolation range.
Definition cosmo.cpp:747

COSMOLOGY::Omo
double Omo
Omega matter, dark matter + baryons.
Definition cosmo.h:316

COSMOLOGY::haloNumberInBufferedCone
double haloNumberInBufferedCone(double mass, double z1, double z2, double fov, double buffer, int type, double alpha=0.0)
The number of halos in a buffered cone between redshift z1 and z2.
Definition cosmo.cpp:994

COSMOLOGY::haloMassInBufferedCone
double haloMassInBufferedCone(double mass, double z1, double z2, double fov, double buffer, int type, double alpha=0.0)
Total mass contained in halos in a buffered cone between redshift z1 and z2.
Definition cosmo.cpp:1025

COSMOLOGY::TopHatVarianceR
double TopHatVarianceR(double R, double z)
Variance within a spherical top-hat filter of size R (Mpc), , at redshift z.
Definition powerCDM.cpp:188

COSMOLOGY::DRradius2
double DRradius2(double zo, double z)
Comoving Dyer-Roeder angular size distance for lambda=0 and pfrac = 1 (all matter in particles)
Definition cosmo.cpp:525

COSMOLOGY::power_linear
double power_linear(double k, double z)
Linear power spectrum P(k,z)/a^2.
Definition powerCDM.cpp:133

COSMOLOGY::sethubble
void sethubble(double ht)
accesser functions:
Definition cosmo.h:192

COSMOLOGY::setOmega_baryon
void setOmega_baryon(double Omega_baryon)
Omega baryon, renormalizes P(k) to keep sig8 fixed.
Definition cosmo.h:211

COSMOLOGY::getZfromDeltaC
double getZfromDeltaC(double dc)
Return the redshift given .
Definition cosmo.cpp:1124

COSMOLOGY::powerEHv2
double powerEHv2(double k)
This is the power spectrum from Eisinstein & Hu with BAO but no neutrinos
Definition powerCDM.cpp:235

COSMOLOGY::Omb
double Omb
Omega baryon.
Definition cosmo.h:320

COSMOLOGY::d_coorDist_dw
double d_coorDist_dw(double zo, double z) const
derivative of coordinate distance with respect to w. Useful for Fisher matrix calculations
Definition cosmo.cpp:676

COSMOLOGY::CorrelationFunction
double CorrelationFunction(double radius, double redshift, double k_max=100, double k_min=1.0e-3)
Dark matter correlation function.
Definition powerCDM.cpp:246

COSMOLOGY::totalMassDensityinHalos
double totalMassDensityinHalos(int t, double alpha, double m_min, double z, double z1, double z2)
The halo total surface mass density in haloes with mass larger than m_min (solar masses) in the redsh...
Definition cosmo.cpp:915

COSMOLOGY::n
double n
Primordial spectral index.
Definition cosmo.h:314

COSMOLOGY::getTimefromDeltaC
double getTimefromDeltaC(double dc)
Return the time from the Big Bang in Gyr given .
Definition cosmo.cpp:1151

COSMOLOGY::ww1
double ww1
Dark energy equation of state parameter p/rho = ww + ww_1 (1+z)
Definition cosmo.h:328

COSMOLOGY::powerCDMz
double powerCDMz(double k, double z)
powerCDM.c calculates the nonlinear P(k,z)/a(r)^2
Definition powerCDM.cpp:10

COSMOLOGY::TopHatVarianceM
double TopHatVarianceM(double M, double z)
Variance within a spherical top-hat filter at mass scale M at redshift z.
Definition powerCDM.cpp:207

COSMOLOGY::haloNumberDensity
double haloNumberDensity(double m, double z, double a, int t, double alpha=0.0)
The cumulative comoving number density of haloes with mass larger than m (in solar masses) at redshif...
Definition cosmo.cpp:874

NFW_Utility
Class for calculating properties of NFW halo profile.
Definition cosmo.h:451

NFW_Utility::NFW_M
double NFW_M(double x, double cons, double M200)
Mass within a radius x in Msun.
Definition nfw.cpp:53

NFW_Utility::match_nfw
void match_nfw(float Vmax, float R_half, float mass, float *cons, float *Rsize)
Returns the concentration and radius of an NFW halo with the mass, half mass radius and Vmax provided...
Definition nfw.cpp:110

NFW_Utility::NFW_rho
double NFW_rho(double cons, double x)
The density of an NFW profile in units of the critical density.
Definition nfw.cpp:101

NFW_Utility::NFW_Vr
double NFW_Vr(double x, double cons, double M200, double R200)
Circular velocity in km/s.
Definition nfw.cpp:42

NFW_Utility::NFW_Concentration
double NFW_Concentration(double Vmax, double M200, double R200)
Concentration of NFW given Vmax in km/s.
Definition nfw.cpp:69

NFW_Utility::NFW_Vmax
double NFW_Vmax(double cons, double M200, double R200)
Maximum circular velocity in km/s.
Definition nfw.cpp:32

NFW_Utility::NFW_V200
double NFW_V200(double M200, double R200)
Circular velocity at R200 in km/s.
Definition nfw.cpp:25

NFW_Utility::NFW_deltac
double NFW_deltac(double cons)
central over-density of nfw halo
Definition nfw.cpp:63