00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011 class Cosmology {
00012
00013 private:
00014 double omega_m;
00015 double omega_x;
00016 double omega_k;
00017 double w_x;
00018 double h;
00020 public:
00021
00026 Cosmology(double o_m=0.260, double omega_x=0.740, \
00027 double omega_k=0.0, double w_x=-1.0, double h=0.72);
00031 ~Cosmology();
00032
00037 double D_C(double z1, double z2);
00038
00043 double D_M(double z1, double z2);
00044
00049 double ang_dist(double z);
00050
00055 double ang_dist_z1_z2(double z1, double z2);
00056
00061 double lum_dist(double z);
00062
00068 double unit_comoving_volume(double z, void *params);
00069
00074 double comoving_volume_z1_z2(double z1, double z2);
00075
00080 double comoving_volume_0_z(double z);
00081
00085 double lookback_time(double z);
00086
00090 double age(double z);
00091
00095 double age_now(void);
00096 };
00097
00098 double E_z(double x, void *params);
00099 double time_integral(double x, void *params);
00100
00106 double schechter(double L, double phi_star, double alpha, double L_star);
00107
00113 double number_schechter(double L, double phi_star, double alpha, double L_star);
00114
00120 double total_number_schechter(double phi_star, double alpha);
00121
00127 double luminosity_schechter(double L, double phi_star, double alpha, double L_star);
00128
00134 double total_luminosity_schechter(double phi_star, double alpha, double L_star);
00135
00142 double double_schechter(double L, double phi_star1, double alpha1, \
00143 double L_star1, double phi_star2, double alpha2, double L_star2);
00144
00150 double number_double_schechter(double L, double phi_star1, double alpha1, double L_star1, \
00151 double phi_star2, double alpha2, double L_star2);
00152
00158 double luminosity_double_schechter(double L, double phi_star1, double alpha1, double L_star1, \
00159 double phi_star2, double alpha2, double L_star2);
00160
00166 double total_luminosity_double_schechter(double phi_star1, double alpha1, double L_star1, \
00167 double phi_star2, double alpha2, double L_star2);