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

added rng allocation for joggling

Browse source
This commit is contained in:
P.M. Sutter 2013-05-06 09:41:55 -05:00
parent ab13d013e0
commit 726ad2d01e
1 changed files with 10 additions and 8 deletions
Download patch file Download diff file Expand all files Collapse all files

18
c_tools/mock/generateMock.cpp
View file

@ -154,15 +154,17 @@ void metricTransform(SimuData *data, int axis, bool reshift, bool pecvel, double
// slightly perturb particle positions
void joggleParticles(SimuData *data) {
gsl_rng *rng;
gsl_rng_set(rng, 31415926535);
cout << "Joggling particle positions..." << endl;
gsl_rng *myRng = gsl_rng_alloc(gsl_rng_taus);
int seed = 314159;
gsl_rng_set(myRng, seed);
for (uint32_t i = 0; i < data->NumPart; i++) {
data->Pos[0][i] += 1.e-3*gsl_rng_uniform(rng);
data->Pos[1][i] += 1.e-3*gsl_rng_uniform(rng);
data->Pos[2][i] += 1.e-3*gsl_rng_uniform(rng);
data->Pos[0][i] -= 1.e-3*gsl_rng_uniform(rng);
data->Pos[1][i] -= 1.e-3*gsl_rng_uniform(rng);
data->Pos[2][i] -= 1.e-3*gsl_rng_uniform(rng);
data->Pos[0][i] += 1.e-3*gsl_rng_uniform(myRng);
data->Pos[1][i] += 1.e-3*gsl_rng_uniform(myRng);
data->Pos[2][i] += 1.e-3*gsl_rng_uniform(myRng);
data->Pos[0][i] -= 1.e-3*gsl_rng_uniform(myRng);
data->Pos[1][i] -= 1.e-3*gsl_rng_uniform(myRng);
data->Pos[2][i] -= 1.e-3*gsl_rng_uniform(myRng);
}
} // end joggleParticles