From d9c7898df45e647da01df1f9a2852e7bffd2f501 Mon Sep 17 00:00:00 2001 From: Yu Feng Date: Fri, 25 May 2018 17:22:54 -0700 Subject: [PATCH] using non-transposed ParticleMesh object. The non-transposed PM object is way faster than the transposed PM object in filling gaussian fields. O(Nmesh) vs O(Nmesh^2 / Np**0.5) This depends on https://github.com/rainwoodman/pmesh/pull/33 --- nbodykit/source/catalog/lognormal.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/nbodykit/source/catalog/lognormal.py b/nbodykit/source/catalog/lognormal.py index c2c8ed5d0..b48e80385 100644 --- a/nbodykit/source/catalog/lognormal.py +++ b/nbodykit/source/catalog/lognormal.py @@ -130,7 +130,7 @@ def _makesource(self, BoxSize, Nmesh): # the particle mesh for gridding purposes _Nmesh = numpy.empty(3, dtype='i8') _Nmesh[:] = Nmesh - pm = ParticleMesh(BoxSize=BoxSize, Nmesh=_Nmesh, dtype='f4', comm=self.comm) + pm = ParticleMesh(BoxSize=BoxSize, Nmesh=_Nmesh, dtype='f4', comm=self.comm, transposed=False) # growth rate to do RSD in the Zel'dovich approx f = self.cosmo.scale_independent_growth_rate(self.attrs['redshift'])