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

removed unused code to fill grids with random numbers

Browse source
This commit is contained in:
Oliver Hahn 2019-11-15 11:54:16 +01:00
parent 675ba19a95
commit 07b430f25c
2 changed files with 0 additions and 36 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

10
include/grid_fft.hh
View file

@ -240,16 +240,6 @@ public:
std::array<size_t,3> get_k3(const size_t i, const size_t j, const size_t k) const
{
return bdistributed? std::array<size_t,3>({j,i+local_1_start_,k}) : std::array<size_t,3>({i,j,k});
// vec3<size_t> kk;
// if( bdistributed ){
// kk[0] = j;
// kk[1] = i + local_1_start_;
// }else{
// kk[0] = i;
// kk[1] = j;
// }
// kk[2] = k;
// return kk;
}
data_t get_cic( const vec3<real_t>& v ) const{

26
src/grid_fft.cc
View file

@ -2,32 +2,6 @@
#include <grid_fft.hh>
#include <thread>
#include <gsl/gsl_rng.h>
#include <gsl/gsl_randist.h>
template <typename data_t,bool bdistributed>
void Grid_FFT<data_t,bdistributed>::FillRandomReal(unsigned long int seed)
{
gsl_rng *RNG = gsl_rng_alloc(gsl_rng_mt19937);
if( bdistributed ){
seed += 17321 * CONFIG::MPI_task_rank;
}
gsl_rng_set(RNG, seed);
for (size_t i = 0; i < sizes_[0]; ++i)
{
for (size_t j = 0; j < sizes_[1]; ++j)
{
for (size_t k = 0; k < sizes_[2]; ++k)
{
this->relem(i, j, k) = gsl_ran_ugaussian_ratio_method(RNG);
}
}
}
gsl_rng_free(RNG);
}
template <typename data_t,bool bdistributed>
void Grid_FFT<data_t,bdistributed>::Setup(void)
{