-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathgamma_t.h
88 lines (77 loc) · 1.64 KB
/
gamma_t.h
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
struct info_for_hankel
{
int max_sum;
double *zero;
double *phi;
double *phidash;
double *weights;
double *besselterm;
};
struct cosmo_params
{
double Om;
double Omb;
double w;
double ns;
double h;
double sigma8;
double redshift;
double growth_factor;
};
struct HOD_params
{
float Mcut;
float M1;
float sigma;
float alpha;
float kappa;
};
struct survey_info
{
int nbins_s;
double *z_s;
double *chi_s;
double *n_z_s;
double z_mean_s;
int nbins_l;
double *z_l;
double *chi_l;
double *n_z_l;
double z_mean_l;
};
struct powerspectra
{
int nk;
double *k;
double *Pk;
};
struct info_for_source_dist_integrand
{
int npts;
double *integrand;
double chi;
double z_l;
};
struct info_for_lens_dist_integrand
{
double *theta; // tracks which value of theta we are up to
int npts; // number of bins for the lens distribution
int nchi; // total number of points evaluated for the radial distance
int ntheta; // total number of points evaluated for the angular scale
double last; // last value of chi at which source distribution is non-zero
double z_l;
double *chi;
double *chi_table; // this is for the power spectrum only
double **pk_table; // holds results of int dk k P(k) J2
double *f_chi_table; // hols results of collective factor n_l*f_K/a
double *n_l; // lens distribution
gsl_spline *spline_f_K; // spline for n_l*f_K/a
gsl_interp_accel *acc_f_K;
gsl_spline *spline_pkJ; // spline for int dk k Pk J2
gsl_interp_accel *acc_pkJ;
};
struct info_for_sigma_crit_integrand
{
struct cosmo_params cosmo;
struct survey_info survey;
};