Add support for passing an array of scale factors into halo_model_base.__call__

[mattkwiecien]

In the cluster count likelihood, we evaluate the halo mass function within the integrand and therefore need this evaluation to be fast. We vectorized the code, but cannot fully because the halo mass function evaluation can only accept a float value for the scale factor. We would like this function call to also support passing an array of scale factor values.

Specifically, when we invoke

CCL/pyccl/halos/halo_model_base.py

Line 180 in 29d4697

def __call__(self, cosmo, M, a):

we can pass an array of mass M, but not an array of scale factors a. It would be great to have this added into CCL, and have the for loop over the scale factors be implemented in the C code.

I think the for loop that needs to be modified is here (?):

CCL/pyccl/ccl_sigM.i

Line 22 in 29d4697

void sigM_vec(ccl_cosmology * cosmo, double a,