void StateSet::init_to_gaussian_noise(RNG& rng) { for (size_t n=0; n<N; n++) for (size_t y=0; y<datalayout.sizey; y++) for (size_t x=0; x<datalayout.sizex; x++) data(n,x,y) = comp(rng.gaussian_rand(), rng.gaussian_rand()); }