forgot to include quasirandom as an arg

fastpm · Apr 23, 2024 · 2c06681 · 2c06681
1 parent 9c2da13
commit 2c06681
Show file tree

Hide file tree

Showing 3 changed files with 9 additions and 4 deletions.
diff --git a/api/fastpm/thermalvelocity.h b/api/fastpm/thermalvelocity.h
@@ -11,6 +11,7 @@ typedef struct FastPMncdmInitData{
     double m_ncdm_sum; /* total eV mass of all ncdm */
     double z;          /* initialization redshift of ncdm species */
 
+    int quasirandom; /* whether velocity sampling is quasirandom or pseudorandom */
     int n_shells;    /* number of shells for splitting velocity magnitude distribution */
     FastPMncdmSphereScheme ncdm_sphere_scheme;
     int n_side;      /* n_side for healpix, also n_fib for fibonacci */
@@ -24,8 +25,8 @@ typedef struct FastPMncdmInitData{
 } FastPMncdmInitData;
 
 FastPMncdmInitData *
-fastpm_ncdm_init_create(double BoxSize, FastPMCosmology * c, double z, int n_shells, int n_side,
-                        int lvk, FastPMncdmSphereScheme ncdm_sphere_scheme);
+fastpm_ncdm_init_create(double BoxSize, FastPMCosmology * c, double z, int quasirandom, int n_shells,
+                        int n_side, int lvk, FastPMncdmSphereScheme ncdm_sphere_scheme);
 
 void
 fastpm_ncdm_init_free(FastPMncdmInitData* nid);

diff --git a/libfastpm/thermalvelocity.c b/libfastpm/thermalvelocity.c
@@ -356,6 +356,7 @@ fastpm_ncdm_init_create(
     double BoxSize,
     FastPMCosmology * c,
     double z,
+    int quasirandom,
     int n_shells,
     int n_side,
     int lvk,
@@ -375,7 +376,8 @@ fastpm_ncdm_init_create(
     nid->Omega_ncdm = Omega_ncdmTimesHubbleEaSq(1, c);
     nid->z = z;
     fastpm_info("ncdm reference redshift = %g\n", z);
-
+
+    nid->quasirandom = quasirandom;
     if (nid->quasirandom){
         nid->n_shells = n_shells;
         nid->n_side = n_side;

diff --git a/src/fastpm.c b/src/fastpm.c
@@ -740,7 +740,9 @@ prepare_ncdm(FastPMSolver * fastpm, RunData * prr, double a0, MPI_Comm comm)
     // init the nid
     FastPMncdmInitData* nid = fastpm_ncdm_init_create(
             boxsize,
-            fastpm->cosmology, 1 / CONF(prr->lua, time_step)[0] - 1, n_shell, n_side, lvk,
+            fastpm->cosmology, 1 / CONF(prr->lua, time_step)[0] - 1,
+            CONF(prr->lua, ncdm_quasirandom),
+            n_shell, n_side, lvk,
             CONF(prr->lua, ncdm_sphere_scheme));
 
     size_t total_np_ncdm_sites = nc_ncdm * nc_ncdm * nc_ncdm;