Fast Methods for Cosmological Simulations
FastSim serves as a tool for quick N-body simulations in modified gravity.
Classes | Macros | Functions | Variables
ccl_test_distances_class_allz.c File Reference
#include "ccl.h"
#include "ctest.h"
#include <stdio.h>
#include <math.h>
Include dependency graph for ccl_test_distances_class_allz.c:

Go to the source code of this file.

Classes

struct  distances_class_data
 

Macros

#define DISTANCES_TOLERANCE   1.0e-6
 
#define N_MODEL   10
 
#define N_Z   10
 

Functions

static void read_benchmark_file (const char *filename, double z[10], double benchmark[10][10])
 
void __attribute__ ((weak))
 
static void compare_distances (int model, struct distances_class_data *data)
 
void __ctest_distances_class_model_1_run (struct distances_class_data *data)
 
void __ctest_distances_class_model_2_run (struct distances_class_data *data)
 
void __ctest_distances_class_model_3_run (struct distances_class_data *data)
 
void __ctest_distances_class_model_4_run (struct distances_class_data *data)
 
void __ctest_distances_class_model_5_run (struct distances_class_data *data)
 
void __ctest_distances_class_model_7_run (struct distances_class_data *data)
 
void __ctest_distances_class_model_8_run (struct distances_class_data *data)
 
void __ctest_distances_class_model_9_run (struct distances_class_data *data)
 
void __ctest_distances_class_model_10_run (struct distances_class_data *data)
 
void __ctest_distances_class_model_11_run (struct distances_class_data *data)
 

Variables

static struct distances_class_data __ctest_distances_class_data
 

Macro Definition Documentation

#define DISTANCES_TOLERANCE   1.0e-6

Definition at line 7 of file ccl_test_distances_class_allz.c.

Referenced by compare_distances().

#define N_MODEL   10

Definition at line 10 of file ccl_test_distances_class_allz.c.

Referenced by __attribute__().

#define N_Z   10

Definition at line 12 of file ccl_test_distances_class_allz.c.

Referenced by compare_distances(), and read_benchmark_file().

Function Documentation

void __attribute__ ( (weak)  )

Definition at line 66 of file ccl_test_distances_class_allz.c.

References ccl_mnu_list, ccl_test_distances::mnu, N_MODEL, ccl_test_distances::Neff, and read_benchmark_file().

66  {
67  // Values that are the same for all 5 models
68  data->Omega_c = 0.25;
69  data->Omega_b = 0.05;
70  data->h = 0.7;
71  data->A_s = 2.1e-9;
72  data->n_s = 0.96;
73  data->mnu_type = ccl_mnu_list;
74 
75  // Values that are different for the different models
76  double Omega_k[N_MODEL] = { 0.0, 0.0, 0.0, 0.05, -0.05,
77  0.0, 0.0, 0.0, 0.05, -0.05 };
78  double w_0[N_MODEL] = { -1.0, -0.9, -0.9, -0.9, -0.9,
79  -1.0, -0.9, -0.9, -0.9, -0.9 };
80  double w_a[N_MODEL] = { 0.0, 0.0, 0.1, 0.1, 0.1,
81  0.0, 0.0, 0.1, 0.1, 0.1 };
82 
83  double mnu[N_MODEL][3] = { {0.0, 0.0, 0.0},
84  {0.0, 0.0, 0.0},
85  {0.0, 0.0, 0.0},
86  {0.0, 0.0, 0.0},
87  {0.0, 0.0, 0.0},
88  {0.04, 0.0, 0.0},
89  {0.05, 0.01, 0.0},
90  {0.03, 0.02, 0.04},
91  {0.05, 0.0, 0.0},
92  {0.03, 0.02, 0.0} };
93 
94  double Neff[N_MODEL] = {3.046, 3.046, 3.046, 3.046, 3.046,
95  3.013, 3.026, 3.040, 3.013, 3.026};
96 
97 
98  // Fill in the values from these constant arrays.
99  for(int i=0; i<N_MODEL; i++) {
100  data->Omega_k[i] = Omega_k[i];
101  data->w_0[i] = w_0[i];
102  data->w_a[i] = w_a[i];
103  data->Neff[i] = Neff[i];
104  for(int j=0; j<3; j++) {
105  data->mnu[i][j] = mnu[i][j];
106  }
107  }
108 
109  // The file of benchmark data.
110  read_benchmark_file("./tests/benchmark/chi_class_allz.txt", data->z_chi, data->chi_benchmark);
111  read_benchmark_file("./tests/benchmark/dm_class_allz.txt", data->z_dm, data->dm_benchmark);
112 }
ccl_test_distances.mnu
list mnu
Definition: ccl_test_distances.py:32
ccl_test_distances.Neff
int Neff
Definition: ccl_test_distances.py:25
ccl_mnu_list
Definition: ccl_core.h:143
read_benchmark_file
static void read_benchmark_file(const char *filename, double z[10], double benchmark[10][10])
Definition: ccl_test_distances_class_allz.c:36
N_MODEL
#define N_MODEL
Definition: ccl_test_distances_class_allz.c:10
void __ctest_distances_class_model_10_run ( struct distances_class_data *  data)

