This PR mainly contains improvements to the HMF tables used in the halo sampler.
For the inverse CDF tables used in the Sampling:
I replaced the interpolation on the CDF with a rootfind to get the inverse, this is more accurate so we require fewer bins.
The output is no longer log-mass, but the ratio of mass to condtition mass, this is slightly faster and makes the table interpolation a little more accurate.
Both tables for grid and halo conditions are now in log-probability, this should allow longer timesteps in the halo sampler (see #376).
Below the minimum log-probability I have added an extrapolation function, which assumes exponential decay in p(M)
The Sheth-Tormen conditional mass function has been fixed and optimized (fixing #370), The Delos+24 conditional mass function is available for the default case but does not perform well in the sampler due to the small timestep.
Alternate sampling methods (Sheth+1999 Partition and Parkinson+08 binary split) have been fixed (#374), they are less conforming to the given CMF but as far as I can tell they are working.
Old global parameters have been removed and many sampler related global parameters have been moved to the input structures.
This PR mainly contains improvements to the HMF tables used in the halo sampler.
For the inverse CDF tables used in the Sampling:
The Sheth-Tormen conditional mass function has been fixed and optimized (fixing #370), The Delos+24 conditional mass function is available for the default case but does not perform well in the sampler due to the small timestep.
Alternate sampling methods (Sheth+1999 Partition and Parkinson+08 binary split) have been fixed (#374), they are less conforming to the given CMF but as far as I can tell they are working.
Old global parameters have been removed and many sampler related global parameters have been moved to the input structures.