IA> 0.5 is breaking the CCL code. Possibly, a way of circumventing this gsl breaking is changing the CCL precision parameters, e.g.,
pyccl.gsl_params["INTEGRATION_EPSREL"] = <somevalue>
pyccl.gsl_params["INTEGRATION_LIMBER_GAUSS_KRONROD_POINTS"] *=<somevalue>
I am opening this discussion if and how we should implement this modification of CCL precision through TJPCov configuration.