MNT: Handle trapz for numpy>=2 (#1971)

trapz has been deprecated in favor of the newly available trapezoid
function. This wraps the new function and avoids a DeprecationWarning on
numpy>=2.

pint/facets/numpy/numpy_func.py

@@ -741,22 +741,25 @@ def _base_unit_if_needed(a):
  raise OffsetUnitCalculusError(a.units)
 
 
+# Can remove trapz wrapping when we only support numpy>=2
 @implements("trapz", "function")
+@implements("trapezoid", "function")
 def _trapz(y, x=None, dx=1.0, **kwargs):
+ trapezoid = np.trapezoid if hasattr(np, "trapezoid") else np.trapz
  y = _base_unit_if_needed(y)
  units = y.units
  if x is not None:
  if hasattr(x, "units"):
  x = _base_unit_if_needed(x)
  units *= x.units
  x = x._magnitude
- ret = np.trapz(y._magnitude, x, **kwargs)
+ ret = trapezoid(y._magnitude, x, **kwargs)
  else:
  if hasattr(dx, "units"):
  dx = _base_unit_if_needed(dx)
  units *= dx.units
  dx = dx._magnitude
- ret = np.trapz(y._magnitude, dx=dx, **kwargs)
+ ret = trapezoid(y._magnitude, dx=dx, **kwargs)
 
  return y.units._REGISTRY.Quantity(ret, units)