21cmfast / 21cmFAST

Official repository for 21cmFAST: a code for generating fast simulations of the cosmological 21cm signal
MIT License
56 stars 37 forks source link

Interpolating Discrete Halofields #395

Open daviesje opened 1 month ago

daviesje commented 1 month ago

Currently, the implementation of the halo quantities on the lightcone is the same as all other fields. The gridded quantities in HaloBox are interpolated to the lightcone slice redshift.

An improved approach would be to call perturb_halo_field at the slice redshift, gridding these moved halos, and then taking the correct slice. which would avoid the case of a single halo being spread across two cells.

It would be far too slow to call determine_halo_list at each lightcone slice so there may be some difficulty in determining the halo masses and rng between the two snapshots, since interpolating between two halo catalogues is not easy.