possibility to call trapezoid for numpy >= 2.0

mflike/theoryforge.py

@@ -72,6 +72,11 @@
 from itertools import product
 
 import numpy as np
+if np.__version__ >= "2.0.0":
+  from numpy import trapezoid as trapz
+else:
+  from numpy import trapz
+
 from cobaya.log import LoggedError
 
 
@@ -214,7 +219,7 @@ def _bandpass_construction(self, **params):
                     tranb = _cmb2bb(nub)
                     # normalization integral to be evaluated at the shifted freqs
                     # in order to have cmb component calibrated to 1
-                    tranb_norm = np.trapezoid(_cmb2bb(nub), nub)
+                    tranb_norm = np.trapz(_cmb2bb(nub), nub)
                     self.bandint_freqs.append([nub, tranb / tranb_norm])
                 # in case we don't want to do band integration, e.g. when we have multifreq bandpass in sacc file
                 if self.bandint_nsteps == 1:
@@ -229,7 +234,7 @@ def _bandpass_construction(self, **params):
                     data_are_monofreq = True
                     self.bandint_freqs.append(nub[0])
                 else:
-                    trans_norm = np.trapezoid(bp * _cmb2bb(nub), nub)
+                    trans_norm = np.trapz(bp * _cmb2bb(nub), nub)
                     trans = bp / trans_norm * _cmb2bb(nub)
                     self.bandint_freqs.append([nub, trans])