Definition at line 199 of file ccl_test_distances_class_allz.c.

References compare_distances().

199  {
200  int model = 8;
201  compare_distances(model, data);
202 }
compare_distances
static void compare_distances(int model, struct distances_class_data *data)
Definition: ccl_test_distances_class_allz.c:114
void __ctest_distances_class_model_11_run ( struct distances_class_data *  data)

Definition at line 204 of file ccl_test_distances_class_allz.c.

References compare_distances().

204  {
205  int model = 9;
206  compare_distances(model, data);
207 }
compare_distances
static void compare_distances(int model, struct distances_class_data *data)
Definition: ccl_test_distances_class_allz.c:114
void __ctest_distances_class_model_1_run ( struct distances_class_data *  data)

Definition at line 159 of file ccl_test_distances_class_allz.c.

References compare_distances().

159  {
160  int model = 0;
161  compare_distances(model, data);
162 }
compare_distances
static void compare_distances(int model, struct distances_class_data *data)
Definition: ccl_test_distances_class_allz.c:114
void __ctest_distances_class_model_2_run ( struct distances_class_data *  data)

Definition at line 164 of file ccl_test_distances_class_allz.c.

References compare_distances().

164  {
165  int model = 1;
166  compare_distances(model, data);
167 }
compare_distances
static void compare_distances(int model, struct distances_class_data *data)
Definition: ccl_test_distances_class_allz.c:114
void __ctest_distances_class_model_3_run ( struct distances_class_data *  data)

Definition at line 169 of file ccl_test_distances_class_allz.c.

References compare_distances().

169  {
170  int model = 2;
171  compare_distances(model, data);
172 }
compare_distances
static void compare_distances(int model, struct distances_class_data *data)
Definition: ccl_test_distances_class_allz.c:114
void __ctest_distances_class_model_4_run ( struct distances_class_data *  data)

Definition at line 174 of file ccl_test_distances_class_allz.c.

References compare_distances().

174  {
175  int model = 3;
176  compare_distances(model, data);
177 }
compare_distances
static void compare_distances(int model, struct distances_class_data *data)
Definition: ccl_test_distances_class_allz.c:114
void __ctest_distances_class_model_5_run ( struct distances_class_data *  data)

Definition at line 179 of file ccl_test_distances_class_allz.c.

References compare_distances().

179  {
180  int model = 4;
181  compare_distances(model, data);
182 }
compare_distances
static void compare_distances(int model, struct distances_class_data *data)
Definition: ccl_test_distances_class_allz.c:114
void __ctest_distances_class_model_7_run ( struct distances_class_data *  data)

Definition at line 184 of file ccl_test_distances_class_allz.c.

References compare_distances().

184  {
185  int model = 5;
186  compare_distances(model, data);
187 }
compare_distances
static void compare_distances(int model, struct distances_class_data *data)
Definition: ccl_test_distances_class_allz.c:114
void __ctest_distances_class_model_8_run ( struct distances_class_data *  data)

Definition at line 189 of file ccl_test_distances_class_allz.c.

References compare_distances().

189  {
190  int model = 6;
191  compare_distances(model, data);
192 }
compare_distances
static void compare_distances(int model, struct distances_class_data *data)
Definition: ccl_test_distances_class_allz.c:114
void __ctest_distances_class_model_9_run ( struct distances_class_data *  data)

Definition at line 194 of file ccl_test_distances_class_allz.c.

References compare_distances().

194  {
195  int model = 7;
196  compare_distances(model, data);
197 }
compare_distances
static void compare_distances(int model, struct distances_class_data *data)
Definition: ccl_test_distances_class_allz.c:114
static void compare_distances ( int  model,
struct distances_class_data *  data 
)
static

Definition at line 114 of file ccl_test_distances_class_allz.c.

