Interpolating Discrete Halofields

[daviesje]

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.