From 54a419c3e03d1dede03cc195c7801a987b7d7620 Mon Sep 17 00:00:00 2001 From: Sandy Yuan Date: Tue, 16 Jan 2024 12:17:43 -0800 Subject: [PATCH] fix reseed bug --- abacusnbody/hod/abacus_hod.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/abacusnbody/hod/abacus_hod.py b/abacusnbody/hod/abacus_hod.py index 0571cfb0..1e36336e 100644 --- a/abacusnbody/hod/abacus_hod.py +++ b/abacusnbody/hod/abacus_hod.py @@ -567,7 +567,7 @@ def run_hod(self, tracers = None, want_rsd = True, want_nfw = False, NFW_draw = if len(self.halo_data['hveldev'].shape) == 1: self.halo_data['hveldev'] = r20*self.halo_data['hsigma3d']/np.sqrt(3) else: - self.halo_data['hveldev'] = r2*self.halo_data['hsigma3d']/np.sqrt(3) + self.halo_data['hveldev'] = r2*np.repeat(self.halo_data['hsigma3d'], 3).reshape((-1, 3))/np.sqrt(3) self.particle_data['prandoms'] = r3 print("gen randoms took, ", time.time() - start)