References distances_class_data::A_s, ASSERT_DBL_NEAR_TOL, ASSERT_EQUAL, ASSERT_NOT_NULL, ccl_comoving_radial_distance(), ccl_cosmology_create(), ccl_cosmology_free(), ccl_distance_modulus(), ccl_parameters_create(), distances_class_data::chi_benchmark, cl_cmbl_bm::cosmo, default_config, DISTANCES_TOLERANCE, distances_class_data::dm_benchmark, distances_class_data::h, distances_class_data::mnu, distances_class_data::mnu_type, distances_class_data::n_s, N_Z, distances_class_data::Neff, distances_class_data::Omega_b, distances_class_data::Omega_c, distances_class_data::Omega_k, halomod_bm::params, ccl_cosmology::status_message, distances_class_data::w_0, distances_class_data::w_a, distances_class_data::z_chi, and distances_class_data::z_dm.

Referenced by __ctest_distances_class_model_10_run(), __ctest_distances_class_model_11_run(), __ctest_distances_class_model_1_run(), __ctest_distances_class_model_2_run(), __ctest_distances_class_model_3_run(), __ctest_distances_class_model_4_run(), __ctest_distances_class_model_5_run(), __ctest_distances_class_model_7_run(), __ctest_distances_class_model_8_run(), and __ctest_distances_class_model_9_run().

115 {
116  int status=0;
117  // Make the parameter set from the input data
118  // Values of some parameters depend on the model index
119 
120  // The arrays of massive neutrions are different lengths, so we need to have an if-statement here to deal with that.
121 
122  ccl_parameters params;
123 
124  params = ccl_parameters_create(data->Omega_c, data->Omega_b, data->Omega_k[model],
125  data->Neff[model], data->mnu[model], data->mnu_type,
126  data->w_0[model], data->w_a[model],
127  data->h, data->A_s, data->n_s,
128  -1,-1,-1,-1,NULL,NULL, &status);
129  ASSERT_EQUAL(0, status);
130 
131  // Make a cosmology object from the parameters with the default configuration
132  ccl_cosmology * cosmo = ccl_cosmology_create(params, default_config);
133  ASSERT_NOT_NULL(cosmo);
134 
135  // Compare to benchmark data
136  for(int i=0; i<N_Z; i++) {
137  // Check comoving radial distance
138  double a = 1/(1.+data->z_chi[i]);
139  double chi_ccl = ccl_comoving_radial_distance(cosmo, a, &status);
140  if(status) printf("%s\n",cosmo->status_message);
141 
142  double absolute_tolerance = DISTANCES_TOLERANCE*data->chi_benchmark[i][model];
143  if (fabs(absolute_tolerance)<1e-12) absolute_tolerance = 1e-12;
144  ASSERT_DBL_NEAR_TOL(data->chi_benchmark[i][model], chi_ccl, absolute_tolerance);
145 
146  // Check distance modulus
147  a = 1/(1.+data->z_dm[i]);
148  double dm_ccl = ccl_distance_modulus(cosmo, a, &status);
149  if(status) printf("%s\n",cosmo->status_message);
150 
151  absolute_tolerance = DISTANCES_TOLERANCE*data->dm_benchmark[i][model];
152  if (fabs(absolute_tolerance)<1e-12) absolute_tolerance = 1e-12;
153  ASSERT_DBL_NEAR_TOL(data->dm_benchmark[i][model], dm_ccl, absolute_tolerance);
154  }
155 
156  ccl_cosmology_free(cosmo);
157 }
N_Z
#define N_Z
Definition: ccl_test_distances_class_allz.c:12
distances_class_data::Neff
double Neff[10]
Definition: ccl_test_distances_class_allz.c:20
distances_class_data::w_a
double w_a[10]
Definition: ccl_test_distances_class_allz.c:25
ccl_parameters
Definition: ccl_core.h:16
distances_class_data::Omega_k
double Omega_k[10]
Definition: ccl_test_distances_class_allz.c:23
distances_class_data::n_s
double n_s
Definition: ccl_test_distances_class_allz.c:19
pz_params::status
int * status
Definition: ccl_redshifts.c:41
distances_class_data::h
double h
Definition: ccl_test_distances_class_allz.c:17
distances_class_data::A_s
double A_s
Definition: ccl_test_distances_class_allz.c:18
distances_class_data::Omega_b
double Omega_b
Definition: ccl_test_distances_class_allz.c:16
distances_class_data::z_dm
double z_dm[10]
Definition: ccl_test_distances_class_allz.c:28
ccl_cosmology_create
ccl_cosmology * ccl_cosmology_create(ccl_parameters params, ccl_configuration config)
Definition: ccl_core.c:173
distances_class_data::w_0
double w_0[10]
Definition: ccl_test_distances_class_allz.c:24
distances_class_data::chi_benchmark
double chi_benchmark[10][10]
Definition: ccl_test_distances_class_allz.c:30
ccl_cosmology_free
void ccl_cosmology_free(ccl_cosmology *cosmo)
Definition: ccl_core.c:829
halomod_bm.params
dictionary params
Definition: halomod_bm.py:27
distances_class_data::mnu_type
ccl_mnu_convention mnu_type
Definition: ccl_test_distances_class_allz.c:22
ASSERT_DBL_NEAR_TOL
#define ASSERT_DBL_NEAR_TOL(exp, real, tol)
Definition: ctest.h:152
distances_class_data::Omega_c
double Omega_c
Definition: ccl_test_distances_class_allz.c:15
DISTANCES_TOLERANCE
#define DISTANCES_TOLERANCE
Definition: ccl_test_distances_class_allz.c:7
distances_class_data::z_chi
double z_chi[10]
Definition: ccl_test_distances_class_allz.c:27
default_config
const ccl_configuration default_config
Definition: ccl_core.c:21
cl_cmbl_bm.cosmo
cosmo
Definition: cl_cmbl_bm.py:10
ccl_cosmology
Definition: ccl_core.h:122
ASSERT_EQUAL
#define ASSERT_EQUAL(exp, real)
Definition: ctest.h:121
ASSERT_NOT_NULL
#define ASSERT_NOT_NULL(real)
Definition: ctest.h:139
ccl_distance_modulus
double ccl_distance_modulus(ccl_cosmology *cosmo, double a, int *status)
Definition: ccl_background.c:871
ccl_parameters_create
ccl_parameters ccl_parameters_create(double Omega_c, double Omega_b, double Omega_k, double Neff, double *mnu, ccl_mnu_convention mnu_type, double w0, double wa, double h, double norm_pk, double n_s, double bcm_log10Mc, double bcm_etab, double bcm_ks, int nz_mgrowth, double *zarr_mgrowth, double *dfarr_mgrowth, int *status)
Definition: ccl_core.c:294
ccl_comoving_radial_distance
double ccl_comoving_radial_distance(ccl_cosmology *cosmo, double a, int *status)
Definition: ccl_background.c:752
distances_class_data::mnu
double mnu[10][3]
Definition: ccl_test_distances_class_allz.c:21
ccl_cosmology::status_message
char status_message[500]
Definition: ccl_core.h:136
distances_class_data::dm_benchmark
double dm_benchmark[10][10]
Definition: ccl_test_distances_class_allz.c:31
static void read_benchmark_file ( const char *  filename,
double  z[10],
double  benchmark[10][10] 
)
static

