def compute_integral(self, x_s, y_s): if len(x_s) == 0: return np.zeros((y_s.shape[0],)) * np.nan closer = np.nanargmin(abs(x_s - self.limits[0])) return y_s[:, closer]