Cosmology.h

Go to the documentation of this file.

//
// File:        Cosmology.h
// Package:     C++ implementation of cosmology part for PyGL
// Copyright:   the GNU General Public License
// Description: Cosmology class and other routines supporting PyGL
//


class Cosmology {

        private:
                double omega_m; 
                double omega_x; 
                double omega_k; 
                double w_x; 
                double h; 
        public:

                Cosmology(double o_m=0.260, double omega_x=0.740, \
                double omega_k=0.0, double w_x=-1.0, double h=0.72);
                ~Cosmology();

                double D_C(double z1, double z2);

                double D_M(double z1, double z2);

                double ang_dist(double z);

                double ang_dist_z1_z2(double z1, double z2);

                double lum_dist(double z);

                double unit_comoving_volume(double z, void *params);

                double comoving_volume_z1_z2(double z1, double z2);

                double comoving_volume_0_z(double z);

                double lookback_time(double z);

                double age(double z);

                double age_now(void);

    double conformal_time(double z);
};

double E_z(double x, void *params);
double time_integral(double x, void *params);

double schechter(double L, double phi_star, double alpha, double L_star);

double number_schechter(double L, double phi_star, double alpha, double L_star);

double total_number_schechter(double phi_star, double alpha);

double luminosity_schechter(double L, double phi_star, double alpha, double L_star);

double total_luminosity_schechter(double phi_star, double alpha, double L_star);

double double_schechter(double L, double phi_star1, double alpha1, \
double L_star1, double phi_star2, double alpha2, double L_star2);

double number_double_schechter(double L, double phi_star1, double alpha1, double L_star1, \
double phi_star2, double alpha2, double L_star2);

double luminosity_double_schechter(double L, double phi_star1, double alpha1, double L_star1, \
double phi_star2, double alpha2, double L_star2);

double total_luminosity_double_schechter(double phi_star1, double alpha1, double L_star1, \
double phi_star2, double alpha2, double L_star2);