Definition at line 36 of file ccl_test_distances_class_allz.c.

References ASSERT_EQUAL, ASSERT_NOT_NULL, Catch::Matchers::Vector::Detail::count(), N_Z, and z.

Referenced by __attribute__().

37 {
38  //Distances are in Mpc
39  FILE * f = fopen(filename, "r");
40  ASSERT_NOT_NULL(f);
41 
42  // Ignore header line
43  fscanf(f, "%*[^\n]\n", NULL);
44 
45  double val;
46  // Read the file
47  for(int i_z=0; i_z<N_Z; i_z++) {
48  for(int i_model=0; i_model<N_MODEL+1; i_model++) {
49  int count = fscanf(f, "%le", &val);
50  // Check that all the value was read successfully
51  ASSERT_EQUAL(1, count);
52  if(i_model == 0) {
53  // The first column holds the redshift values
54  z[i_z] = val;
55  }
56  else {
57  benchmark[i_z][i_model-1] = val;
58  }
59  }
60  }
61  fclose(f);
62 }
N_Z
#define N_Z
Definition: ccl_test_distances_class_allz.c:12
Catch::Matchers::Vector::Detail::count
size_t count(InputIterator first, InputIterator last, T const &item)
Definition: catch.hpp:2747
z
static double z[8]
Definition: ccl_emu17_params.h:104517
ASSERT_EQUAL
#define ASSERT_EQUAL(exp, real)
Definition: ctest.h:121
ASSERT_NOT_NULL
#define ASSERT_NOT_NULL(real)
Definition: ctest.h:139
N_MODEL
#define N_MODEL
Definition: ccl_test_distances_class_allz.c:10

Variable Documentation

static struct distances_class_data __ctest_distances_class_data
static

Definition at line 159 of file ccl_test_distances_class_allz.c.