allow older scipy too

number_density/gaussian_window/gaussian_window.py

@@ -1,7 +1,10 @@
 from cosmosis.datablock import option_section
 from cosmosis.datablock import names
 import numpy as np
-import scipy.integrate
+try:
+    from scipy.integrate import trapezoid
+except ImportError:
+    from scipy.integrate import trapz as trapezoid
 
 def setup(options):
     z = options[option_section, "z"]
@@ -43,7 +46,7 @@ def execute(block, config):
         nz_bin = gaussian(z, mu[i - 1], sigma[i - 1])
         # the bin may not quite go to zero before we get to the
         # edges so normalize it
-        norm = scipy.integrate.trapezoid(nz_bin, z)
+        norm = trapezoid(nz_bin, z)
         nz_bin /= norm
         # Save n(z)
         block[section, "BIN_%d" % i] = nz_bin

number_density/nz_multirank/nz_gz.py

@@ -12,7 +12,10 @@
 
 def nz_to_gchi(z, nz, cosmo=cosmology.Planck15):
     from numpy import apply_along_axis, gradient, multiply, newaxis
-    from scipy.integrate import cumulative_trapezoid
+    try:
+        from scipy.integrate import cumulative_trapezoid
+    except ImportError:
+        from scipy.integrate import cumtrapz as cumulative_trapezoid
     chi = cosmo.comoving_distance(z).value
     dchi = apply_along_axis(gradient, -1, chi)
     dz = apply_along_axis(gradient, -1, z)

structure/baccoemu/baccoemu_vendored/matter_powerspectrum.py

@@ -1578,7 +1578,10 @@ def _smeared_bao_pk(k_lin=None, pk_lin=None, k_emu=None, pk_lin_emu=None, pk_nw=
     :return: smeared BAO pk computed at k_emu
     :rtype: array_like
     """
-    from scipy.integrate import trapezoid
+    try:
+        from scipy.integrate import trapezoid
+    except ImportError:
+        from scipy.integrate import trapz as trapezoid
 
     if grid is None:
         sigma_star_2 = trapezoid(k_lin * pk_lin, x=np.log(k_lin)) / (3 * np.pi**2)