With just density, the l=2,4 multipoles are extremely small, and the GSL
integrator would abort the computation, so now I just turn it off for both
`compute_corrfunc` and `compute_multipole`.