From 0b197e8e755885b98ae1f65d3b03b4ca687dbe8e Mon Sep 17 00:00:00 2001 From: Matteo Bachetti Date: Mon, 2 Sep 2024 13:04:33 +0200 Subject: [PATCH] Fix trapezoid call --- stingray/multitaper.py | 8 +++++++- 1 file changed, 7 insertions(+), 1 deletion(-) diff --git a/stingray/multitaper.py b/stingray/multitaper.py index 6de5bfab7..3b34da6fa 100644 --- a/stingray/multitaper.py +++ b/stingray/multitaper.py @@ -6,6 +6,12 @@ import warnings import numpy as np + +try: + from numpy import trapezoid +except ImportError: + from numpy import trapz as trapezoid + import scipy.optimize import scipy.stats from scipy import signal, interpolate @@ -460,7 +466,7 @@ def _get_adaptive_psd(self, freq_response, eigvals, max_iter=150): psd_est = self.psd_from_freq_response(freq_response, sqrt_eigvals[:, np.newaxis]) - var = np.trapz(psd_est, dx=np.pi / n_freqs) / (2 * np.pi) + var = trapezoid(psd_est, dx=np.pi / n_freqs) / (2 * np.pi) del psd_est psd = np.empty(n_freqs)