-
Notifications
You must be signed in to change notification settings - Fork 6
/
confirm_param.c
117 lines (96 loc) · 4.02 KB
/
confirm_param.c
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
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
//
// Check necessary parameters are initialized and print the parameters
//
#include <stdio.h>
#include "parameters.h"
#include "msg.h"
void confirm_parameters(Parameters* const param, const char filename[])
{
msg_printf(normal, "parameter file %s read.\n", filename);
if(param->nc <= 0)
msg_abort(11001, "nc is not given is the parameter file: %d\n", param->nc);
else
msg_printf(normal, "nc = %d\n", param->nc);
msg_printf(normal, "pm_nc_factor = %d\n", param->pm_nc_factor);
if(param->boxsize <= 0.0)
msg_abort(11002, "boxsize is not given in the parameter file: %f\n",
param->boxsize);
else
msg_printf(normal, "boxsize = %lg\n", param->boxsize);
msg_printf(normal, "a_final = %lg\n", param->a_final);
if(param->ntimestep <= 0)
msg_abort(11003, "ntimestep is not given in the parameter file: %d\n",
param->ntimestep);
else
msg_printf(normal, "ntimestep = %d\n", param->ntimestep);
if(param->n_zout == 0)
msg_abort(11004, "output_redshifts are not given in the parameter file.\n");
else {
msg_printf(normal, "output_redshift = {%lg", param->zout[0]);
for(int i=1; i<param->n_zout; i++)
msg_printf(normal, ", %lg", param->zout[i]);
msg_printf(normal, "}\n");
}
msg_printf(normal, "random_seed = %d\n", param->random_seed);
msg_printf(normal, "nrealization = %d\n", param->nrealization);
if(param->omega_m <= 0.0)
msg_abort(11005, "omega_m is not given in the parameter file.");
else
msg_printf(normal, "omega_m = %lg\n", param->omega_m);
if(param->h <= 0.0)
msg_abort(11006, "h is not given in the parameter file");
else
msg_printf(normal, "h = %lg\n", param->h);
if(param->sigma8 <= 0.0)
msg_abort(11007, "sigma8 is not given in the parameter file");
else
msg_printf(normal, "sigma8 (check) = %lg\n", param->sigma8);
msg_printf(normal, "np_alloc_factor = %lg\n", param->np_alloc_factor);
msg_printf(normal, "loglevel = %d\n", param->loglevel);
if(param->power_spectrum_filename == 0)
msg_abort(11008, "powerspectrum is not given in the parameter file\n");
else if(param->strlen_power_spectrum_filename <= 1)
msg_abort(11018, "Error: powerspectrum filename is empty\n");
else
msg_printf(normal, "powerspectrum = \"%s\"\n",
param->power_spectrum_filename);
if(param->fof_filename == 0)
msg_printf(normal, "No FoF halo finding. fof_filename is not given.\n");
else if(param->strlen_fof_filename <= 1)
msg_abort(11019, "Error: fof_filename is empty\n");
else {
msg_printf(normal, "fof_filename = \"%s\"\n", param->fof_filename);
msg_printf(normal, "linking_factor = %lg\n", param->fof_linking_factor);
}
if(param->snapshot_filename == 0)
msg_printf(normal, "No N-body particle output (snapshot)\n");
else if(param->strlen_snapshot_filename <= 1)
msg_abort(11020, "Error: snapshot filename is empty.\n");
else {
msg_printf(normal, "snapshot = \"%s\"\n", param->snapshot_filename);
msg_printf(normal, "write_longid = %d\n", param->write_longid);
}
if(param->subsample_filename == 0)
msg_printf(normal, "No subsampled particle output (subsample)\n");
else if(param->strlen_subsample_filename <= 1)
msg_abort(10021, "Error: subsample filename is empty.\n");
else {
msg_printf(normal, "subsample = \"%s\"\n", param->subsample_filename);
msg_printf(normal, "subsample_factor = %lg\n", param->subsample_factor);
}
if(param->cgrid_filename == 0)
msg_printf(normal, "No coarse density grid output (coarse_grid)\n");
else if(param->strlen_cgrid_filename <= 1)
msg_abort(10022, "coarse_grid filename is empty.\n");
else {
msg_printf(normal, "coarse_grid = \"%s\"\n", param->cgrid_filename);
msg_printf(normal, "coarse_grid_nc = %d\n", param->cgrid_nc);
}
if(param->init_filename) {
if(param->strlen_init_filename <= 1)
msg_abort(10023, "init filename is empty.\n");
msg_printf(normal, "inital= \"%s\"\n", param->init_filename);
if(param->snapshot_filename == 0)
msg_printf(normal, "write_longid= %d\n", param->write_longid);
}
}