mirror/MUSIC
1
0
Fork 0
mirror of https://github.com/cosmo-sims/MUSIC.git synced 2024-09-19 17:03:46 +02:00
Code Releases Activity

fixed another bug with the hash-key random numbers, when coarsening of noise is needed

Browse source
This commit is contained in:
Oliver Hahn 2012-09-10 16:26:33 +02:00
parent 6e8234c49e
commit 449dbfc431
1 changed files with 6 additions and 3 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

9
random.cc
View file

@ -370,6 +370,7 @@ random_numbers<T>::random_numbers( /*const*/ random_numbers <T>& rc, bool kdegra
#endif
#endif
rnums_.push_back( new Meshvar<T>( res_, 0, 0, 0 ) );
cubemap_[0] = 0; // map all to single array
#pragma omp parallel for reduction(+:sum,sum2,count)
for( int i=0; i<nxc; i++ )
@ -412,7 +413,8 @@ random_numbers<T>::random_numbers( /*const*/ random_numbers <T>& rc, bool kdegra
//... use restriction to get consistent random numbers on coarser grid
mg_straight gop;
rnums_.push_back( new Meshvar<T>( res_, 0, 0, 0 ) );
rnums_.push_back( new Meshvar<T>( res_, 0, 0, 0 ) );
cubemap_[0] = 0; // map all to single array
gop.restrict( *rc.rnums_[0], *rnums_[0] );
#pragma omp parallel for reduction(+:sum,sum2,count)
@ -435,6 +437,7 @@ random_numbers<T>::random_numbers( /*const*/ random_numbers <T>& rc, bool kdegra
baseseed_ = -2;
rnums_.push_back( new Meshvar<T>( res_, 0, 0, 0 ) );
cubemap_[0] = 0;
double fac = 1.0/sqrt(8);
#pragma omp parallel for reduction(+:sum,sum2,count)
@ -723,8 +726,8 @@ double random_numbers<T>::fill_cube( int i, int j, int k)
size_t cubeidx = it->second;
if( rnums_[cubeidx] == NULL )
rnums_[cubeidx] = new Meshvar<T>( cubesize_, 0, 0, 0 );
rnums_[cubeidx] = new Meshvar<T>( cubesize_, 0, 0, 0 );
double mean = 0.0;
for( int ii=0; ii<(int)cubesize_; ++ii )