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

Reshift zero

Browse source
This commit is contained in:
Guilhem Lavaux 2011-02-25 16:08:42 -05:00
parent 3f3ecfa8a2
commit 18a0b5d6cc
2 changed files with 7 additions and 4 deletions
Download patch file Download diff file Expand all files Collapse all files

7
mytools/generateMock.cpp
View file

@ -100,7 +100,7 @@ Interpolate make_cosmological_redshift(double OM, double OL, double z0, double z
return buildFromVector(pairs); return buildFromVector(pairs);
} }
void metricTransform(SimuData *data, int axis) void metricTransform(SimuData *data, int axis, bool reshift)
{ {
int x0, x1, x2; int x0, x1, x2;
@ -121,6 +121,7 @@ void metricTransform(SimuData *data, int axis)
Interpolate z_vs_D = make_cosmological_redshift(data->Omega_M, data->Omega_Lambda, 0., 2.0); // Redshift 2 should be sufficient ? Interpolate z_vs_D = make_cosmological_redshift(data->Omega_M, data->Omega_Lambda, 0., 2.0); // Redshift 2 should be sufficient ?
double z0 = 1/data->time - 1; double z0 = 1/data->time - 1;
double z_base = reshift ? z0 : 0;
TotalExpansion e_computer; TotalExpansion e_computer;
double baseComovingDistance; double baseComovingDistance;
@ -141,7 +142,7 @@ void metricTransform(SimuData *data, int axis)
double reduced_red = (z + baseComovingDistance)*100./LIGHT_SPEED; double reduced_red = (z + baseComovingDistance)*100./LIGHT_SPEED;
// Distorted redshift // Distorted redshift
z = z_vs_D.compute(reduced_red)*LIGHT_SPEED/100.; z = (z_vs_D.compute(reduced_red)-z_base)*LIGHT_SPEED/100.;
// Add peculiar velocity // Add peculiar velocity
z += v/100; z += v/100;
} }
@ -320,7 +321,7 @@ int main(int argc, char **argv)
cout << "Omega_M = " << simu->Omega_M << endl; cout << "Omega_M = " << simu->Omega_M << endl;
cout << "Omega_Lambda = " << simu->Omega_Lambda << endl; cout << "Omega_Lambda = " << simu->Omega_Lambda << endl;
metricTransform(simu, args_info.axis_arg); metricTransform(simu, args_info.axis_arg, args_info.preReShift_flag);
makeBox(simu, simuOut, args_info); makeBox(simu, simuOut, args_info);
delete simu; delete simu;

2
mytools/generateMock.ggo
View file

@ -18,3 +18,5 @@ option "rangeY_min" - "Minimum range in Y for making the box" double require
option "rangeY_max" - "Maximum range in Y for making the box" double required option "rangeY_max" - "Maximum range in Y for making the box" double required
option "rangeZ_min" - "Minimum range in Z for making the box (after distortion)" double required option "rangeZ_min" - "Minimum range in Z for making the box (after distortion)" double required
option "rangeZ_max" - "Maximum range in Z for making the box (after distortion)" double required option "rangeZ_max" - "Maximum range in Z for making the box (after distortion)" double required
option "preReShift" - "Reshift the zero of the Z axis" flag off