forked from rajeshmondal18/ReionYuga
-
Notifications
You must be signed in to change notification settings - Fork 0
/
nbody.h
68 lines (59 loc) · 2.81 KB
/
nbody.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
/* in allotarrays.c */
float **allocate_float_2d(long N1,int N2);
float ***allocate_fftwf_3d(long N1,long N2,long N3);
float Hf(float aa),Df(float aa),ff(float aa);
void Setting_Up_Memory_For_Ro(float av);
void cic(float** rra);
void Get_Phi(int i);
void Update_v(float aa,float delta_aa,float **rra,float **vva);
void Update_x(float aa,float delta_aa,float **rra,float **vva);
void calpow_k(int f_flag, float kmin, float kmax, int Nbin,double* power, double* powerk, double* kmode, long *no);
void calpow(int f_flag,int Nbin,double* power, double* powerk, double* kmode,long *no);
void Zel_move_gradphi(float av,float **rra,float **vva);
void grad_phi(int ix);
int write_output(char *fname,long int seed,int output_flag,float **rra,float **vva,float vaa);
int read_output(char *fname, int read_flag,long int *seed,int *output_flag,int *in_flag,float **rra,float **vva,float *aa);
void read_fof(char *fname, int read_flag,int *output_flag, long *totcluster, float **halo, float *aa);
//Functions required to generate the power spectrum by Hu and Eisentine
void TFset_parameters(float omega0hh, float f_baryon, float Tcmb);
float TFfit_onek(float k, float *tf_baryon, float *tf_cdm);
float Pk(float kk); // power spectrum
float sigma_func(float kk); // for calculating sigma 8
//*************************************************************************
float simp(float (*func)(float),float a,float b,int N); // for integration
void Setting_Up_Memory_For_ionz();
void delta_fill(long*); // fills value of phases for Fourier modes
//*************************************************************************
void cic_vmass(float ***ro_dum, float **data,int xin,int yin,int zin,int min);
void smooth(float ***ro_dum,float Radii);
//void calpow_mom(float ***ro_dum,int Nbin,double* power, double* kmode, double* power_P2,double* power_P4);
void calpow_mom(float ***ro_dum,int Nbin,double* power, double* kmode, double* power_P2,double* power_P4, double* no);
void calpow_mom_k(float ***ro_dum,int Nbin,float kmin,float kmax,double* power, double* kmode, double* power_P2,double* power_P4, double *no);
void density_2_mass(float ***ro_dum, float **data,int xin,int yin,int zin,int min);
//*************************************************************************
typedef struct
{
long npart[6];
double mass[6];
double time;
double redshift;
int flag_sfr;
int flag_feedback;
long npartTotal[6];
int flag_cooling;
int num_files;
double BoxSize;
double Omega0;
double OmegaLambda;
double HubbleParam;
double Omegab;
double sigma_8_present;
long Nx;
long Ny;
long Nz;
float LL;
int output_flag;
int in_flag;
long int seed;
char fill[256- 6*4- 6*8- 2*8- 2*4- 6*4- 2*4 - 6*8 -6*4 -8]; /* fills to 256 Bytes */
} io_header;