guilhem_lavaux/vide_public
0
0
Fork 0
mirror of https://bitbucket.org/cosmicvoids/vide_public.git synced 2025-07-04 23:31:12 +00:00
Code Issues Projects Releases Packages Wiki Activity Actions

Made a more generic version for building table

Browse source
This commit is contained in:
Guilhem Lavaux 2011-02-22 15:35:25 -06:00
parent 9f3be29781
commit 6d338a022a
1 changed files with 4 additions and 4 deletions
Download patch file Download diff file Expand all files Collapse all files

8
mytools/generateMock.cpp
View file

@ -75,7 +75,7 @@ struct TotalExpansion
} }
}; };
Interpolate make_cosmological_redshift(double OM, double OL, double z0, double z1) Interpolate make_cosmological_redshift(double OM, double OL, double z0, double z1, int N = 1000)
{ {
TotalExpansion e_computer; TotalExpansion e_computer;
double D_tilde, Q, Qprime; double D_tilde, Q, Qprime;
@ -84,12 +84,12 @@ Interpolate make_cosmological_redshift(double OM, double OL, double z0, double z
e_computer.Omega_M = OM; e_computer.Omega_M = OM;
e_computer.Omega_L = OL; e_computer.Omega_L = OL;
pairs.resize(100); pairs.resize(N);
ofstream f("comoving_distance.txt"); ofstream f("comoving_distance.txt");
for (int i = 0; i < 100; i++) for (int i = 0; i < N; i++)
{ {
double z = z0 + (z1-z0)/100*i; double z = z0 + (z1-z0)/N*i;
pairs[i].second = z; pairs[i].second = z;
pairs[i].first = gslIntegrate(e_computer, 0, z, 1e-3); pairs[i].first = gslIntegrate(e_computer, 0, z, 1e-3);